diff options
| author | dos-reis <gdr@axiomatics.org> | 2010-05-27 22:31:23 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2010-05-27 22:31:23 +0000 |
| commit | 307e4ad0a2676ef94b42121e79e71888f4778883 (patch) | |
| tree | 15cc67e1c57f7b0efd639bd14472bf3e559c08e3 /src/input/tutchap67.input.pamphlet | |
| parent | 8f5399e77de0741651b0d058b1aa7f57f43aedb0 (diff) | |
| download | open-axiom-307e4ad0a2676ef94b42121e79e71888f4778883.tar.gz | |
cleanup
Diffstat (limited to 'src/input/tutchap67.input.pamphlet')
| -rw-r--r-- | src/input/tutchap67.input.pamphlet | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/src/input/tutchap67.input.pamphlet b/src/input/tutchap67.input.pamphlet index f730cf8f..1d5c1f56 100644 --- a/src/input/tutchap67.input.pamphlet +++ b/src/input/tutchap67.input.pamphlet @@ -127,8 +127,8 @@ solve(matD,[5,6,7,9]) hilbert3 :: Matrix DoubleFloat -- continuing the previous session % * inverse % matrix [[1/(i+j) for i in 1..11] for j in 1..11]::Matrix DoubleFloat; -badUnit := % * inverse %; -diagEls := set [%(i,i) for i in 1..11]; +badUnit := % * inverse % +diagEls := set [%(i,i) for i in 1..11] min diagEls max diagEls offDiags := empty()$Set DoubleFloat @@ -137,14 +137,14 @@ for i in 1..11 repeat _ offDiags := union(offDiags,badUnit(i,j)) min offDiags max offDiags -hilbert11 := matrix [[1/(i+j) for i in 1..11] for j in 1..11]; +hilbert11 := matrix [[1/(i+j) for i in 1..11] for j in 1..11] % * inverse % detHilbert3 := determinant hilbert3 detHilbert11 := determinant hilbert11 % :: DoubleFloat determinant(hilbert11::Matrix DoubleFloat) -test3 := hilbert3 :: Matrix Polynomial Fraction Integer; -test3(1,1) := (1 + eps)/2; +test3 := hilbert3 :: Matrix Polynomial Fraction Integer +test3(1,1) := (1 + eps)/2 determinant test3 (% - detHilbert3)/detHilbert3 for i in 1..3 repeat for j in 1..3 repeat _ @@ -152,8 +152,8 @@ for i in 1..3 repeat for j in 1..3 repeat _ test3 (determinant test3 - detHilbert3)/detHilbert3 error3 := matrix [[eps[i,j] for i in 1..3] for j in 1..3] -test3 := hilbert3 + t*error3; -detErr := (determinant test3 - detHilbert3)/detHilbert3; +test3 := hilbert3 + t*error3 +detErr := (determinant test3 - detHilbert3)/detHilbert3 detErrReduced := coefficient(%,'t,1) coefficient(detErr,'t,0) epses := variables detErrReduced |