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