A new Fortran 90 program to compute regular and irregular associated Legendre functions Original Research Article

We present a modern Fortran 90 code to compute the regular image and irregular image associated Legendre functions for all image (on the cut) and image and integer degree (l) and order (m). The code applies either forward or backward recursion in (l) and (m) in the stable direction, starting with analytically known values for forward recursion and considering both a Wronskian based and a modified Millerʹs method for backward recursion. While some Fortran 77 codes existed for computing the functions off the cut, no Fortran 90 code was available for accurately computing the functions for all real values of x different from image where the irregular functions are not defined.