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