diff options
Diffstat (limited to 'src/algebra/irexpand.spad.pamphlet')
| -rw-r--r-- | src/algebra/irexpand.spad.pamphlet | 6 |
1 files changed, 2 insertions, 4 deletions
diff --git a/src/algebra/irexpand.spad.pamphlet b/src/algebra/irexpand.spad.pamphlet index c01605cb..47fc84a2 100644 --- a/src/algebra/irexpand.spad.pamphlet +++ b/src/algebra/irexpand.spad.pamphlet @@ -152,8 +152,7 @@ IntegrationResultToFunction(R, F): Exports == Implementation where lg2func lg == zero?(d := degree(p := lg.coeff)) => error "poly has degree 0" --- one? d => [linear(p, lg.logand)] - (d = 1) => [linear(p, lg.logand)] + one? d => [linear(p, lg.logand)] d = 2 => quadratic(p, lg.logand) odd? d and ((r := retractIfCan(reductum p)@Union(F, "failed")) case F) => @@ -213,8 +212,7 @@ IntegrationResultToFunction(R, F): Exports == Implementation where -- s = 1 if b > 0, -1 if b < 0, 0 if the sign of b cannot be determined insqrt y == rec := froot(y, 2)$PolynomialRoots(IndexedExponents K, K, R, P, F) --- one?(rec.exponent) => [[rec.coef * rec.radicand, 1], 1] - ((rec.exponent) = 1) => [[rec.coef * rec.radicand, 1], 1] + one?(rec.exponent) => [[rec.coef * rec.radicand, 1], 1] rec.exponent ~=2 => error "Should not happen" [[rec.coef, rec.radicand], ((s := sign(rec.radicand)) case "failed" => 0; s::Z)] |