dsolve(diff(y(x),x) = (x^3*y(x)^3+6*x^2*y(x)^2+12*x*y(x)+8+2*x)/x^3,y(x))