u1:x>=0;

u2:y>=0;

u3:z>=0;

u4:x+y<=100;

u5:2*x+y<=150;

u6:x+y+z<=200;

u7:y+z<=120;

NB:[u1,u2,u3,u4,u5,u6,u7];

zf:x+2*y+z;

load("simplex");

maximize_lp(zf,NB);