I would like to add that, since I need to evaluate this function very close to theta=n*pi, I cannot just use the result from a Legendre polynomial routine and divide it by sin(theta)^m.
Il giorno venerdì 21 giugno 2013 20:25:51 UTC+2, Guido Walter Pettinari ha scritto: > Dear all, > > > > I would be grateful if somebody could point me to a numerical routine (possibly in C) to compute the following special function: > > > > P_lm (cos(theta)) / (sin(theta)^m) , > > > > where P_lm is an associated Legendre polynomial. The function can be also expresses as > > > > P_lm (x) / (1-x^2)^(m/2) . > > > > Since (1-x^2)^(m/2) is a prefactor in P_lm, I guess that the ratio is even easier to estimate than the P_lm's. I checked in GSL and alglib but alas I could not find this function. > > > > Thank you for your attention, > > Guido