オイラーの和公式

数学において、オイラーの和公式（オイラーのわこうしき、オイラー・マクローリンの公式、英: Euler–Maclaurin formula）は1735年頃オイラーとマクローリンにより独立に発見された級数の和を与える公式である。この公式は収束の遅い無限級数の和を求めるときに便利であるが、 f ( x ) {\displaystyle f(x)} が多項式であるような場合を除き、 m → ∞ {\displaystyle m\to \infty } とすればベルヌーイ数が急速に大きくなって発散する。従って、漸近展開のように発散する前の適当なところで打ち切らなければならない。また、この公式は台形公式による数値積分の誤差を示すものと考えることもできる。 ∑ j = 0 n − 1 f ( j ) = ∫ 0 n f ( x ) d x + ∑ k = 1 m B k k ! ( f ( k − 1 ) ( n ) − f ( k − 1 ) ( 0 ) ) + R m {\displaystyle \sum _{j=0}^{n-1}f(j)=\int _{0}^{n}f(x)dx+\sum _{k=1}^{m}{\frac {B_{k}}{k!}}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)+R_{m}} ∑ j = 1 n − 1 f ( j ) + 1 2 ( f ( 0 ) + f ( n ) ) = ∫ 0 n f ( x ) d x + ∑ k = 1 m B 2 k ( 2 k ) ! ( f ( 2 k − 1 ) ( n ) − f ( 2 k − 1 ) ( 0 ) ) + R 2 m + 1 {\displaystyle \sum _{j=1}^{n-1}f(j)+{\frac {1}{2}}\left(f(0)+f(n)\right)=\int _{0}^{n}f(x)dx+\sum _{k=1}^{m}{\frac {B_{2k}}{(2k)!}}\left(f^{(2k-1)}(n)-f^{(2k-1)}(0)\right)+R_{2m+1}} R m = ( − 1 ) m + 1 ∫ 0 n B m ( x − ⌊ x ⌋ ) m ! f ( m ) ( x ) d x {\displaystyle R_{m}=(-1)^{m+1}\int _{0}^{n}{\frac {B_{m}(x-\lfloor {x}\rfloor )}{m!}}f^{(m)}(x)dx} 但し、 B n {\displaystyle B_{n}} はベルヌーイ数、 B n ( x ) {\displaystyle B_{n}(x)} はベルヌーイ多項式である。 B 1 = − 1 2 , B 2 = 1 6 , B 3 = 0 , B 4 = − 1 30 , B 5 = 0 , B 6 = 1 42 , B 7 = 0 , B 8 = − 1 30 , B 9 = 0 , B 10 = 5 66 , ⋯ {\displaystyle B_{1}=-{\frac {1}{2}},B_{2}={\frac {1}{6}},B_{3}=0,B_{4}=-{\frac {1}{30}},B_{5}=0,B_{6}={\frac {1}{42}},B_{7}=0,B_{8}=-{\frac {1}{30}},B_{9}=0,B_{10}={\frac {5}{66}},\cdots } B 0 ( x ) = 1 , B 1 ( x ) = x − 1 2 , B 2 ( x ) = x 2 − x + 1 6 , B 3 ( x ) = x 3 − 3 2 x 2 + 1 2 x , B 4 ( x ) = x 4 − 2 x 3 + x 2 − 1 30 , … {\displaystyle B_{0}(x)=1,B_{1}(x)=x-{\frac {1}{2}},B_{2}(x)=x^{2}-x+{\frac {1}{6}},B_{3}(x)=x^{3}-{\frac {3}{2}}x^{2}+{\frac {1}{2}}x,B_{4}(x)=x^{4}-2x^{3}+x^{2}-{\frac {1}{30}},\dots } なお、 f ( k ) {\displaystyle f^{(k)}} は導関数、 ⌊ x ⌋ {\displaystyle \lfloor {x}\rfloor } は床関数を表す。 ダルブーの公式（英語版）はこれの一般化である。 ベルヌーイ多項式の性質（若しくは定義）により ∫ 0 1 B k − 1 ( x ) ( k − 1 ) ! f ( k − 1 ) ( x ) d x = ∫ 0 1 ( B k ( x ) k ! ) ′ f ( k − 1 ) ( x ) = [ B k ( x ) k ! f ( k − 1 ) ( x ) ] 0 1 − ∫ 0 1 B k ( x ) k ! f ( k ) ( x ) d x {\displaystyle \int _{0}^{1}{\frac {B_{k-1}(x)}{(k-1)!}}f^{(k-1)}(x)dx=\int _{0}^{1}\left({\frac {B_{k}(x)}{k!}}\right)'f^{(k-1)}(x)=\left[{\frac {B_{k}(x)}{k!}}f^{(k-1)}(x)\right]_{0}^{1}-\int _{0}^{1}{\frac {B_{k}(x)}{k!}}f^{(k)}(x)dx} である。有限回の部分積分を繰り返して ∫ 0 1 f ( x ) d x = ∫ 0 1 B 0 ( x ) f ( x ) d x = ∑ k = 1 m [ ( − 1 ) k − 1 B k ( x ) k ! f ( k − 1 ) ( x ) ] 0 1 + ( − 1 ) m ∫ 0 1 B m ( x ) m ! f ( m ) ( x ) d x {\displaystyle \int _{0}^{1}f(x)dx=\int _{0}^{1}B_{0}(x)f(x)dx=\sum _{k=1}^{m}\left[(-1)^{k-1}{\frac {B_{k}(x)}{k!}}f^{(k-1)}(x)\right]_{0}^{1}+(-1)^{m}\int _{0}^{1}{\frac {B_{m}(x)}{m!}}f^{(m)}(x)dx} となるが、これは f ( x ) {\displaystyle f(x)} を f ( j + x ) {\displaystyle f(j+x)} に置き換えても成り立つから ∫ 0 n f ( x ) d x = ∑ j = 0 n − 1 ∫ 0 1 f ( j + x ) d x = ∑ j = 1 n − 1 ∑ k = 1 m [ ( − 1 ) k − 1 B k ( x ) k ! f ( k − 1 ) ( x ) ] 0 1 + ( − 1 ) m ∫ 0 n B m ( x − ⌊ x ⌋ ) m ! f ( m ) ( x ) d x {\displaystyle {\begin{aligned}\int _{0}^{n}f(x)dx&=\sum _{j=0}^{n-1}\int _{0}^{1}f(j+x)dx\\&=\sum _{j=1}^{n-1}\sum _{k=1}^{m}\left[(-1)^{k-1}{\frac {B_{k}(x)}{k!}}f^{(k-1)}(x)\right]_{0}^{1}+(-1)^{m}\int _{0}^{n}{\frac {B_{m}(x-\lfloor {x}\rfloor )}{m!}}f^{(m)}(x)dx\end{aligned}}} である。 B 1 ( 0 ) = − 1 2 , {\displaystyle B_{1}(0)=-\textstyle {\frac {1}{2}},} B 1 ( 1 ) = 1 2 , B 2 k ( 0 ) = B 2 k ( 1 ) = B 2 k , {\displaystyle B_{1}(1)=\textstyle {\frac {1}{2}},B_{2k}(0)=B_{2k}(1)=B_{2k},} B 2 k + 1 ( 0 ) = B 2 k + 1 ( 1 ) = B 2 k + 1 = 0 {\displaystyle B_{2k+1}(0)=B_{2k+1}(1)=B_{2k+1}=0} を代入すれば ∫ 0 n f ( x ) d x = ∑ j = 0 n − 1 f ( j ) − 1 2 f ( 0 ) + 1 2 f ( n ) − ∑ k = 2 m ( − 1 ) k B k k ! ( f ( k − 1 ) ( n ) − f ( k − 1 ) ( 0 ) ) − R m {\displaystyle \int _{0}^{n}f(x)dx=\sum _{j=0}^{n-1}f(j)-{\frac {1}{2}}f(0)+{\frac {1}{2}}f(n)-\sum _{k=2}^{m}(-1)^{k}{\frac {B_{k}}{k!}}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)-R_{m}} R m = ( − 1 ) m + 1 ∫ 0 n B m ( x − ⌊ x ⌋ ) m ! f ( m ) ( x ) d x {\displaystyle R_{m}=(-1)^{m+1}\int _{0}^{n}{\frac {B_{m}(x-\lfloor {x}\rfloor )}{m!}}f^{(m)}(x)dx} を得る。移項して形式を整えると ∑ j = 0 n − 1 f ( j ) = ∫ x = 0 n f ( x ) d x + ∑ k = 1 m B k k ! ( f ( k − 1 ) ( n ) − f ( k − 1 ) ( 0 ) ) + R m {\displaystyle \sum _{j=0}^{n-1}f(j)=\int _{x=0}^{n}f(x)dx+\sum _{k=1}^{m}{\frac {B_{k}}{k!}}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)+R_{m}} となる。或いは ∑ j = 1 n − 1 f ( j ) + 1 2 ( f ( 0 ) + f ( n ) ) = ∫ 0 n f ( x ) d x + ∑ k = 2 2 m + 1 B k k ! ( f ( k − 1 ) ( n ) − f ( k − 1 ) ( 0 ) ) + R 2 m + 1 = ∫ 0 n f ( x ) d x + ∑ k = 1 m B 2 k ( 2 k ) ! ( f ( 2 k − 1 ) ( n ) − f ( 2 k − 1 ) ( 0 ) ) + R 2 m + 1 {\displaystyle {\begin{aligned}\sum _{j=1}^{n-1}f(j)+{\frac {1}{2}}\left(f(0)+f(n)\right)&=\int _{0}^{n}f(x)dx+\sum _{k=2}^{2m+1}{\frac {B_{k}}{k!}}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)+R_{2m+1}\\&=\int _{0}^{n}f(x)dx+\sum _{k=1}^{m}{\frac {B_{2k}}{(2k)!}}\left(f^{(2k-1)}(n)-f^{(2k-1)}(0)\right)+R_{2m+1}\\\end{aligned}}} となる。 アーベル・プラナの公式 レオンハルト・オイラー コリン・マクローリン ^ Springer Online Reference Works: Euler–MacLaurin formula