Numerical solution for quasinormal modes for potentials with exponential tails

The numerical determination of quasinormal mode frequencies is notoriously difficult for long-ranged potentials. Of particular interest are potentials with exponential tails, for which direct numerical integration fails. The reason for such a failure is analysed, and an effective numerical algorithm is developed.