load("simplex");

u1:x>=0;

u2:y>=0;

u3:x<=100;

u4:y<=150;

u5:-4*x+y<=0;

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

zf:-40*x+70*y+43500;

maximize_lp(zf,NB);