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