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