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