dsolve((g1(x)*y(x)+g0(x))*diff(y(x),x)-f1(x)*y(x)-f2(x)*y(x)^2-f3(x)*y(x)^3-f0(x) = 0,y(x))