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