U Zʹj2@sFddlmZmZeddZeddZededdZedfd d Zedgd 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<Zedhd=d>Z edid?d@Z!djdBdCZ"edkdDdEZ#edldFdGZ$dHdIZ%edJdKZ&edLdMZ'eifdNdOZ(edmdQdRZ)eifdSdTZ*edndUdVZ+dWdXZ,dYdZZ-d[d\Z.d]ifd^d_Z/edod`daZ0edpdbdcZ1ddS)q)defun defun_wrappedcCs |d|S)zCComputes the Bessel function `J_0(x)`. See :func:`~mpmath.besselj`.besseljctxxr V/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mpmath/functions/bessel.pyj0sr cCs |d|S)zDComputes the Bessel function `J_1(x)`. See :func:`~mpmath.besselj`.rrrr r r j1sr rc  s&ttkrd}n&}|r8t|r`dkr`dj |f|S|r$|rdkrtj}z8jd7_fddt dD}W5|_X| d 9}n$fd d }j |gf|}n|sr|rrt d krrt d krrz WStk rpYnXssj}n(dkr}nj}ndj}zTjtd t j7_jdddfdd }j |gf|}W5|_X| }|S)NTrc3s:|]2}d||d|VqdS)rN)Zbinomialr).0k)rdnzr r %szbesselj..rrcsjjjdddd}d||dd||dg}djg|d|d||gg||dd|ddg||dg|fg}|S)NprecпTexact?rrfmulrpirrrBTMrrr r h+s" Lzbesselj..h rrcsJjjtdjddd}g|gg|dgg|dg|fgS)NrrTrr)Zfnegrmaxr)rr!)r%rwr r r&Gs&)typeintconvertisint_rermagrfsumrangempf hypercombabsZ_besseljNotImplementedErroronereinfminr) rrr derivativekwargsZn_isintorigvr&r )r%rrrr+rr r sX        ( rc  s|}s|r t|s0d|S|rFd|S|}|dkrfj|S|dkrzd|Sj|S|rƈ|}fdd}j|||gf|}n"fdd}j||gf|}|S)Nrrcsjjjdddd}d||dd||d|dg}djg|d|d||g|dg||dd|ddg||fg}|S)Nr?Trrrrrr r$r r r&fs"&Hzbesseli..hcsPjddd}j||tdjd}|g|gg|dgg|dg|fgS)NrTrrrr)rr*r)rr+r!r$r r r&ms)r. ValueErrorr/r9nanr:r1r5) rrrr<r=r!rr&r?r r$r besseliPs.         rCc KsJ|s|r t|s |j ||S||r8|j||S||}|d}||rx|dkrl|j ||Sd||S|dkrt||dr|j||S|j||S|j d7_ | |\}}||j kr|j } |j d9_ || 7}n|dkr|j |8_ | |\} } |j |||f|| |j | ||f|| S)Nrrrr')rAr:imrBr9r/r-floorninfrZ nint_distanceeps cospi_sinpir) rrrr<r=r!qmrr&cossinr r r besselyts<       rMc sXs jS}|dkr*fdd}nj|7_fdd}j||gf|S)Nrcsfdd}dg| |dg|gggd|g|f}dg|| dg| gggd|g|f}||fS)Nrrr )rr!T1T2rr r r&s (*zbesselk..hc sBjd gdddggg|dd|ggddfgS)Nrrrr)rexprrrr r r&s " )r:r1rr5)rrrr=r%r&r rTr besselks  rUcKs&|j||f||j|j||f|SNrjrMrrr r=r r r hankel1srZcKs&|j||f||j|j||f|SrVrWrYr r r hankel2sr[cKs|dkr<||dkr|S||dkr2|j|S|j|S|jd|dd}d|}|||||j||dd||f|S)NrrQTrrrr)r9r:rBrrRZhyp1f1)rrrJrr=r yr r r whitms  r]cKs|dkr>t||}|dkr"|S|dkr4|j|S|j|S|jd|dd}d|}|||||j||dd||f|S)NrrrQTrrr)r6r9r:rBrrRhyperu)rrrJrr=gr r\r r r whitws  r`c  s |\}}|\}}s^|dkrTd|g||dgSjSd||}|\}}zXj} zDjd7_jdd||f||gdjd} | |WWS| _XWnjk rYnXfdd} j| ||gf|S) Nrr'rrr)maxtermscs||}j|gddgg||d|g|g|gf}j |gddd|gg|d|g||dgd|gf}||fS)Nrrr)sinpir)abr+rNrOrTr r r&s ,>zhyperu..h) _convert_paramr.r9Z gammaprodr:rZhypsum NoConvergencer5) rrcrdrr=ZatypeZbtypeZbbZbbtyper>r?r&r rTr r^s(   $ r^c s4|}fdd}j||gf|S)NcsJddjg|ddgg|dgdgd|dgdd fgSNrrrr?sqrtrrSrTr r r&szstruveh..hr.r5rrrr=r&r rTr struvehs  rmc s4|}fdd}j||gf|S)NcsHddjg|ddgg|dgdgd|dgddfgSrgrirSrTr r r&szstruvel..hrkrlr rTr struvels  rnc s:|d}fdd}j||gf|S)Nrcsj}||}|d}||||d|d|f\}}}}|\}} dkrd|| g|g} } dkr|| g| g} } jdd} | ddgg||gdg||g| f} | dgg||gdg||g| f}| |fS)Nr)rrrZmult)mpq_1_2rHsquare_exp_arg)r?rdurJa1a2b1b2csAr"r+rNrOrwhichrr r r&s$ z_anger..hrer.r5)rr{r?rr=r&r rzr _angers r}cKst|d||f|SNrr}rr?rr=r r r angerjsrcKst|d||f|SNrrrr r r weberesrc sH|d}|d}fdd}j|||gf|S)Nrc sfj}jdd}||d||dgdd|dgggdg|||d|||dg|ffS)Nrrorrr)rprq)rrr?rdr+rTr r r&"s ,zlommels1..hr|rrrr?rr=r&r rTr lommels1s  rc sH|d}|d}fdd}j|||gf|S)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)Nrrorrr)rr)rrr?rdr+rNrOT3rTr r r&4s PHH zlommels2..hr|rr rTr lommels2)s  rc s4|}fdd}j||gf|S)Nc sdd }d|\}}|dgd|gg|dggdd|dd|dg|f}|dgd|dgg|dggdd|dd|dg|f}||fS)Nrrrrhr)rHrr!rKrLrNrOrTr r r&Ks :>zber..hrkrlr rTr berFs  rc s4|}fdd}j||gf|S)Nc sdd }d|\}}|dgd|dgg|dggdd|dd|dg|f}|dgd|gg|dggdd|dd|dg|f}||fS)Nr?rrrhrr)rrrTr r r&Xs >:zbei..hrkrlr rTr beiSs  rc s4|}fdd}j||gf|S)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) Nrr@rrr)rrrhr rr!Zcos1Zsin1Zcos2Zsin2rNrOrZT4rTr r r&esBH@Dzker..hrkrlr rTr ker`s   rc s4|}fdd}j||gf|S)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) Nrrr@rrr)rhrrrrTr r r&usD>@Hzkei..hrkrlr rTr keips   rcsjfdd}|S)NcsJ|j}|j}|d\}}||kr*| S||f|<|dSdS)N)rrr)Z_misc_const_cacherget)rcacherpr?fnamer r f_wrappedszc_memo..f_wrapped)__name__)rrr rr c_memos rcCs"d|d||ddS)Nr rr)Zcbrtgammar4rr r r _airyai_C1srcCs"d|d||ddS)Nrr)rrrr r r _airyai_C2srcCs$d|dd||ddS)Nrr)rnthrootrr4rr r r _airybi_C1srcCs |dd||ddS)Nr)rrrrr r r _airybi_C2srcCs:|j}z$|dd|dd|j}W5||_X| S)Nr)2/3r)rpowerrr)rrr?r r r _airybi_n2_infs $rcCs|dkr|dkr|S|j}|j}z|jd7_||d||d|||j}|dkr||d|d|9}||dd}n.|t|d|d|9}||dd}W5||_X| |StdS) NZrr'rr)rrz1/6)mpq_1_3rrrrrbr6r7)rrrntyper{r!rr?r r r _airyderiv_0s ( rc s|r|\}}nd}sڈr|r|dkr|dkrjkrddddSjkrdddS|dkrjkrSjkrd| S|sjksʈjkrdStdrtdtd nd|r|dkr,fd d }j |gf|SdkrFt ||dSfd d }j ||gf|} r |r|}|Sn fd d }j |gf|SdS) Nrrrrr)zessential singularity of Ai(z)rhcs dkrj7_d}d|}d|d}j8_| djd}|gdgggdd gg|ffSj7_dd }j8_td }t}|gddggggjg|f}|gdggggj g|f}||fSdS) Nrrhrrr)rr)rr)rrr) r0rrRrjrrrrmpq_5_3rr+r!rrCC1C2rNrOr 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+q13q23q43rsrtrurvb3rNrOrr r r&s0   $csdkrj7_d}d|}d|d}j8_|djd}|gdgggdd gg|ffSj7_dd }j8_t}t}|gdggggjg|f}|gdggggj g|f}||fSdS) Nrrhrrr)rr)rr)rr) r0rrRrjrrrrrrrrr r r&s  & )r.reisnormalr:r4rFrAr*r-r1r5r _is_real_typer/r0rrr<r=rrr&r?r rr airyaisF          rc s|r|\}}nd}sr|r|dkrjkrJSjkr|dkrddS|dkrttS|dkrd| S|sjkrSjkrdStdrtdtd nd|rl|dkrfdd }j |gf|Sdkr"t ||dSfd d }j ||gf|} rf |rf|}|Sn fd d }j |gf|SdS) Nrrrrrzessential singularity of Bi(z)rhcs|j7_dd}j8_td}t}|gddggggjg|f}|gdggggjg|f}||fS)Nr)rrrr)rrrrrr+rrrNrOrr 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)rrrrZmpq_1_6Zmpq_5_6)rr+rrrZq16Zq56rsrtrurvrrNrOrr r r&Hs4    csxj7_dd}j8_t}t}|gdggggjg|f}|gdggggjg|f}||fS)Nr)rr)rrrrrrrr r r&[s )r.rerr:rFrrAr*r-r1r5rrr/r0rr rr airybisJ           rFc sdd}dd}t|}|dkr(td|dkr8td|d kr|rtfd d |d jd |d d Sj|d jd |dd S|dkr|dkr|rfdd |d jd |dd Sj|d jd |d d S|dkr|dkr|r~d jd |d ddj}dd ||}fdd |Sd jd |dddj}dd ||}j|SdS)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 1rcs |dSr)rrPrr r qrz_airy_zero..r)rFcs |dSrrrPrr r rvrTy?cs |dSrrrPrr r r}r) r-rAfindrootrrrZln2Zexpjpir4) rr{rr<complexrr#rrxr rr _airy_zerofs6**$$rcCst|d||dS)NrFr)rrr<r r r airyaizerosrcCst|d|||Srr)rrr<rr r r airybizerosrcsprTjkr:dkr.dSdkr:SjkrLdStdrrtdtdnddrt zdkr0dkrt j ddkrއfdd }j |gjd d WSdkr0t  d j ddkr0fd d }j |gjd d WSWnjk rJYnXfdd }j |gfS)Nrrzessential singularityrhr<r)g+?c s,jgddgggdddggddffSNr)rr))rr)rrr)rr rTr r r&sz_scorer..hT)raZ force_seriesrc s.j gddgggdddggddffSrrr rTr r r&scsjfd}dj}dkr4|d9}|d9}j7_dd}j8_|gdgggggdf}|gddgggdgjjg|f}||fS)Nr)rrrrrr)rrrrr)ryr"r+rNrOrrr=r{rr r r&s  $)r.isinfr:rFrAr*r-r1rr7r6argrr5rrf)rrr{r=r&r rr _scorers6     $. rcKst||d|Sr~rrrr=r r r scorergisrcKst||d|Srrrr r r scorerhisrcCs||f|kr4|||fd|jkr4|||fd S|d|d}|d||j|}|d||j|}d|||j |||d|}||s||s||}|j|f|||f<|S)Nrrr)rloggammarXrRrrDr9)rleta_cacheZG3ZG1ZG2r?r r r coulombcs"* rTc  sbfdd}j|||gf|}|r^|s^|s^s^dkr^|}|S)Nc szzj}j|dd}j|ddd}||}||gd|ddgggd|||gd|dg|f}Wn*tk rdgdgggggdf}YnX|fS)NTrrrrrr)rXrrrRrA)rrjwZjwzZjwz2rrNrr+rr r r&s  . zcoulombf..hr)r5rDr9 rrrrr+Zchopr=r&r?r rr coulombfs "  rcsdf|kr2|fdjkr2|fdSfdd}|d}j|f|f<|S)Nrrcsv d}j}d|dd|dd||dd||dd jgS)Nryy?r)rXrr)l2Zjetarrrr r termss  z_coulomb_chi..terms)rZsum_accurately)rrrrrr?r rr _coulomb_chis " rc  sv|s|}fdd}j|||gf|}|rr|sr|srsrdkrr|}|S)Nc sT|dr*dgdgggggdf}|fS| d}z||}j}|}|}||}||} |} d|} ||| |gdd|dddgggd|||gd|dg| f}| | | gdd|ddgggd|||gd|dg| f} || fWStk rNdgdgggggdf}|fYSXdS)Nrrrrr)r/rrXrLrKrrRrA) rrrNrchirrxrwrrrrr rOrr r r&s6           zcoulombg..hr)Z_imr0r5rr rr coulombgs  "  rcCsd|d}|dkr6|s6d|d|d|jd}|dkr`|s`d|d|d|jd}|dkr|rd|d|d|jd}|dkr|rd|d|d|jd}|sl|}|d d|}d|dd|ddd|d} d |dd |dd |d d d|d} d|dd|dd|dd|ddd|d} |r0|}|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| |dkr| | |7} n t| |}qT|t| dkrt| 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)rr3lenr6)rkindprimer?rJrrrds1s2Zs3Zs4Zs5rrxerrir r r mcmahons>     (4@,8D rcs||dkrtd|d}gg|||fddDfddt|dD}t||krn|S|d}q dS)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 )rrr r Jsz)generalized_bisection..cs8g|]0}||ddkr||dfqS)rrr )rr)pointssignsr r rKsrN)rAZlinspacer3r)rrrcrdrNZ ok_intervalsr )rrrrr generalized_bisection;s rcCs|j||dddS)NZillinoisF)Zsolververify)r)rrabr r r find_in_intervalQsrg{Gz?csj}t||d}zx|_t|}t|}dkrVtd|dkrftd|dkrvtd|dkr|rfdd } nfd d } |d kr|rfd d } nfd d } |dkr>|r>|dkr>dkrjWSdkr>d dd } t| | dd | fWbS|||f|krnt| ||||fW2St |||\} } | |krt| | || |fWS|dkr|sd} |dkr|rd} |d kr|sd} |d kr|rd} |d} t ||| \}} | |krt ||| d\}}t | | d||| }t |D]\}}|||||df<qVt| ||dWS| d } qW5|_XdS)Nr'rzv cannot be negativerzm cannot be less than 1rz prime should lie between 0 and 1csj|ddSNr)r<rr rr?r r rcrzbessel_zero..cs |SrVrrrr r rdrrcsj|ddSrrMrrr r rfrcs |SrVrrrr r rgrg333333@g?g?g@r) rr*r1r4r-rAzerorjrrr enumerate)rrrr?rJZisoltolZ_interval_cacherZworkprecrr!rlowrr1r2Zerr2Z intervalsrrr rr bessel_zeroTsd     r cCst|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 rr rr?rJr<r r r besseljzerosprcCst|d||| S)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 rrrr r r besselyzerosYrN)r)r)r)r)r)F)r)rF)rT)rT)r)r)2Z functionsrrr r rrCrMrUrZr[r]r`r^rmrnr}rrrrrrrrrrrrrrrrrrrrrrrrrrrrrrr rrr r r r s   B # #                       [ K   +     "'5 q