[SciPy-dev] Numerical Recipes (was tagging 0.7rc1 this weekend?)

[Alan G Isaac]

On 12/16/2008 1:59 AM Jarrod Millman apparently wrote:
> 3.  scipy/special/orthogonal.py
> 
> Unresolved.


>> The coefficients can be found as:
>>
>> A_n = -a2n / a3n
>>
>> B_n = ( a4n / a3n sqrt(h_n-1 / h_n))**2
>>
>>     where
>>             h_n = int_a^b w(x) f_n(x)^2
>> assume:
>> P_0(x) = 1
>> P_-1(x) == 0
>>
>> See Numerical Recipies in C, page 156 and
>> Abramowitz and Stegun p. 774, 782



Travis O. is listed as writing this orthogonal.py in 2000
and doing bug fixes in 2003.  Note that the
reference to Numerical Recipes is one of two
references documenting an algebraic recursion
relation and is not referring to an implementation.

http://www.fizyka.umk.pl/nrbook/c4-5.pdf

http://www-lsp.ujf-grenoble.fr/recherche/a3t2/a3t2a2/bahram/biblio/abramowitz_and_stegun/page_774.htm

I do not see any problem with the code.
I think this one should also be considered resolved.

Alan