o l h2@sFddlmZmZeddZeddZededdZeded d Zeded d Zed dZeddZ eddZ eddZ eddZ eddZ eddZeddZddZedd Zed!d"Zed#d$Zed%d&Zed'd(Zed)d*Zed+d,Zed-d.Zd/d0Zed1d2Zed3d4Zed5d6Zed7d8Zd9d:Zd;d<Zeded=d>Z eded?d@Z!dfdBdCZ"ededDdEZ#edgdFdGZ$dHdIZ%edJdKZ&edLdMZ'eifdNdOZ(edhdQdRZ)eifdSdTZ*edhdUdVZ+dWdXZ,dYdZZ-d[d\Z.d]ifd^d_Z/eded`daZ0ededbdcZ1ddS)i)defun defun_wrappedcC |d|S)zCComputes the Bessel function `J_0(x)`. See :func:`~mpmath.besselj`.besseljctxxr k/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/mpmath/functions/bessel.pyj0 r cCr)zDComputes the Bessel function `J_1(x)`. See :func:`~mpmath.besselj`.rrrr r r j1rrrc  s,ttur d}n}|rt|r2dkr2dj |fi|S|r|rdkrtj}zjd7_fddt dD}W|_n|_w| d 9}|Sfd d }j |gfi|}|S|s|rt d krt d krz WStyYnwsڈsƈj}nKdkr҈}n?j}n7j}z/jtd t j7_jdddfdd }j |gfi|}W|_n|_w| }|S)NTrc3s<|]}d||d|VqdS)rN)binomialr).0k)r dnzr r %s4zbesselj..rrcsjjjdddd}d||dd||dg}djg|d|d||gg||dd|ddg||dg|fg}|S)NprecпTexact?rrfmulrpirrrBTMr rr r h+s" Lzbesselj..h rrcsJjjtdjddd}g|gg|dgg|dg|fgS)NrrTrr)fnegr!maxr)rr$)r(r wr r r)Gs&$)typeintconvertisint_rermagrfsumrangempf hypercombabs_besseljNotImplementedErroronereinfminr!) r rr derivativekwargsn_isintorigvr)r )r(r rrr/rr r s^          %   rc  s|}sD|rt|sd|S|r#d|S|}|dkr3j|S|dkr=d|Sj|S|rf|}fdd}j|||gfi|}|Sfdd}j||gfi|}|S)Nrrcsjjjdddd}d||dd||d|dg}djg|d|d||g|dg||dd|ddg||fg}|S)Nr?Trrrrr r#r'r r r)fs"&Hzbesseli..hcsPjddd}j||tdjd}|g|gg|dgg|dg|fgS)NrTrrrr)r!r.r)rr/r$r'r r r)ms$)r2 ValueErrorr3r>nanr?r5r9) r rrrArBr$rr)rEr r'r besseliPs0         rIc KsP|sW|rt|s|j ||S||r|j||S||}|d}||r<|dkr6|j ||Sd||S|dkrPt||drP|j||S|j||S|j d7_ | |\}}||j kr{|j } |j d9_ || 7}n |dkr|j |8_ | |\} } |j |||fi|| |j | ||fi|| S)Nrrrr*)rGr?imrHr>r3r1floorninfr nint_distanceeps cospi_sinpir) r rrrArBr$qmrr)cossinr r r besselyts<      rTc s\sjS}|dkrfdd}nj|7_fdd}j||gfi|S)Nrcsfdd}dg| |dg|gggd|g|f}dg|| dg| gggd|g|f}||fS)Nrrr )rr$T1T2rr r r)s (*zbesselk..hc s@jd ggdgg|dd|ggddfgS)Nr)rrrr)r"exprr rr r r)s )r?r5rr9)r rrrBr(r)r r[r besselks  r\cKs.|j||fi||j|j||fi|SNrjrTr rr rBr r r hankel1.racKs.|j||fi||j|j||fi|Sr]r^r`r r r hankel2rbrccKs|dkr||dkr |S||dkr|j|S|j|S|jd|dd}d|}|||||j||dd||fi|S)NrrXTrrrr)r>r?rHr!rYhyp1f1)r rrQrrBr yr r r whitms  4rfcKs|dkrt||}|dkr|S|dkr|j|S|j|S|jd|dd}d|}|||||j||dd||fi|S)NrrrXTrrr)r:r>r?rHr!rYhyperu)r rrQrrBgr rer r r whitws  4ric  s |\}}|\}}s/|dkr*d|g||dgSjSd||}|\}}z,j} z$jd7_jdd||f||gdjd} | |W| _WS| _wjyrYnwfdd} j| ||gfi|S) Nrr*rrr)maxtermscs||}j|gddgg||d|g|g|gf}j |gddd|gg|d|g||dgd|gf}||fS)Nrrr)sinpir")abr/rUrVr[r r r)s ,>zhyperu..h) _convert_paramr2r> gammaprodr?rhypsum NoConvergencer9) r rlrmrrBatypebtypebbbbtyperDrEr)r r[r rgs*   $ rgc 8|}fdd}j||gfi|S)NcsJddjg|ddgg|dgdgd|dgdd fgSNrrrr?sqrtr"rZr[r r r)sJzstruveh..hr2r9r rrrBr)r r[r struveh  r}c rv)NcsHddjg|ddgg|dgdgd|dgddfgSrwryrZr[r r r)sHzstruvel..hr{r|r r[r struvelr~rc s>|d}fdd}j||gfi|S)Nrcsj}||}|d}||||d|d|f\}}}}|\}} dkr2|| g|g} } dkrA|| g| g} } jdd} | ddgg||gdg||g| f} | dgg||gdg||g| f}| |fS)Nr,rrrmult)mpq_1_2rOsquare_exp_arg)rErmurQa1a2b1b2csAr%r/rUrVr whichrr r r)s$ z_anger..hrnr2r9)r rrErrBr)r rr _angers rcKt|d||fi|SNrrr rErrBr r r angerjrcKrNrrrr r r webererrc L|d}|d}fdd}j|||gfi|S)Nrc sfj}jdd}||d||dgdd|dgggdg|||d|||dg|ffS)Nrrrrr,rr)rrErmr/r[r r r)"s , zlommels1..hrr rrErrBr)r r[r lommels1s  rc r)Nrc sj}jdd}||d||dgdd|dgggdg|||d|||dg|f}dg||d| g||||dg|||dggd|g|f}dg||d|g| |||dg|d||ggd|g|f}|||fS)Nrrrrr,rr)rrErmr/rUrVT3r[r r r)4s PHH zlommels2..hrrr r[r lommels2)s  rc rv)Nc sdd }d|\}}|dgd|gg|dggdd|dd|dg|f}|dgd|dgg|dggdd|dd|dg|f}||fS)Nrrrrxr,rOrr$rRrSrUrVr[r r r)Ks :>zber..hr{r|r r[r berF  rc rv)Nc sdd }d|\}}|dgd|dgg|dggdd|dd|dg|f}|dgd|gg|dggdd|dd|dg|f}||fS)Nr?rrrxrr,rrr[r r r)Xs >:zbei..hr{r|r r[r beiSrrc rv)Nc sLdd }d|\}}d|\}}dd|g| d|dg| gggddd|d|dg|f}d| g| dd|dg| dgggddd|d|dg|f}dd|g|d| dg|gggddd|dd|g|f}d| g|dd|dg|dgggddd|dd|g|f} |||| fS) NrrFrrr,rrrxr rr$cos1sin1cos2sin2rUrVrT4r[r r r)esBH@D zker..hr{r|r r[r ker`   rc rv)Nc sHdd }d|\}}d|\}}| dgd|dd|g|dgggddd|dd|g|f}| dgd|d| g|gggddd|dd|g|f}| dgd| d|g| gggdd|dd|dg|f}| dgd| d|dg| dgggdd|dd|dg|f} |||| fS) NrrrFrrr,rxrrrr[r r r)usD>@H zkei..hr{r|r r[r keiprrcsjfdd}|S)NcsF|j}|j}|d\}}||kr| S||f|<|dS)N)rrr)_misc_const_cacherget)r cacherprEfnamer r f_wrappeds zc_memo..f_wrapped)__name__)rrr rr c_memos rcCs"d|d||ddS)Nr rr,cbrtgammar8r r r r _airyai_C1"rcCs"d|d||ddS)Nrr,rrrr r r _airyai_C2rrcCs$d|dd||ddS)Nrr,rnthrootrr8rr r r _airybi_C1s$rcCs |dd||ddS)Nr,rrrrr r r _airybi_C2s rcCs>|j}z|dd|dd|j}W||_| S||_w)Nr,2/3r)rpowerrr")r rrEr r r _airybi_n2_infs "rcCs|dkrk|dkr |S|j}|j}zQ|jd7_||d||d|||j}|dkrF||d|d|9}||dd}n|t|d|d|9}||dd}W||_n||_w| |St) NZrr*rr,rrz1/6)mpq_1_3rrrr"rkr:r<)r rrntyperr$rrEr r r _airyderiv_0s"( rc s|r|\}}nd}smrm|rY|dkrY|dkrBjkr2dddSjkrBdddS|dkrYjkrMSjkrYd| S|s`jksejkridStdr|tdtd nd|r|dkrfd d }j |gfi|Sdkrt ||dSfd d }j ||gfi|} rĈ |rĈ|}|Sfd d }j |gfi|S) Nrrrrr,zessential singularity of Ai(z)rxcsdkrEj7_d}d|}d|d}j8_| djd}|gdgggdd gg|ffSj7_dd }j8_td }t}|gddggggjg|f}|gdggggj g|f}||fS) Nrrxrrr,rr)rr)rrr) r4rrYrzr"rrrmpq_5_3rr/r$rCC1C2rUrVr extraprecrr r r)s(  zairyai..hc sj7_dd}j8_jjj}}}|}d}d||}d||}d||} dg||| g|g||| g||g||| g|f} |}d||}d||}d||} d g||| dg|g||| g||g||| g|f} | | fSNr,rrrr)rrmpq_2_3mpq_4_3) rr/q13q23q43rrrrb3rUrVrr r r)s ,($csdkrDj7_d}d|}d|d}j8_|djd}|gdgggdd gg|ffSj7_dd }j8_t}t}|gdggggjg|f}|gdggggj g|f}||fS) Nrrxrrr,rr)rr)rr) r4rrYrzr"rrrrrrrr r r)s& )r2rnisnormalr?r8rLrGr.r1r5r9r _is_real_typer3r4r rrArBrrr)rEr rr airyaisF        rc s|r|\}}nd}s[r[|rE|dkrEjkr%SjkrE|dkr2dS|dkr:tS|dkrEd| S|sWjkrNSjkrWdStdrjtdtd nd|r|dkrfdd }j |gfi|Sdkrt ||dSfd d }j ||gfi|} r |r|}|Sfd d }j |gfi|S) Nrrrrrzessential singularity of Bi(z)rxcs|j7_dd}j8_td}t}|gddggggjg|f}|gdggggjg|f}||fS)Nr,rrrr)rrrrrr/rrrUrVrr r r);s  zairybi..hcs j7_dd}j8_jjj}}}j}j}|}d}d||} d||} d||} dg||| g|g| | | g||g| | | g|f} |}d||} d||} d||} dg||d|g|g| | | g||g| | | g|f} | | fSr)rrrrmpq_1_6mpq_5_6)rr/rrrq16q56rrrrrrUrVrr r r)Hs ,( csxj7_dd}j8_t}t}|gdggggjg|f}|gdggggjg|f}||fS)Nr,rr)rrrrrrrr r r)[s )r2rnrr?rLrrGr.r1r5r9rrr3r4rr rr airybisJ         rFc sdd}dd}t|}|dkrtd|dvrtd|d krO|r:fd d |d jd |d d Sj|d jd |dd S|dkr|dkr|rqfdd |d jd |dd Sj|d jd |d d S|dkr|dkr|rd jd |d ddj}dd ||}fdd |Sd jd |dddj}dd ||}j|SdSdS)NcSs|ddd|ddS)NUUUUUU?rrr0r tr r r Uhz_airy_zero..UcSs|ddd|ddS)Nrrrrrr rr r r r&irz_airy_zero..Trzk cannot be less than 1rrz%Derivative should lie between 0 and 1rc |dSr)rrWrr r q z_airy_zero..r,rFcrrrrWrr r rvrTy?crrrrWrr r r}r) r1rGfindrootr"rrln2expjpir8) r rrrAcomplexrr&rrr rr _airy_zerofs8**$$rcCst|d||dS)NrFr)r rrAr r r airyaizerorcCst|d|||Srr)r rrArr r r airybizerorrcsjr*jkrdkrdSdkrSjkr&dStdr9tdtdnddrBt zTdkrdkrnt j ddkrnfdd }j |gjd d WSdkrt  d j ddkrfd d }j |gjd d WSWn jyYnwfdd }j |gfiS)Nrrzessential singularityrxrAr,g+?c s*jgddggggdgddffSNr))rr,)rr,rrr,r"r r[r r r)s*z_scorer..hT)rj force_seriesrc s,j gddggggdgddffSrrr r[r r r)s,csjfid}dj}dkr|d9}|d9}j7_dd}j8_|gdgggggdf}|gddgggdgjjg|f}||fS)Nr,rrrrrr)rr"rrr)rr%r/rUrVr rrBrrr r r)s  $)r2isinfr?rLrGr.r1r5rr<r:argr"r9rrq)r rrrBr)r rr _scorers6      $* rcKt||d|Srrr rrBr r r scorergircKrrrrr r r scorerhirr cCs||f|vr|||fd|jkr|||fd S|d|d}|d||j|}|d||j|}d|||j |||d|}||s_||s_||}|j|f|||f<|S)Nrrr)rloggammar_rYr"rJr>)r leta_cacheG3G1G2rEr r r coulombcs"* rTc  sffdd}j|||gfi|}|r1|s1|s1s1dkr1|}|S)Nc sz?j}j|dd}j|ddd}||}||gd|ddgggd|||gd|dg|f}W|fStyUdgdgggggdf}Y|fSw)NTrrrrrr)r_r!rrYrG)r r jwjwzjwz2rrUr r/rr r r)s  . zcoulombf..hr)r9rJr> r r r rr/choprBr)rEr rr coulombfs  " rcsdf|vr|fdjkr|fdSfdd}|d}j|f|f<|S)Nrrcsv d}j}d|dd|dd||dd||dd jgS)Nryy?r)r_r r")l2jetar r r r r termss  z_coulomb_chi..terms)rsum_accurately)r r r r rrEr rr _coulomb_chis " rc  sz|s |}fdd}j|||gfi|}|r;|s;|s;s;dkr;|}|S)Nc sL|drdgdgggggdf}|fS| d}zt||}j}|}|}||}||} |} d|} ||| |gdd|dddgggd|||gd|dg| f}| | | gdd|ddgggd|||gd|dg| f} || fWStydgdgggggdf}|fYSw)Nrrrrr)r3rr_rSrRrrYrG) r r rUrchirrrrrrr rVrr r r)s.          zcoulombg..hr)_imr4r9rr rr coulombgs  " r!cCsd|d}|dkr|sd|d|d|jd}|dkr0|s0d|d|d|jd}|dkrE|rEd|d|d|jd}|dkrZ|rZd|d|d|jd}|s|}|d d|}d|dd|ddd|d} d |dd |dd |d d d|d} d|dd|dd|dd|ddd|d} |r|}|d d|}dd|dd|ddd|d} d d |dd|dd|dd d|d} dd|dd|dd|dd|ddd|d} ||| | | g} |} d}tdt| D] }t| |t| |dkrB| | |7} q)t| |}q)|t| dkrYt| d}| |fS) aj Computes an estimate for the location of the Bessel function zero j_{v,m}, y_{v,m}, j'_{v,m} or y'_{v,m} using McMahon's asymptotic expansion (Abramowitz & Stegun 9.5.12-13, DLMF 20.21(vi)). Returns (r,err) where r is the estimated location of the root and err is a positive number estimating the error of the asymptotic expansion. rrrr,rriSiirrii%iXiO2iu_iRrii i i,i il"qiQYgr)r"r7lenr:)r kindprimerErQrrms1s2s3s4s5rrerrir r r mcmahons6 ****(4@,8D r1csz|dkrtd|d}gg |||fddDfddt|dD}t||kr8|S|d}q)z Given f known to have exactly n simple roots within [a,b], return a list of n intervals isolating the roots and having opposite signs at the endpoints. TODO: this can be optimized, e.g. by reusing evaluation points. rzn cannot be less than 1csg|] }|qSr )sign)rr )r rr r Jsz)generalized_bisection..cs8g|]}||ddkr||dfqS)rrr )rr0)pointssignsr r r3Ksr)rGlinspacer7r')r rrlrmrN ok_intervalsr )r rr4r5r generalized_bisection;s r9cCs|j||dddS)NillinoisF)solververify)r)r rabr r r find_in_intervalQsr>g{Gz?csj}t||d}z3|_t|}t|}dkr+td|dkr3td|dvr;td|dkrP|rIfdd } nfd d } |d kre|r^fd d } nfd d } |dkr|r|dkrdkrzjW|_Sdkrd dd } t| | dd | fW|_S|||f|vrt| ||||fW|_St |||\} } | |krt| | || |fW|_S|dkr|sd} |dkr|rd} |d kr|sd} |d kr|rd} |d} t ||| \}} | |krAt ||| d\}}t | | d||| }t |D]\}}|||||df<q#t| ||dW|_S| d } q|_w)Nr*rzv cannot be negativerzm cannot be less than 1rz prime should lie between 0 and 1cj|ddSNr)rArr r rEr r rczbessel_zero..c |Sr]rrArBr r rdrrcr?r@rTrArBr r rfrCcrDr]rErArBr r rgrg333333@g?g?g@r) rr.r5r8r1rGzerorzr>r1r9 enumerate)r r(r)rErQisoltol_interval_cacherworkprecrr$r/lowrr1r2err2 intervalsrr=r rBr bessel_zeroTsf   rPcCt|d||| S)a For a real order `\nu \ge 0` and a positive integer `m`, returns `j_{\nu,m}`, the `m`-th positive zero of the Bessel function of the first kind `J_{\nu}(z)` (see :func:`~mpmath.besselj`). Alternatively, with *derivative=1*, gives the first nonnegative simple zero `j'_{\nu,m}` of `J'_{\nu}(z)`. The indexing convention is that used by Abramowitz & Stegun and the DLMF. Note the special case `j'_{0,1} = 0`, while all other zeros are positive. In effect, only simple zeros are counted (all zeros of Bessel functions are simple except possibly `z = 0`) and `j_{\nu,m}` becomes a monotonic function of both `\nu` and `m`. The zeros are interlaced according to the inequalities .. math :: j'_{\nu,k} < j_{\nu,k} < j'_{\nu,k+1} j_{\nu,1} < j_{\nu+1,2} < j_{\nu,2} < j_{\nu+1,2} < j_{\nu,3} < \cdots **Examples** Initial zeros of the Bessel functions `J_0(z), J_1(z), J_2(z)`:: >>> from mpmath import * >>> mp.dps = 25; mp.pretty = True >>> besseljzero(0,1); besseljzero(0,2); besseljzero(0,3) 2.404825557695772768621632 5.520078110286310649596604 8.653727912911012216954199 >>> besseljzero(1,1); besseljzero(1,2); besseljzero(1,3) 3.831705970207512315614436 7.01558666981561875353705 10.17346813506272207718571 >>> besseljzero(2,1); besseljzero(2,2); besseljzero(2,3) 5.135622301840682556301402 8.417244140399864857783614 11.61984117214905942709415 Initial zeros of `J'_0(z), J'_1(z), J'_2(z)`:: 0.0 3.831705970207512315614436 7.01558666981561875353705 >>> besseljzero(1,1,1); besseljzero(1,2,1); besseljzero(1,3,1) 1.84118378134065930264363 5.331442773525032636884016 8.536316366346285834358961 >>> besseljzero(2,1,1); besseljzero(2,2,1); besseljzero(2,3,1) 3.054236928227140322755932 6.706133194158459146634394 9.969467823087595793179143 Zeros with large index:: >>> besseljzero(0,100); besseljzero(0,1000); besseljzero(0,10000) 313.3742660775278447196902 3140.807295225078628895545 31415.14114171350798533666 >>> besseljzero(5,100); besseljzero(5,1000); besseljzero(5,10000) 321.1893195676003157339222 3148.657306813047523500494 31422.9947255486291798943 >>> besseljzero(0,100,1); besseljzero(0,1000,1); besseljzero(0,10000,1) 311.8018681873704508125112 3139.236339643802482833973 31413.57032947022399485808 Zeros of functions with large order:: >>> besseljzero(50,1) 57.11689916011917411936228 >>> besseljzero(50,2) 62.80769876483536093435393 >>> besseljzero(50,100) 388.6936600656058834640981 >>> besseljzero(50,1,1) 52.99764038731665010944037 >>> besseljzero(50,2,1) 60.02631933279942589882363 >>> besseljzero(50,100,1) 387.1083151608726181086283 Zeros of functions with fractional order:: >>> besseljzero(0.5,1); besseljzero(1.5,1); besseljzero(2.25,4) 3.141592653589793238462643 4.493409457909064175307881 15.15657692957458622921634 Both `J_{\nu}(z)` and `J'_{\nu}(z)` can be expressed as infinite products over their zeros:: >>> v,z = 2, mpf(1) >>> (z/2)**v/gamma(v+1) * \ ... nprod(lambda k: 1-(z/besseljzero(v,k))**2, [1,inf]) ... 0.1149034849319004804696469 >>> besselj(v,z) 0.1149034849319004804696469 >>> (z/2)**(v-1)/2/gamma(v) * \ ... nprod(lambda k: 1-(z/besseljzero(v,k,1))**2, [1,inf]) ... 0.2102436158811325550203884 >>> besselj(v,z,1) 0.2102436158811325550203884 rrPr rErQrAr r r besseljzerosprTcCrQ)a For a real order `\nu \ge 0` and a positive integer `m`, returns `y_{\nu,m}`, the `m`-th positive zero of the Bessel function of the second kind `Y_{\nu}(z)` (see :func:`~mpmath.bessely`). Alternatively, with *derivative=1*, gives the first positive zero `y'_{\nu,m}` of `Y'_{\nu}(z)`. The zeros are interlaced according to the inequalities .. math :: y_{\nu,k} < y'_{\nu,k} < y_{\nu,k+1} y_{\nu,1} < y_{\nu+1,2} < y_{\nu,2} < y_{\nu+1,2} < y_{\nu,3} < \cdots **Examples** Initial zeros of the Bessel functions `Y_0(z), Y_1(z), Y_2(z)`:: >>> from mpmath import * >>> mp.dps = 25; mp.pretty = True >>> besselyzero(0,1); besselyzero(0,2); besselyzero(0,3) 0.8935769662791675215848871 3.957678419314857868375677 7.086051060301772697623625 >>> besselyzero(1,1); besselyzero(1,2); besselyzero(1,3) 2.197141326031017035149034 5.429681040794135132772005 8.596005868331168926429606 >>> besselyzero(2,1); besselyzero(2,2); besselyzero(2,3) 3.384241767149593472701426 6.793807513268267538291167 10.02347797936003797850539 Initial zeros of `Y'_0(z), Y'_1(z), Y'_2(z)`:: >>> besselyzero(0,1,1); besselyzero(0,2,1); besselyzero(0,3,1) 2.197141326031017035149034 5.429681040794135132772005 8.596005868331168926429606 >>> besselyzero(1,1,1); besselyzero(1,2,1); besselyzero(1,3,1) 3.683022856585177699898967 6.941499953654175655751944 10.12340465543661307978775 >>> besselyzero(2,1,1); besselyzero(2,2,1); besselyzero(2,3,1) 5.002582931446063945200176 8.350724701413079526349714 11.57419546521764654624265 Zeros with large index:: >>> besselyzero(0,100); besselyzero(0,1000); besselyzero(0,10000) 311.8034717601871549333419 3139.236498918198006794026 31413.57034538691205229188 >>> besselyzero(5,100); besselyzero(5,1000); besselyzero(5,10000) 319.6183338562782156235062 3147.086508524556404473186 31421.42392920214673402828 >>> besselyzero(0,100,1); besselyzero(0,1000,1); besselyzero(0,10000,1) 313.3726705426359345050449 3140.807136030340213610065 31415.14112579761578220175 Zeros of functions with large order:: >>> besselyzero(50,1) 53.50285882040036394680237 >>> besselyzero(50,2) 60.11244442774058114686022 >>> besselyzero(50,100) 387.1096509824943957706835 >>> besselyzero(50,1,1) 56.96290427516751320063605 >>> besselyzero(50,2,1) 62.74888166945933944036623 >>> besselyzero(50,100,1) 388.6923300548309258355475 Zeros of functions with fractional order:: >>> besselyzero(0.5,1); besselyzero(1.5,1); besselyzero(2.25,4) 1.570796326794896619231322 2.798386045783887136720249 13.56721208770735123376018 rrRrSr r r besselyzerosYrUN)r)F)rF)rT)2 functionsrrr rrrIrTr\rarcrfrirgr}rrrrrrrrrrrrrrrrrrrrrrrrr rrrr!r1r9r>rPrTrUr r r r s   B # #                      [  K  +    "'5 q