U :qLe|-@sjddlmZmZeddZeddZeddZedd Zed;d d Zed dZeddZ edddZ ed?ddZ ed@ddZedAddZedBdd Zed!d"ZedCd#d$Zed%d&Zed'd(Zed)d*Zed+d,Zed-d.Zed/d0Zed1d2Zed3d4Zed5d6Zed7d8Zed9d:ZdS)D)defun defun_wrappedcCsL||d}d||j||dd|}||sH|||j}|S)Nrrrsquare_exp_argsqrtpihyp1f1_re_imjctxzz2vrt/home/p21-0144/sympy/latex2sympy2solve-back-end/sympyEq/lib/python3.8/site-packages/mpmath/functions/expintegrals.py _erf_complexs  " rcCs|||dkrL||}|j|dd}||||j|dd|}nd||}||sxd| ||j }|S)NrTexactrr) rer fnegexpr r hyperurrrr)rrrZnz2rrrr _erfc_complex s & rcCs|||}||r6z ||WStk r4YnX||rr|jsrzt|||jWStk rpYnX||SN) convert _is_real_type_erfNotImplementedError_is_complex_typeimagtyperealrrrrrrerfs   r*cCs|||}||r6z ||WStk r4YnX||rr|jsrzt|||jWStk rpYnX||Sr ) r!r"_erfcr$r%r&r'r(rr)rrrerfc's   r,FcCsb|jdd}|r6|j|||d}|j|j||d}n|j|||d}|dkr^|j||dd}|S)N)precrTr)r/fmulfdivone)rrmult reciprocalr/rrrrr 6sr cCsR|s|S||}d||j||dd|}||sN|||j}|S)Nrrrr rrrrerfiBs " r5cs}|ks"|dks"|dkr,dS|s8SdkrFjSdkrTjStdkrvddd}nHdjtdd}||d}j d 7_ fd d |S) Nrrz*erfinv(x) is defined only for -1 <= x <= 1g?goe1?rgϷ?r cs|Sr )r*)trxrr]zerfinv..) r bad_domaininfninfabslnr signr r/findroot)rr9xreaurr8rerfinvLs     (rFcCs>||}|||d d|d||d|jS)Nr)r!rr r )rr9musigmarrrnpdf_s rJcCsD||||d}|dkr.|| dSd||dSdS)NrrGr)r r,r*)rr9rHrIrDrrrncdfdsrKc Cs||krd}n|sV|dkr0|dkr0|||}q||||d||d||}n||\}}|dkr||j kr|j } |jd9_|| 7}n|dkr|j|8_||||d||d|} ||||d||d|} | | |}|r||||}|S)NrGrr)betahyp2f1 nint_distancer/eps) rrDbx1x2 regularizedrmdhs1s2rrrbetaincls(&     rZNcCs t|}||}|dkr&|j}d}n||}||jk}|dkrN|j}d}n||}||jk}|s|s|r||dkr|jS||dkr|jS|jS||S||kr|jS||||kr||||| S|r|r| |||| S|r| |||S|r| |||SdS)NFrG) boolr!zeror=rr2nangammagammainc_gamma3 _upper_gamma _lower_gamma)rrrDrQrTZlower_modifiedZupper_modifiedrrrr_s<      r_csP|rt|jS|g|jddfdd}||gS)NTrcs2|gd|dggdgd|gf}|fS)Nrrr)rT1GrQrZnegbrrrWs,z_lower_gamma..h)isnpintr'r=r hypercomb)rrrQrTrWrrerrbs   rbcsH|rzl|rf|r*t|jWSj}z*jd7_||WWS|_Xn|WSWntk rYnX|dkrdkr|dS|dkr̈dksdkr|dSjdd |g|z&fd d }j ||gdd WSj k rBfd d } ||gYSXdS)Nr6rrrGryy?Trcs4|d}gd|ggd| ggdfgSNrrc)rrrfrDrZnegarrrWsz_upper_gamma..h) force_seriescsRgd|dg|gggdf} |gd|dggdgd|gf}||fS)NrrGrrc)rrdT2rkrrrWs.) isintrgr'r\r/_gamma_upper_intr^r$rrh NoConvergence)rrrDrTorigrWrrkrras2      rac Cs||}|r|r|jSz|jd7_|j|||d}|j|||d}||}||t||||dkr||W>S|s|j|d||d}|j|d||d}||}|WSW5|jd8_XtdS)N)rTirG)rgr\r/r_magmaxr$) rrrDrQrTZpolerdrmRrrrr`s" $ r`cCs||r8||r8z|||WStk r6YnX||sL||rT||S||jkrfd|S|dkr||dkrt||jS|j|dS|dkr| | |S|dkr| | |d|dS||d| d||S)NrrGrr) rnr" _expint_intr$isnanr=rr'r2rr_)rnrrrrexpints$ rycCsT|r.|dkr|jS|||||jS|s6|S|dkrD|jS|||S)Nrr)r\eir@ln2r>)rroffsetrrrlisr}cCs0z ||WStk r*||YSXdSr )_eir$ _ei_genericr)rrrrzs rzcCs||jkr|S||jkr|jS||dkrzr|j|}|||jddgg||jdd|}||}|dkr||j |j 7}|dkr||j |j 8}|WS|j k rYnX|| dddd||j }||r|d||||j|7}n||t|7}|S)NrT)maxtermsrlrGr?)r=r>r\rsr2rhyperr/rr rrphyp2f2eulerlogr?)rrrjrimrrrr&s4     $rcCs2z ||WStk r,|d|YSXdSri)_e1r$ryr)rrre1Bs rcCs0z ||WStk r*||YSXdSr )_cir$ _ci_genericr)rrrciIs rcCs||r.||jkr|jS||jkr.|jdS|j|j|dd}|j|dd}d||||}| |}| |}|dkr|dkr||jd7}|dkr||jd8}|dkr|dkr||jd7}|dkr||jd8}| |r|dkr| |}|S)N?TrrrG?) isinfr=r\r>r r0rrrzrrr")rrjznjzrzrealzimagrrrrPs.      rcCs0z ||WStk r*||YSXdSr )_sir$ _si_genericr)rrrsids rcCs||r2||jkrd|jS||jkr2d|jS||dkr|j|j|dd}|j|dd}d||||}| |}|dkr|d|j8}|dkr|d|j7}| |r| |}|S|| dd d d ||SdS) NrrTryrGrrgп) rr=r r>rsr0rrrzrr"hyp1f2)rrrrrrrrrrks$      rcCs|j|dd}d||||}||}||}|dkrR||jd7}n.|dkrj||jd8}n|dkr||jd7}|S)NTrrrGrr)rrzrrr )rrnzrrrrrrchis  rcCs||dkrn|j|dd}d||||}||}|dkrT|d|j8}|dkrj|d|j7}|S||dddd ||SdS) NrTrrrGrrrg?)rsrrzrr r)rrrrrrrrshis rc Cs^||jkr|dS||jkr(|dS|j|dd|ddd|jd |d d S) Nrrr)rr-r)r-rr-)r=mpfr>r rr)rrrfresnelss     rc CsP||jkr|dS||jkr(|dS||ddd|jd |ddS) Nrr)rr-r)r-rr-r)r=rr>rr r)rrrfresnelcs     r)rF)rGr)rGr)rGrF)rGNF)F)F)F)F) functionsrrrrr*r,r r5rFrJrKrZr_rbrar`ryr}rzrrrrrrrrrrrrrrsl           (  &