> with(LinearAlgebra):
> A := Matrix(3,5,[5,2,1,0,0, 2,1,0,1,0,  1,1,0,0,1]);
                                          [5    2    1    0    0]
                                          [                     ]
                                     A := [2    1    0    1    0]
                                          [                     ]
                                          [1    1    0    0    1]

> b := Matrix(3,1,[1000, 450,400]);
                                                  [1000]
                                                  [    ]
                                             b := [ 450]
                                                  [    ]
                                                  [ 400]

> c := Matrix(1,5,[80,50,0,0,0]);
                                    c := [80    50    0    0    0]

> z:=0;
                                                z := 0

> read lppivot:  # Suppress output with colon instead of semicolon
> dict(x3,x4,x5);

     z   =     0.00  +  80.00 x1  +  50.00 x2 
    ---------------------------------------------
     x3  =  1000.00  -   5.00 x1  -   2.00 x2 
     x4  =   450.00  -   2.00 x1  -   1.00 x2 
     x5  =   400.00  -   1.00 x1  -   1.00 x2 

> ratios(1);

    row  1 (x3 basic):  Upper bound = 200.0000 
    row  2 (x4 basic):  Upper bound = 225.0000 
    row  3 (x5 basic):  Upper bound = 400.0000 

> lppivot(3,1);

     z   = 16000.00  +  18.00 x2  -  16.00 x3 
    ---------------------------------------------
     x1  =   200.00  -   0.40 x2  -   0.20 x3 
     x4  =    50.00  -   0.20 x2  +   0.40 x3 
     x5  =   200.00  -   0.60 x2  +   0.20 x3 

> ratios(2);

    row  1 (x1 basic):  Upper bound = 500.0000 
    row  2 (x4 basic):  Upper bound = 250.0000 
    row  3 (x5 basic):  Upper bound = 333.3333 

> lppivot(4,2);

     z   = 20500.00  +  20.00 x3  -  90.00 x4 
    ---------------------------------------------
     x1  =   100.00  -   1.00 x3  +   2.00 x4 
     x2  =   250.00  +   2.00 x3  -   5.00 x4 
     x5  =    50.00  -   1.00 x3  +   3.00 x4 

> ratios(3);

    row  1 (x1 basic):  Upper bound = 100.0000 
    row  2 (x2 basic):  No upper bound 
    row  3 (x5 basic):  Upper bound =  50.0000 

> lppivot(5,3);

     z   = 21500.00  -  30.00 x4  -  20.00 x5 
    ---------------------------------------------
     x1  =    50.00  -   1.00 x4  +   1.00 x5 
     x2  =   350.00  +   1.00 x4  -   2.00 x5 
     x3  =    50.00  +   3.00 x4  -   1.00 x5 

>