normalized ALF (Assotiated Legendre Function)

[Ben Abbott]

On Tuesday, February 12, 2008, at 10:02AM, "Marco Caliari" <marco.caliari at univr.it> 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")(end-1:end)
>ans =
>
>   -3.3248e-53
>           Inf
>octave:2> legendre(151,-0.9,"norm")(end-1:end)
>ans =
>
>   -3.3248e-53
>    9.2660e-55
>
>Marco
>

I'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

-------------- next part --------------
A non-text attachment was scrubbed...
Name: legendre.m
Type: application/octet-stream
Size: 5370 bytes
Desc: not available
Url : https://www.cae.wisc.edu/pipermail/help-octave/attachments/20080212/d5d83577/attachment.obj