[GiNaC-devel] series() errors with divide by zero
J Sage
jlhsage at gmail.com
Mon Jun 27 07:09:46 CEST 2011
...
> series(sqrt(x),x,0,3)
> = series(sqrt(y^2),y,0,3*2)
> = y
> = x^(1/2).
>
> Note that the member functions degree() and ldegree() currently return int,
> so this would have to be generalized somehow, when implementing Puiseux
> series directly in GiNaC.
>
> Bye
> -richy.
> --
> Richard B. Kreckel
> <http://www.ginac.de/~kreckel/>
>
Actually, ginac cannot do a series on x^n or on sqrt(x^2). Here's the ginsh
output:
> series(x^n,x,1);
(0^n)+Order(x)
> series(x^n,x,2);
power::eval(): division by zero
> series(sqrt(x^2),x,3);
power::eval(): division by zero
>
John
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.cebix.net/pipermail/ginac-devel/attachments/20110626/843e31b4/attachment.html>
More information about the GiNaC-devel
mailing list