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