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