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