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

Alan G Isaac aisaac@american....
Tue Dec 16 09:34:36 CST 2008


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



More information about the Scipy-dev mailing list