c=[1.5  2]';
A=[1  1;
   1  3;
   2  5];
b=[550  1000  2000]';
lb=[0  0]';
ub=[400  inf]';
ctype="UUU";
vartype="CC";
s=-1; 
param.msglev=1;
param.itlim=100;
[xmin,fmin,status,extra]=glpk(c,A,b,lb,ub,ctype,vartype,s, param);
