%auto %default_mode maxima
partfrac ( 1/(x^2*(x^2 + 1)), x);
f : integrate(x/(x^3 + 1), x);
g : diff(f, x)
ratsimp(g)