diff options
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) == |