є <Zc@s╪ddlmZmZedДГZedДГZeddДГZeddДГZeddДГZedДГZed ДГZ ed ДГZ edДГZedДГZed ДГZ edДГZedДГZdДZedДГZedДГZedДГZedДГZedДГZedДГZedДГZedДГZdДZedДГZedДГZedДГZedДГZdДZdДZedd ДГZ edd!ДГZ!e"d"ДZ#edd#ДГZ$ede"d$ДГZ%d%ДZ&ed&ДГZ'ed'ДГZ(eid(ДГZ)ede*d)ДГZ+eid*ДГZ,ede*d+ДГZ-d,ДZ.d-ДZ/d.ДZ0d/id0ДZ1edd1ДГZ2edd2ДГZ3d3S(4i(tdefunt defun_wrappedcCs|jd|ГS(sCComputes the Bessel function `J_0(x)`. See :func:`~mpmath.besselj`.i(tbesselj(tctxtx((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytj0scCs|jd|ГS(sDComputes the Bessel function `J_1(x)`. See :func:`~mpmath.besselj`.i(R(RR((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytj1sic s╠tИГtkrt}n<ИjИГЙИjИГ}|rWtИjИГГЙn|rИИdkrИdИИjИИ||НSИjИГЙИjИГЙ|rЛИj|ГЙИjИГrXИdkrXtИГЙИj}zEИjd7_Иj ЗЗЗЗfdЖt ИdГDГГ}Wd|И_X|ИjdГИ9}q╚ЗЗЗfdЖ}Иj|ИИg|Н}n=|rч|rчt ИГdkrчt ИГd krчyИjИИГSWqчtk rуqчXnИs=ИsИjИИ}q┴ИjИГdkr)ИИ}q┴ИjИИ}nДИj}znИjtd t ИГИjГ7_ИjИddtГЙЗЗЗfd Ж}Иj|Иg|Н}Wd|И_X| }|S(Nii ic3sE|];}d|ИjИ|ГИjd|ИИИГVqdS(i iN(tbinomialR(t.0tk(Rtdtntz(s6lib/python2.7/site-packages/mpmath/functions/bessel.pys %siics╣ИjИjИИdИjИГddtГ}d||dd||dg}dИjИg|d|d||gg||dd|ddg||dg|fg}|S(Ntprecg╨┐texactgр?ii(tfmulR tTruetpi(RR trtBtT(tMRR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyth+s1$`i iigр?Rc sgИjИjИИdtdИjИГГdtГ}Иg|gg|dgg|dg|fgS(NR iRi(tfnegRtmaxR R(RR(RRtw(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRGs7(ttypetintRtconverttisintt_reRtmagR tfsumtrangetmpft hypercombtabst_besseljtNotImplementedErrortonetretinftminR( RRRt derivativetkwargstn_isinttorigtvR((RRR RRRs6lib/python2.7/site-packages/mpmath/functions/bessel.pyR sV 1 % c s<Иj|Г}ИjИГЙИs┤|r3tВn|sEd|ИSИj|Гr`d|ИSИj|Г}|dkrКИj|ИS|dkrвd|ИSИj|ИSnИjИГЙ|rИj|Г}ЗЗЗfdЖ}Иj|||g|Н}n-ЗЗЗfdЖ}Иj||g|Н}|S(Niics╝ИjИjИИdИjИГddtГ}d||dd||d|dg}dИjИg|d|d||g|dg||dd|ddg||fg}|S(NR g╨?Rgр?ii(RR RR(RR RRR(RRR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRfs1+\cspИjИddtГ}Иj||dtdИjИГГ}|g|gg|dgg|dg|fgS(Ngр?RR ii(RRRR (RRR(RRR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRms((Rt ValueErrorRR(tnanR)RR#( RRRR+R,RR RR/((RRRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytbesseliPs. cKsд|sт|rtВn|s+|j||S|j|ГrI|j||S|j|Г}|d}|j|ГrЬ|dkrН|j||Sd||Sn|dkr╨t|j|ГГdr╨|j||S|j||Sn|j d7_ |j |Г\}}||j kr<|j } |j d9_ || 7}n|dkrZ|j |8_ n|j|Г\} }|j ||||Н| |j ||||Н|S(Ngр?iii (R0R)timR1R(RRtfloortninfR t nint_distancetepstcospi_sinpiR(RRRR+R,RtqtmR Rtcostsin((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytbesselyts6 % csqИs ИjSИjИГ}|dkr:ЗfdЖ}n!Иj|7_ЗЗfdЖ}Иj||g|НS(NicsЕИdd}Иdg||dg|gggd|g|f}Иdg||dg|gggd|g|f}||fS(Nii((RRtT1tT2(R(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRЯs67c sUИjdИИjИГgdddggg|dd|ggddИfgS(Nigр?gр┐ii (Rtexp(R(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRкs/(R)RR R#(RRRR,RR((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytbesselkШs cKs-|j|||Н|j|j|||НS(N(RtjR=(RRRR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pythankel1пscKs-|j|||Н|j|j|||НS(N(RRBR=(RRRR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pythankel2│scKsм|dkrS|j|Гdkr%|S|j|ГdkrE|j|S|j|Sn|jd|dtГ}d|}|j|Г|||j||dd|||НS(Nigр┐Rgр?ii(R(R)R1RRR@thyp1f1(RR R:RR,Rty((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytwhitm╖s cKsп|dkrVt|j|ГГ}|dkr1|S|dkrH|j|S|j|Sn|jd|dtГ}d|}|j|Г|||j||dd|||НS(Nigр?gр┐Rii(R$R(R)R1RRR@thyperu(RR R:RR,tgRRF((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytwhitw┼s csVИj|Г\}}Иj|Г\}}ИjИГЙИsДИj|ГdkrvИjd|g||dgГSИjИSnd||}Иj|Г\}}ymИj} zSИjd7_Иjdd||f||gdИdИjГ} | И|SWd| И_XWnИjk r*nXЗЗfdЖ}Иj|||g|НS(Nii iii tmaxtermscsиИj|Г}Иj|gddgg||d|g|g|gИf}Иj|Иgddd|gg|d|g||dgd|gИf}||fS(Nii i(tsinpiR(tatbRR>R?(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRщs>Q( t_convert_paramRR(t gammaprodR)R thypsumt NoConvergenceR#(RRMRNRR,tatypetbtypetbbtbbtypeR.R/R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pyRH╙s(" 4csFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Ncs]ИddИjИjГg|ddgg|dgdgd|dgИddfgS(Nigр?ii g°?(tsqrtR(R(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRїs(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytstruvehЁscsFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Ncs\ИddИjИjГg|ddgg|dgdgd|dgИddfgS(Nigр?ii g°?(RWR(R(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR■s(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytstruvel∙scsMИj|Гd}ИjИГЙЗЗЗfdЖ}Иj||g|НS(Nics/Иj}||}|d}||||d|d|f\}}}}Иj|Г\}} ИdkrЙ|И| g|g} }nИdkr│|И|g| g} }nИjИddГ}| ddgg||gdg||g|f} |dgg||gdg||g|f}| |fS(Niiitmultg╨┐(tmpq_1_2R8tsquare_exp_arg(R/RNtuR:ta1ta2tb1tb2tctstARRR>R?(RtwhichR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRs .0-(RORR#(RReR/RR,R((RReRs6lib/python2.7/site-packages/mpmath/functions/bessel.pyt_angerscKst|d|||НS(Ni(Rf(RR/RR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytangerjscKst|d|||НS(Ni(Rf(RR/RR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytweberescs`Иj|Гd}Иj|Гd}ИjИГЙЗЗfdЖ}Иj|||g|НS(Nic sБИj}ИjИddГ}||d||dИgdd|dgggdg|||d|||dg|ffS(NRZg╨┐ii i(R[R\(R]R/RNR(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR"s 8(RORR#(RR]R/RR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytlommels1s cs`Иj|Гd}Иj|Гd}ИjИГЙЗЗfdЖ}Иj|||g|НS(Nic s=Иj}ИjИddГ}||d||dИgdd|dgggdg|||d|||dg|f}dИg||d|g||||dg|||dggd|g|f}dИg||d|g||||dg|d||ggd|g|f}|||fS(NRZg╨┐ii ii(R[R\(R]R/RNRR>R?tT3(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR4s bXX (RORR#(RR]R/RR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytlommels2)s csFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Nc s╠Иdd}Иjd|Г\}}|Иdgd|gg|dggdd|dd|dg|f}|Иdgd|dgg|dggdd|dd|dg|f}||fS(Nigш┐iigр?g°?i(R8(RRR;R<R>R?(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRKs KO(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytberFscsFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Nc s╠Иdd}Иjd|Г\}}|Иdgd|dgg|dggdd|dd|dg|f}|Иdgd|gg|dggdd|dd|dg|f}||fS(Nigш?iig°?gр?i(R8(RRR;R<R>R?(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRXs OK(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytbeiSscsFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Nc sдИdd}Иjd|Г\}}Иjd|Г\}}dИd|g|d|dg|gggddd|d|dg|f}dИ|g|dd|dg|dgggddd|d|dg|f}dИd|g|d|dg|gggddd|dd|g|f}dИ|g|dd|dg|dgggddd|dd|g|f} |||| fS( Nig╨?gш?iiigр?g°?(R8( RRtcos1tsin1tcos2tsin2R>R?RjtT4(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyResSXRV(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytker`s csFИj|Г}ИjИГЙЗЗfdЖ}Иj||g|НS(Nc sЮИdd}Иjd|Г\}}Иjd|Г\}}|dИgd|dd|g|dgggddd|dd|g|f}|dИgd|d|g|gggddd|dd|g|f}|dИgd|d|g|gggdd|dd|dg|f}|dИgd|d|dg|dgggdd|dd|dg|f} |||| fS( Nigш?g╨?iiig°?gр?(R8( RRRnRoRpRqR>R?RjRr(RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRusVOPX(RR#(RRRR,R((RRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytkeips csИjЙЗЗfdЖ}|S(Ncsa|j}|j}|jИdГ\}}||kr;| S|И|Гf|И<|ИdSdS(Ni ii(i i(t_misc_const_cacheR tget(RtcacheR tpR/(tftname(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt f_wrappedГs (t__name__(RyR{((RyRzs6lib/python2.7/site-packages/mpmath/functions/bessel.pytc_memoБs cCs+d|jdГ|j|jdГdГS(Nii ii(tcbrttgammaR"(R((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt _airyai_C1ОscCs+d|jdГ|j|jdГdГS(Ni ii(R~RR"(R((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt _airyai_C2ТscCs.d|jddГ|j|jdГdГS(Niiii(tnthrootRR"(R((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt _airybi_C1ЦscCs*|jddГ|j|jdГdГS(Niii(RВRR"(R((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt _airybi_C2ЪscCsI|j}z.|jddГ|jdГd|j}Wd||_X| S(Nis2/3i(R tpowerRR(RR R/((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt_airybi_n2_infЮs . cCs|dkr|dkr|S|j}|j}z─|jd7_|j|d|Г|jd||Г|j}|dkr╢||jd|d|Г9}||jddГ:}n;|t|jd|d|ГГ9}||jddГ:}Wd||_X| |StВdS( NtZii iiis2/3s1/6(tmpq_1_3R RRЕRRLR$R&(RRRtntypeReRR R/((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt_airyderiv_0иs 2% csTИjИГЙ|r-Иj|Г\}}nd}ИjИГr7Иr7|r√|dkr√|dkr║ИИjkrПИjdГddИSИИjkr║ИjdГddИSn|dkr√ИИjkr┘ИSИИjkr°d|ИSq√n|rИИjks ИИjkr(dИStdГВnИrbtdtdИj ИГГГЙndЙ|r(|dkrвЗЗЗfd Ж}Иj |g|НSИdkr─tИИ||dГSЗЗЗfd Ж}Иj ||g|Н}ИjИГr!Иj |Гr!Иj|Г}n|Sn(ЗЗЗfdЖ}Иj |g|НSdS(NiRЗi iii■ sessential singularity of Ai(z)g°?cs[ИjИГdkr▒ИjИ7_Иd}d|}d|d}ИjИ8_Иj|ГdИjИjГИjИdГ}|gdgggd dgg|ffSИjИ7_Иdd}ИjИ8_tИГd}tИГ}|ИgddggggИjg|f}|gdggggИj g|f}||fSdS(Nig°?gш┐i■ iiii iii gр?(i i(ii( RR R@RWRRВRАRБtmpq_5_3RИ(RRR]tCtC1tC2R>R?(Rt extraprecR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR▌s 4(-'csNИjИ7_Иdd}ИjИ8_ИjИjИj}}}|}d}d||}d||}d||} dИg|||g|g||| g||g||| g|f} |}d||}d||}d||} dИИg|||dg|g||| g||g||| g|f}| |fS(Nii iii(R RИtmpq_2_3tmpq_4_3(RRtq13tq23tq43R^R_R`Ratb3R>R?(RRПR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRЇs$)0csTИjИГdkr░ИjИ7_Иd}d|}d|d}ИjИ8_Иj|ГdИjИjГИjИdГ}|gdgggddgg|ffSИjИ7_Иdd }ИjИ8_tИГ}tИГ}|gdggggИjg|f}И|gdggggИj g|f}||fSdS( Nig°?gш┐i■ iiiiii (ii(ii( RR R@RWRRВRАRБRРRС(RRR]RМRНRОR>R?(RRПR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRs 3('+(RROtisnormalR)R"R5R0RRRR#RКt _is_real_typeRR(RRR+R,RRЙRR/((RRПRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytairyai╛sF%%cs%ИjИГЙ|r-Иj|Г\}}nd}ИjИГrИr|r╞|dkr╞ИИjkrnИSИИjkr╞|dkrСdИS|dkrзtИГS|dkr├d|ИSq╞n|s∙ИИjkr▀ИSИИjkr∙dИSntdГВnИr3tdtdИj ИГГГЙndЙ|r∙|dkrsЗЗЗfdЖ}Иj |g|НSИdkrХtИИ||dГSЗЗЗfd Ж}Иj ||g|Н}ИjИГrЄИj |ГrЄИj|Г}n|Sn(ЗЗЗfd Ж}Иj |g|НSdS(NiRЗi ii■ sessential singularity of Bi(z)g°?csжИjИ7_Иdd}ИjИ8_tИГd}tИГ}|ИgddggggИjg|f}|gdggggИjg|f}||fS(Nii gр?ii(R RГRДRЛRИ(RRНRОR>R?(RRПR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR;s-'cs\ИjИ7_Иdd}ИjИ8_ИjИjИj}}}Иj}Иj}|}d}d||} d||} d||}dИg|||g|g| | |g||g| | |g|f}|}d||} d||} d||}dИg||d|g|g| | |g||g| | |g|f} || fS(Nii iii(R RИRРRСtmpq_1_6tmpq_5_6(RRRТRУRФtq16tq56R^R_R`RaRХR>R?(RRПR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRHs( ),csаИjИ7_Иdd}ИjИ8_tИГ}tИГ}|gdggggИjg|f}И|gdggggИjg|f}||fS(Nii i(R RГRДRРRС(RRНRОR>R?(RRПR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR[s'+(RRORЦR)R5RЖR0RRRR#RКRЧRR(RRR+R,RRЙRR/((RRПRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytairybisJ % c s>dД}dД}t|Г}|dkr9tdГВn|dkrTtdГВn|dkr╬|rЭИjЗfdЖ|dИjd |dd ГГSИjИj|dИjd |dd ГГS|dkrT|tkrT|r#ИjЗfdЖ|dИjd |dd ГГSИjИj|dИjd |dd ГГS|dkr:|tkr:|r┘dИjd |dd dИj}Иj Иj dГdГ||Г}ИjЗfd Ж|ГSdИjd |dd dИj}Иj Иj dГdГ||Г}ИjИj|ГSdS(NcSs |dddd|ddS(Nig@iii0((tt((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytUhscSs |dddd|ddS(Nig@iii0((RЮ((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRisisk cannot be less than 1is%Derivative should lie between 0 and 1csИj|dГS(Ni(RШ(R(R(s6lib/python2.7/site-packages/mpmath/functions/bessel.pytqsiiicsИj|dГS(Ni(RЭ(R(R(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRаvsyш?csИj|dГS(Ni(RЭ(R(R(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRа}s(ii(RR0tfindrootRRШtFalseRЭRtln2texpjpiR"( RReR R+tcomplexRЯRRЮRc((Rs6lib/python2.7/site-packages/mpmath/functions/bessel.pyt _airy_zerofs2 %1%1(&(&cCst|d||tГS(Ni(RжRв(RR R+((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt airyaizeroВscCst|d|||ГS(Ni(Rж(RR R+Rе((s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt airybizeroЖscsяИjИГЙИjИГrzИИjkrTИdkrAdИSИdkrTИSnИИjkrkdИStdГВnИrеtdtdИjИГГГЙndЙИjdГr├t ВnyчИjИГdkrйИdkr>t ИjИГГИjddkr>ЗЗfdЖ}Иj |gd Иjd tГSИdkrйt ИjИГГdИjddkrйЗЗfdЖ}Иj |gd Иjd tГSnWnИjk r└nXЗЗЗЗЗfd Ж}Иj |gИНS(Niisessential singularityg°?R+ig+З┘╬ўя?c s<ИjИgddgggdddggdИdffS(Ni iiii (ii(ii(R((RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRЮsRKtforce_seriesic s=ИjИgddgggdddggdИdffS(Ni iiii (ii(ii(R((RR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRвscs╙ИjИИНd}dИj}ИdkrF|d9}|d9}nИjИ7_Иdd}ИjИ8_|gdgggggdf}|ИgddgggdgИjИjg|f}||fS(Nii■ iii i i(RЭRR RСRЛ(RdRRR>R?(RRПR,ReR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyRзs !6(RtisinfR)R5R0RRRRvR&R$targRR#R RRR(RRReR,R((RRПR,ReRs6lib/python2.7/site-packages/mpmath/functions/bessel.pyt_scorerКs6% 27&cKst||d|ГS(Ni(Rм(RRR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytscorergi╡scKst||d|ГS(Ni(Rм(RRR,((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytscorerhi╣scCs||f|krB|||fd|jkrB|||fd S|jd|dГ}|jd||j|Г}|jd||j|Г}d||j|j|||d|Г}|j|Гp▀|j|ГsЇ|j|Г}n|j|f|||f<|S(Niii(R tloggammaRBR@RR3R((Rtltetat_cachetG3tG1tG2R/((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytcoulombc╜s//c sСЗЗЗfdЖ}Иj|||g|Н}|rНИj|ГrНИj|ГrНИjИГrНИjИГdkrНИj|Г}n|S(NcsуyзИjИ}Иj|ИdtГ}Иj|ddtГ}Иj||Г}|ИИj|Гgd|ddgggd|||gd|dg|f}Wn2tk r█dgdgggggdf}nX|fS(NRi■ iiii (RBRRR╢R@R0(R░R▒tjwtjwztjwz2RМR>(RRR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR╨s = %i(R#R3R(( RR░R▒RRtchopR,RR/((RRRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytcoulombf╩s6csЕИИf|krA|ИИfdИjkrA|ИИfdSЗЗЗfdЖ}Иj|dГ}Иj|f|ИИf<|S(NiicsЛИd}ИjИ}ИjdИ|ГdИjdИ|ГdИjd||ГdИjd||ГdИdИjgS(Niyр┐yр?gр?(RBRпR(tl2tjeta(RR▒R░(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyttermsхs (R tsum_accurately(RR░R▒R▓R╛R/((RR▒R░s6lib/python2.7/site-packages/mpmath/functions/bessel.pyt_coulomb_chiсs/c s▓Иj|Гs!Иj|Г}nЗЗЗfdЖ}Иj|||g|Н}|rоИj|ГrоИj|ГrоИjИГrоИjИГdkrоИj|Г}n|S(Nc s╝Иj|dГr;dgdgggggdf}|fS|d}y9Иj||Г}ИjИ}Иj|Г}Иj|Г}Иj||Г}Иj||Г} Иj|ИГ} d|И}||И| |gdd|dddgggd|||gd|dg|f}|| И| gdd|ddgggd|||gd|dg|f}||fSWn6tk r╖dgdgggggdf}|fSXdS(Niii ii■ (RR└RBR<R;R╢R@R0( R░R▒R>R╝tchiR╖RcRbRНRОR]RR?(RRR(s6lib/python2.7/site-packages/mpmath/functions/bessel.pyR°s*! .))) !i(t_imRR#( RR░R▒RRR║R,RR/((RRRs6lib/python2.7/site-packages/mpmath/functions/bessel.pytcoulombgёs6cCs5d|d}|dkrE|rEd|d|d|jd}n|dkr||r|d|d|d|jd}n|dkr▓|r▓d|d|d|jd}n|dkrш|rшd|d|d|jd}n|sм|}|dd|}d|dd|ddd|d} d |dd |dd|dd d|d} d|dd|dd|dd|ddd|d}n|r||}|dd|}dd|dd|ddd|d} d d |dd|dd|dd d|d} dd|dd|dd|dd|ddd|d}n||| | |g}|} d}xbtdt|ГГD]K}t||Гt||dГkrю| ||7} q│t||Г}q│W|t|Гdkr+t|dГ}n| |fS( sj 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. iiiiii№ iiiр iSi╓i├iii└ i%i XiO2iu╚_iiiRi ii▀i╤ i,Жiil"qi╗QYgi (RR!tlenR$(RtkindtprimeR/R:R]RNts1ts2ts3ts4ts5R╛Rcterrti((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytmcmahons> $$$$*6E.:I$cCsъ|dkrtdГВn|d}g}g}x▓|j|||Г}g|D]}|j||ГГ^qP}gt|dГD]:} || || ddkrВ|| || df^qВ} t| Г|kr╪| S|d}q4WdS(sю 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. isn cannot be less than 1i iN(R0tlinspacetsignR!R─(RRyRMRNRtNtpointstsignsRR═tok_intervals((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytgeneralized_bisection;s (:cCs|j||dddtГS(Ntsolvertillinoistverify(RбRв(RRytab((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytfind_in_intervalQsg{оGсzД?csДИj}t|ИjИГИj|ГГd}zB|И_ИjИГЙt|Г}t|Г}ИdkrtdГВn|dkrЪtdГВn|dkr╡tdГВn|dkrё|r▄ЗЗfdЖ} qёЗЗfdЖ} n|d kr-|rЗЗfd Ж} q-ЗЗfdЖ} n|dkrо|rо|dkrоИdkr^ИjSИdkrоd ИjИdИИd Г} tИ| | dd | fГSn||И|f|krцtИ| |||И|fГSt И||И|Г\} }||kr.tИ| | || |fГS|dkrJ|rJd}n|dkre|red }n|d krБ|rБd}n|d krЬ|rЬd}n|d} x╔t И||И| Г\}}||krdt И||И| dГ\}}t И| |d||| Г}x4t|ГD]&\}}||||И|df>> 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 i(Rх(RR/R:R+((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytbesseljzeroЙspcCst|d|||Г S(s╞ 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 i(Rх(RR/R:R+((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytbesselyzero√sYN(4t functionsRRRRRR2R=RARCRDRGRJRHRXRYRfRgRhRiRkRlRmRsRtR}RАRБRГRДRЖRКRШRЭRвRжRзRиRмRнRоR╢RR╗R└R├R╬R╒R┌RхRцRч(((s6lib/python2.7/site-packages/mpmath/functions/bessel.pytsxB## [K +" ' 5q