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