diff options
Diffstat (limited to 'src/algebra/intaf.spad.pamphlet')
| -rw-r--r-- | src/algebra/intaf.spad.pamphlet | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/algebra/intaf.spad.pamphlet b/src/algebra/intaf.spad.pamphlet index 23f73b17..8a4e2b94 100644 --- a/src/algebra/intaf.spad.pamphlet +++ b/src/algebra/intaf.spad.pamphlet @@ -389,7 +389,7 @@ PureAlgebraicIntegration(R, F, L): Exports == Implementation where linearInXIfCan(x, y) == a := b := 0$UP p := clearDenominator lift(minPoly y, x) - while p ^= 0 repeat + while p ~= 0 repeat degree(q := numer leadingCoefficient p) > 1 => return "failed" a := a + monomial(coefficient(q, 1), d := degree p) b := b - monomial(coefficient(q, 0), d) @@ -401,7 +401,7 @@ PureAlgebraicIntegration(R, F, L): Exports == Implementation where prootintegrate(f, x, y) == modulus := lift(p := minPoly y, x) rf := reductum(ff := univariate(f, x, y, p)) - ((r := retractIfCan(rf)@Union(RF,"failed")) case RF) and rf ^= 0 => + ((r := retractIfCan(rf)@Union(RF,"failed")) case RF) and rf ~= 0 => -- in this case, ff := lc(ff) y^i + r so we integrate both terms -- separately to gain time map(#1(x::F), integrate(r::RF)) + @@ -481,7 +481,7 @@ PureAlgebraicIntegration(R, F, L): Exports == Implementation where -- returns either "failed" or r(u, z) chvarIfCan(p, d, u, u1) == ans:UPUP := 0 - while p ^= 0 repeat + while p ~= 0 repeat (v := composite(u1 * leadingCoefficient(p) / d ** degree(p), u)) case "failed" => return "failed" ans := ans + monomial(v::RF, degree p) |