``` from sympy import * f = x**-3 * (1 + x**4)**-Rational(1, 2) F = integrate(f, x, meijerg=True) print(F) ``` gives the result `-sqrt(1 + x**(-4))/2` which does not equal to the integrand when differentiated ``` print(F.diff(x).xreplace({x: I})) print(f.xreplace({x: I})) ``` ``` -sqrt(2)*I/2 sqrt(2)*I/2 ```