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