[SciPy-dev] Numerical Recipes (was tagging 0.7rc1 this weekend?)
Alan G Isaac
Tue Dec 16 09:34:36 CST 2008
On 12/16/2008 1:59 AM Jarrod Millman apparently wrote:
> 3. scipy/special/orthogonal.py
>> The coefficients can be found as:
>> A_n = -a2n / a3n
>> B_n = ( a4n / a3n sqrt(h_n-1 / h_n))**2
>> h_n = int_a^b w(x) f_n(x)^2
>> 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.
I do not see any problem with the code.
I think this one should also be considered resolved.
More information about the Scipy-dev