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