Otj^dZddlmZddlmZddlmZmZmZddl m Z ddl m Z m Z mZddlmZddlmZdd lmZdd lmZdd lmZmZmZmZdd lmZmZmZm Z dd l!m"Z"m#Z#m$Z$ddl%m&Z&m'Z'ddl(m)Z)ddl*m+Z+ddl,m-Z-GddeZ.GddeZ/GddeZ0GddeZ1GddeZ2GddeZ3edZ4dS) z$ Riemann zeta and related function. )Add)cacheit)ArgumentIndexError expand_mulDefinedFunction) fuzzy_not)piIInteger)Eq)S)Dummy)sympify) bernoulli factorialgenocchiharmonic)re unpolarifyAbs polar_lift)log exp_polarexp)ceilingfloor)sqrt) Piecewise)Polyc2eZdZdZdZd dZdZdZdZdS) lerchphia^ Lerch transcendent (Lerch phi function). Explanation =========== For $\operatorname{Re}(a) > 0$, $|z| < 1$ and $s \in \mathbb{C}$, the Lerch transcendent is defined as .. math :: \Phi(z, s, a) = \sum_{n=0}^\infty \frac{z^n}{(n + a)^s}, where the standard branch of the argument is used for $n + a$, and by analytic continuation for other values of the parameters. A commonly used related function is the Lerch zeta function, defined by .. math:: L(q, s, a) = \Phi(e^{2\pi i q}, s, a). **Analytic Continuation and Branching Behavior** It can be shown that .. math:: \Phi(z, s, a) = z\Phi(z, s, a+1) + a^{-s}. This provides the analytic continuation to $\operatorname{Re}(a) \le 0$. Assume now $\operatorname{Re}(a) > 0$. The integral representation .. math:: \Phi_0(z, s, a) = \int_0^\infty \frac{t^{s-1} e^{-at}}{1 - ze^{-t}} \frac{\mathrm{d}t}{\Gamma(s)} provides an analytic continuation to $\mathbb{C} - [1, \infty)$. Finally, for $x \in (1, \infty)$ we find .. math:: \lim_{\epsilon \to 0^+} \Phi_0(x + i\epsilon, s, a) -\lim_{\epsilon \to 0^+} \Phi_0(x - i\epsilon, s, a) = \frac{2\pi i \log^{s-1}{x}}{x^a \Gamma(s)}, using the standard branch for both $\log{x}$ and $\log{\log{x}}$ (a branch of $\log{\log{x}}$ is needed to evaluate $\log{x}^{s-1}$). This concludes the analytic continuation. The Lerch transcendent is thus branched at $z \in \{0, 1, \infty\}$ and $a \in \mathbb{Z}_{\le 0}$. For fixed $z, a$ outside these branch points, it is an entire function of $s$. Examples ======== The Lerch transcendent is a fairly general function, for this reason it does not automatically evaluate to simpler functions. Use ``expand_func()`` to achieve this. If $z=1$, the Lerch transcendent reduces to the Hurwitz zeta function: >>> from sympy import lerchphi, expand_func >>> from sympy.abc import z, s, a >>> expand_func(lerchphi(1, s, a)) zeta(s, a) More generally, if $z$ is a root of unity, the Lerch transcendent reduces to a sum of Hurwitz zeta functions: >>> expand_func(lerchphi(-1, s, a)) zeta(s, a/2)/2**s - zeta(s, a/2 + 1/2)/2**s If $a=1$, the Lerch transcendent reduces to the polylogarithm: >>> expand_func(lerchphi(z, s, 1)) polylog(s, z)/z More generally, if $a$ is rational, the Lerch transcendent reduces to a sum of polylogarithms: >>> from sympy import S >>> expand_func(lerchphi(z, s, S(1)/2)) 2**(s - 1)*(polylog(s, sqrt(z))/sqrt(z) - polylog(s, sqrt(z)*exp_polar(I*pi))/sqrt(z)) >>> expand_func(lerchphi(z, s, S(3)/2)) -2**s/z + 2**(s - 1)*(polylog(s, sqrt(z))/sqrt(z) - polylog(s, sqrt(z)*exp_polar(I*pi))/sqrt(z))/z The derivatives with respect to $z$ and $a$ can be computed in closed form: >>> lerchphi(z, s, a).diff(z) (-a*lerchphi(z, s, a) + lerchphi(z, s - 1, a))/z >>> lerchphi(z, s, a).diff(a) -s*lerchphi(z, s + 1, a) See Also ======== polylog, zeta References ========== .. [1] Bateman, H.; Erdelyi, A. (1953), Higher Transcendental Functions, Vol. I, New York: McGraw-Hill. Section 1.11. .. [2] https://dlmf.nist.gov/25.14 .. [3] https://en.wikipedia.org/wiki/Lerch_transcendent c |j\dkrtSjrdkrtd}t |z z|dd|z z }t j}tD]"}|||zz }|| |z}#| |Sj rt j}t j }dkrKtkrdzz z}tfdtD}nHdkrBt dzz z}tfdtD}t jjg\}t%dt&zt(zz } dz z} t+| } g} tD][} t-| | z| zt/t,r jd i|| | | z| z|zz \||dz zzt| zzSt/t4r"jdt&t(zz j sdt(t( fvrdkrt ddg\nt(krt ddg\n_t( krt ddg\n>jddt&zt(zz }t |j|jg\tfd tDSt7S) Nrtc4g|]}|z z |zzz Sr&.0kanszs d/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/sympy/functions/special/zeta_functions.py z.lerchphi._eval_expand_func..s0EEEqQQZKQ 2EEEc>g|]}dz |z z|z dz zz Sr#r&r's r.r/z.lerchphi._eval_expand_func..s6LLLaAA NAEAI>9LLLr0cg|]G}tdtztz|zzz zz t|zz zHS)r3)rr r zeta)r(r)r*pqr,s r.r/z.lerchphi._eval_expand_func..se,,,QrT!VAXaZ\**1a4/QQ 0B0BB,,,r0r&)argsr7 is_Integerrrr Zeroreversed all_coeffsdiffsubs is_RationalOnerrranger8r9rr r rpolylog isinstance_eval_expand_funcappendrr!)selfhintsr$startrescaddmulmzetrootup_zetaddargsr)argr*r+r8r9r,r-s @@@@@@r.rFzlerchphi._eval_expand_funcs)1a 661::  < "AFFc Aa!er]A&&Aq1uIE&Callnn-- ( (qw%**Q--88Aq>> ! = 6 &C%C1uu!HH66FAQ1"gEEEEEEEE!HHEEEFa1"IIMQdLLLLLLL588LLLMac13Z==DAqAbDF1H%%Cqs8D__FG1XX 6 6AsAvd{++a))5++44e44Aq&!)D.1!445555QQZW 55 5 a   -16!9bd#3"@ -A"aRSQSDTDTBww!Qyy11a!Qyy11qb"azz11fQi2a(#%((1,,,,,,,"'((,,,- -1a   r0r#c|j\}}}|dkr| t||dz|zS|dkr,t||dz ||t|||zz |z St)Nr#)r:r!r)rHargindexr-r,r*s r.fdiffzlerchphi.fdiffsu)1a q==2hq!a%+++ + ]]QAq))Ahq!Q.?.?,??B B$ $r0c\|}||r|S|SN)rFhas)rHtargetrKs r._eval_rewrite_helperzlerchphi._eval_rewrite_helpers/$$&& 776?? JKr0c 6|tSrZ)r]r7rHr-r,r*kwargss r._eval_rewrite_as_zetazlerchphi._eval_rewrite_as_zetas((...r0c 6|tSrZ)r]rDr_s r._eval_rewrite_as_polylogz!lerchphi._eval_rewrite_as_polylogs((111r0Nr2) __name__ __module__ __qualname____doc__rFrXr]rarcr&r0r.r!r!srggR?!?!?!B%%%%///22222r0r!cTeZdZdZedZd dZdZdZdZ d fd Z xZ S) rDa Polylogarithm function. Explanation =========== For $|z| < 1$ and $s \in \mathbb{C}$, the polylogarithm is defined by .. math:: \operatorname{Li}_s(z) = \sum_{n=1}^\infty \frac{z^n}{n^s}, where the standard branch of the argument is used for $n$. It admits an analytic continuation which is branched at $z=1$ (notably not on the sheet of initial definition), $z=0$ and $z=\infty$. The name polylogarithm comes from the fact that for $s=1$, the polylogarithm is related to the ordinary logarithm (see examples), and that .. math:: \operatorname{Li}_{s+1}(z) = \int_0^z \frac{\operatorname{Li}_s(t)}{t} \mathrm{d}t. The polylogarithm is a special case of the Lerch transcendent: .. math:: \operatorname{Li}_{s}(z) = z \Phi(z, s, 1). Examples ======== For $z \in \{0, 1, -1\}$, the polylogarithm is automatically expressed using other functions: >>> from sympy import polylog >>> from sympy.abc import s >>> polylog(s, 0) 0 >>> polylog(s, 1) zeta(s) >>> polylog(s, -1) -dirichlet_eta(s) If $s$ is a negative integer, $0$ or $1$, the polylogarithm can be expressed using elementary functions. This can be done using ``expand_func()``: >>> from sympy import expand_func >>> from sympy.abc import z >>> expand_func(polylog(1, z)) -log(1 - z) >>> expand_func(polylog(0, z)) z/(1 - z) The derivative with respect to $z$ can be computed in closed form: >>> polylog(s, z).diff(z) polylog(s - 1, z)/z The polylogarithm can be expressed in terms of the lerch transcendent: >>> from sympy import lerchphi >>> polylog(s, z).rewrite(lerchphi) z*lerchphi(z, s, 1) See Also ======== zeta, lerchphi c|jru|tjurt|S|tjurt | S|tjur tjS|dkrt}||vr||S|jr tjS| tj}|rt|S|dur>|tjur|d|z z S|tjur |d|z dzz S|jr|d|z z S| ttr<|s!t|tjkdkr||t|SdSdS)Nr3Fr#T) is_numberr rBr7 NegativeOne dirichlet_etar< _dilogtableis_zeroequalsr[rrrr)clsr,r- dilogtablezones r.evalz polylog.eval%sf ; )AEzzAwwam##%a((((afv a(]]  ??%a=( 9 6Mxx  !77N U]] AF{{!a%y am##!a%!|#y !!a%y  55J ' ' )T )c!ffo$5N5N3q*Q--(( ( ) )5N5Nr0r#c\|j\}}|dkrt|dz ||z St)Nr3r#)r:rDr)rHrWr,r-s r.rXz polylog.fdiffLs5y1 q==1q5!$$Q& &  r0c *|t||dzSNr#r!)rHr,r-r`s r._eval_rewrite_as_lerchphiz!polylog._eval_rewrite_as_lerchphiRs!Q""""r0c L|j\}}|dkrtd|z  S|jrk|dkretd}|d|z z }t | D]}|||z}t |||St||S)Nr#ru) r:rr;rrCr?rr@rD)rHrIr,r-rzrJ_s r.rFzpolylog._eval_expand_funcUsy1 66AJJ;  < 0AFFc Aq1uIEA2YY ( (%**Q--e$$))!Q// /q!}}r0c2|jd}|jrdSdS)Nr#T)r:rn)rHr-s r. _eval_is_zerozpolylog._eval_is_zeroas& IaL 9 4  r0rcddlm}|j\}}||d}|tjur.||dt|jrdnd}|j r | |\} } n#ttf$r|cYSwxYw| j rt|| z } |||z|} |||||} | tjur| S| }|g}t%d| D]"}|| z}||||zz #t)|| zSt+t,|||||S)Nr)Order-+)dirr3)sympy.series.orderrr:r@r NaNlimitr is_negativernleadterm ValueErrorNotImplementedError is_positiver _eval_nseriesremoveOr<rCrGrsuperrD)rHxr+logxcdirrnur-z0r{rnewnortermr,r) __class__s r.rzpolylog._eval_nseriesfs,,,,,, A VVAq\\ ;;A"T((*>#G33CHHB : # A33 34      #qu~~E!Q$NNOOAq$55==??;;HFq$))AAIDHHT!R%Z((((Aw{"Wd##11!QdCCCs,BBBr2)r) rdrerfrg classmethodrsrXrxrFr}r __classcell__rs@r.rDrDsCCJ$)$)[$)L!!!! ###    DDDDDDDDDDr0rDcfeZdZdZed dZd dZd dZd dZdZ d Z d d Z fd Z xZ S)r7a Hurwitz zeta function (or Riemann zeta function). Explanation =========== For $\operatorname{Re}(a) > 0$ and $\operatorname{Re}(s) > 1$, this function is defined as .. math:: \zeta(s, a) = \sum_{n=0}^\infty \frac{1}{(n + a)^s}, where the standard choice of argument for $n + a$ is used. For fixed $a$ not a nonpositive integer the Hurwitz zeta function admits a meromorphic continuation to all of $\mathbb{C}$; it is an unbranched function with a simple pole at $s = 1$. The Hurwitz zeta function is a special case of the Lerch transcendent: .. math:: \zeta(s, a) = \Phi(1, s, a). This formula defines an analytic continuation for all possible values of $s$ and $a$ (also $\operatorname{Re}(a) < 0$), see the documentation of :class:`lerchphi` for a description of the branching behavior. If no value is passed for $a$ a default value of $a = 1$ is assumed, yielding the Riemann zeta function. Examples ======== For $a = 1$ the Hurwitz zeta function reduces to the famous Riemann zeta function: .. math:: \zeta(s, 1) = \zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s}. >>> from sympy import zeta >>> from sympy.abc import s >>> zeta(s, 1) zeta(s) >>> zeta(s) zeta(s) The Riemann zeta function can also be expressed using the Dirichlet eta function: >>> from sympy import dirichlet_eta >>> zeta(s).rewrite(dirichlet_eta) dirichlet_eta(s)/(1 - 2**(1 - s)) The Riemann zeta function at nonnegative even and negative integer values is related to the Bernoulli numbers and polynomials: >>> zeta(2) pi**2/6 >>> zeta(4) pi**4/90 >>> zeta(0) -1/2 >>> zeta(-1) -1/12 >>> zeta(-4) 0 The specific formulae are: .. math:: \zeta(2n) = -\frac{(2\pi i)^{2n} B_{2n}}{2(2n)!} .. math:: \zeta(-n,a) = -\frac{B_{n+1}(a)}{n+1} No closed-form expressions are known at positive odd integers, but numerical evaluation is possible: >>> zeta(3).n() 1.20205690315959 The derivative of $\zeta(s, a)$ with respect to $a$ can be computed: >>> from sympy.abc import a >>> zeta(s, a).diff(a) -s*zeta(s + 1, a) However the derivative with respect to $s$ has no useful closed form expression: >>> zeta(s, a).diff(s) Derivative(zeta(s, a), s) The Hurwitz zeta function can be expressed in terms of the Lerch transcendent, :class:`~.lerchphi`: >>> from sympy import lerchphi >>> zeta(s, a).rewrite(lerchphi) lerchphi(1, s, a) See Also ======== dirichlet_eta, lerchphi, polylog References ========== .. [1] https://dlmf.nist.gov/25.11 .. [2] https://en.wikipedia.org/wiki/Hurwitz_zeta_function Nc|tjur ||S|tjus|tjur tjS|tjur tjS|tjur tjS|tjur tjS|j}| tj}|r |jrtd|z ||dz z S|tjurF|r@|j r;dtztz|z t|zdt|zz SdSdS|r-|jr&|j r||t|dz |z S|jr%|jr |jdus |jdurtjSdSdSdS)Nr#r3F)r rBrComplexInfinityInfinityr<r;is_nonpositiveris_evenr r rrr is_integer)rpr,r*sints r.rsz zeta.evals ::3q66M !%ZZ1::5L !%ZZ$ $ !*__5L !*__6M| 9A  A$ QqS!$$!, , !%ZZ F  F2a! |ill2a ! nEE F F F F  al q} 3q66HQqS!,,, , \ a. &&!*:e*C*C5L    *C*Cr0r#c |dkrN|jrG|jr@|jr9dtztz|z t |zdt |zz St d|z ||dz z S)Nr#r3)ris_nonnegativerr r rrrHr,r*r`s r._eval_rewrite_as_bernoullizzeta._eval_rewrite_as_bernoullisi 66al6q'76AI6rT!VaK<)A,,.!IaLL.A A1a  AaC((r0c b|dkr|S|jd}t|ddd|z zz z S)Nr#rr3)r:rlrs r._eval_rewrite_as_dirichlet_etaz#zeta._eval_rewrite_as_dirichlet_etas: 66K IaLQQQZ00r0c $td||Srvrwrs r.rxzzeta._eval_rewrite_as_lerchphis1a   r0cFt|jddz jS)Nrr#)rr:rn)rHs r._eval_is_finitezzeta._eval_is_finites$)A,*3444r0c .|jd}t|jdkr |jdn tj}|jrO|jr#t |t|dz |z S|jr|jdus |jdur tj S|S)Nrr#F) r:lenr rBrrr7rrr)rHrIr,r*s r.rFzzeta._eval_expand_func"s IaL NNQ..DIaLLAE < } 2Aww!A#q!1!111 Q\U%:%:$--u  r0ct|jdkr |j\}}n |jdz\}}|dkr| t|dz|zSt)Nr3r2r#)rr:r7r)rHrWr,r*s r.rXz zeta.fdiff-sZ ty>>Q  9DAqq9t#DAq q==2d1q5!nn$ $$ $r0cRt|jdkr |j\}}n|jtjfz\}} ||\}}n#t $r|cYSwxYw|jr|jst tt| |||S)Nr3)rr) rr:r rBrrrrrr7_eval_as_leading_term) rHrrrr,r*rLers r.rzzeta._eval_as_leading_term7s ty>>Q  9DAqq9x'DAq ::a==DAqq"   KKK  = & &% %T4  66qt$6OOOsA A&%A&rZr2)rdrerfrgrrsrrrxrrFrXrrrs@r.r7r7shhT[4)))) 1111 !!!!555   %%%%PPPPPPPPPr0r7cLeZdZdZeddZd dZejfdZ dZ dS) rla Dirichlet eta function. Explanation =========== For $\operatorname{Re}(s) > 0$ and $0 < x \le 1$, this function is defined as .. math:: \eta(s, a) = \sum_{n=0}^\infty \frac{(-1)^n}{(n+a)^s}. It admits a unique analytic continuation to all of $\mathbb{C}$ for any fixed $a$ not a nonpositive integer. It is an entire, unbranched function. It can be expressed using the Hurwitz zeta function as .. math:: \eta(s, a) = \zeta(s,a) - 2^{1-s} \zeta\left(s, \frac{a+1}{2}\right) and using the generalized Genocchi function as .. math:: \eta(s, a) = \frac{G(1-s, a)}{2(s-1)}. In both cases the limiting value of $\log2 - \psi(a) + \psi\left(\frac{a+1}{2}\right)$ is used when $s = 1$. Examples ======== >>> from sympy import dirichlet_eta, zeta >>> from sympy.abc import s >>> dirichlet_eta(s).rewrite(zeta) Piecewise((log(2), Eq(s, 1)), ((1 - 2**(1 - s))*zeta(s), True)) See Also ======== zeta References ========== .. [1] https://en.wikipedia.org/wiki/Dirichlet_eta_function .. [2] Peter Luschny, "An introduction to the Bernoulli function", https://arxiv.org/abs/2009.06743 Nc|tjur ||S|N|dkrtdSt|}|tsddd|z zz |zSdS|dkr3ddlm}td||z ||dzdz zSt||}t||dzdz }|ts(|ts|dd|z z|zz SdSdS)Nr#r3rdigamma)r rBrr7r['sympy.functions.special.gamma_functionsr)rpr,r*r-rz1z2s r.rszdirichlet_eta.evalws) ::3q66M 9Avv1vv QA55;; *A!H )) F !VV G G G G G Gq66GGAJJ&!A#q)9)99 9 !QZZ !ac1W  vvd|| &BFF4LL &AaC2 % % & & & &r0r#c ddlm}|dkrHttdt |dfddd|z zz t |zdfSttd||z ||dzdz zt |dft ||dd|z zt ||dzdz zz dfS)Nrrr#r3T)rrrrr r7rHr,r*r`rs r.raz#dirichlet_eta._eval_rewrite_as_zetasCCCCCC 66c!ffbAhh/1q1Q3x<4772JD1QRR R#a&&771::-1a0@0@@"Q((Kaa!A#ha!A#q)9)9994@BB Br0c ddlm}ttd||z ||dzdz zt |dft d|z |d|dz zz dfS)Nrrr3r#T)rrrrr rrs r._eval_rewrite_as_genocchiz'dirichlet_eta._eval_rewrite_as_genocchisCCCCCC#a&&771::-1a0@0@@"Q((K!A#q!!Q!A#Y/688 8r0ctd|jDr-|t|SdS)Nc3$K|] }|jV dSrZ)rj)r(is r. z,dirichlet_eta._eval_evalf..s$..qq{......r0)allr:rewriter7 _eval_evalf)rHprecs r.rzdirichlet_eta._eval_evalfsM ..DI... . . 8<<%%11$77 7 8 8r0rZr2) rdrerfrgrrsrar rBrrr&r0r.rlrlHs,,\&&&[&$BBBB./U8888 88888r0rlc.eZdZdZedZdZdS) riemann_xia Riemann Xi function. Examples ======== The Riemann Xi function is closely related to the Riemann zeta function. The zeros of Riemann Xi function are precisely the non-trivial zeros of the zeta function. >>> from sympy import riemann_xi, zeta >>> from sympy.abc import s >>> riemann_xi(s).rewrite(zeta) s*(s - 1)*gamma(s/2)*zeta(s)/(2*pi**(s/2)) References ========== .. [1] https://en.wikipedia.org/wiki/Riemann_Xi_function cddlm}t|}|tjtjfvr tjSt|ts+||dz z||dz z|zdt|dz zzz SdSNr)gammar#r3) rrr7r r<rBHalfrEr )rpr,rr-s r.rszriemann_xi.evalsAAAAAA GG   6M!T"" 8a!e9UU1Q3ZZ')1R!A#Y;7 7 8 8r0c ~ddlm}||dz z||dz zt|zdt|dz zzz Sr)rrr7r )rHr,r`rs r.raz riemann_xi._eval_rewrite_as_zetasOAAAAAA!a%yqs#DGG+QrAaCy[99r0N)rdrerfrgrrsrar&r0r.rrsH.88[8:::::r0rc*eZdZdZeddZdS) stieltjesa Represents Stieltjes constants, $\gamma_{k}$ that occur in Laurent Series expansion of the Riemann zeta function. Examples ======== >>> from sympy import stieltjes >>> from sympy.abc import n, m >>> stieltjes(n) stieltjes(n) The zero'th stieltjes constant: >>> stieltjes(0) EulerGamma >>> stieltjes(0, 1) EulerGamma For generalized stieltjes constants: >>> stieltjes(n, m) stieltjes(n, m) Constants are only defined for integers >= 0: >>> stieltjes(-1) zoo References ========== .. [1] https://en.wikipedia.org/wiki/Stieltjes_constants Nc|Ct|}|tjur tjS|jr|jr tjS|jr]|tjur tjS|dkr tjS|js tjS|tjur|dvr tjS|j r tjS|j r|dvr tjS|j dkr tjSdS)NrrvF) rr rr;rr is_Numberr< EulerGammais_extended_negativernr)rpr+r*s r.rszstieltjes.evals = AAEzzu | ) 0 )(( ; $AEzzu Q((\ $((afi|# ! %$ $ 9 i<  <5 $ $ ! r0rZ)rdrerfrgrrsr&r0r.rrs?""H%%%[%%%r0rctjtdzdz tddzdz z t dtdzdz t tztdzz t ddz dz tdz dz tt ddz dz dzdz zt ddz dz tdz dz tt ddzdz dzz dt dz dz tdzdz tt ddz dz dzz t ddz dz tdzdz tt ddz dz dzz t t tjztdzd z z t t tjztdzd z z dt z tdzd z t tjzz tt zdz tdzz dt ztdzd z t tjzztt zdz tdzzdt z dz tddz d z tt ztdzd z zdtdzzd z zt tjzz i S) Nr3 r5r# rV0`)r rr rr r rCatalanr&r0r.rmrms; Ab3q6619Q;& RU1WqtCFF{* q''A+qBE6"9sDGGAIq='9'91'> q''A+qBE6"9sDGGAIq='9'91'<< T!WWa"a%(S$q''!)Q%7%7%:: a1a"a%(S$q''!)Q%7%7%:: AaiK"a%( " aR \BE"H $ AAb1QY;&AaA6 AAb1QY;&AaA6 Q SVVQYJqL2a4A;q=01RU72:=!) K  r0N)5rgsympy.core.addrsympy.core.cachersympy.core.functionrrrsympy.core.logicrsympy.core.numbersr r r sympy.core.relationalr sympy.core.singletonr sympy.core.symbolrsympy.core.sympifyr%sympy.functions.combinatorial.numbersrrrr$sympy.functions.elementary.complexesrrrr&sympy.functions.elementary.exponentialrrr#sympy.functions.elementary.integersrr(sympy.functions.elementary.miscellaneousr$sympy.functions.elementary.piecewisersympy.polys.polytoolsrr!rDr7rlrrrmr&r0r.rs**$$$$$$OOOOOOOOOO&&&&&&----------$$$$$$""""""######&&&&&&ZZZZZZZZZZZZPPPPPPPPPPPPFFFFFFFFFF>>>>>>>>999999::::::&&&&&&22222222LeDeDeDeDeDoeDeDeDXzPzPzPzPzP?zPzPzPzP8P8P8P8P8OP8P8P8f$:$:$:$:$:$:$:$:N?%?%?%?%?%?%?%?%D       r0