A New High-Performing Method for Solving the Homogeneous Teukolsky Equation

The numerical waveforms for the extreme mass-ratio inspirals (EMRIs) require a huge amount of homogeneous solutions of the Teukolsky equation in the frequency domain. The calculation accuracy and efficiency of the homogeneous solutions are the key performance bottleneck in waveform generation. In this paper, we propose a new numerical method based on the analytical series expansion which is most efficient for computing the homogeneous solutions with very high accuracy and a wider frequency range compared with the existing methods. Our new method is definitely useful for constructing the waveform templates of EMRIs.