u1:x>=0;

u2:y>=0;

u3:z>=0;

u4:w>=0;

u5:x+y+z+w<=400;

u6:x+w>=230;

u7:x+y<=200;

u8:z+w<=180;

u9:y+z+2*w<=300;

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

load(simplex);

zf:3*x+y+z+4*w;

maximize_lp(zf,NB);