DSolve[Derivative[2][y][x] == -(((2*a*x^2 + n*(1 + n)*(-1 + x^2) + (a - n)*(1 + a + n)*x^2*(-1 + x^2))*y[x])/(x^2*(-1 + x^2))) + (2*x*Derivative[1][y][x])/(-1 + x^2),y[x],x]