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