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