dsolve(diff(diff(y(x),x),x) = -2/x*(x^2-1)/(x-1)^2*diff(y(x),x)-(-2*x^2+2*x+2)/x^2/(x-1)^2*y(x),y(x))