ex

// Test and Demonstrate the use of lpsolve

>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 = [-1, 2]';
>A = [2, 1;-4, 4];
>b = [5, 5];
>e = -[1, 1];
>xint = [1, 2];
>{v,x} = lpsolve(f,A,b,e,[],[],xint);
>v
3
>x
                  1 
                  2 
>f = [50, 100];
>A = [10, 5;4, 10; 1, 1.5];
>b = [2500, 2000, 450];
>e = [-1, -1, -1];
>{v,x} = lpsolve(f,A,b,e);
>v
21875
>x
              187.5 
                125 
>f = -[40, 36];
>vub = [8, 10];
>A = [5, 3];
>b = [45];
>e = 1;
>{v,x} = lpsolve(f,A,b,e,[],vub);
>v
-380
>x
                  8 
      1.66666666667 
>f = [10, 6, 4];
>A = [1, 1, 1;10, 4, 5;2, 2, 6];
>b = [100, 600, 300];
>e = [-1, -1, -1];
>xint = [2];
>{v,x} = lpsolve(f,A,b,e,[],[],xint);
>v
732
>x
                 33 
                 67 
                  0 
>f = -[3, -7, -12];
>a = [-3, 6, 8;6, -3, 7;-6, 3, 3];
>b = [12, 8, 5];
>e = [-1, -1, -1];
>xint = [1, 2, 3];
>{v,x} = lpsolve(f,a,b,e,[],[],xint);
>v
15
>x
                  2 
                  3 
                  0 
>f = -[2, 3, 7, 7];
>a = [1, 1, -2, -5;-1, 2, 1, 4];
>b = [2, -3];
>e = [1, 1];
>xint = [1, 2, 3, 4];
>vub = [1, 1, 1, 1];
>{v,x} = lpsolve(f,a,b,e,[],vub,xint);
>v
-5
>x
                  1 
                  1 
                  0 
                  0 
>f = -[1, 2, 3, 7, 8, 8];
>a = [5, -3, 2, -3, -1, 2; -1, 0, 2, 1, 3, -3;1, 2, -1, 0, 5, -1];
>b = [-5, -1, 3];
>e = [1, 1, 1];
>xint = [1, 2, 3, 4, 5, 6];
>vub = [1, 1, 1, 1, 1, 1];
>{v,x} = lpsolve(f,a,b,e,[],vub,xint);
>v
-3
>x
                  1 
                  1 
                  0 
                  0 
                  0 
                  0 
>f=[8, 15];
>a = [10, 21;2, 1];
>b = [156, 22];
>e = [-1, -1];
>{v,x} = lpsolve(f,a,b,e);
>v
119.625
>x
             9.5625 
              2.875 
>f=[3, 13];
>a = [2, 9;11, -8];
>b = [40, 82];
>e = [-1, -1];
>{v,x} = lpsolve(f,a,b,e);
>v
58.8
>x
                9.2 
                2.4 
>f=[592, 381, 273, 55, 48, 37, 23];
>a = [3534, 2356, 1767, 589, 528, 451, 304];
>b = [119567];
>e = [-1];
>xint = [1, 2, 3, 4, 5, 6, 7];
>vub = [];
>{v,x} = lpsolve(f,a,b,e,[],vub,xint);
>v
19979
>x
                 32 
                  2 
                  1 
                  0 
                  0 
                  0 
                  0