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