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