example3
>load lpsolve
LPSOLVE package (ported by Peter Notebaert), and ilpsolve() routine.
The full documentation is available on
http://lpsolve.sourceforge.net/5.5/Euler.htm
>X = 0.1:0.05:125;
>Y1 = (15000. - 120*X)/210;
>plot2d(0|X,0|Y1,filled=true);
>insimg;

>X = 0.1:0.05:38;
>Y1 = (15000. - 120*X)/210;
>Y2 = max((4000 - 110.*X)./30, 0);
>Y3 = max(75 - X, 0.);
>Ytop = min(min(Y1, Y2), Y3);
>plot2d(0|X,0|Ytop,filled=true);
>insimg;

>X = 0.1:0.05:38;
>Y1 = (15000. - 120*X)/210;
>Y2 = max((4000 - 110.*X)./30, 0);
>Y3 = max(75 - X, 0.);
>Ytop = min(min(Y1, Y2), Y3);
>plot2d(0|X,0|Ytop,filled=true);
>n=(1000:1000:9000)';
>plot2d("(n-143*x)/60", add=1, color=2);
>insimg;

>x = [1, 1; 110, 30] \ [75; 4000]
21.875
53.125
>P = [143, 60] . x
6315.625