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