dsolve(diff(y(x),x) = ((x^2+1)^(3/2)*x^2+(x^2+1)^(3/2)+y(x)^2*(x^2+1)^(3/2)+x^2*y(x)^3+y(x)^3)*x/(x^2+1)^3,y(x))