example6

>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
>f = [110*1.3, 30*2.0, 125*1.56, 75*1.8, 95*.95, 100*2.25, 50*1.35];
>A = [120,210,150.75,115,186,140,85; 110,30,125,75,95,100,50; 1,1,1,1,1,1,1; 1,-1,0,0,0,0,0; 0,0,1,0,-2,0,0; 0,0,0,-1,0,-1,1];
>b = [55000; 40000; 400; 0; 0; 0];
>lp = lpmaker(f, A, b, [-1,-1,-1,-1,-1,-1], ...
>  [10,10,10,10,20,20,20],[100,1.0e30,50,1.0e30,1.0e30,250,1.0e30],[],1,0);
>solvestat = eulpsolve("solve", lp);
>obj = eulpsolve("get_objective", lp)
75398.0434783
>x = eulpsolve("get_variables", lp)
                 10 
                 10 
                 40 
       45.652173913 
                 20 
                250 
                 20 
>eulpsolve("delete_lp", lp);