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