using Base64 using PyPlot using QuadGK using SpecialFunctions using SymPy const e = ℯ linspace(start, stop, length) = range(start, stop=stop, length=length) Euler変換¶函数 $f(x)$ に作用する差分作用素 $T$ を $Tf(x)=f(x+1)$ と定めると, 数列 $f(x+n)$ に関する交代級数は形式的に次のように表わされる: $$ \sum_{n=0}^\infty (-1)^n f(x+n)=\sum_{n=0}^\infty (-T)^nf(x)=\frac{1}{1+T}\,f(x). $$そして, $$ \frac{1}{1+T} = \frac{1}{2-(1-T)} = \sum_{n=0}^\infty