From 307e4ad0a2676ef94b42121e79e71888f4778883 Mon Sep 17 00:00:00 2001 From: dos-reis Date: Thu, 27 May 2010 22:31:23 +0000 Subject: cleanup --- src/input/tutchap67.input.pamphlet | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) (limited to 'src/input/tutchap67.input.pamphlet') 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 -- cgit v1.2.3