dsolve(diff(diff(y(x),x),x) = -b*x/(x^2-1)/a*diff(y(x),x)-(c*x^2+d*x+e)/a/(x^2-1)^2*y(x),y(x))