c=[8 6 10 9 9 12 13 7 14 9 16 5]; A=[1 1 1 1 0 0 0 0 0 0 0 0; 0 0 0 0 1 1 1 1 0 0 0 0; 0 0 0 0 0 0 0 0 1 1 1 1; -1 0 0 0 -1 0 0 0 -1 0 0 0; 0 -1 0 0 0 -1 0 0 0 -1 0 0; 0 0 -1 0 0 0 -1 0 0 0 -1 0; 0 0 0 -1 0 0 0 -1 0 0 0 -1]; b=[35 50 40 -45 -20 -30 -30]; lb=[0 0 0 0 0 0 0 0 0 0 0 0]; [X,z,exitflag,output,lambda]=linprog(c,A,b,[],[],lb);