Hi,<br><br>There is an error in ginac when taking a series of a function at<br>a singular point, for instance in ginsh,<br><br>> series(sqrt(x),x,3);<br>power::eval(): division by zero<br><br>from ginac-1.6 file power.cpp line 408. <br>
<br>I instrumented power.cpp with some printf functions and found<br>the offending dispatch to be in line 665:<br><br>         coeff = deriv.subs(r, subs_options::no_pattern);<br><br>while evaluating the third term in the expansion, presumably with<br>
a 1/x factor.  I'm new to ginac so its taking me a while to learn the<br>code -- next step is to try and get it running in gdb -- so I'd appreciate<br>any input from the list about how ginac is handling and/or is <br>
supposed to handle this kind of error.  Other software, like maxima,<br>handles the problem fine:<br><br>(%o6)                      maxima(sqrt(x), x, 0, 3)<br>(%i7) taylor(sqrt(x),x,0,3);<br><br>Thanks in advance,<br><br>
John<br><br><i><br></i><br>