アイゼンシュタイン級数の収束が遅すぎる！
んで高速化しようというこのプロジェクトの結論。
まず定義。
$$G_{2k}(\tau) := \sum_{(m,n)\in \mathbb{Z}^2 \backslash (0,0)} \frac{1}{(m\tau + n)^{2k}}$$
で、高速版。(ついでに上記は上半平面上でしか収束しないのにこっち実軸以外全平面で収束するバグ)
$$G_{2k}(\tau)=2\zeta(2k)-\frac{2\pi^{2k}}{(2k-1)!}\sum^{\infty}_{n=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z)\right|_{z=\pi n \tau} \right)$$
よーし頑張ってこー
$$ \begin{align} &G_{2k}(\tau) := \sum_{(m,n)\in \mathbb{Z}^2 \backslash (0,0)} \frac{1}{(m\tau + n)^{2k}} \\ &\qquad\quad \ = 2\sum^{\infty}_{m=1} \sum^{\infty}_{n=1} \left( \frac{1}{(m\tau + n)^{2k}} + \frac{1}{(m\tau - n)^{2k}} \right)+2\sum^{\infty}_{m=1}\frac{1}{(m\tau)^{2k}}+2\sum^{\infty}_{n=1}\frac{1}{n^{2k}} \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( \psi^{(2k-1)}(1+m\tau) + \psi^{(2k-1)}(1-m\tau) \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( \psi^{(2k-1)}(1+m\tau) - \psi^{(2k-1)}(m\tau) - \pi^{2k} \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\frac{1}{(2k-1)!} \sum^{\infty}_{m=1} \left( -\frac{(2k-1)!}{(\tau m)^{2k}} - \pi^{2k} \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = - 2\frac{1}{\tau^{2k}}\zeta(2k)- \frac{2\pi^{2k}}{(2k-1)!} \sum^{\infty}_{m=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right)+2\frac{1}{\tau^{2k}}\zeta(2k)+2\zeta(2k) \\ &\qquad\quad \ = 2\zeta(2k) - \frac{2\pi^{2k}}{(2k-1)!} \sum^{\infty}_{m=1} \left( \left. \frac{d^{2k-1}}{dz^{2k-1}}\cot(z) \right|_{z=\pi m \tau} \right) \end{align} $$
ほぼゼータの定義とディガンマ関数の級数展開、相反公式をただぶっこむだけだね。
あってるよな...?(自分を信じたら負けだと思ってる)
DomeainColoringでもするのだ。とりあえず$k=2,3$でやろうではないか。
$$G_4(\tau) = \frac{\pi^4}{45}-\frac{2\pi^4}{3}\sum^{\infty}_{n=1}\left( \frac{2}{\sin^2(\pi n\tau)} - \frac{3}{\sin^4(\pi n\tau)} \right)$$
$$G_6(\tau) = \frac{2\pi^6}{945}+\frac{2\pi^6}{15}\sum^{\infty}_{n=1}\left( \frac{2}{\sin^2(\pi n\tau)} - \frac{15}{\sin^4(\pi n\tau)} + \frac{15}{\sin^6(\pi n\tau)} \right)$$
計算ミスが怖いけど、まぁグラフ書かせるかぁー
G4の画像
お、いい感じ。
G4の拡大図
思ったより緻密だねぇ?
G6の画像
うむうむ、よろしい。
ちなみに、一緒にこうやって遊んでたやつが、「微分のところ一般化できないかな」なんていってこんな式を得ていた。
$$ \begin{align} &G_{2k}(\tau)=(-1)^{k+1}\frac{(2\pi)^{2k}}{(2k)!}\left( B_{2k} -2k\sum^{\infty}_{n=1}\frac{A_{2k-1}(q^n)}{(1-q^n)^{2k}} \right) \\ &ただしq=e^{2\pi i \tau} \end{align} $$
なお、$A_{n}(x)$はオイラー多項式、$B_{n}$はおなじみベルヌーイ数ですぞ。
要望があれば導出かくかも。でも収束遅いしまた下半平面で収束してないし、作った本人もじゃあ何ってところまで行ってなかったから多分そんな要望は来ない。たぶん。
ではまた。
この記事を高評価した人
この記事に送られたバッジ
投稿者
