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