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