diff options
Diffstat (limited to 'src/algebra/fparfrac.spad.pamphlet')
| -rw-r--r-- | src/algebra/fparfrac.spad.pamphlet | 9 |
1 files changed, 3 insertions, 6 deletions
diff --git a/src/algebra/fparfrac.spad.pamphlet b/src/algebra/fparfrac.spad.pamphlet index 9afe1d78..7d7c48a1 100644 --- a/src/algebra/fparfrac.spad.pamphlet +++ b/src/algebra/fparfrac.spad.pamphlet @@ -106,8 +106,7 @@ FullPartialFractionExpansion(F, UP): Exports == Implementation where -- returns Trace_k^k(a) (h(a) / (x - a)^n) where d(a) = 0 REC2RF(d, h, n) == --- one?(m := degree d) => - ((m := degree d) = 1) => + one?(m := degree d) => a := - (leadingCoefficient reductum d) / (leadingCoefficient d) h(a)::UP / (x - a::UP)**n dd := UP2SUP d @@ -174,16 +173,14 @@ FullPartialFractionExpansion(F, UP): Exports == Implementation where ans output(n, d, h) == --- one? degree d => - (degree d) = 1 => + one? degree d => a := - leadingCoefficient(reductum d) / leadingCoefficient(d) h(a)::O / outputexp((x - a::UP)::O, n) sum(outputForm(makeSUP h, alpha) / outputexp(xx - alpha, n), outputForm(makeSUP d, alpha) = zr) outputexp(f, n) == --- one? n => f - (n = 1) => f + one? n => f f ** (n::O) @ |