c=[-200 -400];
A=[1/40 1/60;
    1/50 1/50];
b=[1 1];
lb=[0,0];
options=optimset('Display','iter');
[X,z,exitflag,output,lambda]=linprog(c,A,b,[],[],lb,[],[],options);