In[]:=
ClearAll
["`*"];
In[]:=
hilb[l_,th_]:=
1/2
(th/Sin[th])
*BesselJ[0,(l+1/2)*th];
In[]:=
th=π/7;​​hilb[1000,th]//N​​​​
Out[]=
-0.0203653
In[]:=
Plot[LegendreP[l,th]//N,{l,0,1000}]
Out[]=
200
400
600
800
1000
-0.10
-0.05
0.05
0.10
​
In[]:=
LegendreP[1000,th]//N
General
:-
24272717090968164314746593613950077377666444318575025881796960704551665021551621974598214051672375
17198153866121549618259314463156589972989898676883651386288636902717108677272125334352538871267328
is too small to represent as a normalized machine number; precision may be lost.
General
:
521979410233208052427947069829530667633454389637574240636209030823525725315770998172986373695057875
337083815775982372517882563477869163470602014066884671712572832932553300745336565533097618768396288
is too small to represent as a normalized machine number; precision may be lost.
General
:-
86458888716007267027620404644199146431252024983573182725000248470995358120187059796684817561821725
51615959290697300791800767532548715656435933403885997309877150427972241766296615972555728739106816
is too small to represent as a normalized machine number; precision may be lost.
General
:Further output of General::munfl will be suppressed during this calculation.
Out[]=
-3.86639×
170
10