Octave glpk help

[elektroneg]

Hi!
I try to write some simple code for flux balance analysis, using glpk as a
solver and it gives me an error: 214 
"No dual feasible solution (LP presolver)".
Here's example that I tried to solve:

a = [1 -1 0 0; 0 1 -1 1; 1 0 0 0];
c = [0 0 1 0]';
b = [0 0 1]';
vartype = "CCCC";
ctype = "SSS";
sense = -1;
ub = [];
lb = [];
[xopt, fopt, status] = glpk(c,a,b,lb,ub,ctype,vartype,s)

The solution should be x3 = 3.. Please tell me what am I doin wrong? Thank
you!

--
View this message in context: 
http://octave.1599824.n4.nabble.com/Octave-glpk-help-tp4630979.html
Sent from the Octave - General mailing list archive at Nabble.com.