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