diff --git a/books/bookvol10.4.pamphlet b/books/bookvol10.4.pamphlet
index 050bf59..f5526b9 100644
--- a/books/bookvol10.4.pamphlet
+++ b/books/bookvol10.4.pamphlet
@@ -102867,11 +102867,1284 @@ NagLapack(): Exports == Implementation where
 )set message auto off
 )clear all
 
---S 1 of 1
+--S 1 of 130
 )show NagMatrixOperationsPackage
+--R 
+--R NagMatrixOperationsPackage  is a package constructor
+--R Abbreviation for NagMatrixOperationsPackage is NAGF01 
+--R This constructor is exposed in this frame.
+--R Issue )edit bookvol10.4.pamphlet to see algebra source code for NAGF01 
+--R
+--R------------------------------- Operations --------------------------------
+--R f01brf : (Integer,Integer,Integer,Integer,DoubleFloat,Boolean,Boolean,List(Boolean),Matrix(DoubleFloat),Matrix(Integer),Matrix(Integer),Integer) -> Result
+--R f01bsf : (Integer,Integer,Integer,Matrix(Integer),Matrix(Integer),Matrix(Integer),Matrix(Integer),Boolean,DoubleFloat,Boolean,Matrix(Integer),Matrix(DoubleFloat),Integer) -> Result
+--R f01maf : (Integer,Integer,Integer,Integer,List(Boolean),Matrix(DoubleFloat),Matrix(Integer),Matrix(Integer),DoubleFloat,DoubleFloat,Integer) -> Result
+--R f01mcf : (Integer,Matrix(DoubleFloat),Integer,Matrix(Integer),Integer) -> Result
+--R f01qcf : (Integer,Integer,Integer,Matrix(DoubleFloat),Integer) -> Result
+--R f01qdf : (String,String,Integer,Integer,Matrix(DoubleFloat),Integer,Matrix(DoubleFloat),Integer,Integer,Matrix(DoubleFloat),Integer) -> Result
+--R f01qef : (String,Integer,Integer,Integer,Integer,Matrix(DoubleFloat),Matrix(DoubleFloat),Integer) -> Result
+--R f01rcf : (Integer,Integer,Integer,Matrix(Complex(DoubleFloat)),Integer) -> Result
+--R f01rdf : (String,String,Integer,Integer,Matrix(Complex(DoubleFloat)),Integer,Matrix(Complex(DoubleFloat)),Integer,Integer,Matrix(Complex(DoubleFloat)),Integer) -> Result
+--R f01ref : (String,Integer,Integer,Integer,Integer,Matrix(Complex(DoubleFloat)),Matrix(Complex(DoubleFloat)),Integer) -> Result
+--R
 --E 1
 
-)spool
+)clear all
+
+--S 2 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 2
+
+--S 3 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 3
+
+--S 4 of 130
+n:=6
+--R 
+--R
+--R   (3)  6
+--R                                                        Type: PositiveInteger
+--E 4
+
+--S 5 of 130
+nz:=15
+--R 
+--R
+--R   (4)  15
+--R                                                        Type: PositiveInteger
+--E 5
+
+--S 6 of 130
+licn:=150
+--R 
+--R
+--R   (5)  150
+--R                                                        Type: PositiveInteger
+--E 6
+
+--S 7 of 130
+lirn:=75
+--R 
+--R
+--R   (6)  75
+--R                                                        Type: PositiveInteger
+--E 7
+
+--S 8 of 130
+pivot:=0.1
+--R 
+--R
+--R   (7)  0.1
+--R                                                                  Type: Float
+--E 8
+
+--S 9 of 130
+lblock:=true
+--R 
+--R
+--R   (8)  true
+--R                                                                Type: Boolean
+--E 9
+
+--S 10 of 130
+grow:=true
+--R 
+--R
+--R   (9)  true
+--R                                                                Type: Boolean
+--E 10
+
+--S 11 of 130
+abort:=[true ,true ,false ,true ]
+--R 
+--R
+--R   (10)  [true,true,false,true]
+--R                                                          Type: List(Boolean)
+--E 11
+
+--S 12 of 130
+a:Matrix SF:=
+ [[5 ,2 ,-1 ,2 ,3 ,-2 ,1 ,1 ,-1 ,-1 ,2 ,-3 ,-1 ,-1 ,6_
+ ,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 ,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 ,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 ,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 ,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     [5., 2., - 1., 2., 3., - 2., 1., 1., - 1., - 1., 2., - 3., - 1., - 1., 6.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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 12
+
+--S 13 of 130
+irn:Matrix Integer:=
+ [[1 ,2 ,2 ,2 ,3 ,4 ,4 ,4 ,5 ,5 ,5 ,5 ,6 ,6 ,6_
+  ,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 ,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 ,0 ,0 ,0 ,0 ]]
+--R 
+--R
+--R   (12)
+--R   [
+--R     [1, 2, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+--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,
+--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,
+--R      0, 0, 0]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 13
+
+--S 14 of 130
+icn:Matrix Integer:=
+ [[1 ,2 ,3 ,4 ,3 ,1 ,4 ,5 ,1 ,4 ,5 ,6 ,1 ,2 ,6_
+  ,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 ,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 ,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 ,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 ,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   (13)
+--R   [
+--R     [1, 2, 3, 4, 3, 1, 4, 5, 1, 4, 5, 6, 1, 2, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+--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,
+--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,
+--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,
+--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,
+--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,
+--R      0, 0, 0, 0, 0, 0]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 14
+
+--S 15 of 130
+ifail:= 110
+--R 
+--R
+--R   (14)  110
+--R                                                        Type: PositiveInteger
+--E 15
+
+--S 16 of 130
+-- result:=f01brf(n,nz,licn,lirn,pivot,lblock,grow,abort,a,irn,icn,ifail)
+--E 16
+
+)clear all
+
+--S 17 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 17
+
+--S 18 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 18
+
+--S 19 of 130
+n := 6
+--R 
+--R
+--R   (3)  6
+--R                                                        Type: PositiveInteger
+--E 19
+
+--S 20 of 130
+nz := 15
+--R 
+--R
+--R   (4)  15
+--R                                                        Type: PositiveInteger
+--E 20
+
+--S 21 of 130
+nzmax := 50
+--R 
+--R
+--R   (5)  50
+--R                                                        Type: PositiveInteger
+--E 21
+
+--S 22 of 130
+licn := 3*nzmax
+--R 
+--R
+--R   (6)  150
+--R                                                        Type: PositiveInteger
+--E 22
+
+--S 23 of 130
+ta1:List SF := new(licn-nz,0.0);
+--R 
+--R
+--R                                                      Type: List(DoubleFloat)
+--E 23
+
+--S 24 of 130
+ta2:List SF := [5.,2.,-1.,2.,3.,-2.,1.,1.,-1.,-1.,2.,-3.,-1.,-1.,6.] ;
+--R 
+--R
+--R                                                      Type: List(DoubleFloat)
+--E 24
+
+--S 25 of 130
+ta3 := concat(ta2,ta1);
+--R 
+--R
+--R                                                      Type: List(DoubleFloat)
+--E 25
+
+--S 26 of 130
+a :Matrix SF := matrix [ta3] 
+--R 
+--R
+--R   (10)
+--R   [
+--R     [5., 2., - 1., 2., 3., - 2., 1., 1., - 1., - 1., 2., - 3., - 1., - 1., 6.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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 26
+
+--S 27 of 130
+lirn := (3*nzmax/2)::Integer
+--R 
+--R
+--R   (11)  75
+--R                                                                Type: Integer
+--E 27
+
+--S 28 of 130
+tr1:List Integer := new(lirn-nz,0);
+--R 
+--R
+--R                                                          Type: List(Integer)
+--E 28
+
+--S 29 of 130
+tr2 := [1,2,2,2,3,4,4,4,5,5,5,5,6,6,6];
+--R 
+--R
+--R                                                  Type: List(PositiveInteger)
+--E 29
+
+--S 30 of 130
+tr3 := concat(tr2,tr1);
+--R 
+--R
+--R                                                          Type: List(Integer)
+--E 30
+
+--S 31 of 130
+irn := matrix[tr3]
+--R 
+--R
+--R   (15)
+--R   [
+--R     [1, 2, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+--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,
+--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,
+--R      0, 0, 0]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 31
+
+--S 32 of 130
+tc1:List Integer := new(licn-nz,0);
+--R 
+--R
+--R                                                          Type: List(Integer)
+--E 32
+
+--S 33 of 130
+tc2 := [1,2,3,4,3,1,4,5,1,4,5,6,1,2,6];
+--R 
+--R
+--R                                                  Type: List(PositiveInteger)
+--E 33
+
+--S 34 of 130
+tc3 := concat(tc2,tc1);
+--R 
+--R
+--R                                                          Type: List(Integer)
+--E 34
+
+--S 35 of 130
+icn :=  matrix [tc3] 
+--R 
+--R
+--R   (19)
+--R   [
+--R     [1, 2, 3, 4, 3, 1, 4, 5, 1, 4, 5, 6, 1, 2, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0,
+--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,
+--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,
+--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,
+--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,
+--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,
+--R      0, 0, 0, 0, 0, 0]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 35
+
+--S 36 of 130
+pivot :SF:= 0.1
+--R 
+--R
+--R   (20)  0.10000000000000001
+--R                                                            Type: DoubleFloat
+--E 36
+
+--S 37 of 130
+lblock := true
+--R 
+--R
+--R   (21)  true
+--R                                                                Type: Boolean
+--E 37
+
+--S 38 of 130
+grow := true
+--R 
+--R
+--R   (22)  true
+--R                                                                Type: Boolean
+--E 38
+
+--S 39 of 130
+abort := [true,true,false,true]
+--R 
+--R
+--R   (23)  [true,true,false,true]
+--R                                                          Type: List(Boolean)
+--E 39
+
+--S 40 of 130
+ifail := 110
+--R 
+--R
+--R   (24)  110
+--R                                                        Type: PositiveInteger
+--E 40
+
+--S 41 of 130
+--f01brfResult:=f01brf(n,nz,licn,lirn,pivot,lblock,grow,abort,a,irn,icn,ifail);
+--E 41
+
+--S 42 of 130
+eta :SF:= .1 -- now run f01bsf
+--R 
+--R
+--R   (25)  0.10000000000000001
+--R                                                            Type: DoubleFloat
+--E 42
+
+--S 43 of 130
+abort := true
+--R 
+--R
+--R   (26)  true
+--R                                                                Type: Boolean
+--E 43
+
+--S 44 of 130
+-- keep := f01brfResult.'ikeep;
+--E 44
+
+--S 45 of 130
+-- disp := f01brfResult.'idisp;
+--E 45
+
+--S 46 of 130
+-- cn :=  f01brfResult.'icn;
+--E 46
+
+--S 47 of 130
+ta2:=[10.,12.,-3.,-1.,15.,-2.,10.,-1.,-1.,-5.,1.,-1.,-1.,-2.,6.];
+--R 
+--R
+--R                                                      Type: List(DoubleFloat)
+--E 47
+
+--S 48 of 130
+a:= matrix [concat (ta2,ta1)]
+--R 
+--R
+--R   (28)
+--R   [
+--R     [10., 12., - 3., - 1., 15., - 2., 10., - 1., - 1., - 5., 1., - 1., - 1.,
+--R      - 2., 6., 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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0.]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 48
+
+--S 49 of 130
+-- result:=f01bsf(n,nz,licn,matrix [tr2],matrix [tc2],cn,keep,grow,eta,abort,_
+--   matrix [[disp(1,1),disp(1,2)]],a,ifail)
+--E 49
+
+)clear all
+
+--S 50 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 50
+
+--S 51 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 51
+
+--S 52 of 130
+n:=16
+--R 
+--R
+--R   (3)  16
+--R                                                        Type: PositiveInteger
+--E 52
+
+--S 53 of 130
+nz:=40
+--R 
+--R
+--R   (4)  40
+--R                                                        Type: PositiveInteger
+--E 53
+
+--S 54 of 130
+licn:=90
+--R 
+--R
+--R   (5)  90
+--R                                                        Type: PositiveInteger
+--E 54
+
+--S 55 of 130
+lirn:=50
+--R 
+--R
+--R   (6)  50
+--R                                                        Type: PositiveInteger
+--E 55
+
+--S 56 of 130
+abort:=[true ,true ,true ]
+--R 
+--R
+--R   (7)  [true,true,true]
+--R                                                          Type: List(Boolean)
+--E 56
+
+--S 57 of 130
+avals:Matrix SF:= 
+ [[1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0_
+ ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
+ ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
+ ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
+ ,-0.25 ,-0.25 ,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 ,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   (8)
+--R   [
+--R     [1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., - 0.25,
+--R      - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25,
+--R      - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, - 0.25,
+--R      - 0.25, - 0.25, - 0.25, - 0.25, - 0.25, 0., 0., 0., 0., 0., 0., 0., 0.,
+--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., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
+--R      0., 0., 0., 0., 0., 0.]
+--R     ]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 57
+
+--S 58 of 130
+irn:Matrix Integer:=
+ [[1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14_
+ ,15 ,16 ,1 ,2 ,3 ,5 ,6 ,7 ,9 ,10 ,11 ,13 ,14_
+ ,15 ,1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12_
+ ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
+--R 
+--R
+--R   (9)
+--R   [
+--R     [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 1, 2, 3, 5, 6, 7,
+--R      9, 10, 11, 13, 14, 15, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 0, 0, 0, 0,
+--R      0, 0, 0, 0, 0, 0]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 58
+
+--S 59 of 130
+icn:Matrix Integer:=
+ [[1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14_
+ ,15 ,16 ,2 ,3 ,4 ,6 ,7 ,8 ,10 ,11 ,12 ,14 ,15 ,16_
+ ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14 ,15 ,16 ,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 ,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   (10)
+--R   [
+--R     [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 2, 3, 4, 6, 7, 8,
+--R      10, 11, 12, 14, 15, 16, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 0, 0,
+--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,
+--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]
+--R     ]
+--R                                                        Type: Matrix(Integer)
+--E 59
+
+--S 60 of 130
+droptl:= 0.1
+--R 
+--R
+--R   (11)  0.1
+--R                                                                  Type: Float
+--E 60
+
+--S 61 of 130
+densw:= 0.8
+--R 
+--R
+--R   (12)  0.8
+--R                                                                  Type: Float
+--E 61
+
+--S 62 of 130
+ifail:=111
+--R 
+--R
+--R   (13)  111
+--R                                                        Type: PositiveInteger
+--E 62
+
+--S 63 of 130
+-- result:=f01maf(n,nz,licn,lirn,abort,avals,irn,icn,droptl,densw,ifail)
+--E 63
+
+)clear all
+
+--S 64 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 64
+
+--S 65 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 65
+
+--S 66 of 130
+n:=6
+--R 
+--R
+--R   (3)  6
+--R                                                        Type: PositiveInteger
+--E 66
+
+--S 67 of 130
+avals:Matrix SF:=
+ [[1.0 ,2.0 ,5.0 ,3.0 ,13.0 ,16.0 ,5.0 ,14.0_
+ ,18.0 ,8.0 ,55.0 ,24.0 ,17.0 ,77.0 ]]
+--R 
+--R
+--R   (4)  [1.  2.  5.  3.  13.  16.  5.  14.  18.  8.  55.  24.  17.  77.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 67
+
+--S 68 of 130
+lal:= 14
+--R 
+--R
+--R   (5)  14
+--R                                                        Type: PositiveInteger
+--E 68
+
+--S 69 of 130
+nrow:Matrix Integer:= [[1 ,2 ,2 ,1 ,5 ,3 ]]
+--R 
+--R
+--R   (6)  [1  2  2  1  5  3]
+--R                                                        Type: Matrix(Integer)
+--E 69
+
+--S 70 of 130
+-- result:=f01mcf(n,avals,lal,nrow,-1)
+--E 70
+
+)clear all
+
+--S 71 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 71
+
+--S 72 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 72
+
+--S 73 of 130
+m:=5
+--R 
+--R
+--R   (3)  5
+--R                                                        Type: PositiveInteger
+--E 73
+
+--S 74 of 130
+n:=3
+--R 
+--R
+--R   (4)  3
+--R                                                        Type: PositiveInteger
+--E 74
+
+--S 75 of 130
+lda:=5
+--R 
+--R
+--R   (5)  5
+--R                                                        Type: PositiveInteger
+--E 75
+
+--S 76 of 130
+a:Matrix SF:=
+ [[2.0 ,2.5 ,2.5 ],_
+  [2.0 ,2.5 ,2.5 ],_
+  [1.6 ,-0.4 ,2.8 ],_
+  [2.0 ,-0.5 ,0.5 ],_
+  [1.2 ,-0.3 ,-2.9 ]]
+--R 
+--R
+--R        +        2.                   2.5                   2.5         +
+--R        |                                                               |
+--R        |        2.                   2.5                   2.5         |
+--R        |                                                               |
+--R   (6)  |1.6000000000000001  - 0.39999999999999997   2.7999999999999998 |
+--R        |                                                               |
+--R        |        2.                  - 0.5                  0.5         |
+--R        |                                                               |
+--R        +       1.2          - 0.30000000000000004  - 2.9000000000000004+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 76
+
+--S 77 of 130
+-- result:=f01qcf(m,n,lda,a,-1)
+--E 77
+
+)clear all
+
+--S 78 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 78
+
+--S 79 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 79
+
+--S 80 of 130
+trans:="n"
+--R 
+--R
+--R   (3)  "n"
+--R                                                                 Type: String
+--E 80
+
+--S 81 of 130
+wheret:="i"
+--R 
+--R
+--R   (4)  "i"
+--R                                                                 Type: String
+--E 81
+
+--S 82 of 130
+m:=5
+--R 
+--R
+--R   (5)  5
+--R                                                        Type: PositiveInteger
+--E 82
+
+--S 83 of 130
+n:=3
+--R 
+--R
+--R   (6)  3
+--R                                                        Type: PositiveInteger
+--E 83
+
+--S 84 of 130
+a:Matrix SF:=
+ [[2.0 ,2.5 ,2.5 ],_
+  [2.0 ,2.5 ,2.5 ],_
+  [1.6 ,-0.4 ,2.8 ],_
+  [2.0 ,-0.5 ,0.5 ],_
+  [1.2 ,-0.3 ,-2.9 ]]
+--R 
+--R
+--R        +        2.                   2.5                   2.5         +
+--R        |                                                               |
+--R        |        2.                   2.5                   2.5         |
+--R        |                                                               |
+--R   (7)  |1.6000000000000001  - 0.39999999999999997   2.7999999999999998 |
+--R        |                                                               |
+--R        |        2.                  - 0.5                  0.5         |
+--R        |                                                               |
+--R        +       1.2          - 0.30000000000000004  - 2.9000000000000004+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 84
+
+--S 85 of 130
+lda:=5
+--R 
+--R
+--R   (8)  5
+--R                                                        Type: PositiveInteger
+--E 85
+
+--S 86 of 130
+zeta:Matrix SF:= [[0.0 ,0.0 ,0.0 ]]
+--R 
+--R
+--R   (9)  [0.  0.  0.]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 86
+
+--S 87 of 130
+ncolb:=2
+--R 
+--R
+--R   (10)  2
+--R                                                        Type: PositiveInteger
+--E 87
+
+--S 88 of 130
+ldb:=5
+--R 
+--R
+--R   (11)  5
+--R                                                        Type: PositiveInteger
+--E 88
+
+--S 89 of 130
+b:Matrix SF:=
+ [[1.1 ,0.0 ],_
+  [0.9 ,0.0 ],_
+  [0.6 ,1.32 ],_
+  [0.0 ,1.1 ],_
+  [-0.8 ,-0.26 ]]
+--R 
+--R
+--R         + 1.1000000000000001             0.          +
+--R         |                                            |
+--R         | 0.89999999999999991            0.          |
+--R         |                                            |
+--R   (12)  | 0.59999999999999998    1.3199999999999998  |
+--R         |                                            |
+--R         |         0.             1.1000000000000001  |
+--R         |                                            |
+--R         +- 0.79999999999999993  - 0.25999999999999995+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 89
+
+--S 90 of 130
+-- result:=f01qdf(trans,wheret,m,n,a,lda,zeta,ncolb,ldb,b,-1)
+--E 90
+
+)clear all
+
+--S 91 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 91
+
+--S 92 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 92
+
+--S 93 of 130
+wheret:="s"
+--R 
+--R
+--R   (3)  "s"
+--R                                                                 Type: String
+--E 93
+
+--S 94 of 130
+m:=5
+--R 
+--R
+--R   (4)  5
+--R                                                        Type: PositiveInteger
+--E 94
+
+--S 95 of 130
+n:=3
+--R 
+--R
+--R   (5)  3
+--R                                                        Type: PositiveInteger
+--E 95
+
+--S 96 of 130
+ncolq:=5
+--R 
+--R
+--R   (6)  5
+--R                                                        Type: PositiveInteger
+--E 96
+
+--S 97 of 130
+lda:=5
+--R 
+--R
+--R   (7)  5
+--R                                                        Type: PositiveInteger
+--E 97
+
+--S 98 of 130
+zeta:Matrix SF:= [[1.2247 ,1.1547 ,1.2649 ]]
+--R 
+--R
+--R   (8)  [1.2246999999999999  1.1547000000000001  1.2648999999999999]
+--R                                                    Type: Matrix(DoubleFloat)
+--E 98
+
+--S 99 of 130
+a:Matrix SF:=
+ [[-4.0 ,-2.0 ,-3.0 ,0.0 ,0.0 ],_
+  [0.4085 ,-3.0 ,-2.0 ,0.0 ,0.0 ],_
+  [0.3266 ,-0.4619 ,-4.0 ,0.0 ,0.0 ],_
+  [0.4082 ,-0.5774 ,0.0 ,0.0 ,0.0 ],_
+  [0.2449 ,-0.3464 ,-0.6326 ,0.0 ,0.0 ]]
+--R 
+--R
+--R   (9)
+--R   +       - 4.                  - 2.                   - 3.           0.  0.+
+--R   |                                                                         |
+--R   |0.40849999999999997          - 3.                   - 2.           0.  0.|
+--R   |                                                                         |
+--R   |      0.3266         - 0.46189999999999998          - 4.           0.  0.|
+--R   |                                                                         |
+--R   |0.40820000000000001  - 0.57739999999999991           0.            0.  0.|
+--R   |                                                                         |
+--R   +0.24490000000000001  - 0.34640000000000004  - 0.63259999999999994  0.  0.+
+--R                                                    Type: Matrix(DoubleFloat)
+--E 99
+
+--S 100 of 130
+-- result:=f01qef(wheret,m,n,ncolq,lda,zeta,a,-1)
+--E 100
+
+)clear all
+
+--S 101 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 101
+
+--S 102 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 102
+
+--S 103 of 130
+m:=5
+--R 
+--R
+--R   (3)  5
+--R                                                        Type: PositiveInteger
+--E 103
+
+--S 104 of 130
+n:=3
+--R 
+--R
+--R   (4)  3
+--R                                                        Type: PositiveInteger
+--E 104
+
+--S 105 of 130
+lda:=5
+--R 
+--R
+--R   (5)  5
+--R                                                        Type: PositiveInteger
+--E 105
+
+--S 106 of 130
+a:Matrix Complex SF:=
+ [[complex(0,0.5) ,complex(-0.5,1.5),complex(-1.0,1.0)],_
+  [complex(0.4,0.3),complex(0.9,1.3),complex(0.2,1.4) ],_
+  [complex(0.4,0), complex(-0.4,0.4),complex(1.8,0)],_
+  [complex(0.3,-0.4),complex(0.1,0.7),complex(0,0)],_
+  [complex(0,-0.3),complex(0.3,0.3),complex(0,2.4) ]]
+--R 
+--R
+--R   (6)
+--R   [[0.5 %i,- 0.5 + 1.5 %i,- 1. + %i],
+--R
+--R     [0.40000000000000002 + 0.29999999999999999 %i,
+--R      0.89999999999999991 + 1.2999999999999998 %i,
+--R      0.20000000000000001 + 1.3999999999999999 %i]
+--R     ,
+--R
+--R     [0.40000000000000002, - 0.39999999999999997 + 0.40000000000000002 %i,
+--R      1.7999999999999998]
+--R     ,
+--R
+--R     [0.29999999999999999 - 0.39999999999999997 %i,
+--R      0.10000000000000001 + 0.69999999999999996 %i, 0.]
+--R     ,
+--R
+--R     [- 0.30000000000000004 %i, 0.29999999999999999 + 0.29999999999999999 %i,
+--R      2.3999999999999999 %i]
+--R     ]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 106
+
+--S 107 of 130
+-- result:=f01rcf(m,n,lda,a,-1)
+--E 107
+
+)clear all
+
+--S 108 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 108
+
+--S 109 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 109
+
+--S 110 of 130
+trans:="n"
+--R 
+--R
+--R   (3)  "n"
+--R                                                                 Type: String
+--E 110
+
+--S 111 of 130
+wheret:="i"
+--R 
+--R
+--R   (4)  "i"
+--R                                                                 Type: String
+--E 111
+
+--S 112 of 130
+m:=5
+--R 
+--R
+--R   (5)  5
+--R                                                        Type: PositiveInteger
+--E 112
+
+--S 113 of 130
+n:=3
+--R 
+--R
+--R   (6)  3
+--R                                                        Type: PositiveInteger
+--E 113
+
+--S 114 of 130
+a:Matrix Complex SF:=
+ [[      0.5*%i ,-0.5 + 1.5*%i ,-1.0 + 1.0*%i ],_
+  [0.4 + 0.3*%i , 0.9 + 1.3*%i , 0.2 + 1.4*%i ],_
+  [0.4          ,-0.4 + 0.4*%i , 1.8 ],_
+  [0.3 - 0.4*%i , 0.1 + 0.7*%i , 0.0 ],_
+  [     -0.3*%i , 0.3 + 0.3*%i , 2.4 ]]
+--R 
+--R
+--R   (7)
+--R   [[0.5 %i,- 0.5 + 1.5 %i,- 1. + %i],
+--R
+--R     [0.40000000000000002 + 0.29999999999999999 %i,
+--R      0.89999999999999991 + 1.2999999999999998 %i,
+--R      0.20000000000000001 + 1.3999999999999999 %i]
+--R     ,
+--R
+--R     [0.40000000000000002, - 0.39999999999999997 + 0.40000000000000002 %i,
+--R      1.7999999999999998]
+--R     ,
+--R
+--R     [0.29999999999999999 - 0.40000000000000002 %i,
+--R      0.10000000000000001 + 0.69999999999999996 %i, 0.]
+--R     ,
+--R
+--R     [- 0.29999999999999999 %i, 0.29999999999999999 + 0.29999999999999999 %i,
+--R      2.3999999999999999]
+--R     ]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 114
+
+--S 115 of 130
+lda:=5
+--R 
+--R
+--R   (8)  5
+--R                                                        Type: PositiveInteger
+--E 115
+
+--S 116 of 130
+theta:Matrix Complex SF:= [[0.0 ,0.0 ,0.0 ]]
+--R 
+--R
+--R   (9)  [0.  0.  0.]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 116
+
+--S 117 of 130
+ncolb:=2
+--R 
+--R
+--R   (10)  2
+--R                                                        Type: PositiveInteger
+--E 117
+
+--S 118 of 130
+ldb:=5
+--R 
+--R
+--R   (11)  5
+--R                                                        Type: PositiveInteger
+--E 118
+
+--S 119 of 130
+b:Matrix Complex SF :=
+ [[-0.55 + 1.05*%i , 0.45 + 1.05*%i ],_
+  [ 0.49 + 0.93*%i , 1.09 + 0.13*%i ],_
+  [ 0.56 - 0.16*%i , 0.64 + 0.16*%i ],_
+  [ 0.39 + 0.23*%i ,-0.39 - 0.23*%i ],_
+  [ 1.13 + 0.83*%i ,-1.13 + 0.77*%i ]]
+--R 
+--R
+--R   (12)
+--R   [
+--R     [- 0.54999999999999993 + 1.0499999999999998 %i,
+--R      0.44999999999999996 + 1.0499999999999998 %i]
+--R     ,
+--R    [0.48999999999999999 + 0.92999999999999994 %i,1.0899999999999999 + 0.13 %i],
+--R
+--R     [0.56000000000000005 - 0.15999999999999998 %i,
+--R      0.6399999999999999 + 0.15999999999999998 %i]
+--R     ,
+--R
+--R     [0.39000000000000001 + 0.22999999999999998 %i,
+--R      - 0.38999999999999996 - 0.22999999999999998 %i]
+--R     ,
+--R
+--R     [1.1299999999999999 + 0.82999999999999996 %i,
+--R      - 1.1299999999999999 + 0.77000000000000002 %i]
+--R     ]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 119
+
+--S 120 of 130
+-- result:=f01rdf(trans,wheret,m,n,a,lda,theta,ncolb,ldb,b,-1)
+--E 120
+
+)clear all
+
+--S 121 of 130
+showArrayValues true
+--R 
+--R
+--R   (1)  true
+--R                                                                Type: Boolean
+--E 121
+
+--S 122 of 130
+showScalarValues true
+--R 
+--R
+--R   (2)  true
+--R                                                                Type: Boolean
+--E 122
+
+--S 123 of 130
+wheret:="s"
+--R 
+--R
+--R   (3)  "s"
+--R                                                                 Type: String
+--E 123
+
+--S 124 of 130
+m:=5
+--R 
+--R
+--R   (4)  5
+--R                                                        Type: PositiveInteger
+--E 124
+
+--S 125 of 130
+n:=3
+--R 
+--R
+--R   (5)  3
+--R                                                        Type: PositiveInteger
+--E 125
+
+--S 126 of 130
+ncolq:=2
+--R 
+--R
+--R   (6)  2
+--R                                                        Type: PositiveInteger
+--E 126
+
+--S 127 of 130
+lda:=5
+--R 
+--R
+--R   (7)  5
+--R                                                        Type: PositiveInteger
+--E 127
+
+--S 128 of 130
+theta:Matrix Complex SF:= [[1 + 0.5*%i ,1.0954-0.3333*%i ,1.2649-1.1565*%i ]]
+--R 
+--R
+--R   (8)
+--R   [
+--R     [1. + 0.5 %i, 1.0953999999999999 - 0.33329999999999999 %i,
+--R      1.2648999999999999 - 1.1564999999999999 %i]
+--R     ]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 128
+
+--S 129 of 130
+a:Matrix Complex SF :=
+ [[ 1             , 1      +1     *%i , 1     + %i ],_
+  [-0.2  - 0.4*%i ,-2                 ,-1     - %i ],_
+  [-0.32 - 0.16*%i,-0.3505 +0.2629*%i ,-3          ],_
+  [-0.4  + 0.2*%i ,         0.5477*%i , 0.0        ],_
+  [-0.12 + 0.24*%i, 0.1972 +0.2629*%i , 0.6325     ]]
+--R 
+--R
+--R   (9)
+--R   [[1.,1. + %i,1. + %i],
+--R    [- 0.19999999999999998 - 0.40000000000000002 %i,- 2.,- 1. - %i],
+--R
+--R     [- 0.31999999999999995 - 0.15999999999999998 %i,
+--R      - 0.35050000000000003 + 0.26290000000000002 %i, - 3.]
+--R     ,
+--R    [- 0.39999999999999997 + 0.20000000000000001 %i,0.54769999999999996 %i,0.],
+--R
+--R     [- 0.12 + 0.23999999999999999 %i,
+--R      0.19719999999999999 + 0.26290000000000002 %i, 0.63249999999999995]
+--R     ]
+--R                                           Type: Matrix(Complex(DoubleFloat))
+--E 129
+
+--S 130 of 130
+-- result:=f01ref(wheret,m,n,ncolq,lda,theta,a,-1)
+--E 130
+
+)spool
+ 
 )lisp (bye)
 \end{chunk}
 
diff --git a/changelog b/changelog
index 2720fdc..d400896 100644
--- a/changelog
+++ b/changelog
@@ -1,5 +1,19 @@
+20141207 tpd src/axiom-website/patches.html 20141207.02.tpd.patch
+20141207 tpd books/bookvol10.4 NAGF01 absorb src/input test case files
+20141207 tpd src/input/Makefile remove f01 files
+20141207 tpd src/input/f01brf.input absorbed, removed
+20141207 tpd src/input/f01bsf.input absorbed, removed
+20141207 tpd src/input/f01maf.input absorbed, removed
+20141207 tpd src/input/f01mcf.input absorbed, removed
+20141207 tpd src/input/f01qcf.input absorbed, removed
+20141207 tpd src/input/f01qdf.input absorbed, removed
+20141207 tpd src/input/f01qef.input absorbed, removed
+20141207 tpd src/input/f01rcf.input absorbed, removed
+20141207 tpd src/input/f01rdf.input absorbed, removed
+20141207 tpd src/input/f01ref.input absorbed, removed
 20141207 tpd src/axiom-website/patches.html 20141207.01.tpd.patch
 20141207 tpd books/bookvol10.4 D03AGNT absorb src/input test case files
+20141207 tpd src/input/Makefile remove d03 files
 20141207 tpd src/input/d03edf.input absorbed, removed
 20141207 tpd src/input/d03eef.input absorbed, removed
 20141207 tpd src/input/d03faf.input absorbed, removed
@@ -7,6 +21,7 @@
 20141206 tpd books/bookvol10.4 D02AGNT absorb src/input test case files
 20141206 tpd src/axiom-website/patches.html 20141206.05.tpd.patch
 20141206 tpd books/bookvol10.4 D02AGNT absorb src/input test case files
+20141206 tpd src/input/Makefile remove d02 files
 20141206 tpd src/input/d02bbf.input absorbed, removed
 20141206 tpd src/input/d02bhf.input absorbed, removed
 20141206 tpd src/input/d02cjf.input absorbed, removed
@@ -17,6 +32,7 @@
 20141206 tpd src/input/d02raf.input absorbed, removed
 20141206 tpd src/axiom-website/patches.html 20141206.04.tpd.patch
 20141206 tpd books/bookvol10.4 D01AGNT absorb src/input test case files
+20141206 tpd src/input/Makefile remove d01 files
 21041206 tpd src/input/d01ajf.input absorbed, removed
 21041206 tpd src/input/d01akf.input absorbed, removed
 21041206 tpd src/input/d01alf.input absorbed, removed
@@ -31,6 +47,7 @@
 21041206 tpd src/input/d01gbf.input absorbed, removed
 20141206 tpd src/axiom-website/patches.html 20141206.03.tpd.patch
 20141206 tpd books/bookvol10.4 NAGC06 absorb src/input test case files
+20141206 tpd src/input/Makefile remove c06 files
 20141206 tpd src/input/c06eaf.input absorbed, removed
 20141206 tpd src/input/c06ebf.input absorbed, removed
 20141206 tpd src/input/c06ecf.input absorbed, removed
diff --git a/patch b/patch
index c4cb5c4..b08194c 100644
--- a/patch
+++ b/patch
@@ -1,2 +1,2 @@
-books/bookvol10.4 D03AGNT absorb src/input test case files
+books/bookvol10.4 NAGF01 absorb src/input test case files
 
diff --git a/src/axiom-website/patches.html b/src/axiom-website/patches.html
index 348c766..9f20820 100644
--- a/src/axiom-website/patches.html
+++ b/src/axiom-website/patches.html
@@ -4772,6 +4772,8 @@ books/bookvol10.4 D02AGNT absorb src/input test case files<br/>
 books/bookvol10.4 D02AGNT absorb src/input test case files<br/>
 <a href="patches/20141207.01.tpd.patch">20141207.01.tpd.patch</a>
 books/bookvol10.4 D03AGNT absorb src/input test case files<br/>
+<a href="patches/20141207.02.tpd.patch">20141207.02.tpd.patch</a>
+books/bookvol10.4 NAGF01 absorb src/input test case files<br/>
  </body>
 </html>
 
diff --git a/src/input/Makefile.pamphlet b/src/input/Makefile.pamphlet
index e7891ae..007072f 100644
--- a/src/input/Makefile.pamphlet
+++ b/src/input/Makefile.pamphlet
@@ -198,9 +198,7 @@ SKIP= as-eg1.output as-eg2.output    \
       drawpoly.output  drawx.output     dropt.output     \
       egg.output       \
       explot2d.output  explot3d.output  \
-      f01brf.output    f01bsf.output    f01maf.output    \
-      f01mcf.output    f01qcf.output    f01qdf.output    f01qef.output    \
-      f01rcf.output    f01rdf.output    f01ref.output    f02aaf.output    \
+      f02aaf.output    \
       f02abf.output    f02adf.output    f02aef.output    f02aff.output    \
       f02agf.output    f02ajf.output    f02akf.output    f02awf.output    \
       f02axf.output    f02bbf.output    f02bjf.output    f02fjf.output    \
@@ -626,10 +624,7 @@ ASFILES=${OUT}/aseg6.as   ${OUT}/aseg7.as  ${OUT}/ecfact.as \
         ${OUT}/hilbert.as ${OUT}/matops.as ${OUT}/pdecomp0.as \
         ${OUT}/romnum.as
 
-NAGLIB=${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 \
-       ${OUT}/f01ref.input   ${OUT}/f02aaf.input     ${OUT}/f02abf.input \
+NAGLIB=${OUT}/f02aaf.input   ${OUT}/f02abf.input \
        ${OUT}/f02adf.input   ${OUT}/f02aef.input     ${OUT}/f02aff.input \
        ${OUT}/f02agf.input   ${OUT}/f02ajf.input     ${OUT}/f02akf.input \
        ${OUT}/f02awf.input   ${OUT}/f02axf.input     ${OUT}/f02bbf.input \
@@ -1066,11 +1061,6 @@ DOCFILES= \
   ${DOC}/expr1.input.dvi       ${DOC}/expr.input.dvi       \
   ${DOC}/exprode.input.dvi     ${DOC}/exprpoly.input.dvi   \
   ${DOC}/exseries.input.dvi    ${DOC}/exsum.input.dvi      \
-  ${DOC}/f01brf.input.dvi      ${DOC}/f01bsf.input.dvi     \
-  ${DOC}/f01maf.input.dvi      ${DOC}/f01mcf.input.dvi     \
-  ${DOC}/f01qcf.input.dvi      ${DOC}/f01qdf.input.dvi     \
-  ${DOC}/f01qef.input.dvi      ${DOC}/f01rcf.input.dvi     \
-  ${DOC}/f01rdf.input.dvi      ${DOC}/f01ref.input.dvi     \
   ${DOC}/f02aaf.input.dvi      ${DOC}/f02abf.input.dvi     \
   ${DOC}/f02adf.input.dvi      ${DOC}/f02aef.input.dvi     \
   ${DOC}/f02aff.input.dvi      ${DOC}/f02agf.input.dvi     \
diff --git a/src/input/f01brf.input.pamphlet b/src/input/f01brf.input.pamphlet
deleted file mode 100644
index 1ae3936..0000000
--- a/src/input/f01brf.input.pamphlet
+++ /dev/null
@@ -1,65 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01brf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-n:=6
-nz:=15
-licn:=150
-lirn:=75
-pivot:=0.1
-lblock:=true
-grow:=true
-abort:=[true ,true ,false ,true ]
-a:Matrix SF:=
- [[5 ,2 ,-1 ,2 ,3 ,-2 ,1 ,1 ,-1 ,-1 ,2 ,-3 ,-1 ,-1 ,6_
- ,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 ,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 ,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 ,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 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0_
- ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-irn:Matrix Integer:=
- [[1 ,2 ,2 ,2 ,3 ,4 ,4 ,4 ,5 ,5 ,5 ,5 ,6 ,6 ,6_
-  ,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 ,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 ,0 ,0 ,0 ,0 ]]
-icn:Matrix Integer:=
- [[1 ,2 ,3 ,4 ,3 ,1 ,4 ,5 ,1 ,4 ,5 ,6 ,1 ,2 ,6_
-  ,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 ,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 ,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 ,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 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0_
-  ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-ifail:= 110
-result:=f01brf(n,nz,licn,lirn,pivot,lblock,grow,abort,a,irn,icn,ifail)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01bsf.input.pamphlet b/src/input/f01bsf.input.pamphlet
deleted file mode 100644
index 0f242d9..0000000
--- a/src/input/f01bsf.input.pamphlet
+++ /dev/null
@@ -1,60 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01bsf.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 := 6
-nz := 15
-nzmax := 50
-licn := 3*nzmax
-ta1:List SF := new(licn-nz,0.0);
-ta2:List SF := [5.,2.,-1.,2.,3.,-2.,1.,1.,-1.,-1.,2.,-3.,-1.,-1.,6.] ;
-ta3 := concat(ta2,ta1);
-a :Matrix SF := matrix [ta3] 
-lirn := (3*nzmax/2)::Integer
-tr1:List Integer := new(lirn-nz,0);
-tr2 := [1,2,2,2,3,4,4,4,5,5,5,5,6,6,6];
-tr3 := concat(tr2,tr1);
-irn := matrix[tr3]
-tc1:List Integer := new(licn-nz,0);
-tc2 := [1,2,3,4,3,1,4,5,1,4,5,6,1,2,6];
-tc3 := concat(tc2,tc1);
-icn :=  matrix [tc3] 
-pivot :SF:= 0.1
-lblock := true
-grow := true
-abort := [true,true,false,true]
-ifail := 110
-f01brfResult:=f01brf(n,nz,licn,lirn,pivot,lblock,grow,abort,a,irn,icn,ifail);
-eta :SF:= .1 -- now run f01bsf
-abort := true
-keep := f01brfResult.'ikeep;
-disp := f01brfResult.'idisp;
-cn :=  f01brfResult.'icn;
-ta2:=[10.,12.,-3.,-1.,15.,-2.,10.,-1.,-1.,-5.,1.,-1.,-1.,-2.,6.];
-a:= matrix [concat (ta2,ta1)]
-result:=f01bsf(n,nz,licn,matrix [tr2],matrix [tc2],cn,keep,grow,eta,abort,_
-matrix [[disp(1,1),disp(1,2)]],a,ifail)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01maf.input.pamphlet b/src/input/f01maf.input.pamphlet
deleted file mode 100644
index 8446dcc..0000000
--- a/src/input/f01maf.input.pamphlet
+++ /dev/null
@@ -1,57 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01maf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-n:=16
-nz:=40
-licn:=90
-lirn:=50
-abort:=[true ,true ,true ]
-avals:Matrix SF:= 
- [[1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,1.0_
- ,1.0 ,1.0 ,1.0 ,1.0 ,1.0 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
- ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
- ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25 ,-0.25_
- ,-0.25 ,-0.25 ,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 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0_
- ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-irn:Matrix Integer:=
- [[1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14_
- ,15 ,16 ,1 ,2 ,3 ,5 ,6 ,7 ,9 ,10 ,11 ,13 ,14_
- ,15 ,1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12_
- ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-icn:Matrix Integer:=
- [[1 ,2 ,3 ,4 ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14_
- ,15 ,16 ,2 ,3 ,4 ,6 ,7 ,8 ,10 ,11 ,12 ,14 ,15 ,16_
- ,5 ,6 ,7 ,8 ,9 ,10 ,11 ,12 ,13 ,14 ,15 ,16 ,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 ,0 ,0 ,0 ,0 ,0 ,0_
- ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ]]
-droptl:= 0.1
-densw:= 0.8
-ifail:=111
-result:=f01maf(n,nz,licn,lirn,abort,avals,irn,icn,droptl,densw,ifail)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01mcf.input.pamphlet b/src/input/f01mcf.input.pamphlet
deleted file mode 100644
index 4333a44..0000000
--- a/src/input/f01mcf.input.pamphlet
+++ /dev/null
@@ -1,35 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01mcf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-n:=6
-avals:Matrix SF:=
- [[1.0 ,2.0 ,5.0 ,3.0 ,13.0 ,16.0 ,5.0 ,14.0_
- ,18.0 ,8.0 ,55.0 ,24.0 ,17.0 ,77.0 ]]
-lal:= 14
-nrow:Matrix Integer:= [[1 ,2 ,2 ,1 ,5 ,3 ]]
-result:=f01mcf(n,avals,lal,nrow,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01qcf.input.pamphlet b/src/input/f01qcf.input.pamphlet
deleted file mode 100644
index 9e4b1b2..0000000
--- a/src/input/f01qcf.input.pamphlet
+++ /dev/null
@@ -1,38 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01qcf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-m:=5
-n:=3
-lda:=5
-a:Matrix SF:=
- [[2.0 ,2.5 ,2.5 ],_
-  [2.0 ,2.5 ,2.5 ],_
-  [1.6 ,-0.4 ,2.8 ],_
-  [2.0 ,-0.5 ,0.5 ],_
-  [1.2 ,-0.3 ,-2.9 ]]
-result:=f01qcf(m,n,lda,a,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01qdf.input.pamphlet b/src/input/f01qdf.input.pamphlet
deleted file mode 100644
index e3bf9e8..0000000
--- a/src/input/f01qdf.input.pamphlet
+++ /dev/null
@@ -1,50 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01qdf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-trans:="n"
-wheret:="i"
-m:=5
-n:=3
-a:Matrix SF:=
- [[2.0 ,2.5 ,2.5 ],_
-  [2.0 ,2.5 ,2.5 ],_
-  [1.6 ,-0.4 ,2.8 ],_
-  [2.0 ,-0.5 ,0.5 ],_
-  [1.2 ,-0.3 ,-2.9 ]]
-lda:=5
-zeta:Matrix SF:=
- [[0.0 ,0.0 ,0.0 ]]
-ncolb:=2
-ldb:=5
-b:Matrix SF:=
- [[1.1 ,0.0 ],_
-  [0.9 ,0.0 ],_
-  [0.6 ,1.32 ],_
-  [0.0 ,1.1 ],_
-  [-0.8 ,-0.26 ]]
-result:=f01qdf(trans,wheret,m,n,a,lda,zeta,ncolb,ldb,b,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01qef.input.pamphlet b/src/input/f01qef.input.pamphlet
deleted file mode 100644
index fa90202..0000000
--- a/src/input/f01qef.input.pamphlet
+++ /dev/null
@@ -1,42 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01qef.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-wheret:="s"
-m:=5
-n:=3
-ncolq:=5
-lda:=5
-zeta:Matrix SF:=
- [[1.2247 ,1.1547 ,1.2649 ]]
-a:Matrix SF:=
- [[-4.0 ,-2.0 ,-3.0 ,0.0 ,0.0 ],_
-  [0.4085 ,-3.0 ,-2.0 ,0.0 ,0.0 ],_
-  [0.3266 ,-0.4619 ,-4.0 ,0.0 ,0.0 ],_
-  [0.4082 ,-0.5774 ,0.0 ,0.0 ,0.0 ],_
-  [0.2449 ,-0.3464 ,-0.6326 ,0.0 ,0.0 ]]
-result:=f01qef(wheret,m,n,ncolq,lda,zeta,a,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01rcf.input.pamphlet b/src/input/f01rcf.input.pamphlet
deleted file mode 100644
index 4634c77..0000000
--- a/src/input/f01rcf.input.pamphlet
+++ /dev/null
@@ -1,38 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01rcf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-m:=5
-n:=3
-lda:=5
-a:Matrix Complex SF:=
- [[complex(0,0.5) ,complex(-0.5,1.5),complex(-1.0,1.0)],_
-  [complex(0.4,0.3),complex(0.9,1.3),complex(0.2,1.4) ],_
-  [complex(0.4,0), complex(-0.4,0.4),complex(1.8,0)],_
-  [complex(0.3,-0.4),complex(0.1,0.7),complex(0,0)],_
-  [complex(0,-0.3),complex(0.3,0.3),complex(0,2.4) ]]
-result:=f01rcf(m,n,lda,a,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01rdf.input.pamphlet b/src/input/f01rdf.input.pamphlet
deleted file mode 100644
index 5594d04..0000000
--- a/src/input/f01rdf.input.pamphlet
+++ /dev/null
@@ -1,50 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f001rdf.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-trans:="n"
-wheret:="i"
-m:=5
-n:=3
-a:Matrix Complex SF:=
- [[      0.5*%i ,-0.5 + 1.5*%i ,-1.0 + 1.0*%i ],_
-  [0.4 + 0.3*%i , 0.9 + 1.3*%i , 0.2 + 1.4*%i ],_
-  [0.4          ,-0.4 + 0.4*%i , 1.8 ],_
-  [0.3 - 0.4*%i , 0.1 + 0.7*%i , 0.0 ],_
-  [     -0.3*%i , 0.3 + 0.3*%i , 2.4 ]]
-lda:=5
-theta:Matrix Complex SF:=
- [[0.0 ,0.0 ,0.0 ]]
-ncolb:=2
-ldb:=5
-b:Matrix Complex SF :=
- [[-0.55 + 1.05*%i , 0.45 + 1.05*%i ],_
-  [ 0.49 + 0.93*%i , 1.09 + 0.13*%i ],_
-  [ 0.56 - 0.16*%i , 0.64 + 0.16*%i ],_
-  [ 0.39 + 0.23*%i ,-0.39 - 0.23*%i ],_
-  [ 1.13 + 0.83*%i ,-1.13 + 0.77*%i ]]
-result:=f01rdf(trans,wheret,m,n,a,lda,theta,ncolb,ldb,b,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}
diff --git a/src/input/f01ref.input.pamphlet b/src/input/f01ref.input.pamphlet
deleted file mode 100644
index e59a3db..0000000
--- a/src/input/f01ref.input.pamphlet
+++ /dev/null
@@ -1,42 +0,0 @@
-\documentclass{article}
-\usepackage{axiom}
-\setlength{\textwidth}{400pt}
-\begin{document}
-\title{\$SPAD/src/input f01ref.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}
-)cl all
-showArrayValues true
-showScalarValues true
-
-wheret:="s"
-m:=5
-n:=3
-ncolq:=2
-lda:=5
-theta:Matrix Complex SF:=
- [[1 + 0.5*%i ,1.0954-0.3333*%i ,1.2649-1.1565*%i ]]
-a:Matrix Complex SF :=
- [[ 1             , 1      +1     *%i , 1     + %i ],_
-  [-0.2  - 0.4*%i ,-2                 ,-1     - %i ],_
-  [-0.32 - 0.16*%i,-0.3505 +0.2629*%i ,-3          ],_
-  [-0.4  + 0.2*%i ,         0.5477*%i , 0.0        ],_
-  [-0.12 + 0.24*%i, 0.1972 +0.2629*%i , 0.6325     ]]
-result:=f01ref(wheret,m,n,ncolq,lda,theta,a,-1)
-\end{chunk}
-\eject
-\begin{thebibliography}{99}
-\bibitem{1} nothing
-\end{thebibliography}
-\end{document}