MathGroup Archive 2007

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: LegendreP error (bug?) in Mathematica


John Ralston wrote:

> LegendreP[ l, m] and SphericalHarmonicY[ t, p, l, m]  go
> wrong for large index l .  
> 
> For l> 40 or so neither can be used reliably everywhere. 
> 
> To see the breakdown plot the functions.  Not all index cases fail.
> Here's some examples: 
> 
> Plot[LegendreP[ 44, 13, x] , {x, -1, 1} ]
> 
> Plot[LegendreP[ 66, 9, x] , {x, -1, 1} ]
> 
> ListPlot[ Table[ {j,
>      Sqrt[ 4Pi/ (
>       2l + 1) ] Abs[ SphericalHarmonicY[ j, 0, 0,
>                      3.7 ]]}, {j, 1, 55}], PlotJoined -> True]
> 
> Has anyone fixed this?  Does anyone care? 
> I need l ranges above 200. 

<snip>

Seems to work like a charm with Mathematica 6.0.1. See

http://homepages.nyu.edu/~jmg336/mathematica/LegendreP601.pdf

-- 
Jean-Marc


  • Prev by Date: Re: FactorInteger slow in V6.0 vs V5.2
  • Next by Date: Re: LegendreP error (bug?) in Mathematica
  • Previous by thread: Re: LegendreP error (bug?) in Mathematica
  • Next by thread: Re: LegendreP error (bug?) in Mathematica