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