Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project! See Demo

Close

#2471 partfrac might need to gfactor ?

open
nobody
None
5
2012-11-18
2012-09-20
Barton Willis
No

(%i35) partfrac(1/((z+2*%i*w)*(z^2+4*w^2)),z);
expt: undefined: 0 to a negative exponent.
-- an error. To debug this try: debugmode(true);

A workaround is to map gfactor on the the expression before computing the pfd:

(%i41) partfrac(scanmap('gfactor, 1/((z+2*%i*w)*(z^2+4*w^2))),z);
(%o41) 1/(16*w^2*(z+2*%i*w))+%i/(4*w*(z+2*%i*w)^2)-1/(16*w^2*(z-2*%i*w))

Discussion