*h^ S r SSKJr SSKJr SSKJrJrJr SSK J
r
SSK J r J
r
Jr SSKJr SSKJr SS KJr SS
KJr SS KJrJrJrJr SS K J r J r J r J r SS
K!J"r"J#r#J$r$ SSK%J&r&J'r' SSK(J)r) SSK*J+r+ SSK,J-r- " S S\5 r. " S S\5 r/ " S S\5 r0 " S S\5 r1 " S S\5 r2 " S S \5 r3\S 5 r4g ) 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)Polyc : \ rS rSrSrS rS
S jrS rS rS r Sr
g ) lerchphi aB
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 t U R u p#nUS:X a [ X45 $ UR ( a US:: a [ S5 n[ XT- U* - U5 nSSU-
- n[
R n[ UR 5 5 H n XU- -
nXWR U5 - nM UR XR5 $ UR ( Ga [
R n
[
R n US: aN [ U5 n X:X a U S- n XL- nX,* - n [ [ U 5 V
s/ s H oX-
- * XM- U- - PM sn
6 n
ORUS:: aL [ U* 5 S- n XL-
nX,- n [ [ U 5 V
s/ s H oU S-
U
-
- XM-
S-
U- - PM sn
6 n
[ UR UR" /5 u p[% S[& - [( - U - 5 nUSU - - n[+ U5 n/ n[ U 5 HX n
[- X?U
- U- 5 n[/ U[, 5 ( a UR0 " S0 UD6nUR3 UUU
- U- U- - 5 MZ XXS-
- - [ U6 - - $ [/ U[4 5 ( a, UR S [& [( - - R ( d US[( [( * 4; a US:X a [ SS/5 u nnO{U[( :X a [ SS/5 u nnO`U[( * :X a [ SS/5 u nnODUR S S[& - [( - - n[ UR UR" /5 u nn[ [ U5 V
s/ s H= n
[5 S[& - [( - U
- U- U- 5 UU- - [ X=U- U- 5 - PM? sn
6 $ [7 X#U5 $ s sn
f s sn
f s sn
f )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% r5 startrescaddmulnkmzetrootup_zetaddargsr6 args c/home/mande/repo/quber/.venv/lib/python3.13/site-packages/sympy/functions/special/zeta_functions.pyr9 lerchphi._eval_expand_func sn ))a 6: < ! === &&C%%C1u!H6FA"gE!HEHqZK
2HEFa1"IMd58L8aA NAEAI> 98LMacc133Z=DAAbDF1H%Cqs8D _FG1XAAvd{+ a))++4e4Aq&!)D.1!4 45
QQZ W
55 5 a 166!9bd#3"@"@A"aRSQSDTBw!Qy1a!Qy1qb"az1ffQi2a(#%% (1"'(,"*Q QrT!VAXaZ\*1a4/QQ 0BB"*, -
- a E F
M4,s N+- N0AN5c U R u p#nUS:X a U* [ X#S- U5 - $ US:X a" [ X#S-
U5 U[ X#U5 - -
U- $ [ e)N r$ )r* r! r )r; argindexr= r> r? s rM fdifflerchphi.fdiff sb ))a q=2hqa%++ +
] QAq)AhqQ.?,??B B$ $ c V U R 5 nUR U5 ( a U$ U $ N)r9 has)r; targetrA s rM _eval_rewrite_helper lerchphi._eval_rewrite_helper s' $$& 776??JKrT c , U R [ 5 $ rV )rY r+ r; r= r> r? kwargss rM _eval_rewrite_as_zeta lerchphi._eval_rewrite_as_zeta s ((..rT c , U R [ 5 $ rV )rY r7 r\ s rM _eval_rewrite_as_polylog!lerchphi._eval_rewrite_as_polylog s ((11rT r) Nr$ ) __name__
__module__ __qualname____firstlineno____doc__r9 rR rY r^ ra __static_attributes__r) rT rM r! r! s# gR?!B%/2rT r! c \ ^ \ rS rSrSr\S 5 rS
S jrS rS r S r
S U 4S jjr S r U =r
$ ) r7 aR
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 UR ( ay U[ R L a [ U5 $ U[ R L a [ U5 * $ U[ R L a [ R $ US:X a [ 5 nX#; a X2 $ UR ( a [ R $ UR [ R 5 nU( a [ U5 $ USL aR U[ R L a USU-
- $ U[ R L a USU-
S- - $ UR ( a USU-
- $ UR [ [ 5 ( a: U( d [ U5 [ R :* S:X a U " U[ U5 5 $ g g )Nr& Fr$ T) is_numberr
r3 r+ NegativeOne
dirichlet_etar- _dilogtableis_zeroequalsrW r r r r )clsr> r=
dilogtablezones rM eval polylog.eval% s, ;;AEEzAwamm#%a(((affvv
a (]
?%=( 9966M xx 7N
U]
AFF{!a%y amm#!a%!|#yy!a%y
55J ' 'Tc!fo$5Nq*Q-( ( 6O 'rT c Z U R u p#US:X a [ US-
U5 U- $ [ e)Nr& r$ )r* r7 r )r; rQ r> r= s rM rR
polylog.fdiffL s0 yy q=1q5!$Q& & rT c U[ X!S5 - $ Nr$ r! )r; r> r= r] s rM _eval_rewrite_as_lerchphi!polylog._eval_rewrite_as_lerchphiR s !"""rT c 2 U R u p#US:X a [ SU-
5 * $ UR ( aY US:: aS [ S5 nUSU-
- n[ U* 5 H nXER U5 - nM [
U5 R XC5 $ [ X#5 $ )Nr$ r u) r* r r, r r4 r0 r r1 r7 )r; r< r> r= r r@ _s rM r9 polylog._eval_expand_funcU s yy 6AJ; < SSK Jn U R u pgUR US5 nU[ R
L a- UR
US[ U5 R ( a SOSS9nUR ( a UR U5 u pU
R ( a [ X*- 5 n U" X- U5 n UR XX45 R! 5 n
U
[ R" L a U $ U
nU/n[% SU 5 H nX- nUR' UUU- - 5 M [) U6 U - $ [* [, U ]? XX45 $ ! [ [ 4 a U s $ f = f)Nr )Order-+)dirr& )sympy.series.orderr r* r1 r
NaNlimitr is_negativerq leadterm
ValueErrorNotImplementedError is_positiver
_eval_nseriesremoveOr- r4 r: r superr7 )r; xrE logxcdirr nur= z0r r newnortermr> rF __class__s rM r polylog._eval_nseriesf s) ,
VVAq\
;A"T(*>*>3CHB
::
A qu~!$NOOA$5==?; HF q$ AIDHHT!R%Z( ( Aw{"Wd1!CC# 34
s :D= =EEr) rc )r )rd re rf rg rh classmethodrv rR r} r9 r r ri
__classcell__r s @rM r7 r7 s? CJ $) $)L! #
D DrT r7 c t ^ \ rS rSrSr\S S j5 rS
S jrS
S jrS
S jr S r
S r S
S jr U 4S
jr
S rU =r$ )r+ i 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 U[ R L a U " U5 $ U[ R L d U[ R L a [ R $ U[ R L a [ R $ U[ R L a [ R $ U[ R L a [ R
$ UR nUc [ R nU( a& UR ( a [ SU-
U5 US-
- $ U[ R L aI U( aA UR ( a/ S[ - [ - U- * [ U5 - S[ U5 - - $ g g U( a: UR ( a) UR ( a U " U5 [ US-
U5 -
$ UR ( aA UR ( a/ UR SL d UR SL a [ R $ g g g )Nr$ r& F)r
r3 r ComplexInfinityInfinityr- r, is_nonpositiver is_evenr r
r r r
is_integer)rs r> r? sints rM rv zeta.eval s[ :q6M
!%%Z1:55L
!%%Z$$ $
!**_55L
!**_66M|| 9A A$$ QqS!$!, ,
!%%Z 2a! |il2a ! nEE "t
allq}}q6HQqS! ,, ,
\\a . . &!*:*:e*C55L +D /\rT c US:X ab UR ( aQ UR ( a@ UR ( a/ S[ - [ - U- * [ U5 - S[
U5 - - $ [ SU-
U5 US-
- $ )Nr$ r& )r is_nonnegativer r r
r r r; r> r? r] s rM _eval_rewrite_as_bernoulli zeta._eval_rewrite_as_bernoulli sa 6allq'7'7AIIrT!VaK<)A,.!IaL.A A1a AaC((rT c ^ US:w a U $ U R S n[ U5 SSSU-
- -
- $ )Nr$ r r& )r* ro r s rM _eval_rewrite_as_dirichlet_eta#zeta._eval_rewrite_as_dirichlet_eta s7 6K IIaL Q QQZ 00rT c [ SX5 $ r{ r| r s rM r} zeta._eval_rewrite_as_lerchphi s 1 rT c L [ U R S S-
R 5 $ )Nr r$ )r r* rq )r; s rM _eval_is_finitezeta._eval_is_finite s $))A,*3344rT c U R S n[ U R 5 S: a U R S O[ R nUR ( ak UR
( a [
U5 [ US-
U5 -
$ UR ( a. UR SL d UR SL a [ R $ U $ )Nr r$ F)
r* lenr
r3 r r r+ r r r )r; r< r> r? s rM r9 zeta._eval_expand_func" s IIaL NQ .DIIaLAEE <<}}Aw!A#q!111 Q\\U%:$$-uu rT c [ U R 5 S:X a U R u p#OU R S- u p#US:X a U* [ US- U5 - $ [ e)Nr& rc r$ )r r* r+ r )r; rQ r> r? s rM rR
zeta.fdiff- sS tyy>Q 99DAq99t#DA q=2d1q5!n$ $$ $rT c T > [ U R 5 S:X a U R u pEO U R [ R 4- u pE UR U5 u pgUR ( a UR ( d [
e[ [ U ]+ XUS9$ ! [
a U s $ f = f)Nr& )r r ) r r* r
r3 r r r r r r+ _eval_as_leading_term) r; r r r r> r? rB er s rM r zeta._eval_as_leading_term7 s tyy>Q 99DAq99x'DA ::a=DA
== % %T46q$6OO
# K s B B'&B'r) rV rc )rd re rf rg rh r rv r r r} r r9 rR r ri r r s @rM r+ r+ sH hT 4)
1 !5 %P PrT r+ c \ \ rS rSrSr\S S j5 rS
S jr\R 4S jr
S r Sr g) ro iH aE
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 U[ R L a U " U5 $ UcE US:X a [ S5 $ [ U5 nUR [ 5 ( d SSSU-
- -
U- $ g US:X a) SSKJn [ S5 U" U5 -
U" US- S- 5 - $ [ X5 n[ XS- S- 5 nUR [ 5 ( d) UR [ 5 ( d USSU-
- U- -
$ g g )Nr$ r& r digamma)r
r3 r r+ rW 'sympy.functions.special.gamma_functionsr )rs r> r? r= r z1z2s rM rv dirichlet_eta.evalw s :q6M 9Av1v
QA55;;A!H ))
!V Gq6GAJ&!A#q)99 9
!Z
!c1W
vvd||BFF4LLAaC2
% % %1|rT c
H SSK Jn US:X a8 [ [ S5 [ US5 4SSSU-
- -
[ U5 - S45 $ [ [ S5 U" U5 -
U" US- S- 5 - [ US5 4[ X5 SSU-
- [ XS- S- 5 - -
S45 $ )Nr r r$ r& T)r r r r r r+ r; r> r? r] r s rM r^ #dirichlet_eta._eval_rewrite_as_zeta s C 6 c!fbAh /1q1Q3x<472JD1QR R#a&71:-1a0@@"Q(Kaa!A#haA#q)9 994@B BrT c SSK Jn [ [ S5 U" U5 -
U" US- S- 5 - [ US5 4[ SU-
U5 SUS-
- - S45 $ )Nr r r& r$ T)r r r r r r r s rM _eval_rewrite_as_genocchi'dirichlet_eta._eval_rewrite_as_genocchi s^ C#a&71:-1a0@@"Q(K!A#q!Q!A#Y/68 8rT c [ S U R 5 5 ( a$ U R [ 5 R U5 $ g )Nc 3 8 # U H oR v M g 7frV )rm ).0is rM ,dirichlet_eta._eval_evalf.. s .Iq{{Is )allr* rewriter+ _eval_evalf)r; precs rM r dirichlet_eta._eval_evalf s6 .DII. . .<<%11$7 7 /rT r) rV rc )
rd re rf rg rh r rv r^ r
r3 r r ri r) rT rM ro ro H s5 ,\ & &$B ./UU 8
8rT ro c . \ rS rSrSr\S 5 rS rSrg)
riemann_xii a
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 SSK Jn [ U5 nU[ R [ R
4; a [ R $ [ U[ 5 ( d& XS-
- U" US- 5 - U- S[ US- - - - $ g Nr )gammar$ r& ) r r r+ r
r- r3 Halfr8 r )rs r> r r= s rM rv riemann_xi.eval sl A G 66M!T""!e9U1Q3Z')1R!A#Y;7 7 #rT c l SSK Jn XS-
- U" US- 5 - [ U5 - S[ US- - - - $ r )r r r+ r )r; r> r] r s rM r^ riemann_xi._eval_rewrite_as_zeta s8 Aa%yqs#DG+QrAaCy[99rT r) N) rd re rf rg rh r rv r^ ri r) rT rM r r s . 8 8:rT r c , \ rS rSrSr\SS j5 rSrg) stieltjesi a$
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 Ub` [ U5 nU[ R L a [ R $ UR ( a! UR ( a [ R
$ UR ( a U[ R L a [ R $ US: a [ R
$ UR ( d [ R
$ U[ R L a US; a [ R $ UR ( a [ R
$ UR ( a US; a [ R $ UR S:X a [ R
$ g )Nr r{ F) r r
r r, r r is_Numberr-
EulerGammais_extended_negativerq r )rs rE r? s rM rv stieltjes.eval s =
AAEEzuu || 0 0((( ;;AEEzuu Q(((\\(((affi ||# ! !$$ $ 99i << <<5 $$ $ !rT r) rV )rd re rf rg rh r rv ri r) rT rM r r s "H % %rT r c [ R [ S- S- [ S5 S- S- -
[ S5 [ S- S- [
[ - [ S5 - -
[
S5 S-
* S- [ S- * S- [ [
S5 S-
S- 5 S- S- - [
S5 S- * S- [ S- * S- [ [
S5 S- S- 5 S- -
S[
S5 -
S- [ S- S- [ [
S5 S-
S- 5 S- -
[
S5 S-
S- [ S- S- [ [
S5 S-
S- 5 S- -
[
[
[ R - [ S- S - -
[
* [
* [ R - [ S- S - -
S[
-
[ S- S
- [
[ R - -
[ [
- S- [ S5 - -
S[
- [ S- S
- [
[ R - - [ [
- S- [ S5 - - S[
-
S- [ S5 S- * S - [ [
- [ S5 - S - - S[ S- - S - - [
[ R - -
0 $ )
Nr& r( r$
rP 0 ` )r
r r r r r
r Catalanr) rT rM rp rp 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 rT N)5rh sympy.core.addr sympy.core.cacher sympy.core.functionr r r sympy.core.logicr 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! r7 r+ ro r r rp r) rT rM r s * $ O O & - - $ " # & Z Z P P F F > 9 : &2 2LeDo eDXzP? zPzP8O P8f$: $:N?% ?%D
rT