<br><br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">...<br>
  series(sqrt(x),x,0,3)<br>
    = series(sqrt(y^2),y,0,3*2)<br>
    = y<br>
    = x^(1/2).<br>
<br>
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.<br>
<br>
Bye<br>
  -richy.<br><font color="#888888">
-- <br>
Richard B. Kreckel<br>
<<a href="http://www.ginac.de/%7Ekreckel/" target="_blank">http://www.ginac.de/~kreckel/</a>></font><br></blockquote></div><br>Actually, ginac cannot do a series on x^n or on sqrt(x^2).  Here's the ginsh output:<br>
<br>> series(x^n,x,1);<br>(0^n)+Order(x)<br>> series(x^n,x,2);<br>power::eval(): division by zero<br>> series(sqrt(x^2),x,3);<br>power::eval(): division by zero<br>> <br><br>John<br>