diff options
author | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
---|---|---|
committer | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
commit | 001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532 (patch) | |
tree | da9e2fe5d81ff4cd7709d12e44b8c3e348b8a8e3 /src/algebra/rderf.spad.pamphlet | |
parent | a7bab9a6c2070d05e2dbd256ce455079c8ced385 (diff) | |
download | open-axiom-001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532.tar.gz |
Replace `^=' with `~='.
Diffstat (limited to 'src/algebra/rderf.spad.pamphlet')
-rw-r--r-- | src/algebra/rderf.spad.pamphlet | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/algebra/rderf.spad.pamphlet b/src/algebra/rderf.spad.pamphlet index 33dcd1dd..f8d8d0b3 100644 --- a/src/algebra/rderf.spad.pamphlet +++ b/src/algebra/rderf.spad.pamphlet @@ -117,7 +117,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where q:UP := 0 db := (degree bb)::Z lb := leadingCoefficient bb - while cc ^= 0 repeat + while cc ~= 0 repeat d < 0 or (n := (degree cc)::Z - db) < 0 or n > d => return [q, true] r := monomial((leadingCoefficient cc) / lb, n::N) cc := cc - bb * r - derivation r @@ -130,7 +130,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where -- dtm1 = degree(Dt) - 1 SPDEnocancel2(bb, cc, d, dtm1, lt, derivation) == q:UP := 0 - while cc ^= 0 repeat + while cc ~= 0 repeat d < 0 or (n := (degree cc)::Z - dtm1) < 0 or n > d => return [[q, true]] if n > 0 then r := monomial((leadingCoefficient cc) / (n * lt), n::N) @@ -139,7 +139,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where q := q + r else -- n = 0 so solution must have degree 0 db:N := (zero? bb => 0; degree bb); - db ^= degree(cc) => return [[q, true]] + db ~= degree(cc) => return [[q, true]] zero? db => return [[bb, cc, 0, 1, q]] r := leadingCoefficient(cc) / leadingCoefficient(bb) cc := cc - r * bb - derivation(r::UP) @@ -162,7 +162,7 @@ TranscendentalRischDE(F, UP): Exports == Implementation where v := polyRDE(u.a, bb, cc, n, differentiate).ans [v.ans / u.t, v.nosol] --- return an a bound on the degree of a solution of A P'+ B P = C,A ^= 0 +-- return an a bound on the degree of a solution of A P'+ B P = C,A ~= 0 -- cancellation at infinity is possible -- base case: F' = 0 getBound(a, b, dc) == |