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