g ddlmZddlmZmZmZmZddlmZddl m Z m Z m Z ddl mZmZddlmZmZmZmZddlmZddlmZdd lmZmZmZdd lmZmZdd l m!Z!m"Z"dd l#m$Z$m%Z%dd l&m'Z'ddl(m)Z)m*Z*m+Z+ddl,m-Z-m.Z.ddl/m0Z0m1Z1m2Z2ddl3m4Z4ddl5m6Z6m7Z7ddl8m9Z9dZ:Gdde Z;Gdde Z<Gdde Z=Gdde Z>Gdde Z?Gdd e Z@Gd!d"e ZAGd#d$e ZBd%S)&)prod)AddSDummy expand_func)Expr)FunctionArgumentIndexError PoleError) fuzzy_and fuzzy_not)RationalpiooIPowzeta)erferfcEi)re unpolarify)explog)ceilingfloor)sqrt)sincoscot) bernoulliharmonic) factorialrfRisingFactorial)as_int)mpworkprec) prec_to_dpscL t|ddS#t$rYdSwxYw)NF)strictT)r( ValueErrorns c/var/www/html/ai-engine/env/lib/python3.11/site-packages/sympy/functions/special/gamma_functions.pyintliker2sAqt uus  ##ceZdZdZdZejfZddZe dZ dZ dZ dZ d Zdd Zd Zdfd ZddZxZS)gammaa The gamma function .. math:: \Gamma(x) := \int^{\infty}_{0} t^{x-1} e^{-t} \mathrm{d}t. Explanation =========== The ``gamma`` function implements the function which passes through the values of the factorial function (i.e., $\Gamma(n) = (n - 1)!$ when n is an integer). More generally, $\Gamma(z)$ is defined in the whole complex plane except at the negative integers where there are simple poles. Examples ======== >>> from sympy import S, I, pi, gamma >>> from sympy.abc import x Several special values are known: >>> gamma(1) 1 >>> gamma(4) 6 >>> gamma(S(3)/2) sqrt(pi)/2 The ``gamma`` function obeys the mirror symmetry: >>> from sympy import conjugate >>> conjugate(gamma(x)) gamma(conjugate(x)) Differentiation with respect to $x$ is supported: >>> from sympy import diff >>> diff(gamma(x), x) gamma(x)*polygamma(0, x) Series expansion is also supported: >>> from sympy import series >>> series(gamma(x), x, 0, 3) 1/x - EulerGamma + x*(EulerGamma**2/2 + pi**2/12) + x**2*(-EulerGamma*pi**2/12 - zeta(3)/3 - EulerGamma**3/6) + O(x**3) We can numerically evaluate the ``gamma`` function to arbitrary precision on the whole complex plane: >>> gamma(pi).evalf(40) 2.288037795340032417959588909060233922890 >>> gamma(1+I).evalf(20) 0.49801566811835604271 - 0.15494982830181068512*I See Also ======== lowergamma: Lower incomplete gamma function. uppergamma: Upper incomplete gamma function. polygamma: Polygamma function. loggamma: Log Gamma function. digamma: Digamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Gamma_function .. [2] https://dlmf.nist.gov/5 .. [3] https://mathworld.wolfram.com/GammaFunction.html .. [4] https://functions.wolfram.com/GammaBetaErf/Gamma/ Tc|dkr<||jdtd|jdzSt||Nr5r)funcargs polygammar selfargindexs r1fdiffz gamma.fdiffrsH q==99TYq\**9Q ! +E+EE E$T844 4c~|jr0|tjur tjS|turtSt |r%|jrt |dz StjS|jr|j dkrt|j |j z}|jr|tj }}n)|dzx}}|dzdkr tj }n tj }|ttdd|zdz}|jr|t!t"zd|zz Sd|zt!t"z|z SdSdSdS)Nr5r) is_NumberrNaNrr2 is_positiver%ComplexInfinity is_RationalqabspOne NegativeOnerrangerr)clsargr0kcoeffs r1evalz gamma.evalxs@ = 5ae||u   5?-$S1W---,, 55A::CE ce+A2#$ae5 !A Aq5A::$%EEE$%MET%1Q3"2"2333E5$T"XX~144 !tDHH}u44; 5 5 5 5:r?c t|jd}|jrt|j|jkrt d}|j|jz}|j||jzz }|||z|t||jS|j rq| \}}|r%|jdkrt|}||z f|z}|}|j |ddi}||t||zS|j|jS)Nrxr5reevalF)r9rGrIrJrHrr8_eval_expand_funcsubsris_Add as_coeff_addr _new_rawargsr') r<hintsrOrTr0rJrQtailintparts r1rVzgamma._eval_expand_funcs0il ? X35zzCE!!#JJESUNEAceGOyyQ''99;;@@HQPSPUDVDVWWW : @**,,KE4 A,,)D0#3#T8%88D99T???4#?#?? ?ty$)$$r?cf||jdSNr)r8r9 conjugater<s r1_eval_conjugatezgamma._eval_conjugates&yy1//11222r?c|jd}|jr |jrdSt|r|dkrdS|js|jrdSdS)NrFT)r9is_nonpositive is_integerr2rE is_nonintegerr<rTs r1 _eval_is_realzgamma._eval_is_realsc IaL    5 1:: !q&&5 = AO 4  r?ch|jd}|jrdS|jrt|jSdS)NrT)r9rErfris_evenrgs r1_eval_is_positivezgamma._eval_is_positives? IaL = $4 _ $88# # $ $r?Nc :tt|SN)rloggamma)r<zlimitvarkwargss r1_eval_rewrite_as_tractablez gamma._eval_rewrite_as_tractables8A;;r?c &t|dz SNr5r%r<rorqs r1_eval_rewrite_as_factorialz gamma._eval_rewrite_as_factorialsQr?rc`|jd|d}|jr|dks#t|||S|jd|z }||dzt |jd| dzz |||SNrr5)r9limit is_Integersuper _eval_nseriesr8r&)r<rTr0logxcdirx0t __class__s r1r}zgamma._eval_nseriess Yq\  1 % %  5"''77((At44 4 IaL2  !a%  DIaL2#'!:!::II!QPTUUUr?cD|jd}||d}|jrM|jrF| }tj|z||dzz }|||z|z S|js||Stry) r9rWrerdrrLr8as_leading_term is_infiniter )r<rTr~rrOrr0ress r1_eval_as_leading_termzgamma._eval_as_leading_termsil XXa^^ = !R. !A-"499QU#3#33Ca00333 3 !99R== kkr?r5rm)rr_)__name__ __module__ __qualname____doc__ unbranchedrrF_singularitiesr> classmethodrRrVrbrhrkrrrwr}r __classcell__rs@r1r4r4"sJJXJ')N5555 55[5@%%%(333$$$       VVVVVV        r?r4cdeZdZdZd dZedZdZdZdZ fdZ d Z d Z d Z xZS) lowergammaa The lower incomplete gamma function. Explanation =========== It can be defined as the meromorphic continuation of .. math:: \gamma(s, x) := \int_0^x t^{s-1} e^{-t} \mathrm{d}t = \Gamma(s) - \Gamma(s, x). This can be shown to be the same as .. math:: \gamma(s, x) = \frac{x^s}{s} {}_1F_1\left({s \atop s+1} \middle| -x\right), where ${}_1F_1$ is the (confluent) hypergeometric function. Examples ======== >>> from sympy import lowergamma, S >>> from sympy.abc import s, x >>> lowergamma(s, x) lowergamma(s, x) >>> lowergamma(3, x) -2*(x**2/2 + x + 1)*exp(-x) + 2 >>> lowergamma(-S(1)/2, x) -2*sqrt(pi)*erf(sqrt(x)) - 2*exp(-x)/sqrt(x) See Also ======== gamma: Gamma function. uppergamma: Upper incomplete gamma function. polygamma: Polygamma function. loggamma: Log Gamma function. digamma: Digamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Incomplete_gamma_function#Lower_incomplete_gamma_function .. [2] Abramowitz, Milton; Stegun, Irene A., eds. (1965), Chapter 6, Section 5, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables .. [3] https://dlmf.nist.gov/8 .. [4] https://functions.wolfram.com/GammaBetaErf/Gamma2/ .. [5] https://functions.wolfram.com/GammaBetaErf/Gamma3/ rAcdddlm}|dkr0|j\}}tt | ||dz zzS|dkr_|j\}}t |t |zt|t||zz |gddgdd|gg|z St||Nr)meijergrAr5) sympy.functions.special.hyperrr9rrr4digammar uppergammar r<r=raros r1r>zlowergamma.fdiffs999999 q==9DAq 1 ~&&q1q5z1 1 ]]9DAq88GAJJ&Q 1a0@0@)@@'"q!fq!QiQ778 8%T844 4r?c Ftjur tjS\}}jr-jr&t }|krt |SnjrVjrO|dkrHdtztz|ztj  zzt z t |zSn<|dkr6tdtztz|zzt |zSj rtjurtjt z Stjur1t!tt#t!zSjs dzjrdz }|jrۉjrSt|t t|zt'fdt)Dzz St+t tjt!tz t t'fdt)dtjzDzz zSjstj tjz ztzt#t!zt+dz z t t'fdt)dt-ddz DzzSjr tjSdS)NrrAr5c:g|]}|zt|z Sru.0rPrTs r1 z#lowergamma.eval..Ks.LuLuLughQRSVV_`aVbVbMbLuLuLur?cng|]1}|tjz zttj|zz 2Sr)rHalfr4rs r1rz#lowergamma.eval..MsZXYXYXY~XY\]`a`f\fXghmnontwxnxhyhyXyXYXYXYr?clg|]0}|zdz ztzt|zz 1Srr4)rrPrrTs r1rz#lowergamma.eval..PsdpdpdpNOdehilmhmpqhqdrsxyzs{s{d{}BCDGHCH}I}IeIdpdpdpr?rB)rZeroextract_branch_factorrerErrrdrrrLr%rrCrKrrrr{rrMr4ris_zero)rNrrTnxr0bs `` r1rRzlowergamma.eval#s" ;;6M''))A < 5AM 5ABQww!!R((( \ 5a. 5AvvtAvax  33IqbMMAJqRTDUDUUU !VVqtAvaxz??:a#4#44 4 ; qAEzzusA2ww&afBxxDGG ,, q!A#!1 qE=[|[(||c1"gg ! .DsLuLuLuLulqrsltltLuLuLuGv.vvv$QxxAFA)>)>tBxx)G#qb''RUXYXYXYXYDIJKMNQRQWMWDXDXXYXYXYSZKZ*Z [[|q=16A:6r9#d1gg,,FuQQRU||SVY[\Z\V]V]^adpdpdpdpdpSXYZ\defhi\j\jmn\nSoSodpdpdp_qWqqq 9 6M  r?cftd|jDr|jd|}|jd|}t|5t j|d|}dddn #1swxYwYt j||S|S)Nc3$K|] }|jV dSrm is_numberrrTs r1 z)lowergamma._eval_evalf..V$..qq{......r?rr5)allr9 _to_mpmathr*r)gammaincr _from_mpmathr<precrrors r1 _eval_evalfzlowergamma._eval_evalfUs ..DI... . .  ! ''--A ! ''--A$ + +k!Q** + + + + + + + + + + + + + + +$S$// /Ks.BBBc|jd}|tjtjfvrE||jd|SdSr7r9rrNegativeInfinityr8r`rgs r1rbzlowergamma._eval_conjugate_X IaL QVQ/0 0 099TYq\3355q{{}}EE E 1 0r?c|j\}}t||||||g}|s|S|||}|jrt|j|jgS|||}t|j|jt|jgSrm) r9r _eval_is_meromorphicrWrerE is_finiter r)r<rTrsro args_meromz0s0s r1rzlowergamma._eval_is_meromorphicds y1 6 6q! < < " "1a ( ( *++    VVAq\\ < <amR\:;; ; VVAq\\", i 6K6KLMMMr?cv ddlm}|j\ |dturt |s_ zt z}t fdt|dz D}| z | zz}||z|zSt ||||S)Nr)Oc3JK|]}|zt|dzz VdS)r5N)r&)rrPrros r1rz+lowergamma._eval_aseries..zs8CC1a41a!e ,CCCCCCr?r5) sympy.series.orderrr9rhasrsumrMr| _eval_aseries) r<r0args0rTr~rrQsum_exprorrors @@r1rzlowergamma._eval_aseriesus((((((y1 8r>>!%%((>qDaRLECCCCCeAEllCCCCCH!Q$qA2w,A>A% %ww$$Qq$777r?c Bt|t||z Srm)r4rr<rrTrqs r1_eval_rewrite_as_uppergammaz&lowergamma._eval_rewrite_as_uppergammaQxx*Q****r?c ddlm}|jr |jr|S|t |S)Nrexpint)'sympy.functions.special.error_functionsrrerdrewriterr<rrTrqrs r1_eval_rewrite_as_expintz"lowergamma._eval_rewrite_as_expintsMBBBBBB < A, K||J''//777r?c2|jd}|jrdSdS)Nr5T)r9rrgs r1 _eval_is_zerozlowergamma._eval_is_zeros& IaL 9 4  r?rA)rrrrr>rrRrrbrrrrrrrs@r1rrs44n 5 5 5 5//[/bFFF NNN"88888+++888 r?rcTeZdZdZd dZdZedZdZdZ dZ d Z d Z d S) ra The upper incomplete gamma function. Explanation =========== It can be defined as the meromorphic continuation of .. math:: \Gamma(s, x) := \int_x^\infty t^{s-1} e^{-t} \mathrm{d}t = \Gamma(s) - \gamma(s, x). where $\gamma(s, x)$ is the lower incomplete gamma function, :class:`lowergamma`. This can be shown to be the same as .. math:: \Gamma(s, x) = \Gamma(s) - \frac{x^s}{s} {}_1F_1\left({s \atop s+1} \middle| -x\right), where ${}_1F_1$ is the (confluent) hypergeometric function. The upper incomplete gamma function is also essentially equivalent to the generalized exponential integral: .. math:: \operatorname{E}_{n}(x) = \int_{1}^{\infty}{\frac{e^{-xt}}{t^n} \, dt} = x^{n-1}\Gamma(1-n,x). Examples ======== >>> from sympy import uppergamma, S >>> from sympy.abc import s, x >>> uppergamma(s, x) uppergamma(s, x) >>> uppergamma(3, x) 2*(x**2/2 + x + 1)*exp(-x) >>> uppergamma(-S(1)/2, x) -2*sqrt(pi)*erfc(sqrt(x)) + 2*exp(-x)/sqrt(x) >>> uppergamma(-2, x) expint(3, x)/x**2 See Also ======== gamma: Gamma function. lowergamma: Lower incomplete gamma function. polygamma: Polygamma function. loggamma: Log Gamma function. digamma: Digamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Incomplete_gamma_function#Upper_incomplete_gamma_function .. [2] Abramowitz, Milton; Stegun, Irene A., eds. (1965), Chapter 6, Section 5, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables .. [3] https://dlmf.nist.gov/8 .. [4] https://functions.wolfram.com/GammaBetaErf/Gamma2/ .. [5] https://functions.wolfram.com/GammaBetaErf/Gamma3/ .. [6] https://en.wikipedia.org/wiki/Exponential_integral#Relation_with_other_functions rAc&ddlm}|dkr1|j\}}tt |  ||dz zzS|dkr?|j\}}t ||t |z|gddgdd|gg|zSt||r)rrr9rrrrr rs r1r>zuppergamma.fdiffs999999 q==9DAqA'''AE 2 2 ]]9DAqa##CFF*WWR!Q!QBPQ-R-RR R$T844 4r?cztd|jDr|jd|}|jd|}t|5t j||tj}dddn #1swxYwYtj||S|S)Nc3$K|] }|jV dSrmrrs r1rz)uppergamma._eval_evalf..rr?rr5) rr9rr*r)rinfrrrs r1rzuppergamma._eval_evalfs ..DI... . . 0 ! ''--A ! ''--A$ 0 0k!Q// 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0$S$// / s.!BB"Bc hddlm}jrYtjur tjSt ur tjSjr#tj rtS \}}j r-j r&t}|krt|Sn҉j rVjrO|dkrHdt zt"z|ztj zzt' z t|zSnu|dkrotdt)dt zt"z|zzz zt)dt zt"z|zzt|zzSjrYtjurj rt+  Stjurt) Stjur1t1t t3t1zSjs dzjrĉdz }|j rj rCt) t'|zt7fdt9DzStt3t1ztjtddz z zt) zt1zt7fdt9tjz DzzS|jr#|| t|dzzzSjstjtjz zt zt3t1ztdz z zt) zt7fd t9tjz Dzz Sjrj rt+  Sjr#tj rtSdSdS) Nrrr5rAc:g|]}|zt|z SrrurrPros r1rz#uppergamma.eval..s?>Q>Q>QBC?@dYq\\>Q>Q>Q>Qr?rBc~g|]9}ttj |z  |zztdz z :Sr)r4rrrrPrros r1rz#uppergamma.eval.. sY(D(D(D,-).qvgk(:(:qb1W(DuQqSzz(Q(D(D(Dr?cfg|]-}|ztzt|zdzz .Srrrs r1rz#uppergamma.eval..sQ5Q5Q5Q9:67TE!HH_uQqSQRU||5S5Q5Q5Qr?)rrrCrrDrrrrrEr4rrerrrdrrrLr%rrrKrrrr{rrM)rNrrorrr0rs `` r1rRzuppergamma.evalskBBBBBB ; $AEzzu bv  $a55$$ 88O''))A < VAM VABBww!!R((( \ Va. VAvv"uQwqy!!44Yr]]BZPQSUEVEVVV !VV88QQrT!VAXaZ01C"Qq OOJqRTDUDU4UU U ; SAF{{q}{A2waeA2wwafBxxT!WW -- S!A#!1 SE= @|F"A2ww15>Q>Q>Q>QGLQxx>Q>Q>Q9R RR!&a4Q== 8 ! AaDDF ;c1"gg EQ O"%(D(D(D(D(D16q16z1B1B(D(D(D#E!E!EF\@!61"a==AQ)???|SMAFQJ7"CAFQJ>O>O5Q5Q5Q0RRRS 9  rFF7N 9 A* 88O    r?c|jd}|tjtjfvrE||jd|SdSr7rr<ros r1rbzuppergamma._eval_conjugaterr?c:t|||Srm)rr)r<rTrs r1rzuppergamma._eval_is_meromorphic"s..tQ:::r?c Bt|t||z Srm)r4rrs r1_eval_rewrite_as_lowergammaz&uppergamma._eval_rewrite_as_lowergamma%rr?c \tt|t||z Srm)rrnrrs r1rrz%uppergamma._eval_rewrite_as_tractable(s%8A;;*Q"2"222r?c 8ddlm}|d|z |||zzS)Nrrr5)rrrs r1rz"uppergamma._eval_rewrite_as_expint+s3BBBBBBva!eQ1$$r?Nr) rrrrr>rrrRrbrrrrrrr?r1rrs>>B 5 5 5 566[6pFFF ;;;+++333%%%%%r?rcxeZdZdZedZdZdZdZdZ dZ dZ d Z dd Z ddZfdZdZxZS)r:a The function ``polygamma(n, z)`` returns ``log(gamma(z)).diff(n + 1)``. Explanation =========== It is a meromorphic function on $\mathbb{C}$ and defined as the $(n+1)$-th derivative of the logarithm of the gamma function: .. math:: \psi^{(n)} (z) := \frac{\mathrm{d}^{n+1}}{\mathrm{d} z^{n+1}} \log\Gamma(z). For `n` not a nonnegative integer the generalization by Espinosa and Moll [5]_ is used: .. math:: \psi(s,z) = \frac{\zeta'(s+1, z) + (\gamma + \psi(-s)) \zeta(s+1, z)} {\Gamma(-s)} Examples ======== Several special values are known: >>> from sympy import S, polygamma >>> polygamma(0, 1) -EulerGamma >>> polygamma(0, 1/S(2)) -2*log(2) - EulerGamma >>> polygamma(0, 1/S(3)) -log(3) - sqrt(3)*pi/6 - EulerGamma - log(sqrt(3)) >>> polygamma(0, 1/S(4)) -pi/2 - log(4) - log(2) - EulerGamma >>> polygamma(0, 2) 1 - EulerGamma >>> polygamma(0, 23) 19093197/5173168 - EulerGamma >>> from sympy import oo, I >>> polygamma(0, oo) oo >>> polygamma(0, -oo) oo >>> polygamma(0, I*oo) oo >>> polygamma(0, -I*oo) oo Differentiation with respect to $x$ is supported: >>> from sympy import Symbol, diff >>> x = Symbol("x") >>> diff(polygamma(0, x), x) polygamma(1, x) >>> diff(polygamma(0, x), x, 2) polygamma(2, x) >>> diff(polygamma(0, x), x, 3) polygamma(3, x) >>> diff(polygamma(1, x), x) polygamma(2, x) >>> diff(polygamma(1, x), x, 2) polygamma(3, x) >>> diff(polygamma(2, x), x) polygamma(3, x) >>> diff(polygamma(2, x), x, 2) polygamma(4, x) >>> n = Symbol("n") >>> diff(polygamma(n, x), x) polygamma(n + 1, x) >>> diff(polygamma(n, x), x, 2) polygamma(n + 2, x) We can rewrite ``polygamma`` functions in terms of harmonic numbers: >>> from sympy import harmonic >>> polygamma(0, x).rewrite(harmonic) harmonic(x - 1) - EulerGamma >>> polygamma(2, x).rewrite(harmonic) 2*harmonic(x - 1, 3) - 2*zeta(3) >>> ni = Symbol("n", integer=True) >>> polygamma(ni, x).rewrite(harmonic) (-1)**(n + 1)*(-harmonic(x - 1, n + 1) + zeta(n + 1))*factorial(n) See Also ======== gamma: Gamma function. lowergamma: Lower incomplete gamma function. uppergamma: Upper incomplete gamma function. loggamma: Log Gamma function. digamma: Digamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Polygamma_function .. [2] https://mathworld.wolfram.com/PolygammaFunction.html .. [3] https://functions.wolfram.com/GammaBetaErf/PolyGamma/ .. [4] https://functions.wolfram.com/GammaBetaErf/PolyGamma2/ .. [5] O. Espinosa and V. Moll, "A generalized polygamma function", *Integral Transforms and Special Functions* (2004), 101-115. c6|tjus|tjur tjS|tur|jrtn tjS|jr|jr tjS|tjur*t|tdtzdz z S|jr|t us)| ttt fvrtS|jrt|dz tjz S|jrF|\}}|dkr+t%t'tj|dSdSdS|jr|jrt-|}||krt'||S|jr6tj|dzzt/|zt1|dz|zS|tjurEtj|dzzt/|zd|dzzdz zt1|dzzSdSdSdS)NrAr5F)evaluate)rrDrrrr{rdrFrLrnrrextract_multiplicativelyrr$ EulerGammarGas_numer_denomrr:reis_nonnegativerr%rr)rNr0rorJrHnzs r1rRzpolygamma.evals ::ae5L "WW.22 . \ Va. V$ $ !-  A;;QrTQ. . Y VRCxx155a88R"IEE  M!}}q|33 M''))166&yU'K'K'KLLL  M M6 \ Va. VABBww B'''| V}qs+ill:T!A#q\\IIaf}qs+ill:a!A#hqjIDQRSTQTIIUU V V V V r?cV|jdjr|jdjrdSdSdS)Nrr5T)r9rEras r1rhzpolygamma._eval_is_reals< 9Q< #  ! (@ 4    r?c|jd}t|j|jg}t|jt |gSrt)r9r is_negativere is_complexr )r<rois_negative_integers r1_eval_is_complexzpolygamma._eval_is_complexsB IaL' (EFF!, 2E(F(FGHHHr?cp|j\}}|jr |jr |jrdS|jr |jrdSdSdSdSNTF)r9rEis_oddis_realrjr<r0ros r1rkzpolygamma._eval_is_positivesey1 = x AI ty Q] u       r?cp|j\}}|jr |jr |jrdS|jr |jrdSdSdSdSr)r9rErjrrrs r1_eval_is_negativezpolygamma._eval_is_negativesey1 = y Q] tx AI u       r?c  |j\ jr\jrT jr jdjrdz dkr5t  fdt dt dzD}n2t  fdt t  D }t z tj ztz|zzSn j r \ jrej r^ fdt t D}dkrt |z tzSt |dzzz S z dkr jr \ tj t$t' t$z z zdz z t z t fdt d Dz} dkr9t)  z |t fdt DzS dkr.sB%I%I%I*+&)E1&&&&%I%I%Ir?c6g|]}t|zSrrrs r1rz/polygamma._eval_expand_func..sB&C&C&C*+'*E1'&'&&C&C&Cr?c Rg|]#}tt|z$Sr)r:r)rrrQr0ros r1rz/polygamma._eval_expand_func..sM???'(&aX5."."*"##???r?rAc g|]M}td|ztzzz tdt|tzz zzNSr)r!rrr )rrPrJrHs r1rz/polygamma._eval_expand_func..sRZZZ#a!ebj1nq())CCB OO0C,D,DDZZZr?c g|] }d|zz  SrrrrPrs r1rz/polygamma._eval_expand_func..s!%E%E%Eqa26l%E%E%Er?c&g|] }ddz |z z Srrrs r1rz/polygamma._eval_expand_func..s&%I%I%I1a26A:&6%I%I%Ir?Fr)r9r{rrXrrMintr:rrLr%is_Mul as_two_termsrErrGrrrr"rrnrerrdiffrWrr4) r<r[r\part_1rdztrQr r0rJrHrors @@@@@@@r1rVzpolygamma._eval_expand_funcsy1 < A, x q #Xa%Aqyy"%I%I%I%I%I/4QE Q/G/G%I%I%I J!$&C&C&C&C&C/4S%[[/A/A&C&C&C!D D$QE 22Q]A5EiPQll5RSW5WWWX >>++q#9(99??????,1#e**,=,=???DAvv"Dz%/#e**<<"Dz%!a%.88U  66am6##%%DAql]R#a"fqj//%9A%==AFZZZZZeTUWXkkZZZJ\\F1uu!HHU%E%E%E%EE!HH%E%E%E FFFQ!a%LLU%I%I%I%Ia%I%I%I JJJ 77A;;QrTQ. . <5 A$4$=$=c Aq!**//!$$))!QqS11C1<'1"++5ac1EEPQrR RAr?c |jr=|jr8tj|dzzt |zt |dz|zSdSdSrt)rerErrLr%rr<r0rorqs r1_eval_rewrite_as_zetazpolygamma._eval_rewrite_as_zeta s\ < FAM F=1q5))A,,6tAE1~~E E F F F Fr?c |jrr|jrt|dz tjz Stj|dzzt |zt|dzt|dz |dzz zSdSrt)rerr$rrrLr%rrs r1_eval_rewrite_as_harmonicz#polygamma._eval_rewrite_as_harmonics < ^y ^A55}qs+ill:d1Q3ii(STUVSVXYZ[X[J\J\>\]]  ^ ^r?Nrcddlm}fd|jD\}}||}|dkrB|dz r*|t n|}||zS|||S)NrOrderc:g|]}|Sr)r)rrrTs r1rz3polygamma._eval_as_leading_term..s'888!!!$$888r?r5)rrr9containsrgetnr8)r<rTr~rrr0rors ` r1rzpolygamma._eval_as_leading_terms,,,,,,8888di8881 E!QKK 66ajj1oo6!\3q666tD6688d? "99Q?? "r?rAcx|dkr%|jdd\}}t|dz|St||NrAr5)r9r:r )r<r=r0ros r1r>zpolygamma.fdiff sB q==9RaR=DAqQUA&& &$T844 4r?cTddlm}|dtks$|jdjr|jdjs$t ||||S|jd|jd}|dkrtddzz z }d}|dkr|dz |}nOt|dzdz} fdtd| D} |t| z}|d|zz |}| ||||zSt|} | || zdzz z} t|dzdz} td| D]L} | d| z|zdz zd| z|zdz zd| zd| zdz zz } | td| z| zd| zzz z } M|dd| zzz |}|dkr|dz |}n|dkr|ddzz |}| |||z}ddz |zz|z |||S)Nrrr5rAcRg|]#}td|zd|zd|zzzz $Srr#rs r1rz+polygamma._eval_aseries..8s8JJJYqs^^qs1qs8|4JJJr?r)rrrr9r{rr|rrrrMrr}r4r#)r<r0rrTr~rNrrmlface0rPrors @r1rzpolygamma._eval_aseries's,,,,,, 8r>>1( -1Yq\-H 77((E1d;; ; IaL IaL 66AAaC AA1uuE!A#qMMQUQJ''JJJJeAqkkJJJS!W E!AqD&!$$??1a..2 2((Cquac{"BQ ##A1a[[ 2 21Q37Q;'1q15!A#!aIi!nnS(QqS11aAaCj!$$AAvvE!A#qMMaE!AqD&!$$  At,,q0A"Q$NQ&55aDAA Ar?ctd|jDsdS|jd|dz}|jd|dz}tj|r|dkr t jSt|dz5tj|r|dkrtj||}notj |dz|}tj |dz|d}|tj tj | z|zztj | z}dddn #1swxYwYtj||S)Nc3$K|] }|jV dSrmr)rrs r1rz(polygamma._eval_evalf..Qs$2211;222222r?r r5)rr9rr)isintrrFr*r:reulerrrgammarr)r<rrrorztrs r1rzpolygamma._eval_evalfPs22 22222  F IaL # #DG , , IaL # #DG , , 8A;; %166$ $ d2g   O Ox{{ OqAvvl1a((WQqS!__gac1a((bhQB72==A2N  O O O O O O O O O O O O O O O d+++sB E  EEr_r)rrrrrrRrhrrkrrVrrrr>rrrrs@r1r:r:4shhTVV[V:III 333jFFF^^^####5555'B'B'B'B'BR,,,,,,,r?r:cdeZdZdZedZdZd fd ZfdZdZ d Z d Z dd Z xZ S)rna The ``loggamma`` function implements the logarithm of the gamma function (i.e., $\log\Gamma(x)$). Examples ======== Several special values are known. For numerical integral arguments we have: >>> from sympy import loggamma >>> loggamma(-2) oo >>> loggamma(0) oo >>> loggamma(1) 0 >>> loggamma(2) 0 >>> loggamma(3) log(2) And for symbolic values: >>> from sympy import Symbol >>> n = Symbol("n", integer=True, positive=True) >>> loggamma(n) log(gamma(n)) >>> loggamma(-n) oo For half-integral values: >>> from sympy import S >>> loggamma(S(5)/2) log(3*sqrt(pi)/4) >>> loggamma(n/2) log(2**(1 - n)*sqrt(pi)*gamma(n)/gamma(n/2 + 1/2)) And general rational arguments: >>> from sympy import expand_func >>> L = loggamma(S(16)/3) >>> expand_func(L).doit() -5*log(3) + loggamma(1/3) + log(4) + log(7) + log(10) + log(13) >>> L = loggamma(S(19)/4) >>> expand_func(L).doit() -4*log(4) + loggamma(3/4) + log(3) + log(7) + log(11) + log(15) >>> L = loggamma(S(23)/7) >>> expand_func(L).doit() -3*log(7) + log(2) + loggamma(2/7) + log(9) + log(16) The ``loggamma`` function has the following limits towards infinity: >>> from sympy import oo >>> loggamma(oo) oo >>> loggamma(-oo) zoo The ``loggamma`` function obeys the mirror symmetry if $x \in \mathbb{C} \setminus \{-\infty, 0\}$: >>> from sympy.abc import x >>> from sympy import conjugate >>> conjugate(loggamma(x)) loggamma(conjugate(x)) Differentiation with respect to $x$ is supported: >>> from sympy import diff >>> diff(loggamma(x), x) polygamma(0, x) Series expansion is also supported: >>> from sympy import series >>> series(loggamma(x), x, 0, 4).cancel() -log(x) - EulerGamma*x + pi**2*x**2/12 - x**3*zeta(3)/3 + O(x**4) We can numerically evaluate the ``loggamma`` function to arbitrary precision on the whole complex plane: >>> from sympy import I >>> loggamma(5).evalf(30) 3.17805383034794561964694160130 >>> loggamma(I).evalf(20) -0.65092319930185633889 - 1.8724366472624298171*I See Also ======== gamma: Gamma function. lowergamma: Lower incomplete gamma function. uppergamma: Upper incomplete gamma function. polygamma: Polygamma function. digamma: Digamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Gamma_function .. [2] https://dlmf.nist.gov/5 .. [3] https://mathworld.wolfram.com/LogGammaFunction.html .. [4] https://functions.wolfram.com/GammaBetaErf/LogGamma/ c|jr2|jrtS|jrt t |Sn|jr~|\}}|jr`|dkrZt ttdd|z zzt |zt |dztj zz S|turtSt|tur tj S|tjur tjSdSr#)rerdrrErr4 is_rationalrrrrrrIrFrD)rNrorJrHs r1rRz loggamma.evals < U %  %588}}$ % ] U##%%DAq} Ua488a!a%j0588;eQUAFN>S>SSTTT 77I VVr\\$ $ ::5L :r?c fddlm}|jd}|jr|\}}||z}|||zz }|jr|jr||krt d}|jrKt||z |t|zz |t|dz |z|z|d|fzS|j r\t||z |t|zz ttz|zz|t||z|z |d| fz S|j rt||z S|S)NrSumrPr5) sympy.concrete.summationsr9r9rGrrErrnrrrrr)r<r[r9rorJrHr0rPs r1rVzloggamma._eval_expand_funcsU111111 IaL = +##%%DAqQAAaCA} + +1q55#JJ=+#AE??Qs1vvX5CQ A S1q\\TUWX[\Z\S]A^A^^^Y+#AE??* r?Nrc|jd|d}|jr&|j|j}||||St |||Sr_)r9rzr_eval_rewrite_as_intractabler}r|)r<rTr0r~rrfrs r1r}zloggamma._eval_nseriessk Yq\  1 % % : /1149=A??1a.. .ww$$Q4000r?c ddlm}|dtkr$t||||S|jd t  tjz z z t dtzdz z} fdtd|D}d}|dkr |d|}n|d |zz |}|t|z ||||zS)NrrrAcjg|]/}td|zd|zd|zdz zd|zdz zzz 0S)rAr5r&rs r1rz*loggamma._eval_aseries..sI P P PqYqs^^qsAaC!G}Q1q\9 : P P Pr?r5) rrrr|rr9rrrrrMrr}) r<r0rrTr~rr(r*rrors @r1rzloggamma._eval_aseriess,,,,,, 8r>>77((E1d;; ; IaL FFAJ ! #c!B$iik 1 P P P PE!QKK P P P  66a AAa1fa  ACG **1a66::r?c :tt|Srm)rr4rvs r1r<z%loggamma._eval_rewrite_as_intractables588}}r?cD|jd}|jrdS|jrdSdS)NrTF)r9rErdrs r1rhzloggamma._eval_is_reals6 IaL = 4   5  r?c|jd}|tjtjfvr'||SdSr_rrs r1rbzloggamma._eval_conjugatesD IaL QVQ/0 0 099Q[[]]++ + 1 0r?r5cd|dkrtd|jdSt||r7)r:r9r r;s r1r>zloggamma.fdiffs1 q==Q ! -- -$T844 4r?r_r)rrrrrrRrVr}rr<rhrbr>rrs@r1rnrnasllZ[&*111111 ; ; ; ; ;,,, 55555555r?rncheZdZdZdZddZdZdZdZdZ e d Z d Z d Z d ZddZd S)rat The ``digamma`` function is the first derivative of the ``loggamma`` function .. math:: \psi(x) := \frac{\mathrm{d}}{\mathrm{d} z} \log\Gamma(z) = \frac{\Gamma'(z)}{\Gamma(z) }. In this case, ``digamma(z) = polygamma(0, z)``. Examples ======== >>> from sympy import digamma >>> digamma(0) zoo >>> from sympy import Symbol >>> z = Symbol('z') >>> digamma(z) polygamma(0, z) To retain ``digamma`` as it is: >>> digamma(0, evaluate=False) digamma(0) >>> digamma(z, evaluate=False) digamma(z) See Also ======== gamma: Gamma function. lowergamma: Lower incomplete gamma function. uppergamma: Upper incomplete gamma function. polygamma: Polygamma function. loggamma: Log Gamma function. trigamma: Trigamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Digamma_function .. [2] https://mathworld.wolfram.com/DigammaFunction.html .. [3] https://functions.wolfram.com/GammaBetaErf/PolyGamma2/ c|jd}t|}td||S)Nrr/r9r+r:evalfr<rronprecs r1rzdigamma._eval_evalfT9 IaLD!!A$$u$---r?r5c`|jd}td|Sr_r9r:r>r<r=ros r1r>z digamma.fdiffY' IaLA$$&&&r?cF|jd}td|jSr_r9r:rrs r1rhzdigamma._eval_is_real] IaLA&&r?cF|jd}td|jSr_r9r:rErs r1rkzdigamma._eval_is_positivea IaLA**r?cF|jd}td|jSr_r9r:rrs r1rzdigamma._eval_is_negativeerTr?c|t}tjg|z}|||||Srm)rr:rrrr<r0rrTr~ as_polygammas r1rzdigamma._eval_aseriesis;||I..  E!))!UAt<<rhrkrrrrRrVrrfrrr?r1rr$s..^... '''''''++++++=== [111...222222r?rcneZdZdZdZddZdZdZdZdZ e d Z d Z d Z d Zd ZddZdS)trigammaa^ The ``trigamma`` function is the second derivative of the ``loggamma`` function .. math:: \psi^{(1)}(z) := \frac{\mathrm{d}^{2}}{\mathrm{d} z^{2}} \log\Gamma(z). In this case, ``trigamma(z) = polygamma(1, z)``. Examples ======== >>> from sympy import trigamma >>> trigamma(0) zoo >>> from sympy import Symbol >>> z = Symbol('z') >>> trigamma(z) polygamma(1, z) To retain ``trigamma`` as it is: >>> trigamma(0, evaluate=False) trigamma(0) >>> trigamma(z, evaluate=False) trigamma(z) See Also ======== gamma: Gamma function. lowergamma: Lower incomplete gamma function. uppergamma: Upper incomplete gamma function. polygamma: Polygamma function. loggamma: Log Gamma function. digamma: Digamma function. sympy.functions.special.beta_functions.beta: Euler Beta function. References ========== .. [1] https://en.wikipedia.org/wiki/Trigamma_function .. [2] https://mathworld.wolfram.com/TrigammaFunction.html .. [3] https://functions.wolfram.com/GammaBetaErf/PolyGamma2/ c|jd}t|}td||S)Nrr5r/rFrHs r1rztrigamma._eval_evalfrJr?r5c`|jd}td|SryrLrMs r1r>ztrigamma.fdiffrNr?cF|jd}td|jSryrPrs r1rhztrigamma._eval_is_realrQr?cF|jd}td|jSryrSrs r1rkztrigamma._eval_is_positiverTr?cF|jd}td|jSryrVrs r1rztrigamma._eval_is_negativerTr?c|t}tjg|z}|||||Srm)rr:rrKrrXs r1rztrigamma._eval_aseriess;||I.. 5 ))!UAt<<rhrkrrrrRrVrrfrrrr?r1rmrms..^... '''''''++++++=== [111///222222r?rmc@eZdZdZdZd dZedZdZdZ dS) multigammaa The multivariate gamma function is a generalization of the gamma function .. math:: \Gamma_p(z) = \pi^{p(p-1)/4}\prod_{k=1}^p \Gamma[z + (1 - k)/2]. In a special case, ``multigamma(x, 1) = gamma(x)``. Examples ======== >>> from sympy import S, multigamma >>> from sympy import Symbol >>> x = Symbol('x') >>> p = Symbol('p', positive=True, integer=True) >>> multigamma(x, p) pi**(p*(p - 1)/4)*Product(gamma(-_k/2 + x + 1/2), (_k, 1, p)) Several special values are known: >>> multigamma(1, 1) 1 >>> multigamma(4, 1) 6 >>> multigamma(S(3)/2, 1) sqrt(pi)/2 Writing ``multigamma`` in terms of the ``gamma`` function: >>> multigamma(x, 1) gamma(x) >>> multigamma(x, 2) sqrt(pi)*gamma(x)*gamma(x - 1/2) >>> multigamma(x, 3) pi**(3/2)*gamma(x)*gamma(x - 1)*gamma(x - 1/2) Parameters ========== p : order or dimension of the multivariate gamma function See Also ======== gamma, lowergamma, uppergamma, polygamma, loggamma, digamma, trigamma, sympy.functions.special.beta_functions.beta References ========== .. [1] https://en.wikipedia.org/wiki/Multivariate_gamma_function TrAc ddlm}|dkrV|j\}}td}||||t d|d|z dz z|d|fzSt ||)Nrr8rArPr5)r:r9r9rr8r:r )r<r=r9rTrJrPs r1r>zmultigamma.fdiff"s111111 q==9DAqc A99Q??33yAQ M'B'BQ1I#N#NN N$T844 4r?cddlm}|jdus |jdurt dt d}t ||dz zdz z|t|d|z dz z|d|fzS) Nr)ProductFz+Order parameter p must be positive integer.rPr5rA) sympy.concrete.productsr~rErer.rrr4doit)rNrTrJr~rPs r1rRzmultigamma.eval+s333333 =E ! !Q\U%:%:JKK K #JJQAYq[!''%QUAI *>*>+,a)#5#556:dff =r?cf|j\}}|||Srm)r9r8r`)r<rTrJs r1rbzmultigamma._eval_conjugate4s)y1yy***r?c|j\}}d|z}|jr ||dz kdurdSt|r ||dz krdS||dz ks|jrdSdS)NrAr5TF)r9rer2rf)r<rTrJys r1rhzmultigamma._eval_is_real8suy1 aC < Q1q5\d225 1:: 1Q<<5 A;;!/;4 ;r?Nr) rrrrrr>rrRrbrhrr?r1r{r{ss77pJ5555==[=+++r?r{N)Cmathr sympy.corerrrrsympy.core.exprrsympy.core.functionr r r sympy.core.logicr r sympy.core.numbersrrrrsympy.core.powerr&sympy.functions.special.zeta_functionsrrrrr$sympy.functions.elementary.complexesrr&sympy.functions.elementary.exponentialrr#sympy.functions.elementary.integersrr(sympy.functions.elementary.miscellaneousr(sympy.functions.elementary.trigonometricr r!r"%sympy.functions.combinatorial.numbersr#r$(sympy.functions.combinatorial.factorialsr%r&r'sympy.utilities.miscr(mpmathr)r*mpmath.libmp.libmpfr+r2r4rrr:rnrrmr{rr?r1rs111111111111 GGGGGGGGGG11111111222222222222 777777AAAAAAAAAA????????;;;;;;;;>>>>>>>>999999BBBBBBBBBBEEEEEEEESSSSSSSSSS''''''++++++uuuuuHuuuxmmmmmmmm`_%_%_%_%_%_%_%_%Lj,j,j,j,j,j,j,j,Z @5@5@5@5@5x@5@5@5FZ2Z2Z2Z2Z2hZ2Z2Z2|]2]2]2]2]2x]2]2]2JYYYYYYYYYYr?