diff options
| author | dos-reis <gdr@axiomatics.org> | 2010-06-13 16:37:27 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2010-06-13 16:37:27 +0000 |
| commit | 701073966d9aba091a99aa3cbef48fe3df226dfe (patch) | |
| tree | ed595066e9e5e307fa4b86c9eb3852145965ac3f /src/algebra/numtheor.spad.pamphlet | |
| parent | 9a4759f38f13d4bf79133846eb5fd2aff58b2235 (diff) | |
| download | open-axiom-701073966d9aba091a99aa3cbef48fe3df226dfe.tar.gz | |
* 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.
Diffstat (limited to 'src/algebra/numtheor.spad.pamphlet')
| -rw-r--r-- | src/algebra/numtheor.spad.pamphlet | 11 |
1 files changed, 8 insertions, 3 deletions
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 |