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