Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

This looks like a bug in giac. The integrand that Sage passes to giac is

(b^2*x^2 + a*b*x*2 + a^2)^(-3/2) * x^4

which is equivalent to the one you defined. This however crashes in giac itself:

0>> integrate((b^2*x^2 + a*b*x*2 + a^2)^(-3/2) * x^4, x)
Warning, integration of abs or sign assumes constant sign by intervals (correct if the argument is real):
Check [abs(b*t_nostep+a)]
Discontinuities at zeroes of b*t_nostep+a were not checked
Segmentation fault (core dumped)

(tested with giac 1.7.0)