sg ^ d Z ddlmZ ddlmZ ddlmZmZmZ ddl m
Z
m Z m Z ddl
mZ ddlmZ ddlmZ dd lmZ dd
lmZmZmZmZ dd lmZm Z m Z m Z dd l m Z m!Z!m"Z" dd
l#m$Z$m%Z% ddl&m'Z' ddl(m)Z) ddl*m+Z+ G d de Z, G d de Z- G d de Z. G d de Z/ G d de Z0 G d d e Z1ed Z2y ) z$ Riemann zeta and related function. )Add)cacheit)ArgumentIndexError
expand_mulFunction)piIInteger)Eq)S)Dummy)sympify) bernoulli factorialgenocchiharmonic)re
unpolarifyAbs
polar_lift)log exp_polarexp)ceilingfloor)sqrt) Piecewise)Polyc 0 e Zd ZdZd ZddZd Zd Zd Zy) 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 \ }}}|dk( r t || S |j r|dk rt d }t ||z | z | }dd|z
z }t
j }t |j D ] } || |z z
}||j | z } |j || S |j rt
j }
t
j } |dkD rPt | } | |k( r| dz } || z }|| z } t t | D
cg c] }
||
| z
z ||
z |z z c}
}
nR|dk rMt | dz } || z
}|| z } t t | D
cg c] }
|| dz
|
z
z ||
z
dz
|z z c}
}
t |j |j" g \ }} t% dt& z t( z | z }|d| z z }t+ | }g }t | D ]S }
t- |||
z |z }t/ |t, r |j0 di |}|j3 |||
z |z |z z U |
| | |dz
z z t | z z S t/ |t4 r'|j d t& t( z z j s|dt( t( fv r|dk( rt ddg \ }}ny|t( k( rt ddg \ }}n_|t( k( rt ddg \ }}nD|j d dt& z t( z z }t |j |j" g \ }}t t | D
cg c]= }
t5 dt& z t( z |
z |z |z ||z z t ||
|z |z z ? c}
S t7 ||| S c c}
w c c}
w c c}
w )N r t ) argszeta
is_Integerr
r r Zeroreversed
all_coeffsdiffsubs is_RationalOner r rangepqr r r r polylog
isinstance_eval_expand_funcappendr r )selfhintszsar# r3 startrescaddmulnkmzetrootup_zetaddargsr4 args Y/var/www/html/venv/lib/python3.12/site-packages/sympy/functions/special/zeta_functions.pyr7 zlerchphi._eval_expand_func s ))1a 61: < ! == &&C%%C1u!H6FAQ1"gE!HEqQQZKQ
2EFa1"IMQd58LaAA NAEAI> 9LMacc133Z=DAqAbDF1H%Cqs8D _FG1X
6AsAvd{+ a)+++4e4Aq&!)D.1!4 45
6
QQZ W
55 5 a 166!9bd#3"@"@A"aRSQSDTBw!Qy1a!Qy1qb"az1ffQi2a(#%% (1"'(, QrT!VAXaZ\*1a4/QQ 0BB , -
- 1a E F
M4,s N' N>ANc | j \ }}}|dk( r| t ||dz | z S |dk( r$t ||dz
| |t ||| z z
|z S t )N r" )r( r r )r9 argindexr; r< r= s rK fdiffzlerchphi.fdiff sj ))1a q=2hq!a%++ +
] QAq)Ahq!Q.?,??B B$ $ c L | j }|j | r|S | S N)r7 has)r9 targetr? s rK _eval_rewrite_helperz lerchphi._eval_rewrite_helper s% $$& 776?JKrP c , | j t S rR )rU r) r9 r; r< r= kwargss rK _eval_rewrite_as_zetaz lerchphi._eval_rewrite_as_zeta s ((..rP c , | j t S rR )rU r5 rW s rK _eval_rewrite_as_polylogz!lerchphi._eval_rewrite_as_polylog s ((11rP Nr" ) __name__
__module__ __qualname____doc__r7 rO rU rY r[ r' rP rK r r s# gR?!B%/2rP r c N e Zd ZdZed ZddZd Zd Zd Z d fd Z
xZ S )
r5 a
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 |j ru|t j u r t | S |t j u r t | S |t j u rt j S |dk( rt }||v r|| S |j rt j S |j t j }|r t | S |du rK|t j u r|d|z
z S |t j u r |d|z
dz z S |j r|d|z
z S |j t t r4|s t | t j k dk( r | |t | S y y )Nr$ Fr" T) is_numberr r1 r) NegativeOne
dirichlet_etar+ _dilogtableis_zeroequalsrS r r r r )clsr< r;
dilogtablezones rK evalz polylog.eval$ s- ;;AEEzAwamm#%a(((affvv
a (]
?%a=( 9966M xx 7N
U]
AFF{!a%y amm#!a%!|#yy!a%y
55J 'Tc!fo$5Nq*Q-( ( 6O 'rP c Z | j \ }}|dk( rt |dz
| |z S t )Nr$ r" )r( r5 r )r9 rN r< r; s rK rO z
polylog.fdiffK s2 yy1 q=1q5!$Q& & rP c " |t ||d z S Nr" r )r9 r< r; rX s rK _eval_rewrite_as_lerchphiz!polylog._eval_rewrite_as_lerchphiQ s !Q"""rP c ( | j \ }}|dk( rt d|z
S |j rX|dk rSt d }|d|z
z }t | D ] }||j | z } t
| j || S t || S )Nr" r u) r( r r* r
r2 r. r r/ r5 )r9 r: r< r; rs r> _s rK r7 zpolylog._eval_expand_funcT s yy1 6AJ; <*>3CHB
::
A 3 qu~!Q$NOOAq$5==?; HF q$ )AAIDHHT!R%Z() Aw{"Wd1!QdCC# 34
s 3D6 6E
E
r\ )r ) r] r^ r_ r` classmethodrl rO rq r7 rw r
__classcell__r s @rK r5 r5 s? CJ $) $)L! #
D DrP r5 c b e Zd ZdZed
d Zd dZd dZd dZd Z d Z
d dZ d fd Z xZ
S )
r) a
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
c |t j u r | | S |t j u s|t j u rt j S |t j u rt j S |t j u rt j S |t j u rt j
S |j }|t j }|r!|j rt d|z
| |dz
z S |t j u r?|r<|j r/dt z t z |z t | z dt | z z S y y |r0|j r$|j r | | t |dz
| z
S |j r:|j r-|j du s|j du rt j S y y y )Nr" r$ F)r r1 r~ ComplexInfinityInfinityr+ r* is_nonpositiver is_evenr r r r r
is_integer)ri r< r= sints rK rl z zeta.eval s[ :q6M
!%%Z1:55L
!%%Z$$ $
!**_55L
!**_66M|| 9A A$$ QqS!$!, ,
!%%Z 2a! |il2a ! nEE "t
allq}}q6HQqS! ,, ,
\\a . . &!*:*:e*C55L +D /\rP c |dk( rS|j rG|j r;|j r/dt z t z |z t | z dt
| z z S t d|z
| |dz
z S )Nr" r$ )r is_nonnegativer r r r r r9 r< r= rX s rK _eval_rewrite_as_bernoulliz zeta._eval_rewrite_as_bernoulli sa 6allq'7'7AIIrT!VaK<)A,.!IaL.A A1a AaC((rP c \ |dk7 r| S | j d }t | ddd|z
z z
z S )Nr" r r$ )r( re r s rK _eval_rewrite_as_dirichlet_etaz#zeta._eval_rewrite_as_dirichlet_eta s7 6K IIaL Q QQZ 00rP c t d|| S ro rp r s rK rq z zeta._eval_rewrite_as_lerchphi s 1a rP c F | j d dz
j }|| S y )Nr r" rv )r9
arg_is_ones rK _eval_is_finitezzeta._eval_is_finite s- iilQ&//
!!> ! "rP c h | j d }t | j dkD r| j d nt j }|j r_|j
rt
| t |dz
| z
S |j r,|j du s|j du rt j S | S )Nr r" F)
r( lenr r1 r r r) r r r~ )r9 r: r< r= s rK r7 zzeta._eval_expand_func# s IIaL NQ .DIIaLAEE <<}}Aw!A#q!111 Q\\U%:$$-uu rP c t | j dk( r| j \ }}n| j dz \ }}|dk( r| t |dz | z S t )Nr$ r\ r" )r r( r) r )r9 rN r< r= s rK rO z
zeta.fdiff. sU tyy>Q 99DAq99t#DAq q=2d1q5!n$ $$ $rP c H t | j dk( r| j \ }}n!| j t j fz \ }} |j | \ }}|j r|j st
t t | + ||| S # t
$ r | cY S w xY w)Nr$ ) r r( r r1 r r r r r r) _eval_as_leading_term) r9 r r r r< r= r@ er s rK r zzeta._eval_as_leading_term8 s tyy>Q 99DAq99x'DAq ::a=DAq
== % %T46q$EE
# K s B B! B!rR r\ )Nr )r] r^ r_ r` r rl r r rq r r7 rO r r r s @rK r) r) sH hT 4)
1 !"
%F FrP r) c N e Zd ZdZedd ZddZej fdZ d Z
y) re a
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 |t j u r | | S |?|dk( r t d S t | }|j t sddd|z
z z
|z S y |dk( r)ddlm} t d || z
||dz dz z S t || }t ||dz dz }|j t s$|j t s|dd|z
z |z z
S y y )Nr" r$ r digamma)r r1 r r) rS 'sympy.functions.special.gamma_functionsr )ri r< r= r; r z1z2s rK rl zdirichlet_eta.evalx s :q6M 9Av1v
QA55;A!H ))
!V Gq6GAJ&!A#q)99 9
!QZ
!ac1W
vvd|BFF4LAaC2
% % %1|rP c
J ddl m} |dk( r8t t d t |d fddd|z
z z
t | z df S t t d || z
||dz dz z t |d ft || dd|z
z t ||dz dz z z
df S )Nr r r" r$ T)r r r r r r) r9 r< r= rX r s rK rY z#dirichlet_eta._eval_rewrite_as_zeta s C 6 c!fbAh /1q1Q3x<472JD1QR R#a&71:-1a0@@"Q(Kaa!A#ha!A#q)9 994@B BrP c ddl m} t t d || z
||dz dz z t |d ft d|z
| d|dz
z z df S )Nr r r$ r" T)r r r r r r r s rK _eval_rewrite_as_genocchiz'dirichlet_eta._eval_rewrite_as_genocchi s^ C#a&71:-1a0@@"Q(K!A#q!Q!A#Y/68 8rP c t d | j D r$| j t j | S y )Nc 3 4 K | ] }|j y wrR )rc ).0is rK z,dirichlet_eta._eval_evalf.. s .qq{{.s )allr( rewriter) _eval_evalf)r9 precs rK r zdirichlet_eta._eval_evalf s3 .DII. .<<%11$7 7 /rP rR r\ ) r] r^ r_ r` r rl rY r r1 r r r' rP rK re re I s5 ,\ & &$B ./UU 8
8rP re c & e Zd ZdZed Zd Zy)
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
c ddl m} t | }|t j t j
fv rt j S t |t s'||dz
z ||dz z |z dt |dz z z z S y Nr )gammar" r$ ) r r r) r r+ r1 Halfr6 r )ri r< r r; s rK rl zriemann_xi.eval sk A G 66M!T"a!e9U1Q3Z')1R!A#Y;7 7 #rP c n ddl m} ||dz
z ||dz z t | z dt |dz z z z S r )r r r) r )r9 r< rX r s rK rY z riemann_xi._eval_rewrite_as_zeta s: A!a%yqs#DG+QrAaCy[99rP N)r] r^ r_ r` r rl rY r' rP rK r r s . 8 8:rP r c " e Zd ZdZedd Zy) 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 r |Ut | }|t j u rt j S |j r |j rt j
S |j ry|t j u rt j S |dk rt j
S |j st j
S |t j u r|dv rt j S |j rt j
S |j r|dv rt j S |j dk( rt j
S y )Nr ro F) r r r~ r* r r is_Numberr+
EulerGammais_extended_negativerg r )ri rC r= s rK rl zstieltjes.eval s =
AAEEzuu || 0 0((( ;;AEEzuu Q(((\\(((affi ||# ! !$$ $ 99i << <<5 $$ $ !rP rR )r] r^ r_ r` r rl r' rP rK r r s "H % %rP r c t j t dz dz t d dz dz z
t d t dz dz t
t z t d z z
t
d dz
dz t dz dz t t
d dz
dz dz dz z t
d dz dz t dz dz t t
d dz dz dz z
dt
d z
dz t dz dz t t
d dz
dz dz z
t
d dz
dz t dz dz t t
d dz
dz dz z
t
t
t j z t dz d z z
t
t
t j z t dz d z z
dt
z
t dz d
z t
t j z z
t t
z dz t d z z
dt
z t dz d
z t
t j z z t t
z dz t d z z dt
z
dz t d dz d z t t
z t d z d z z dt dz z d z z t
t j z z
i S )
Nr$ r& r"
rM 0 ` )r r r r r
r r Catalanr' rP rK rf rf s
Ab3q619Q;&
RU1WqtCF{*
q'A+qBE6"9sDGAIq='91'>
q'A+qBE6"9sDGAIq='91'<<
T!Wa"a%(S$q'!)Q %7%::
a1a"a%(S$q'!)Q %7%:: AaiiK"a%( "
aR \BE"H
$ AAb1QYY;&AaA 6 AAb1QYY;&AaA 6
Q SVQYJqL2a4A;q=01RU72:=!)) K rP N)3r` sympy.core.addr sympy.core.cacher sympy.core.functionr r r sympy.core.numbersr r r
sympy.core.relationalr sympy.core.singletonr sympy.core.symbolr
sympy.core.sympifyr %sympy.functions.combinatorial.numbersr r r r $sympy.functions.elementary.complexesr r r r &sympy.functions.elementary.exponentialr r r #sympy.functions.elementary.integersr r (sympy.functions.elementary.miscellaneousr $sympy.functions.elementary.piecewiser sympy.polys.polytoolsr r r5 r) re r r rf r' rP rK r s * $ H H - - $ " # & Z Z P P F F > 9 : &2x 2LeDh eDX|F8 |F~P8H P8f$: $:N?% ?%D
rP