diff options
Diffstat (limited to 'src/algebra/fparfrac.spad.pamphlet')
| -rw-r--r-- | src/algebra/fparfrac.spad.pamphlet | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/algebra/fparfrac.spad.pamphlet b/src/algebra/fparfrac.spad.pamphlet index b769ba05..f6495e21 100644 --- a/src/algebra/fparfrac.spad.pamphlet +++ b/src/algebra/fparfrac.spad.pamphlet @@ -141,7 +141,7 @@ FullPartialFractionExpansion(F, UP): Exports == Implementation where lvar := concat(makeVariable(u, m), lvar) lval := concat(inv((m+1)::F) * q, lval) qq := q0 quo gcd(pp, q0) -- new center - if (degree(qq) > 0) then + if positive? degree(qq) then ans := concat([(n + 1 - m)::N, qq, (pp * bm * cn * cm) rem qq], ans) cm := (c * cm) rem q0 -- cm = c**m modulo q now um := u1 * um -- um = u**m now |