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