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