27 Jun
2011
27 Jun
'11
7:09 a.m.
...
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/>
series(x^n,x,1); (0^n)+Order(x) series(x^n,x,2);
series(sqrt(x^2),x,3);
Actually, ginac cannot do a series on x^n or on sqrt(x^2). Here's the ginsh output: power::eval(): division by zero power::eval(): division by zero
John