$$\newcommand{acoloneqq}[0]{\ &\hspace-2pt\coloneqq}
\newcommand{ar}[1]{\operatorname{ar{#1}}}
\newcommand{arc}[1]{\operatorname{arc{#1}}}
\newcommand{asupplement}[1]{&\hspace{#1}\textsf}
\newcommand{beginend}[2]{{\begin{#1}#2\end{#1}}}
\newcommand{bm}[0]{\boldsymbol}
\newcommand{bscolor}[1]{\color{var(--bs-#1)}}
\newcommand{bscolorbox}[1]{\colorbox{var(--bs-#1)}}
\newcommand{bsrowcolor}[1]{\rowcolor{var(--bs-#1)}}
\newcommand{C}[0]{\mathbb C}
\newcommand{Defarrow}[0]{\xLeftrightarrow{\textrm{def}}}
\newcommand{DefOp}[1]{\DeclareMathOperator{#1}{#1}}
\newcommand{fqty}[0]{\!\qty}
\newcommand{hen}[1]{{(\textrm{{#1}辺})}}
\newcommand{hygeo}[6]{{{}_{#1}{#2}_{#3}{\qty[\beginend{matrix*}{#4\\ #5}\ ;{#6}]}}}
\newcommand{In}[0]{\in\mathbb}
\newcommand{kfrac}[0]{\mathop{\Large\raise-.8pt{\textrm K}}}
\newcommand{Kfrac}[0]{\mathop{\huge\raise-2.2pt{\textrm K}}}
\newcommand{kome}[0]{\textreferencemark}
\newcommand{leftshiftarrow}[0]{{\substack{\curvearrowright\\ \leftharpoondown}}}
\newcommand{lr}[3]{{\left#1{#2}\right#3}}
\newcommand{lvvr}[3]{\lr{#1}{\negmedspace\lr|{#2}|\negmedspace}{#3}}
\newcommand{N}[0]{\mathbb N}
\newcommand{ot}[0]{\leftarrow}
\newcommand{P}[0]{\mathbb P}
\newcommand{Q}[0]{\mathbb Q}
\newcommand{R}[0]{\mathbb R}
\newcommand{RANGE}[0]{}\newcommand{rangeex}[6][,]{{#2{#3}_{#5}#4#1\cdots#1#2{#3}_{#6}#4}}\newcommand{range}[2][,]{\rangeex[#1]{}{#2}{}}{}
\newcommand{REQUIRE}[0]{}\require{physics}{}
\newcommand{rightshiftarrow}[0]{{\substack{\curvearrowleft\\ \rightharpoondown}}}
\newcommand{rmIm}[0]{\operatorname{Im}}
\newcommand{rmRe}[0]{\operatorname{Re}}
\newcommand{rprod}[0]{\mathop{\prod\!\llap\coprod}}
\newcommand{sahen}[0]{\hen左}
\newcommand{STIRLING}[0]{}\newcommand{stirling}[3][]{{\qty[\beginend{matrix}{{#2}\\ {#3}}{#1}]}}\newcommand{Stirling}[3][]{{\qty{\beginend{matrix}{{#2}\\ {#3}}{#1}}}}{}
\newcommand{uhen}[0]{\hen右}
\newcommand{vbin}[1]{\mathbin{{#1}\!\llap|\ }}
\newcommand{Z}[0]{\mathbb Z}
$$
$
\beginend{array}{l
N\ge1 \\
\beginend{aligned}{
\frac\pi4 &=
\sum_{n=0}^{N-1} \frac{(-1)^n}{2n+1} +
\frac{(-1)^N}2\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2N} \\&=
\sum_{n=0}^{N-1} \frac{(-1)^n}{2n+1} +
\cfrac{\lr.{(-1)^N}/2}
{2N+\cfrac{1^2}
{2N+\cfrac{2^2}
{2N+\cfrac{3^2}
{2N+\ddots}}}}
} \\
N\ge0 \\
\beginend{aligned}{
\ln2 &=
\sum_{n=1}^N \frac{(-1)^{n-1}}n +
(-1)^N\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2N+1} \\&=
\sum_{n=1}^N \frac{(-1)^{n-1}}n +
\cfrac{(-1)^N}
{2N+1+\cfrac{1^2}
{2N+1+\cfrac{2^2}
{2N+1+\cfrac{3^2}
{2N+1+\ddots}}}}
}}$
定義
| |
|---|
| $\delta_{a,b}$ |
クロネッカーのデルタ
|
| $\displaystyle \sum_{n=0}^{-1} f(n) \coloneqq 0$ |
空和
|
| $\displaystyle \Kfrac_{n=0}^\infty \frac{a_n}{b_n}$ |
連分数
|
| $B_n(z),E_n(z)$ |
ベルヌーイ多項式、オイラー多項式
|
一般的な式
以前書いた『
ゼータ関数を加速させよ
』でこのような式を導出しました。
$\displaystyle \sum_{n=0}^\infty \frac{(-1)^n}{z+n} = \sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} + (-1)^N\sum_{n=0}^\infty E_{2n}\cdot(2z+2N-1)^{-2n-1}$
右辺の2つ目の級数は発散級数で、$N$を大きくし分母を大きくした上で
適切な項数の部分和をとることにより一つ目の部分和の誤差が打ち消されるというものです。
この発散級数は
連分数展開
できます。
数値計算によればこの連分数は恐らく収束します。
$N$を大きくしていくことで連分数の収束は急速に速まるため、
数値計算に有用です。
$\displaystyle 4\qty(\frac11-\frac13+\cdots-\frac1{15} +
\cfrac{\frac12}
{16+\cfrac{1^2}
{16+\cfrac{2^2}
{\phantom{16+{}}\cfrac{\small\ddots}
{16+\cfrac{6^2}
{16}}}}}) =
3.141592653\bscolor{secondary}60\cdots$
$N\ge0$
$\beginend{align}{
\sum_{n=0}^\infty \frac{(-1)^n}{z+n} &=
\sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} +
\frac{(-1)^N}2\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2z+2N-1} \\&=
\sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} +
\cfrac{(-1)^N}
{2z+2N-1+\cfrac{1^2}
{2z+2N-1+\cfrac{2^2}
{2z+2N-1+\cfrac{3^2}
{2z+2N-1+\scriptsize\ddots}}}}
}$
オイラーの定数
先程の式の左辺$\displaystyle \sum_{n=0}^\infty \frac{(-1)^n}{z+n}$は$
\dfrac12\qty(\psi\fqty(\dfrac{z+1}2)-\psi\fqty(\dfrac z2))$と等しく、
ディガンマ関数が出現しますが、打ち消し合うためここからオイラーの定数$\gamma$を導くことはできません。
『ゼータ関数を加速させよ』で導いた式により、ディガンマ関数は次のように表されます。
$\displaystyle \psi(z) = \ln(z-a+N) - \sum_{n=0}^{N-1} \frac1{z+n} +
\sum_{n=1}^\infty \frac{(-1)^{n-1}B_n(a)}n(z-a+N)^{-n}$
よって、$\displaystyle \sum_{n=1}^\infty \frac{B_n(a)}nz^n$を連分数展開できればオイラーの定数を効率的に求めることができそうです。
また、別のアプローチとして
指数積分
と
ガウスの連分数
の${}_2F_0$への一般化を用いる方法があります。
$
\beginend{align}{
x&>0\ ? \\
\gamma &=
-\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x -
\cfrac{e^{-x}}
{x+\cfrac1
{1+\cfrac1
{x+\cfrac2
{1+\cfrac2
{x+\cfrac3
{1+\cfrac3
{x+\scriptsize\ddots}}}}}}} \\&=
-\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x -
e^{-x}\Kfrac_{n=0}^\infty \frac{\delta_{n,0}-n^2}{2n+1+x} \\&=
-\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x -
\cfrac{e^{-x}}
{1+x-\cfrac{1^2}
{3+x-\cfrac{2^2}
{5+x-\cfrac{3^2}
{7+x-\ddots}}}}
}$
1つ目の連分数からの
式変形
により2つ目の連分数が得られます。
$x$を大きくすれば連分数は速くなりますがその分級数が遅くなります。
