diff --git a/books/bookvol10.4.pamphlet b/books/bookvol10.4.pamphlet
index 7d44e73..f37da6a 100644
--- a/books/bookvol10.4.pamphlet
+++ b/books/bookvol10.4.pamphlet
@@ -102886,11 +102886,1755 @@ NagMatrixOperationsPackage(): Exports == Implementation where
 )set message auto off
 )clear all
 
---S 1 of 1
+--S 1 of 153
 )show NagOptimisationPackage
+--R 
+--R NagOptimisationPackage  is a package constructor
+--R Abbreviation for NagOptimisationPackage is NAGE04 
+--R This constructor is exposed in this frame.
+--R Issue )edit bookvol10.4.pamphlet to see algebra source code for NAGE04 
+--R
+--R------------------------------- Operations --------------------------------
+--R e04dgf : (Integer,DoubleFloat,DoubleFloat,Integer,DoubleFloat,Boolean,DoubleFloat,DoubleFloat,Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Integer,Union(fn: FileName,fp: Asp49(OBJFUN))) -> Result
+--R e04fdf : (Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Integer,Union(fn: FileName,fp: Asp50(LSFUN1))) -> Result
+--R e04gcf : (Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Integer,Union(fn: FileName,fp: Asp19(LSFUN2))) -> Result
+--R e04jaf : (Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Integer,Union(fn: FileName,fp: Asp24(FUNCT1))) -> Result
+--R e04mbf : (Integer,Integer,Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Boolean,Integer,Integer,Matrix(DoubleFloat),Integer) -> Result
+--R e04naf : (Integer,Integer,Integer,Integer,Integer,Integer,Integer,Integer,DoubleFloat,Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Boolean,Boolean,Boolean,Integer,Integer,Matrix(DoubleFloat),Matrix(Integer),Integer,Union(fn: FileName,fp: Asp20(QPHESS))) -> Result
+--R e04ucf : (Integer,Integer,Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Integer,Integer,Boolean,DoubleFloat,Integer,DoubleFloat,DoubleFloat,Boolean,DoubleFloat,DoubleFloat,DoubleFloat,DoubleFloat,Boolean,Integer,Integer,Integer,Integer,Integer,DoubleFloat,DoubleFloat,DoubleFloat,Integer,Integer,Integer,Integer,Integer,Matrix(Integer),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Matrix(DoubleFloat),Integer,Union(fn: FileName,fp: Asp55(CONFUN)),Union(fn: FileName,fp: Asp49(OBJFUN))) -> Result
+--R e04ycf : (Integer,Integer,Integer,DoubleFloat,Matrix(DoubleFloat),Integer,Matrix(DoubleFloat),Integer) -> Result
+--R
 --E 1
 
-)spool
+--S 2 of 153
+digits 20
+--R 
+--R
+--R   (1)  20
+--R                                                        Type: PositiveInteger
+--E 2
+
+--S 3 of 153
+showArrayValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 3
+
+--S 4 of 153
+showScalarValues true
+--R 
+--R
+--R   (3)  true
+--R                                                                Type: Boolean
+--E 4
+
+--S 5 of 153
+n:=2
+--R 
+--R
+--R   (4)  2
+--R                                                        Type: PositiveInteger
+--E 5
+
+--S 6 of 153
+es:=1.0
+--R 
+--R
+--R   (5)  1.0
+--R                                                                  Type: Float
+--E 6
+
+--S 7 of 153
+fu:=0.4373903597E-14
+--R 
+--R
+--R   (6)  0.4373903597 E -14
+--R                                                                  Type: Float
+--E 7
+
+--S 8 of 153
+it:=50
+--R 
+--R
+--R   (7)  50
+--R                                                        Type: PositiveInteger
+--E 8
+
+--S 9 of 153
+lin:=0.9
+--R 
+--R
+--R   (8)  0.9
+--R                                                                  Type: Float
+--E 9
+
+--S 10 of 153
+list:=true
+--R 
+--R
+--R   (9)  true
+--R                                                                Type: Boolean
+--E 10
+
+--S 11 of 153
+ma:=1.0E+20
+--R 
+--R
+--R   (10)  100 0000000000 00000000.0
+--R                                                                  Type: Float
+--E 11
+
+--S 12 of 153
+op:=3.26E-12
+--R 
+--R
+--R   (11)  0.326 E -11
+--R                                                                  Type: Float
+--E 12
+
+--S 13 of 153
+pr:=10
+--R 
+--R
+--R   (12)  10
+--R                                                        Type: PositiveInteger
+--E 13
+
+--S 14 of 153
+sta:=1
+--R 
+--R
+--R   (13)  1
+--R                                                        Type: PositiveInteger
+--E 14
+
+--S 15 of 153
+sto:=2
+--R 
+--R
+--R   (14)  2
+--R                                                        Type: PositiveInteger
+--E 15
+
+--S 16 of 153
+ve:=0
+--R 
+--R
+--R   (15)  0
+--R                                                     Type: NonNegativeInteger
+--E 16
+
+--S 17 of 153
+x:Matrix SF:= [[-1.0 ,1.0 ]]
+--R 
+--R
+--R   (16)  [- 1.  1.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 17
+
+--S 18 of 153
+ef:Expression Float:=_
+  (exp(X[1])*(4*X[1]^2+2*X[2]^2+4*X[1]*X[2]+2*X[2]+1))::EXPR FLOAT
+--R 
+--R
+--R                                                       X
+--R                2                            2          1
+--R   (17)  (2.0 X   + (4.0 X  + 2.0)X  + 4.0 X   + 1.0)%e
+--R               2          1        2        1
+--R                                                      Type: Expression(Float)
+--E 18
+
+--S 19 of 153
+objfun:ASP49(OBJFUN):= retract ef
+--R 
+--R
+--R   (18)  OBJFUN
+--R                                                          Type: Asp49(OBJFUN)
+--E 19
+
+--S 20 of 153
+-- result:=e04dgf(n,es,fu,it,lin,list,ma,op,pr,sta,sto,ve,x,-1,objfun)
+--E 20
+
+)clear all
+ 
+
+--S 21 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 21
+
+--S 22 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 22
+
+--S 23 of 153
+m:=15
+--R 
+--R
+--R   (3)  15
+--R                                                        Type: PositiveInteger
+--E 23
+
+--S 24 of 153
+n:=3
+--R 
+--R
+--R   (4)  3
+--R                                                        Type: PositiveInteger
+--E 24
+
+--S 25 of 153
+liw:=1
+--R 
+--R
+--R   (5)  1
+--R                                                        Type: PositiveInteger
+--E 25
+
+--S 26 of 153
+lw:=171
+--R 
+--R
+--R   (6)  171
+--R                                                        Type: PositiveInteger
+--E 26
+
+--S 27 of 153
+x:Matrix SF:= [[0.5 ,1.0 ,1.5 ]]
+--R 
+--R
+--R   (7)  [0.5  1.  1.5]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 27
+
+--S 28 of 153
+vef:Vector Expression Float:= vector
+ [(XC[3]+15*XC[2])^(-1)+XC[1]-0.14 ,_
+  2*(2*XC[3]+14*XC[2])^(-1)+XC[1]-0.18 ,_
+  3*(3*XC[3]+13*XC[2])^(-1)+XC[1]-0.22 ,_
+  4*(4*XC[3]+12*XC[2])^(-1)+XC[1]-0.25 ,_
+  5*(5*XC[3]+11*XC[2])^(-1)+XC[1]-0.29 ,_
+  6*(6*XC[3]+10*XC[2])^(-1)+XC[1]-0.32 ,_
+  7*(7*XC[3]+9*XC[2])^(-1)+XC[1]-0.35 ,_
+  8*(8*XC[3]+8*XC[2])^(-1)+XC[1]-0.39 ,_
+  9*(7*XC[3]+7*XC[2])^(-1)+XC[1]-0.37 ,_
+  10*(6*XC[3]+6*XC[2])^(-1)+XC[1]-0.58 ,_
+  11*(5*XC[3]+5*XC[2])^(-1)+XC[1]-0.73 ,_
+  12*(4*XC[3]+4*XC[2])^(-1)+XC[1]-0.96 ,_
+  13*(3*XC[3]+3*XC[2])^(-1)+XC[1]-1.34 ,_
+  14*(2*XC[3]+2*XC[2])^(-1)+XC[1]-2.1 ,_
+  15*(XC[3]+XC[2])^(-1)+XC[1]-4.39 ]
+--R 
+--R
+--R   (8)
+--R    (XC  - 0.14)XC  + (15.0 XC  - 2.1)XC  + 1.0
+--R       1          3           1         2
+--R   [-------------------------------------------,
+--R                   XC  + 15.0 XC
+--R                     3          2
+--R    (XC  - 0.18)XC  + (7.0 XC  - 1.26)XC  + 1.0
+--R       1          3          1          2
+--R    -------------------------------------------,
+--R                   XC  + 7.0 XC
+--R                     3         2
+--R
+--R         (XC  - 0.22)XC
+--R            1          3
+--R       + 
+--R         (4.3333333333 333333333 XC  - 0.9533333333 3333333333)XC  + 1.0
+--R                                   1                             2
+--R    /
+--R       XC  + 4.3333333333 333333333 XC
+--R         3                            2
+--R     ,
+--R    (XC  - 0.25)XC  + (3.0 XC  - 0.75)XC  + 1.0
+--R       1          3          1          2
+--R    -------------------------------------------,
+--R                   XC  + 3.0 XC
+--R                     3         2
+--R    (XC  - 0.29)XC  + (2.2 XC  - 0.638)XC  + 1.0
+--R       1          3          1           2
+--R    --------------------------------------------,
+--R                    XC  + 2.2 XC
+--R                      3         2
+--R
+--R         (XC  - 0.32)XC
+--R            1          3
+--R       + 
+--R         (1.6666666666 666666667 XC  - 0.5333333333 3333333333)XC  + 1.0
+--R                                   1                             2
+--R    /
+--R       XC  + 1.6666666666 666666667 XC
+--R         3                            2
+--R     ,
+--R    (XC  - 0.35)XC  + (1.2857142857 142857143 XC  - 0.45)XC  + 1.0
+--R       1          3                             1          2
+--R    --------------------------------------------------------------,
+--R                   XC  + 1.2857142857 142857143 XC
+--R                     3                            2
+--R    (XC  - 0.39)XC  + (XC  - 0.39)XC  + 1.0
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 0.37)XC  + (XC  - 0.37)XC  + 1.2857142857 142857143
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 0.58)XC  + (XC  - 0.58)XC  + 1.6666666666 666666667
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 0.73)XC  + (XC  - 0.73)XC  + 2.2
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 0.96)XC  + (XC  - 0.96)XC  + 3.0
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 1.34)XC  + (XC  - 1.34)XC  + 4.3333333333 333333333
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 2.1)XC  + (XC  - 2.1)XC  + 7.0
+--R       1         3      1         2
+--R    -------------------------------------,
+--R                  XC  + XC
+--R                    3     2
+--R    (XC  - 4.39)XC  + (XC  - 4.39)XC  + 15.0
+--R       1          3      1          2
+--R    ----------------------------------------]
+--R                    XC  + XC
+--R                      3     2
+--R                                              Type: Vector(Expression(Float))
+--E 28
+
+--S 29 of 153
+lsfun1:ASP50(LSFUN1):= retract vef
+--R 
+--R
+--R   (9)  LSFUN1
+--R                                                          Type: Asp50(LSFUN1)
+--E 29
+
+--S 30 of 153
+-- result:=e04fdf(m,n,liw,lw,x,-1,lsfun1)
+--E 30
+
+)clear all
+ 
+
+--S 31 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 31
+
+--S 32 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 32
+
+--S 33 of 153
+m:=15
+--R 
+--R
+--R   (3)  15
+--R                                                        Type: PositiveInteger
+--E 33
+
+--S 34 of 153
+n:=3
+--R 
+--R
+--R   (4)  3
+--R                                                        Type: PositiveInteger
+--E 34
+
+--S 35 of 153
+liw:=1
+--R 
+--R
+--R   (5)  1
+--R                                                        Type: PositiveInteger
+--E 35
+
+--S 36 of 153
+lw:=177
+--R 
+--R
+--R   (6)  177
+--R                                                        Type: PositiveInteger
+--E 36
+
+--S 37 of 153
+x:Matrix SF:= [[0.5 ,1.0 ,1.5 ]]
+--R 
+--R
+--R   (7)  [0.5  1.  1.5]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 37
+
+--S 38 of 153
+vef:Vector Expression Float:=vector
+ [(XC[3]+15*XC[2])^(-1)+XC[1]-0.14 ,_
+   2*(2*XC[3]+14*XC[2])^(-1)+XC[1]-0.18 ,_
+   3*(3*XC[3]+13*XC[2])^(-1)+XC[1]-0.22 ,_
+   4*(4*XC[3]+12*XC[2])^(-1)+XC[1]-0.25 ,_
+   5*(5*XC[3]+11*XC[2])^(-1)+XC[1]-0.29 ,_
+   6*(6*XC[3]+10*XC[2])^(-1)+XC[1]-0.32 ,_
+   7*(7*XC[3]+9*XC[2])^(-1)+XC[1]-0.35 ,_
+   8*(8*XC[3]+8*XC[2])^(-1)+XC[1]-0.39 ,_
+   9*(7*XC[3]+7*XC[2])^(-1)+XC[1]-0.37 ,_
+   10*(6*XC[3]+6*XC[2])^(-1)+XC[1]-0.58 ,_
+   11*(5*XC[3]+5*XC[2])^(-1)+XC[1]-0.73 ,_
+   12*(4*XC[3]+4*XC[2])^(-1)+XC[1]-0.96 ,_
+   13*(3*XC[3]+3*XC[2])^(-1)+XC[1]-1.34 ,_
+   14*(2*XC[3]+2*XC[2])^(-1)+XC[1]-2.1 ,_
+   15*(XC[3]+XC[2])^(-1)+XC[1]-4.39 ]
+--R 
+--R
+--R   (8)
+--R    (XC  - 0.14)XC  + (15.0 XC  - 2.1)XC  + 1.0
+--R       1          3           1         2
+--R   [-------------------------------------------,
+--R                   XC  + 15.0 XC
+--R                     3          2
+--R    (XC  - 0.18)XC  + (7.0 XC  - 1.26)XC  + 1.0
+--R       1          3          1          2
+--R    -------------------------------------------,
+--R                   XC  + 7.0 XC
+--R                     3         2
+--R
+--R         (XC  - 0.22)XC
+--R            1          3
+--R       + 
+--R         (4.3333333333 333333333 XC  - 0.9533333333 3333333333)XC  + 1.0
+--R                                   1                             2
+--R    /
+--R       XC  + 4.3333333333 333333333 XC
+--R         3                            2
+--R     ,
+--R    (XC  - 0.25)XC  + (3.0 XC  - 0.75)XC  + 1.0
+--R       1          3          1          2
+--R    -------------------------------------------,
+--R                   XC  + 3.0 XC
+--R                     3         2
+--R    (XC  - 0.29)XC  + (2.2 XC  - 0.638)XC  + 1.0
+--R       1          3          1           2
+--R    --------------------------------------------,
+--R                    XC  + 2.2 XC
+--R                      3         2
+--R
+--R         (XC  - 0.32)XC
+--R            1          3
+--R       + 
+--R         (1.6666666666 666666667 XC  - 0.5333333333 3333333333)XC  + 1.0
+--R                                   1                             2
+--R    /
+--R       XC  + 1.6666666666 666666667 XC
+--R         3                            2
+--R     ,
+--R    (XC  - 0.35)XC  + (1.2857142857 142857143 XC  - 0.45)XC  + 1.0
+--R       1          3                             1          2
+--R    --------------------------------------------------------------,
+--R                   XC  + 1.2857142857 142857143 XC
+--R                     3                            2
+--R    (XC  - 0.39)XC  + (XC  - 0.39)XC  + 1.0
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 0.37)XC  + (XC  - 0.37)XC  + 1.2857142857 142857143
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 0.58)XC  + (XC  - 0.58)XC  + 1.6666666666 666666667
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 0.73)XC  + (XC  - 0.73)XC  + 2.2
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 0.96)XC  + (XC  - 0.96)XC  + 3.0
+--R       1          3      1          2
+--R    ---------------------------------------,
+--R                   XC  + XC
+--R                     3     2
+--R    (XC  - 1.34)XC  + (XC  - 1.34)XC  + 4.3333333333 333333333
+--R       1          3      1          2
+--R    ----------------------------------------------------------,
+--R                             XC  + XC
+--R                               3     2
+--R    (XC  - 2.1)XC  + (XC  - 2.1)XC  + 7.0
+--R       1         3      1         2
+--R    -------------------------------------,
+--R                  XC  + XC
+--R                    3     2
+--R    (XC  - 4.39)XC  + (XC  - 4.39)XC  + 15.0
+--R       1          3      1          2
+--R    ----------------------------------------]
+--R                    XC  + XC
+--R                      3     2
+--R                                              Type: Vector(Expression(Float))
+--E 38
+
+--S 39 of 153
+lsfun2:Asp19(LSFUN2):= retract vef
+--R 
+--R
+--R   (9)  LSFUN2
+--R                                                          Type: Asp19(LSFUN2)
+--E 39
+
+--S 40 of 153
+-- result:=e04gcf(m,n,liw,lw,x,-1,lsfun2)
+--E 40
+
+)clear all
+ 
+
+--S 41 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 41
+
+--S 42 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 42
+
+--S 43 of 153
+n:=4
+--R 
+--R
+--R   (3)  4
+--R                                                        Type: PositiveInteger
+--E 43
+
+--S 44 of 153
+ibound:=0
+--R 
+--R
+--R   (4)  0
+--R                                                     Type: NonNegativeInteger
+--E 44
+
+--S 45 of 153
+liw:=6
+--R 
+--R
+--R   (5)  6
+--R                                                        Type: PositiveInteger
+--E 45
+
+--S 46 of 153
+lw:=54
+--R 
+--R
+--R   (6)  54
+--R                                                        Type: PositiveInteger
+--E 46
+
+--S 47 of 153
+bl:Matrix SF:= [[1 ,-2 ,-1.0e-6 ,1 ]]
+--R 
+--R
+--R   (7)  [1.  - 2.  - 1.0000000000000002E-6  1.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 47
+
+--S 48 of 153
+bu:Matrix SF:= [[3 ,0 ,1.0e6 ,3 ]]
+--R 
+--R
+--R   (8)  [3.  0.  1000000.  3.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 48
+
+--S 49 of 153
+x:Matrix SF:= [[3 ,-1 ,0 ,1 ]]
+--R 
+--R
+--R   (9)  [3.  - 1.  0.  1.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 49
+
+--S 50 of 153
+ef:Expression Float:=
+ ((XC[1]+10*XC[2])^2+5*(XC[3]-XC[4])^2+(XC[2]-2*XC[3])^4+_
+   10*(XC[1]-XC[4])^4)::EXPR FLOAT
+--R 
+--R
+--R   (10)
+--R             4              3            2          2
+--R     10.0 XC   - 40.0 XC XC   + (60.0 XC   + 5.0)XC
+--R            4           1  4            1          4
+--R   + 
+--R                           3               4              3
+--R     (- 10.0 XC  - 40.0 XC  )XC  + 16.0 XC   - 32.0 XC XC
+--R               3          1    4          3           2  3
+--R   + 
+--R              2          2          3         4            2
+--R     (24.0 XC   + 5.0)XC   - 8.0 XC  XC  + XC   + 100.0 XC   + 20.0 XC XC
+--R             2          3          2   3     2            2           1  2
+--R   + 
+--R             4      2
+--R     10.0 XC   + XC
+--R            1      1
+--R                                                      Type: Expression(Float)
+--E 50
+
+--S 51 of 153
+funct1:Asp24(FUNCT1):=retract ef
+--R 
+--R
+--R   (11)  FUNCT1
+--R                                                          Type: Asp24(FUNCT1)
+--E 51
+
+--S 52 of 153
+-- result:=e04jaf(n,ibound,liw,lw,bl,bu,x,-1,funct1)
+--E 52
+
+)clear all
+ 
+
+--S 53 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 53
+
+--S 54 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 54
+
+--S 55 of 153
+itmax:=20
+--R 
+--R
+--R   (3)  20
+--R                                                        Type: PositiveInteger
+--E 55
+
+--S 56 of 153
+msglvl:=1
+--R 
+--R
+--R   (4)  1
+--R                                                        Type: PositiveInteger
+--E 56
+
+--S 57 of 153
+n:=7
+--R 
+--R
+--R   (5)  7
+--R                                                        Type: PositiveInteger
+--E 57
+
+--S 58 of 153
+nclin:=7
+--R 
+--R
+--R   (6)  7
+--R                                                        Type: PositiveInteger
+--E 58
+
+--S 59 of 153
+nctotl:=14
+--R 
+--R
+--R   (7)  14
+--R                                                        Type: PositiveInteger
+--E 59
+
+--S 60 of 153
+nrowa:=7
+--R 
+--R
+--R   (8)  7
+--R                                                        Type: PositiveInteger
+--E 60
+
+--S 61 of 153
+a:Matrix SF:=
+ [[1 ,1 ,1 ,1 ,1 ,1 ,1 ],_
+  [0.15 ,0.04 ,0.02 ,0.04 ,0.02 ,0.01 ,0.03 ],_
+  [0.03 ,0.05 ,0.08 ,0.02 ,0.06 ,0.01 ,0 ],_
+  [0.02 ,0.04 ,0.01 ,0.02 ,0.02 ,0 ,0 ],_
+  [0.02 ,0.03 ,0 ,0 ,0.01 ,0 ,0 ],_
+  [0.7 ,0.75 ,0.8 ,0.75 ,0.8 ,0.97 ,0 ],_
+  [0.02 ,0.06 ,0.08 ,0.12 ,0.02 ,0.01 ,0.97 ]]
+--R 
+--R
+--R   (9)
+--R   [[1.,1.,1.,1.,1.,1.,1.],
+--R
+--R     [0.14999999999999999, 3.9999999999999994E-2, 1.9999999999999997E-2,
+--R      3.9999999999999994E-2, 1.9999999999999997E-2, 9.9999999999999985E-3,
+--R      2.9999999999999999E-2]
+--R     ,
+--R
+--R     [2.9999999999999999E-2, 5.0000000000000003E-2, 7.9999999999999988E-2,
+--R      1.9999999999999997E-2, 5.9999999999999998E-2, 9.9999999999999985E-3, 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 3.9999999999999994E-2, 9.9999999999999985E-3,
+--R      1.9999999999999997E-2, 1.9999999999999997E-2, 0., 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 2.9999999999999999E-2, 0., 0.,
+--R      9.9999999999999985E-3, 0., 0.]
+--R     ,
+--R
+--R     [0.69999999999999996, 0.75, 0.80000000000000004, 0.75,
+--R      0.80000000000000004, 0.96999999999999997, 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 5.9999999999999998E-2, 7.9999999999999988E-2,
+--R      0.12, 1.9999999999999997E-2, 9.9999999999999985E-3, 0.96999999999999997]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 61
+
+--S 62 of 153
+bl:Matrix SF:=
+ [[-0.01 ,-0.1 ,-0.01 ,-0.04 ,-0.1 ,-0.01 ,-0.01 ,_
+   -0.13 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,_
+   -0.0992 ,-0.003 ]]
+--R 
+--R
+--R   (10)
+--R   [
+--R     [- 9.9999999999999985E-3, - 9.9999999999999992E-2,
+--R      - 9.9999999999999985E-3, - 3.9999999999999994E-2,
+--R      - 9.9999999999999992E-2, - 9.9999999999999985E-3,
+--R      - 9.9999999999999985E-3, - 0.12999999999999998, - 1.0E+21, - 1.0E+21,
+--R      - 1.0E+21, - 1.0E+21, - 9.920000000000001E-2, - 2.9999999999999996E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 62
+
+--S 63 of 153
+bu:Matrix SF:=
+ [[0.01 ,0.15 ,0.03 ,0.02 ,0.05 ,1.0e+21 ,1.0e+21 ,_
+  -0.13 ,-0.0049 ,-0.0064 ,-0.0037 ,-0.0012 ,1.0e+21 ,0.002 ]]
+--R 
+--R
+--R   (11)
+--R   [
+--R     [9.9999999999999985E-3, 0.14999999999999999, 2.9999999999999999E-2,
+--R      1.9999999999999997E-2, 5.0000000000000003E-2, 1.0E+21, 1.0E+21,
+--R      - 0.12999999999999998, - 4.8999999999999998E-3, - 6.3999999999999994E-3,
+--R      - 3.6999999999999997E-3, - 1.2000000000000001E-3, 1.0E+21, 2.0E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 63
+
+--S 64 of 153
+cvec:Matrix SF:= [[-0.02 ,-0.2 ,-0.2 ,-0.2 ,-0.2 ,0.04 ,0.04 ]]
+--R 
+--R
+--R   (12)
+--R   [
+--R     [- 1.9999999999999997E-2, - 0.19999999999999998, - 0.19999999999999998,
+--R      - 0.19999999999999998, - 0.19999999999999998, 3.9999999999999994E-2,
+--R      3.9999999999999994E-2]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 64
+
+--S 65 of 153
+linobj:=true
+--R 
+--R
+--R   (13)  true
+--R                                                                Type: Boolean
+--E 65
+
+--S 66 of 153
+liwork:=14
+--R 
+--R
+--R   (14)  14
+--R                                                        Type: PositiveInteger
+--E 66
+
+--S 67 of 153
+lwork:=182
+--R 
+--R
+--R   (15)  182
+--R                                                        Type: PositiveInteger
+--E 67
+
+--S 68 of 153
+x:Matrix SF:= [[-0.01 ,-0.03 ,0.0 ,-0.01 ,-0.1 ,0.02 ,0.01 ]]
+--R 
+--R
+--R   (16)
+--R   [
+--R     [- 9.9999999999999985E-3, - 2.9999999999999999E-2, 0.,
+--R      - 9.9999999999999985E-3, - 9.9999999999999992E-2, 1.9999999999999997E-2,
+--R      9.9999999999999985E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 68
+
+--S 69 of 153
+-- result:=e04mbf(itmax,msglvl,n,nclin,nctotl,nrowa,a,bl,bu,_
+--                cvec,linobj,liwork,lwork,x,-1)
+--E 69
+
+)clear all
+ 
+
+--S 70 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 70
+
+--S 71 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 71
+
+--S 72 of 153
+itmax:=20
+--R 
+--R
+--R   (3)  20
+--R                                                        Type: PositiveInteger
+--E 72
+
+--S 73 of 153
+msglvl:=-1
+--R 
+--R
+--R   (4)  - 1
+--R                                                                Type: Integer
+--E 73
+
+--S 74 of 153
+n:=7
+--R 
+--R
+--R   (5)  7
+--R                                                        Type: PositiveInteger
+--E 74
+
+--S 75 of 153
+nclin:=7
+--R 
+--R
+--R   (6)  7
+--R                                                        Type: PositiveInteger
+--E 75
+
+--S 76 of 153
+nctotl:=14
+--R 
+--R
+--R   (7)  14
+--R                                                        Type: PositiveInteger
+--E 76
+
+--S 77 of 153
+nrowa:=7
+--R 
+--R
+--R   (8)  7
+--R                                                        Type: PositiveInteger
+--E 77
+
+--S 78 of 153
+nrowh:=7
+--R 
+--R
+--R   (9)  7
+--R                                                        Type: PositiveInteger
+--E 78
+
+--S 79 of 153
+ncolh:=7
+--R 
+--R
+--R   (10)  7
+--R                                                        Type: PositiveInteger
+--E 79
+
+--S 80 of 153
+bigbnd:=1.0e10
+--R 
+--R
+--R   (11)  100 00000000.0
+--R                                                                  Type: Float
+--E 80
+
+--S 81 of 153
+a:Matrix SF:=
+ [[1 ,1 ,1 ,1 ,1 ,1 ,1 ],_
+  [0.15 ,0.04 ,0.02 ,0.04 ,0.02 ,0.01 ,0.03 ],_
+  [0.03 ,0.05 ,0.08 ,0.02 ,0.06 ,0.01 ,0 ],_
+  [0.02 ,0.04 ,0.01 ,0.02 ,0.02 ,0 ,0 ],_
+  [0.02 ,0.03 ,0 ,0 ,0.01 ,0 ,0 ],_
+  [0.7 ,0.75 ,0.8 ,0.75 ,0.8 ,0.97 ,0 ],_
+  [0.02 ,0.06 ,0.08 ,0.12 ,0.02 ,0.01 ,0.97 ]]
+--R 
+--R
+--R   (12)
+--R   [[1.,1.,1.,1.,1.,1.,1.],
+--R
+--R     [0.14999999999999999, 3.9999999999999994E-2, 1.9999999999999997E-2,
+--R      3.9999999999999994E-2, 1.9999999999999997E-2, 9.9999999999999985E-3,
+--R      2.9999999999999999E-2]
+--R     ,
+--R
+--R     [2.9999999999999999E-2, 5.0000000000000003E-2, 7.9999999999999988E-2,
+--R      1.9999999999999997E-2, 5.9999999999999998E-2, 9.9999999999999985E-3, 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 3.9999999999999994E-2, 9.9999999999999985E-3,
+--R      1.9999999999999997E-2, 1.9999999999999997E-2, 0., 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 2.9999999999999999E-2, 0., 0.,
+--R      9.9999999999999985E-3, 0., 0.]
+--R     ,
+--R
+--R     [0.69999999999999996, 0.75, 0.80000000000000004, 0.75,
+--R      0.80000000000000004, 0.96999999999999997, 0.]
+--R     ,
+--R
+--R     [1.9999999999999997E-2, 5.9999999999999998E-2, 7.9999999999999988E-2,
+--R      0.12, 1.9999999999999997E-2, 9.9999999999999985E-3, 0.96999999999999997]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 81
+
+--S 82 of 153
+bl:Matrix SF:=
+ [[-0.01 ,-0.1 ,-0.01 ,-0.04 ,-0.1 ,-0.01 ,-0.01 ,-0.13 ,_
+   -1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-0.0992 ,-0.003 ]]
+--R 
+--R
+--R   (13)
+--R   [
+--R     [- 9.9999999999999985E-3, - 9.9999999999999992E-2,
+--R      - 9.9999999999999985E-3, - 3.9999999999999994E-2,
+--R      - 9.9999999999999992E-2, - 9.9999999999999985E-3,
+--R      - 9.9999999999999985E-3, - 0.12999999999999998, - 1.0E+21, - 1.0E+21,
+--R      - 1.0E+21, - 1.0E+21, - 9.920000000000001E-2, - 2.9999999999999996E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 82
+
+--S 83 of 153
+bu:Matrix SF:=
+ [[0.01 ,0.15 ,0.03 ,0.02 ,0.05 ,1.0e+21 ,1.0e+21 ,-0.13 ,_
+  -0.0049 ,-0.0064 ,-0.0037 ,-0.0012 ,1.0e+21 ,0.002 ]]
+--R 
+--R
+--R   (14)
+--R   [
+--R     [9.9999999999999985E-3, 0.14999999999999999, 2.9999999999999999E-2,
+--R      1.9999999999999997E-2, 5.0000000000000003E-2, 1.0E+21, 1.0E+21,
+--R      - 0.12999999999999998, - 4.8999999999999998E-3, - 6.3999999999999994E-3,
+--R      - 3.6999999999999997E-3, - 1.2000000000000001E-3, 1.0E+21, 2.0E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 83
+
+--S 84 of 153
+cvec:Matrix SF:= [[-0.02 ,-0.2 ,-0.2 ,-0.2 ,-0.2 ,0.04 ,0.04 ]]
+--R 
+--R
+--R   (15)
+--R   [
+--R     [- 1.9999999999999997E-2, - 0.19999999999999998, - 0.19999999999999998,
+--R      - 0.19999999999999998, - 0.19999999999999998, 3.9999999999999994E-2,
+--R      3.9999999999999994E-2]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 84
+
+--S 85 of 153
+featol:Matrix SF:=
+ [[0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,_
+   0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,_
+   0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ]]
+--R 
+--R
+--R   (16)
+--R   [
+--R     [1.0529999999999999E-8, 1.0529999999999999E-8, 1.0529999999999999E-8,
+--R      1.0529999999999999E-8, 1.0529999999999999E-8, 1.0529999999999999E-8,
+--R      1.0529999999999999E-8, 1.0529999999999999E-8, 1.0529999999999999E-8,
+--R      1.0529999999999999E-8, 1.0529999999999999E-8, 1.0529999999999999E-8,
+--R      1.0529999999999999E-8, 1.0529999999999999E-8]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 85
+
+--S 86 of 153
+hess:Matrix SF :=
+ [[2 ,0 ,0 ,0 ,0 ,0 ,0 ],_
+  [0 ,2 ,0 ,0 ,0 ,0 ,0 ],_
+  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
+  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
+  [0 ,0 ,0 ,0 ,2 ,0 ,0 ],_
+  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ],_
+  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ]]
+--R 
+--R
+--R         +2.  0.  0.  0.  0.   0.    0. +
+--R         |                              |
+--R         |0.  2.  0.  0.  0.   0.    0. |
+--R         |                              |
+--R         |0.  0.  2.  2.  0.   0.    0. |
+--R         |                              |
+--R   (17)  |0.  0.  2.  2.  0.   0.    0. |
+--R         |                              |
+--R         |0.  0.  0.  0.  2.   0.    0. |
+--R         |                              |
+--R         |0.  0.  0.  0.  0.  - 2.  - 2.|
+--R         |                              |
+--R         +0.  0.  0.  0.  0.  - 2.  - 2.+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 86
+
+--S 87 of 153
+cold:=true
+--R 
+--R
+--R   (18)  true
+--R                                                                Type: Boolean
+--E 87
+
+--S 88 of 153
+lpp:=false
+--R 
+--R
+--R   (19)  false
+--R                                                                Type: Boolean
+--E 88
+
+--S 89 of 153
+orthog:=true
+--R 
+--R
+--R   (20)  true
+--R                                                                Type: Boolean
+--E 89
+
+--S 90 of 153
+liwork:=14
+--R 
+--R
+--R   (21)  14
+--R                                                        Type: PositiveInteger
+--E 90
+
+--S 91 of 153
+lwork:=238
+--R 
+--R
+--R   (22)  238
+--R                                                        Type: PositiveInteger
+--E 91
+
+--S 92 of 153
+x:Matrix SF:= [[-0.01 ,-0.03 ,0.0 ,-0.01 ,-0.1 ,0.02 ,0.01 ]]
+--R 
+--R
+--R   (23)
+--R   [
+--R     [- 9.9999999999999985E-3, - 2.9999999999999999E-2, 0.,
+--R      - 9.9999999999999985E-3, - 9.9999999999999992E-2, 1.9999999999999997E-2,
+--R      9.9999999999999985E-3]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 92
+
+--S 93 of 153
+istate:Matrix Integer:= [[0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
+--R 
+--R
+--R   (24)  [0  0  0  0  0  0  0  0  0  0  0  0  0  0]
+--R                                                        Type: Matrix(Integer)
+--E 93
+
+--S 94 of 153
+mei:Matrix Expression Integer:=
+ [[2 ,0 ,0 ,0 ,0 ,0 ,0 ],_
+  [0 ,2 ,0 ,0 ,0 ,0 ,0 ],_
+  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
+  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
+  [0 ,0 ,0 ,0 ,2 ,0 ,0 ],_
+  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ],_
+  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ]]
+--R 
+--R
+--R         +2  0  0  0  0   0    0 +
+--R         |                       |
+--R         |0  2  0  0  0   0    0 |
+--R         |                       |
+--R         |0  0  2  2  0   0    0 |
+--R         |                       |
+--R   (25)  |0  0  2  2  0   0    0 |
+--R         |                       |
+--R         |0  0  0  0  2   0    0 |
+--R         |                       |
+--R         |0  0  0  0  0  - 2  - 2|
+--R         |                       |
+--R         +0  0  0  0  0  - 2  - 2+
+--R                                            Type: Matrix(Expression(Integer))
+--E 94
+
+--S 95 of 153
+qphess:Asp20(QPHESS):= retract mei
+--R 
+--R
+--R   (26)  QPHESS
+--R                                                          Type: Asp20(QPHESS)
+--E 95
+
+--S 96 of 153
+-- result:=e04naf(itmax,msglvl,n,nclin,nctotl,nrowa,nrowh,ncolh,_
+--               bigbnd,a,bl,bu,cvec,featol,_
+--               hess,cold,lpp,orthog,liwork,lwork,x,istate,-1,qphess)
+--E 96
+
+)clear all
+ 
+
+--S 97 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 97
+
+--S 98 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 98
+
+--S 99 of 153
+n:=9
+--R 
+--R
+--R   (3)  9
+--R                                                        Type: PositiveInteger
+--E 99
+
+--S 100 of 153
+nclin:=4
+--R 
+--R
+--R   (4)  4
+--R                                                        Type: PositiveInteger
+--E 100
+
+--S 101 of 153
+ncnln:=14
+--R 
+--R
+--R   (5)  14
+--R                                                        Type: PositiveInteger
+--E 101
+
+--S 102 of 153
+nrowa:=4
+--R 
+--R
+--R   (6)  4
+--R                                                        Type: PositiveInteger
+--E 102
+
+--S 103 of 153
+nrowj:=14
+--R 
+--R
+--R   (7)  14
+--R                                                        Type: PositiveInteger
+--E 103
+
+--S 104 of 153
+nrowr:=9
+--R 
+--R
+--R   (8)  9
+--R                                                        Type: PositiveInteger
+--E 104
+
+--S 105 of 153
+a:Matrix SF:=
+ [[-1.0 ,1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
+  [0.0 ,-1.0 ,1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
+  [0.0 ,0.0 ,1.0 ,-1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
+  [0.0 ,0.0 ,0.0 ,1.0 ,-1.0 ,0.0 ,0.0 ,0.0 ,0.0 ]]
+--R 
+--R
+--R        +- 1.   1.   0.   0.    0.   0.  0.  0.  0.+
+--R        |                                          |
+--R        | 0.   - 1.  1.   0.    0.   0.  0.  0.  0.|
+--R   (9)  |                                          |
+--R        | 0.    0.   1.  - 1.   0.   0.  0.  0.  0.|
+--R        |                                          |
+--R        + 0.    0.   0.   1.   - 1.  0.  0.  0.  0.+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 105
+
+--S 106 of 153
+bl:Matrix SF:=
+ [[0.0 ,-1.0e15 ,-1.0 ,-1.0e15 ,0.0 ,0.0 ,0.0 ,-1.0e15 ,_
+  -1.0e15 ,0.0 ,0.0 ,0.0 ,0.0 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,_
+  -1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,_
+  -1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ]]
+--R 
+--R
+--R   (10)
+--R   [
+--R     [0., - 1.0E+15, - 1., - 1.0E+15, 0., 0., 0., - 1.0E+15, - 1.0E+15, 0., 0.,
+--R      0., 0., - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15,
+--R      - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15, - 1.0E+15,
+--R      - 1.0E+15, - 1.0E+15]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 106
+
+--S 107 of 153
+bu:Matrix SF:=
+ [[1.0e15 ,1.0e15 ,1.0 ,1.0e15 ,1.0e15 ,1.0e15 ,1.0e15 ,_
+  0.0 ,0.0 ,1.0e15 ,1.0e15 ,1.0e15 ,1.0e15 ,0.0 ,0.0 ,_
+  0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ]]
+--R 
+--R
+--R   (11)
+--R   [
+--R     [1.0E+15, 1.0E+15, 1., 1.0E+15, 1.0E+15, 1.0E+15, 1.0E+15, 0., 0.,
+--R      1.0E+15, 1.0E+15, 1.0E+15, 1.0E+15, 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--R      0., 0., 0., 0., 0.]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 107
+
+--S 108 of 153
+liwork:=59
+--R 
+--R
+--R   (12)  59
+--R                                                        Type: PositiveInteger
+--E 108
+
+--S 109 of 153
+lwork:=968
+--R 
+--R
+--R   (13)  968
+--R                                                        Type: PositiveInteger
+--E 109
+
+--S 110 of 153
+sta:=false
+--R 
+--R
+--R   (14)  false
+--R                                                                Type: Boolean
+--E 110
+
+--S 111 of 153
+cra:= 0.01
+--R 
+--R
+--R   (15)  0.01
+--R                                                                  Type: Float
+--E 111
+
+--S 112 of 153
+der:=3
+--R 
+--R
+--R   (16)  3
+--R                                                        Type: PositiveInteger
+--E 112
+
+--S 113 of 153
+fea:=0.1053671201E-7
+--R 
+--R
+--R   (17)  0.1053671201 E -7
+--R                                                                  Type: Float
+--E 113
+
+--S 114 of 153
+fun:=0.4373903510E-14
+--R 
+--R
+--R   (18)  0.437390351 E -14
+--R                                                                  Type: Float
+--E 114
+
+--S 115 of 153
+hes:= true
+--R 
+--R
+--R   (19)  true
+--R                                                                Type: Boolean
+--E 115
+
+--S 116 of 153
+infb:= 1.00E+15
+--R 
+--R
+--R   (20)  10000000 00000000.0
+--R                                                                  Type: Float
+--E 116
+
+--S 117 of 153
+infs:= 1.00E+15
+--R 
+--R
+--R   (21)  10000000 00000000.0
+--R                                                                  Type: Float
+--E 117
+
+--S 118 of 153
+linf:= 0.1053671201E-7
+--R 
+--R
+--R   (22)  0.1053671201 E -7
+--R                                                                  Type: Float
+--E 118
+
+--S 119 of 153
+lint:= 0.9
+--R 
+--R
+--R   (23)  0.9
+--R                                                                  Type: Float
+--E 119
+
+--S 120 of 153
+list:= true
+--R 
+--R
+--R   (24)  true
+--R                                                                Type: Boolean
+--E 120
+
+--S 121 of 153
+maji:= 30
+--R 
+--R
+--R   (25)  30
+--R                                                        Type: PositiveInteger
+--E 121
+
+--S 122 of 153
+majp:= 1
+--R 
+--R
+--R   (26)  1
+--R                                                        Type: PositiveInteger
+--E 122
+
+--S 123 of 153
+mini:= 81
+--R 
+--R
+--R   (27)  81
+--R                                                        Type: PositiveInteger
+--E 123
+
+--S 124 of 153
+minp:= 0
+--R 
+--R
+--R   (28)  0
+--R                                                     Type: NonNegativeInteger
+--E 124
+
+--S 125 of 153
+mon:= -1
+--R 
+--R
+--R   (29)  - 1
+--R                                                                Type: Integer
+--E 125
+
+--S 126 of 153
+nonf:= 1.05E-08
+--R 
+--R
+--R   (30)  0.105 E -7
+--R                                                                  Type: Float
+--E 126
+
+--S 127 of 153
+opt:= 3.26E-08
+--R 
+--R
+--R   (31)  0.326 E -7
+--R                                                                  Type: Float
+--E 127
+
+--S 128 of 153
+ste:= 2.0
+--R 
+--R
+--R   (32)  2.0
+--R                                                                  Type: Float
+--E 128
+
+--S 129 of 153
+stao:= 1
+--R 
+--R
+--R   (33)  1
+--R                                                        Type: PositiveInteger
+--E 129
+
+--S 130 of 153
+stac:= 1
+--R 
+--R
+--R   (34)  1
+--R                                                        Type: PositiveInteger
+--E 130
+
+--S 131 of 153
+stoo:= 9
+--R 
+--R
+--R   (35)  9
+--R                                                        Type: PositiveInteger
+--E 131
+
+--S 132 of 153
+stoc:= 9
+--R 
+--R
+--R   (36)  9
+--R                                                        Type: PositiveInteger
+--E 132
+
+--S 133 of 153
+ve:= 3
+--R 
+--R
+--R   (37)  3
+--R                                                        Type: PositiveInteger
+--E 133
+
+--S 134 of 153
+istate:Matrix Integer:=new(1,27,0)
+--R 
+--R
+--R   (38)
+--R   [[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]]
+--R                                                        Type: Matrix(Integer)
+--E 134
+
+--S 135 of 153
+cjac:Matrix SF:=new(14,9,0.0)
+--R 
+--R
+--R         +0.  0.  0.  0.  0.  0.  0.  0.  0.+
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R   (39)  |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         +0.  0.  0.  0.  0.  0.  0.  0.  0.+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 135
+
+--S 136 of 153
+clamda:Matrix SF:=new(1,27,0.0)
+--R 
+--R
+--R   (40)
+--R   [
+--R     [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--R      0., 0., 0., 0., 0., 0., 0., 0., 0.]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 136
+
+--S 137 of 153
+r:Matrix SF:=new(9,9,0.0)
+--R 
+--R
+--R         +0.  0.  0.  0.  0.  0.  0.  0.  0.+
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R   (41)  |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         |0.  0.  0.  0.  0.  0.  0.  0.  0.|
+--R         |                                  |
+--R         +0.  0.  0.  0.  0.  0.  0.  0.  0.+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 137
+
+--S 138 of 153
+x:Matrix SF:=
+ [[0.1 ,0.125 ,0.666666 ,0.142857 ,0.111111 ,0.2 ,0.25 ,-0.2 ,-0.25 ]]
+--R 
+--R
+--R   (42)
+--R   [
+--R     [0.10000000000000001, 0.125, 0.66666599999999998, 0.14285700000000001,
+--R      0.11111099999999999, 0.20000000000000001, 0.25, - 0.19999999999999998,
+--R      - 0.25]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 138
+
+--S 139 of 153
+vef:Vector Expression Float:=vector
+ [X[1]^2 + X[6]^2 ,_
+ (X[2] - X[1])^2 + (X[7] - X[6])^2 ,_
+ (X[3] - X[1])^2 + X[6]*2 ,_
+ (X[1] - X[4])^2 + (X[6] - X[8])^2 ,_
+ (X[1] - X[5])^2 + (X[6] - X[9])^2 ,_
+ X[2]^2 + X[7]^2 ,_
+ (X[3] - X[2])^2 + X[7]^2 ,_
+ (X[4] - X[2])^2 + (X[8] - X[7])^2 ,_
+ (X[2] - X[5])^2 + (X[7] - X[9])^2 ,_
+ (X[4] - X[3])^2 + X[8]^2 ,_
+ (X[5] - X[3])^2 + X[9]^2 ,_
+ X[4]^2 + X[8]^2 ,_
+ (X[4] - X[5])^2 + (X[9] - X[8])^2 ,_
+ X[5]^2 + X[9]^2 ]
+--R 
+--R
+--R   (43)
+--R      2     2    2                2     2                2
+--R   [X   + X  , X   - 2.0 X X  + X   + X   - 2.0 X X  + X  ,
+--R     6     1    7         6 7    6     2         1 2    1
+--R               2                2    2                2     2                2
+--R    2.0 X  + X   - 2.0 X X  + X  , X   - 2.0 X X  + X   + X   - 2.0 X X  + X  ,
+--R         6    3         1 3    1    8         6 8    6     4         1 4    1
+--R      2                2     2                2    2     2
+--R    X   - 2.0 X X  + X   + X   - 2.0 X X  + X  , X   + X  ,
+--R     9         6 9    6     5         1 5    1    7     2
+--R      2     2                2    2                2     2                2
+--R    X   + X   - 2.0 X X  + X  , X   - 2.0 X X  + X   + X   - 2.0 X X  + X  ,
+--R     7     3         2 3    2    8         7 8    7     4         2 4    2
+--R      2                2     2                2    2     2                2
+--R    X   - 2.0 X X  + X   + X   - 2.0 X X  + X  , X   + X   - 2.0 X X  + X  ,
+--R     9         7 9    7     5         2 5    2    8     4         3 4    3
+--R      2     2                2    2     2
+--R    X   + X   - 2.0 X X  + X  , X   + X  ,
+--R     9     5         3 5    3    8     4
+--R      2                2     2                2    2     2
+--R    X   - 2.0 X X  + X   + X   - 2.0 X X  + X  , X   + X  ]
+--R     9         8 9    8     5         4 5    4    9     5
+--R                                              Type: Vector(Expression(Float))
+--E 139
+
+--S 140 of 153
+confun:Asp55(CONFUN):= retract vef
+--R 
+--R
+--R   (44)  CONFUN
+--R                                                          Type: Asp55(CONFUN)
+--E 140
+
+--S 141 of 153
+ef:Expression Float:=(-X[2]*X[6] + X[1]*X[7] - X[3]*X[7] - _
+                      X[5]*X[8] + X[4]*X[9] + X[3]*X[8])::EXPR FLOAT
+--R 
+--R
+--R   (45)  X X  + (- 1.0 X  + X )X  + (- 1.0 X  + X )X  - 1.0 X X
+--R          4 9           5    3  8           3    1  7        2 6
+--R                                                      Type: Expression(Float)
+--E 141
+
+--S 142 of 153
+objfun:Asp49(OBJFUN):=retract ef
+--R 
+--R
+--R   (46)  OBJFUN
+--R                                                          Type: Asp49(OBJFUN)
+--E 142
+
+--S 143 of 153
+-- result:=e04ucf(n,nclin,ncnln,nrowa,nrowj,nrowr,a,bl,bu,liwork,lwork,_
+-- sta,cra,der,fea,fun,hes,infb,infs,linf,lint,list,maji,majp,_
+-- mini,minp,mon,nonf,opt,ste,stao,stac,stoo,stoc,ve,istate,_
+-- cjac,clamda,r,x,-1,confun,objfun)
+--E 143
+
+)clear all
+ 
+
+--S 144 of 153
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 144
+
+--S 145 of 153
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 145
+
+--S 146 of 153
+job:=0
+--R 
+--R
+--R   (3)  0
+--R                                                     Type: NonNegativeInteger
+--E 146
+
+--S 147 of 153
+m:=15
+--R 
+--R
+--R   (4)  15
+--R                                                        Type: PositiveInteger
+--E 147
+
+--S 148 of 153
+n:=3
+--R 
+--R
+--R   (5)  3
+--R                                                        Type: PositiveInteger
+--E 148
+
+--S 149 of 153
+fsumsq:=0.0082148773065789729
+--R 
+--R
+--R   (6)  0.0082148773 065789729
+--R                                                                  Type: Float
+--E 149
+
+--S 150 of 153
+s:Matrix SF:= [[4.0965034571419325 ,1.5949579400198182 ,0.061258491120317927 ]]
+--R 
+--R
+--R   (7)  [4.0965034571419316  1.5949579400198182  6.1258491120317921E-2]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 150
+
+--S 151 of 153
+lv:=3
+--R 
+--R
+--R   (8)  3
+--R                                                        Type: PositiveInteger
+--E 151
+
+--S 152 of 153
+v:Matrix SF:=
+ [[0.9353959087369782,- 0.25922842715810668, - 0.24048932714625143,_
+   0.35295122094583287, 0.64323460044494585, 0.67946647041872243,_
+   - 0.021445968175122822, - 0.72045115819433947, 0.69317400348358882]]
+--R 
+--R
+--R   (9)
+--R   [
+--R     [0.9353959087369782, - 0.25922842715810668, - 0.2404893271462514,
+--R      0.35295122094583287, 0.64323460044494585, 0.67946647041872232,
+--R      - 2.1445968175122822E-2, - 0.72045115819433958, 0.69317400348358882]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 152
+
+--S 153 of 153
+-- result:=e04ycf(job,m,n,fsumsq,s,lv,v,-1)
+--E 153
+
+)spool
+ 
 )lisp (bye)
 \end{chunk}
 
diff --git a/changelog b/changelog
index 70cc308..c02dc0d 100644
--- a/changelog
+++ b/changelog
@@ -1,3 +1,9 @@
+20141202 tpd src/axiom-website/patches.html 20141202.06.tpd.patch
+20141202 tpd src/input/Makefile remove NAGE04 src/input files from build
+20141202 tpd books/bookvol10.4 NAGE04 absorb src/input test case files
+20141202 tpd src/axiom-website/patches.html 20141202.05.tpd.patch
+20141202 tpd src/input/Makefile remove NAGE02 src/input files from build
+20141202 tpd books/bookvol10.4 NAGE02 absorb src/input test case files
 20141202 tpd src/axiom-website/patches.html 20141202.04.tpd.patch
 20141202 tpd src/input/Makefile remove NAGE01 src/input files from build
 20141202 tpd src/axiom-website/patches.html 20141202.03.tpd.patch
diff --git a/patch b/patch
index 9e24f40..5660b1b 100644
--- a/patch
+++ b/patch
@@ -1 +1,3 @@
-src/input/Makefile remove NAGE01 src/input files from build process
+books/bookvol10.4 NAGE04 absorb src/input test case files
+
+integrate the files into the packages book
diff --git a/src/input/Makefile.pamphlet b/src/input/Makefile.pamphlet
index 5ee4859..6babac1 100644
--- a/src/input/Makefile.pamphlet
+++ b/src/input/Makefile.pamphlet
@@ -206,13 +206,7 @@ SKIP= as-eg1.output as-eg2.output    \
       draw2dsf.output  drawalg.output   drawcfn.output   drawcfun.output  \
       drawcurv.output  drawcx.output    drawex.output    draw.output      \
       drawpoly.output  drawx.output     dropt.output     \
-      e02adf.output    e02aef.output    e02agf.output    \
-      e02ahf.output    e02ajf.output    e02akf.output    e02baf.output    \
-      e02bbf.output    e02bcf.output    e02bdf.output    e02bef.output    \
-      e02daf.output    e02dcf.output    e02ddf.output    e02def.output    \
-      e02dff.output    e02gaf.output    e02zaf.output    e04dgf.output    \
-      e04fdf.output    e04gcf.output    e04jaf.output    e04mbf.output    \
-      e04naf.output    e04ucf.output    e04ycf.output    egg.output       \
+      egg.output       \
       explot2d.output  explot3d.output  \
       f01brf.output    f01bsf.output    f01maf.output    \
       f01mcf.output    f01qcf.output    f01qdf.output    f01qef.output    \
@@ -656,9 +650,6 @@ NAGLIB=${OUT}/c02aff.input   ${OUT}/c02agf.input     ${OUT}/c05adf.input \
        ${OUT}/d02ejf.input   ${OUT}/d02gaf.input     ${OUT}/d02gbf.input \
        ${OUT}/d02kef.input   ${OUT}/d02raf.input     ${OUT}/d03edf.input \
        ${OUT}/d03eef.input   ${OUT}/d03faf.input     \
-       ${OUT}/e04dgf.input     ${OUT}/e04fdf.input \
-       ${OUT}/e04gcf.input   ${OUT}/e04jaf.input     ${OUT}/e04mbf.input \
-       ${OUT}/e04naf.input   ${OUT}/e04ucf.input     ${OUT}/e04ycf.input \
        ${OUT}/f01brf.input   ${OUT}/f01bsf.input     ${OUT}/f01maf.input \
        ${OUT}/f01mcf.input   ${OUT}/f01qcf.input     ${OUT}/f01qdf.input \
        ${OUT}/f01qef.input   ${OUT}/f01rcf.input     ${OUT}/f01rdf.input \
@@ -1097,12 +1088,7 @@ DOCFILES= \
   ${DOC}/drawcurv.input.dvi    ${DOC}/drawcx.input.dvi     \
   ${DOC}/drawex.input.dvi      ${DOC}/draw.input.dvi       \
   ${DOC}/drawpoly.input.dvi    ${DOC}/drawx.input.dvi      \
-  ${DOC}/dropt.input.dvi       \
-  ${DOC}/e04dgf.input.dvi     \
-  ${DOC}/e04fdf.input.dvi      ${DOC}/e04gcf.input.dvi     \
-  ${DOC}/e04jaf.input.dvi      ${DOC}/e04mbf.input.dvi     \
-  ${DOC}/e04naf.input.dvi      ${DOC}/e04ucf.input.dvi     \
-  ${DOC}/e04ycf.input.dvi      ${DOC}/e1.input.dvi         \
+  ${DOC}/dropt.input.dvi       ${DOC}/e1.input.dvi         \
   ${DOC}/ei.input.dvi          ${DOC}/easter.input.dvi     \
   ${DOC}/ecfact.as.dvi         ${DOC}/efi.input.dvi        \
   ${DOC}/egg.input.dvi         ${DOC}/eigen.input.dvi      \
diff --git a/src/input/e04dgf.input.pamphlet b/src/input/e04dgf.input.pamphlet
deleted file mode 100644
index aa109f3..0000000
--- a/src/input/e04dgf.input.pamphlet
+++ /dev/null
@@ -1,46 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04dgf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-n:=2
-es:=1.0
-fu:=0.4373903597E-14
-it:=50
-lin:=0.9
-list:=true
-ma:=1.0E+20
-op:=3.26E-12
-pr:=10
-sta:=1
-sto:=2
-ve:=0
-x:Matrix SF:=
- [[-1.0 ,1.0 ]]
-ef:Expression Float:=_
-  (exp(X[1])*(4*X[1]^2+2*X[2]^2+4*X[1]*X[2]+2*X[2]+1))::EXPR FLOAT
-objfun:ASP49(OBJFUN):= retract ef
-result:=e04dgf(n,es,fu,it,lin,list,ma,op,pr,sta,sto,ve,x,-1,objfun)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04fdf.input.pamphlet b/src/input/e04fdf.input.pamphlet
deleted file mode 100644
index 2baa82f..0000000
--- a/src/input/e04fdf.input.pamphlet
+++ /dev/null
@@ -1,52 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04fdf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-m:=15
-n:=3
-liw:=1
-lw:=171
-x:Matrix SF:=
- [[0.5 ,1.0 ,1.5 ]]
-vef:Vector Expression Float:= vector
- [(XC[3]+15*XC[2])^(-1)+XC[1]-0.14 ,_
-  2*(2*XC[3]+14*XC[2])^(-1)+XC[1]-0.18 ,_
-  3*(3*XC[3]+13*XC[2])^(-1)+XC[1]-0.22 ,_
-  4*(4*XC[3]+12*XC[2])^(-1)+XC[1]-0.25 ,_
-  5*(5*XC[3]+11*XC[2])^(-1)+XC[1]-0.29 ,_
-  6*(6*XC[3]+10*XC[2])^(-1)+XC[1]-0.32 ,_
-  7*(7*XC[3]+9*XC[2])^(-1)+XC[1]-0.35 ,_
-  8*(8*XC[3]+8*XC[2])^(-1)+XC[1]-0.39 ,_
-  9*(7*XC[3]+7*XC[2])^(-1)+XC[1]-0.37 ,_
-  10*(6*XC[3]+6*XC[2])^(-1)+XC[1]-0.58 ,_
-  11*(5*XC[3]+5*XC[2])^(-1)+XC[1]-0.73 ,_
-  12*(4*XC[3]+4*XC[2])^(-1)+XC[1]-0.96 ,_
-  13*(3*XC[3]+3*XC[2])^(-1)+XC[1]-1.34 ,_
-  14*(2*XC[3]+2*XC[2])^(-1)+XC[1]-2.1 ,_
-  15*(XC[3]+XC[2])^(-1)+XC[1]-4.39 ]
-lsfun1:ASP50(LSFUN1):= retract vef
-result:=e04fdf(m,n,liw,lw,x,-1,lsfun1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04gcf.input.pamphlet b/src/input/e04gcf.input.pamphlet
deleted file mode 100644
index 18565d9..0000000
--- a/src/input/e04gcf.input.pamphlet
+++ /dev/null
@@ -1,52 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04gcf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-m:=15
-n:=3
-liw:=1
-lw:=177
-x:Matrix SF:=
- [[0.5 ,1.0 ,1.5 ]]
-vef:Vector Expression Float:=vector
- [(XC[3]+15*XC[2])^(-1)+XC[1]-0.14 ,_
-   2*(2*XC[3]+14*XC[2])^(-1)+XC[1]-0.18 ,_
-   3*(3*XC[3]+13*XC[2])^(-1)+XC[1]-0.22 ,_
-   4*(4*XC[3]+12*XC[2])^(-1)+XC[1]-0.25 ,_
-   5*(5*XC[3]+11*XC[2])^(-1)+XC[1]-0.29 ,_
-   6*(6*XC[3]+10*XC[2])^(-1)+XC[1]-0.32 ,_
-   7*(7*XC[3]+9*XC[2])^(-1)+XC[1]-0.35 ,_
-   8*(8*XC[3]+8*XC[2])^(-1)+XC[1]-0.39 ,_
-   9*(7*XC[3]+7*XC[2])^(-1)+XC[1]-0.37 ,_
-   10*(6*XC[3]+6*XC[2])^(-1)+XC[1]-0.58 ,_
-   11*(5*XC[3]+5*XC[2])^(-1)+XC[1]-0.73 ,_
-   12*(4*XC[3]+4*XC[2])^(-1)+XC[1]-0.96 ,_
-   13*(3*XC[3]+3*XC[2])^(-1)+XC[1]-1.34 ,_
-   14*(2*XC[3]+2*XC[2])^(-1)+XC[1]-2.1 ,_
-   15*(XC[3]+XC[2])^(-1)+XC[1]-4.39 ]
-lsfun2:Asp19(LSFUN2):= retract vef
-result:=e04gcf(m,n,liw,lw,x,-1,lsfun2)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04jaf.input.pamphlet b/src/input/e04jaf.input.pamphlet
deleted file mode 100644
index 2479974..0000000
--- a/src/input/e04jaf.input.pamphlet
+++ /dev/null
@@ -1,43 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04jaf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-n:=4
-ibound:=0
-liw:=6
-lw:=54
-bl:Matrix SF:=
- [[1 ,-2 ,-1.0e-6 ,1 ]]
-bu:Matrix SF:=
- [[3 ,0 ,1.0e6 ,3 ]]
-x:Matrix SF:=
- [[3 ,-1 ,0 ,1 ]]
-ef:Expression Float:=
- ((XC[1]+10*XC[2])^2+5*(XC[3]-XC[4])^2+(XC[2]-2*XC[3])^4+_
-   10*(XC[1]-XC[4])^4)::EXPR FLOAT
-funct1:Asp24(FUNCT1):=retract ef
-result:=e04jaf(n,ibound,liw,lw,bl,bu,x,-1,funct1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04mbf.input.pamphlet b/src/input/e04mbf.input.pamphlet
deleted file mode 100644
index 87690ee..0000000
--- a/src/input/e04mbf.input.pamphlet
+++ /dev/null
@@ -1,58 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04mbf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-itmax:=20
-msglvl:=1
-n:=7
-nclin:=7
-nctotl:=14
-nrowa:=7
-a:Matrix SF:=
- [[1 ,1 ,1 ,1 ,1 ,1 ,1 ],_
-  [0.15 ,0.04 ,0.02 ,0.04 ,0.02 ,0.01 ,0.03 ],_
-  [0.03 ,0.05 ,0.08 ,0.02 ,0.06 ,0.01 ,0 ],_
-  [0.02 ,0.04 ,0.01 ,0.02 ,0.02 ,0 ,0 ],_
-  [0.02 ,0.03 ,0 ,0 ,0.01 ,0 ,0 ],_
-  [0.7 ,0.75 ,0.8 ,0.75 ,0.8 ,0.97 ,0 ],_
-  [0.02 ,0.06 ,0.08 ,0.12 ,0.02 ,0.01 ,0.97 ]]
-bl:Matrix SF:=
- [[-0.01 ,-0.1 ,-0.01 ,-0.04 ,-0.1 ,-0.01 ,-0.01 ,_
-   -0.13 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,_
-   -0.0992 ,-0.003 ]]
-bu:Matrix SF:=
- [[0.01 ,0.15 ,0.03 ,0.02 ,0.05 ,1.0e+21 ,1.0e+21 ,_
-  -0.13 ,-0.0049 ,-0.0064 ,-0.0037 ,-0.0012 ,1.0e+21 ,0.002 ]]
-cvec:Matrix SF:=
- [[-0.02 ,-0.2 ,-0.2 ,-0.2 ,-0.2 ,0.04 ,0.04 ]]
-linobj:=true
-liwork:=14
-lwork:=182
-x:Matrix SF:=
- [[-0.01 ,-0.03 ,0.0 ,-0.01 ,-0.1 ,0.02 ,0.01 ]]
-result:=e04mbf(itmax,msglvl,n,nclin,nctotl,nrowa,a,bl,bu,_
-               cvec,linobj,liwork,lwork,x,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04naf.input.pamphlet b/src/input/e04naf.input.pamphlet
deleted file mode 100644
index 0f80da7..0000000
--- a/src/input/e04naf.input.pamphlet
+++ /dev/null
@@ -1,86 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04naf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-itmax:=20
-msglvl:=-1
-n:=7
-nclin:=7
-nctotl:=14
-nrowa:=7
-nrowh:=7
-ncolh:=7
-bigbnd:=1.0e10
-a:Matrix SF:=
- [[1 ,1 ,1 ,1 ,1 ,1 ,1 ],_
-  [0.15 ,0.04 ,0.02 ,0.04 ,0.02 ,0.01 ,0.03 ],_
-  [0.03 ,0.05 ,0.08 ,0.02 ,0.06 ,0.01 ,0 ],_
-  [0.02 ,0.04 ,0.01 ,0.02 ,0.02 ,0 ,0 ],_
-  [0.02 ,0.03 ,0 ,0 ,0.01 ,0 ,0 ],_
-  [0.7 ,0.75 ,0.8 ,0.75 ,0.8 ,0.97 ,0 ],_
-  [0.02 ,0.06 ,0.08 ,0.12 ,0.02 ,0.01 ,0.97 ]]
-bl:Matrix SF:=
- [[-0.01 ,-0.1 ,-0.01 ,-0.04 ,-0.1 ,-0.01 ,-0.01 ,-0.13 ,_
-   -1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-1.0e+21 ,-0.0992 ,-0.003 ]]
-bu:Matrix SF:=
- [[0.01 ,0.15 ,0.03 ,0.02 ,0.05 ,1.0e+21 ,1.0e+21 ,-0.13 ,_
-  -0.0049 ,-0.0064 ,-0.0037 ,-0.0012 ,1.0e+21 ,0.002 ]]
-cvec:Matrix SF:=
- [[-0.02 ,-0.2 ,-0.2 ,-0.2 ,-0.2 ,0.04 ,0.04 ]]
-featol:Matrix SF:=
- [[0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,_
-   0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,_
-   0.1053e-7 ,0.1053e-7 ,0.1053e-7 ,0.1053e-7 ]]
-hess:Matrix SF :=
- [[2 ,0 ,0 ,0 ,0 ,0 ,0 ],_
-  [0 ,2 ,0 ,0 ,0 ,0 ,0 ],_
-  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
-  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
-  [0 ,0 ,0 ,0 ,2 ,0 ,0 ],_
-  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ],_
-  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ]]
-cold:=true
-lpp:=false
-orthog:=true
-liwork:=14
-lwork:=238
-x:Matrix SF:=
- [[-0.01 ,-0.03 ,0.0 ,-0.01 ,-0.1 ,0.02 ,0.01 ]]
-istate:Matrix Integer:=
- [[0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-mei:Matrix Expression Integer:=
- [[2 ,0 ,0 ,0 ,0 ,0 ,0 ],_
-  [0 ,2 ,0 ,0 ,0 ,0 ,0 ],_
-  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
-  [0 ,0 ,2 ,2 ,0 ,0 ,0 ],_
-  [0 ,0 ,0 ,0 ,2 ,0 ,0 ],_
-  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ],_
-  [0 ,0 ,0 ,0 ,0 ,-2 ,-2 ]]
-qphess:Asp20(QPHESS):= retract mei
-result:=e04naf(itmax,msglvl,n,nclin,nctotl,nrowa,nrowh,ncolh,_
-               bigbnd,a,bl,bu,cvec,featol,_
-hess,cold,lpp,orthog,liwork,lwork,x,istate,-1,qphess)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04ucf.input.pamphlet b/src/input/e04ucf.input.pamphlet
deleted file mode 100644
index efdaf6a..0000000
--- a/src/input/e04ucf.input.pamphlet
+++ /dev/null
@@ -1,103 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04ucf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-n:=9
-nclin:=4
-ncnln:=14
-nrowa:=4
-nrowj:=14
-nrowr:=9
-a:Matrix SF:=
- [[-1.0 ,1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
-  [0.0 ,-1.0 ,1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
-  [0.0 ,0.0 ,1.0 ,-1.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ],_
-  [0.0 ,0.0 ,0.0 ,1.0 ,-1.0 ,0.0 ,0.0 ,0.0 ,0.0 ]]
-bl:Matrix SF:=
- [[0.0 ,-1.0e15 ,-1.0 ,-1.0e15 ,0.0 ,0.0 ,0.0 ,-1.0e15 ,_
-  -1.0e15 ,0.0 ,0.0 ,0.0 ,0.0 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,_
-  -1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,_
-  -1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ,-1.0e15 ]]
-bu:Matrix SF:=
- [[1.0e15 ,1.0e15 ,1.0 ,1.0e15 ,1.0e15 ,1.0e15 ,1.0e15 ,_
-  0.0 ,0.0 ,1.0e15 ,1.0e15 ,1.0e15 ,1.0e15 ,0.0 ,0.0 ,_
-  0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ,0.0 ]]
-liwork:=59
-lwork:=968
-sta:=false
-cra:= 0.01
-der:=3
-fea:=0.1053671201E-7
-fun:=0.4373903510E-14
-hes:= true
-infb:= 1.00E+15
-infs:= 1.00E+15
-linf:= 0.1053671201E-7
-lint:= 0.9
-list:= true
-maji:= 30
-majp:= 1
-mini:= 81
-minp:= 0
-mon:= -1
-nonf:= 1.05E-08
-opt:= 3.26E-08
-ste:= 2.0
-stao:= 1
-stac:= 1
-stoo:= 9
-stoc:= 9
-ve:= 3
-istate:Matrix Integer:=new(1,27,0)
-cjac:Matrix SF:=new(14,9,0.0)
-clamda:Matrix SF:=new(1,27,0.0)
-r:Matrix SF:=new(9,9,0.0)
-x:Matrix SF:=
- [[0.1 ,0.125 ,0.666666 ,0.142857 ,0.111111 ,0.2 ,0.25 ,-0.2 ,-0.25 ]]
-vef:Vector Expression Float:=vector
- [X[1]^2 + X[6]^2 ,_
- (X[2] - X[1])^2 + (X[7] - X[6])^2 ,_
- (X[3] - X[1])^2 + X[6]*2 ,_
- (X[1] - X[4])^2 + (X[6] - X[8])^2 ,_
- (X[1] - X[5])^2 + (X[6] - X[9])^2 ,_
- X[2]^2 + X[7]^2 ,_
- (X[3] - X[2])^2 + X[7]^2 ,_
- (X[4] - X[2])^2 + (X[8] - X[7])^2 ,_
- (X[2] - X[5])^2 + (X[7] - X[9])^2 ,_
- (X[4] - X[3])^2 + X[8]^2 ,_
- (X[5] - X[3])^2 + X[9]^2 ,_
- X[4]^2 + X[8]^2 ,_
- (X[4] - X[5])^2 + (X[9] - X[8])^2 ,_
- X[5]^2 + X[9]^2 ]
-confun:Asp55(CONFUN):= retract vef
-ef:Expression Float:=(-X[2]*X[6] + X[1]*X[7] - X[3]*X[7] - _
-                      X[5]*X[8] + X[4]*X[9] + X[3]*X[8])::EXPR FLOAT
-objfun:Asp49(OBJFUN):=retract ef
-result:=e04ucf(n,nclin,ncnln,nrowa,nrowj,nrowr,a,bl,bu,liwork,lwork,_
-sta,cra,der,fea,fun,hes,infb,infs,linf,lint,list,maji,majp,_
-mini,minp,mon,nonf,opt,ste,stao,stac,stoo,stoc,ve,istate,_
-cjac,clamda,r,x,-1,confun,objfun)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/e04ycf.input.pamphlet b/src/input/e04ycf.input.pamphlet
deleted file mode 100644
index e86d1ad..0000000
--- a/src/input/e04ycf.input.pamphlet
+++ /dev/null
@@ -1,40 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input e04ycf.input}
-\author{The Axiom Team}
-\maketitle
-\begin{abstract}
-\end{abstract}
-\eject
-\tableofcontents
-\eject
-\section{License}
-\begin{chunk}{license}
---Copyright The Numerical Algorithms Group Limited 1994.
-\end{chunk}
-\begin{chunk}{*}
-\getchunk{license}
-)clear all
-showArrayValues true
-showScalarValues true
-
-job:=0
-m:=15
-n:=3
-fsumsq:=0.0082148773065789729
-s:Matrix SF:=
- [[4.0965034571419325 ,1.5949579400198182 ,0.061258491120317927 ]]
-lv:=3
-v:Matrix SF:=
- [[0.9353959087369782,- 0.25922842715810668, - 0.24048932714625143,_
-   0.35295122094583287, 0.64323460044494585, 0.67946647041872243,_
-   - 0.021445968175122822, - 0.72045115819433947, 0.69317400348358882]]
-result:=e04ycf(job,m,n,fsumsq,s,lv,v,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}