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