From 701073966d9aba091a99aa3cbef48fe3df226dfe Mon Sep 17 00:00:00 2001 From: dos-reis Date: Sun, 13 Jun 2010 16:37:27 +0000 Subject: * algebra/array1.spad.pamphlet (IndexedFlexibleArray) [merge!]: Tidy. * algebra/intfact.spad.pamphlet (IntegerRoots) [perfectNthRoot]: Likewise. * algebra/numode.spad.pamphlet (NumericalOrdinaryDifferentialEquations) [rf4a]: Likewise. * algebra/solvelin.spad.pamphlet (LinearSystemMatrixPackage) [aSolution]: Likewise. --- src/algebra/numtheor.spad.pamphlet | 11 ++++++++--- 1 file changed, 8 insertions(+), 3 deletions(-) (limited to 'src/algebra/numtheor.spad.pamphlet') diff --git a/src/algebra/numtheor.spad.pamphlet b/src/algebra/numtheor.spad.pamphlet index e3ed3b72..a4174263 100644 --- a/src/algebra/numtheor.spad.pamphlet +++ b/src/algebra/numtheor.spad.pamphlet @@ -405,8 +405,10 @@ IntegerNumberTheoryFunctions(): Exports == Implementation where even? b and positiveRemainder(a,4) > 1 => error "J(a/b) not defined for b even and a = 2 or 3 (mod 4)" even? b and even? a => 0 - k : Integer - for k in 0.. while even? b repeat b := b quo 2 + k : Integer := 0 + while even? b repeat + b := b quo 2 + k := k + 1 j:I := (odd? k and positiveRemainder(a,8) = 5 => -1; 1) b = 1 => j a := positiveRemainder(a,b) @@ -418,7 +420,10 @@ IntegerNumberTheoryFunctions(): Exports == Implementation where (a,b) := (b rem a,a) else -- J(2*a/b) = J(a/b) (-1) (b**2-1)/8 - for k in 0.. until odd? a repeat a := a quo 2 + k := 0 + until odd? a repeat + a := a quo 2 + k := k + 1 if odd? k and (b+2) rem 8 > 4 then j := -j a = 0 => 0 j -- cgit v1.2.3