dsolve(diff(y(x),x)-a3*y(x)^3-a2*y(x)^2-a1*y(x)-a0 = 0,y(x))