
From:  Ben Abbott 
Subject:  Re: normalized ALF (Assotiated Legendre Function) 
Date:  Tue, 12 Feb 2008 20:42:35 0500 
The method I used detects an underflow for orders greater than 151. octave:104> result_octave1 = legendre (152, [1:0.1:1]); warning: legendre: results may be unstable for high orders. I've attached the script. Ben
legendre.m
Description: Binary data
On Feb 12, 2008, at 1:38 PM, Ben Abbott wrote:
On Tuesday, February 12, 2008, at 10:02AM, "Marco Caliari" <address@hidden > wrote:Hi.The normalized Lagrange functions should (almost) never give NaNs of Infs.The enclosed script fixes a problem in the previous: octave:1> legendreold(151,0.9,"norm")(end1:end) ans = 3.3248e53 Inf octave:2> legendre(151,0.9,"norm")(end1:end) ans = 3.3248e53 9.2660e55 MarcoI've tried to combine the various implementations. I still need to look into adding a warning in then event of an overflow.We should also add more tests to include the current improvements. Ben <legendre.m>_______________________________________________ Helpoctave mailing list address@hidden https://www.cae.wisc.edu/mailman/listinfo/helpoctave
[Prev in Thread]  Current Thread  [Next in Thread] 