線形多段階法

Class of iterative numerical methods for solving differential equations

線形マルチステップ法は、常微分方程式の数値解法に使用されます。概念的には、数値法は初期点から開始し、次に短い時間ステップを踏んで次の解点を探します。このプロセスは、解をマップするために後続のステップに継続されます。シングルステップ法 (オイラー法など) では、前の 1 つの点とその導関数のみを参照して現在の値を決定します。ルンゲ・クッタなどの方法では、より高次の方法を得るためにいくつかの中間ステップ (たとえば、半ステップ) を実行しますが、2 番目のステップを実行する前に以前の情報をすべて破棄します。マルチステップ法では、前のステップからの情報を破棄するのではなく、保持して使用することで効率を上げようとします。したがって、マルチステップ法では、前の複数の点と導関数値を参照します。線形マルチステップ法の場合、前の点と導関数値の線形結合が使用されます。

定義

常微分方程式の数値解析法は、次のような初期値問題の近似解を求める。 $${\displaystyle y'=f(t,y),\quad y(t_{0})=y_{0}.}$$

結果は、離散時間 におけるの値の近似値です。ここで、 は時間ステップ ( と呼ばれることもあります) であり、は整数です。 ${\displaystyle y(t)}$ ${\displaystyle t_{i}}$ $${\displaystyle y_{i}\approx y(t_{i})\quad {\text{where}}\quad t_{i}=t_{0}+ih,}$$ ${\displaystyle h}$ ${\displaystyle \Delta t}$ ${\displaystyle i}$

マルチステップ法は、前のステップの情報を用いて次の値を計算します。特に、線形マルチステップ法は、 と の線形結合を用いて、目的の現在のステップにおけるの値を計算します。したがって、線形マルチステップ法は、という形式の方法です。係数と によって、この方法は決定されます。この方法の設計者は、真の解への良好な近似値を得る必要性と、適用しやすい方法を求める要望のバランスを取りながら、係数を選択します。多くの場合、方法を簡素化するために、多くの係数がゼロになります。 ${\displaystyle s}$ ${\displaystyle y_{i}}$ ${\displaystyle f(t_{i},y_{i})}$ ${\displaystyle y}$ $${\displaystyle {\begin{aligned}&y_{n+s}+a_{s-1}\cdot y_{n+s-1}+a_{s-2}\cdot y_{n+s-2}+\cdots +a_{0}\cdot y_{n}\\&\qquad {}=h\cdot \left(b_{s}\cdot f(t_{n+s},y_{n+s})+b_{s-1}\cdot f(t_{n+s-1},y_{n+s-1})+\cdots +b_{0}\cdot f(t_{n},y_{n})\right)\\&\Leftrightarrow \sum _{j=0}^{s}a_{j}y_{n+j}=h\sum _{j=0}^{s}b_{j}f(t_{n+j},y_{n+j}),\end{aligned}}}$$ ${\displaystyle a_{s}=1}$ ${\displaystyle a_{0},\dotsc ,a_{s-1}}$ ${\displaystyle b_{0},\dotsc ,b_{s}}$

明示的法と暗黙的法を区別することができます。 の場合、式は を直接計算できるため、この方法は「明示的」と呼ばれます。の場合、 の値はの値に依存し、方程式は について解く必要があるため、この方法は「暗黙的」と呼ばれます。暗黙的式を解くには、ニュートン法などの反復法がよく用いられます。 ${\displaystyle b_{s}=0}$ ${\displaystyle y_{n+s}}$ ${\displaystyle b_{s}\neq 0}$ ${\displaystyle y_{n+s}}$ ${\displaystyle f(t_{n+s},y_{n+s})}$ ${\displaystyle y_{n+s}}$

明示的な多段階法を用いて の値を「予測」することがあります。その値は、暗黙的な式で使用され、値を「修正」します。その結果が予測子修正子法です。 ${\displaystyle y_{n+s}}$

例

例として、問題を考えてみましょう 。正確な解は です。 $${\displaystyle y'=f(t,y)=y,\quad y(0)=1.}$$ ${\displaystyle y(t)=e^{t}}$

ワンステップオイラー

単純な数値法はオイラー法です。オイラー法は、1 ステップの退化ケースに対する明示的なマルチステップ法として考えることができます。 $${\displaystyle y_{n+1}=y_{n}+hf(t_{n},y_{n}).}$$

この方法を問題にステップ サイズで適用すると、次の結果が得られます。 ${\displaystyle h={\tfrac {1}{2}}}$ ${\displaystyle y'=y}$ $${\displaystyle {\begin{aligned}y_{1}&=y_{0}+hf(t_{0},y_{0})=1+{\tfrac {1}{2}}\cdot 1=1.5,\\y_{2}&=y_{1}+hf(t_{1},y_{1})=1.5+{\tfrac {1}{2}}\cdot 1.5=2.25,\\y_{3}&=y_{2}+hf(t_{2},y_{2})=2.25+{\tfrac {1}{2}}\cdot 2.25=3.375,\\y_{4}&=y_{3}+hf(t_{3},y_{3})=3.375+{\tfrac {1}{2}}\cdot 3.375=5.0625.\end{aligned}}}$$

ツーステップ・アダムス・バッシュフォース

オイラー法は1 ステップの手法です。単純なマルチステップ手法は 2 ステップの Adams–Bashforth 法です。この手法では、次の値 を計算するためにと の2 つの値が必要です。しかし、初期値問題では という 1 つの値しか提供されません。この問題を解決する 1 つの方法は、オイラー法で計算された を 2 番目の値として使用することです。この選択により、Adams–Bashforth 法では次の式が得られます (4 桁に丸められます)。における正確な解はであるため、2 ステップの Adams–Bashforth 法の方がオイラー法よりも正確です。ステップ サイズが十分に小さい場合は常にこのようになります。 $${\displaystyle y_{n+2}=y_{n+1}+{\tfrac {3}{2}}hf(t_{n+1},y_{n+1})-{\tfrac {1}{2}}hf(t_{n},y_{n}).}$$ ${\displaystyle y_{n+1}}$ ${\displaystyle y_{n}}$ ${\displaystyle y_{n+2}}$ ${\displaystyle y_{0}=1}$ ${\displaystyle y_{1}}$ $${\displaystyle {\begin{aligned}y_{2}&=y_{1}+{\tfrac {3}{2}}hf(t_{1},y_{1})-{\tfrac {1}{2}}hf(t_{0},y_{0})=1.5+{\tfrac {3}{2}}\cdot {\tfrac {1}{2}}\cdot 1.5-{\tfrac {1}{2}}\cdot {\tfrac {1}{2}}\cdot 1=2.375,\\y_{3}&=y_{2}+{\tfrac {3}{2}}hf(t_{2},y_{2})-{\tfrac {1}{2}}hf(t_{1},y_{1})=2.375+{\tfrac {3}{2}}\cdot {\tfrac {1}{2}}\cdot 2.375-{\tfrac {1}{2}}\cdot {\tfrac {1}{2}}\cdot 1.5=3.7812,\\y_{4}&=y_{3}+{\tfrac {3}{2}}hf(t_{3},y_{3})-{\tfrac {1}{2}}hf(t_{2},y_{2})=3.7812+{\tfrac {3}{2}}\cdot {\tfrac {1}{2}}\cdot 3.7812-{\tfrac {1}{2}}\cdot {\tfrac {1}{2}}\cdot 2.375=6.0234.\end{aligned}}}$$ ${\displaystyle t=t_{4}=2}$ ${\displaystyle e^{2}=7.3891\ldots }$

多段階法のファミリー

一般的に使用される線形多段階法には、アダムス・バッシュフォース法、アダムス・モールトン法、および後方微分公式(BDF) の 3 つのファミリーがあります。

アダムズ・バッシュフォース法

アダムズ・バッシュフォース法は陽解法です。係数は とであり、 はsの位数を持つように選択されます（これにより s の位数が一意に決定されます）。 ${\displaystyle a_{s-1}=-1}$ ${\displaystyle a_{s-2}=\cdots =a_{0}=0}$ ${\displaystyle b_{j}}$

s = 1, 2, 3, 4, 5 の場合の Adams–Bashforth 法は次のとおりです (Hairer, Nørsett & Wanner 1993, §III.1; Butcher 2003, p. 103): $${\displaystyle {\begin{aligned}y_{n+1}&=y_{n}+hf(t_{n},y_{n}),\qquad {\text{(This is the Euler method)}}\\y_{n+2}&=y_{n+1}+h\left({\frac {3}{2}}f(t_{n+1},y_{n+1})-{\frac {1}{2}}f(t_{n},y_{n})\right),\\y_{n+3}&=y_{n+2}+h\left({\frac {23}{12}}f(t_{n+2},y_{n+2})-{\frac {16}{12}}f(t_{n+1},y_{n+1})+{\frac {5}{12}}f(t_{n},y_{n})\right),\\y_{n+4}&=y_{n+3}+h\left({\frac {55}{24}}f(t_{n+3},y_{n+3})-{\frac {59}{24}}f(t_{n+2},y_{n+2})+{\frac {37}{24}}f(t_{n+1},y_{n+1})-{\frac {9}{24}}f(t_{n},y_{n})\right),\\y_{n+5}&=y_{n+4}+h\left({\frac {1901}{720}}f(t_{n+4},y_{n+4})-{\frac {2774}{720}}f(t_{n+3},y_{n+3})+{\frac {2616}{720}}f(t_{n+2},y_{n+2})-{\frac {1274}{720}}f(t_{n+1},y_{n+1})+{\frac {251}{720}}f(t_{n},y_{n})\right).\end{aligned}}}$$

係数は次のように決定できます。多項式補間を使用して、次数 の多項式pを求めます。多項式補間のラグランジュの公式は次式となります。多項式pは、解くべき微分方程式の右辺の局所的に良い近似値となるため、代わりに方程式を考えます。この方程式は正確に解くことができ、解はpの積分になります。これは、次式を取ることを示唆しています。 アダムズ・バッシュフォース法は、 pの公式を代入したときに発生します。係数は次式で与えられます 。補間式pで 置き換えると、次数h ^sの誤差が生じるため、sステップのアダムズ・バッシュフォース法は、次数sになります(Iserles 1996、§2.1)。 ${\displaystyle b_{j}}$ ${\displaystyle s-1}$ $${\displaystyle p(t_{n+i})=f(t_{n+i},y_{n+i}),\qquad {\text{for }}i=0,\ldots ,s-1.}$$ $${\displaystyle p(t)=\sum _{j=0}^{s-1}{\frac {(-1)^{s-j-1}f(t_{n+j},y_{n+j})}{j!(s-j-1)!h^{s-1}}}\prod _{i=0 \atop i\neq j}^{s-1}(t-t_{n+i}).}$$ ${\displaystyle y'=f(t,y)}$ ${\displaystyle y'=p(t)}$ $${\displaystyle y_{n+s}=y_{n+s-1}+\int _{t_{n+s-1}}^{t_{n+s}}p(t)\,\mathrm {d} t.}$$ ${\displaystyle b_{j}}$ $${\displaystyle b_{s-j-1}={\frac {(-1)^{j}}{j!(s-j-1)!}}\int _{0}^{1}\prod _{i=0 \atop i\neq j}^{s-1}(u+i)\,\mathrm {d} u,\qquad {\text{for }}j=0,\ldots ,s-1.}$$ ${\displaystyle f(t,y)}$

アダムズ・バッシュフォース法は、フランシス・バッシュフォースによる毛細管現象をモデル化する微分方程式を解くために、ジョン・カウチ・アダムズによって考案されました。バッシュフォース（1883）は自身の理論とアダムズの数値解析法（ゴールドスタイン、1977）を発表しました。

アダムス・モールトン法

Adams–Moulton法は、Adams–Bashforth法と同様に、およびを持つという点で類似しています。この場合も、係数bは可能な限り高い次数が得られるように選択されます。しかし、Adams–Moulton法は暗黙的な手法です。 という制約を取り除くことで、sステップAdams–Moulton法は次数 に到達できますが、sステップAdams–Bashforth法は次数sしかありません。 ${\displaystyle a_{s-1}=-1}$ ${\displaystyle a_{s-2}=\cdots =a_{0}=0}$ ${\displaystyle b_{s}=0}$ ${\displaystyle s+1}$

s = 0, 1, 2, 3, 4 の場合の Adams–Moulton 法がリストされています (Hairer, Nørsett & Wanner 1993, §III.1; Quarteroni, Sacco & Saleri 2000)。最初の 2 つの方法は、それぞれ後退オイラー法と台形法です。 $${\displaystyle {\begin{aligned}y_{n}&=y_{n-1}+hf(t_{n},y_{n}),\\y_{n+1}&=y_{n}+{\frac {1}{2}}h\left(f(t_{n+1},y_{n+1})+f(t_{n},y_{n})\right),\\y_{n+2}&=y_{n+1}+h\left({\frac {5}{12}}f(t_{n+2},y_{n+2})+{\frac {8}{12}}f(t_{n+1},y_{n+1})-{\frac {1}{12}}f(t_{n},y_{n})\right),\\y_{n+3}&=y_{n+2}+h\left({\frac {9}{24}}f(t_{n+3},y_{n+3})+{\frac {19}{24}}f(t_{n+2},y_{n+2})-{\frac {5}{24}}f(t_{n+1},y_{n+1})+{\frac {1}{24}}f(t_{n},y_{n})\right),\\y_{n+4}&=y_{n+3}+h\left({\frac {251}{720}}f(t_{n+4},y_{n+4})+{\frac {646}{720}}f(t_{n+3},y_{n+3})-{\frac {264}{720}}f(t_{n+2},y_{n+2})+{\frac {106}{720}}f(t_{n+1},y_{n+1})-{\frac {19}{720}}f(t_{n},y_{n})\right).\end{aligned}}}$$

アダムス・モールトン法の導出はアダムス・バッシュフォース法と似ているが、補間多項式では上記の点だけでなく、 も用いる。係数は次のように与えられる。 ${\displaystyle t_{n-1},\dots ,t_{n-s}}$ ${\displaystyle t_{n}}$ $${\displaystyle b_{s-j}={\frac {(-1)^{j}}{j!(s-j)!}}\int _{0}^{1}\prod _{i=0 \atop i\neq j}^{s}(u+i-1)\,\mathrm {d} u,\qquad {\text{for }}j=0,\ldots ,s.}$$

アダムズ・モールトン法は、アダムズ・バッシュフォース法と同様に、ジョン・カウチ・アダムズにのみ由来する。フォレスト・レイ・モールトンの名がこれらの方法と結びついたのは、彼がこれらの方法がアダムズ・バッシュフォース法と併用して予測子・修正子のペアとして使えることに気づいたためである（Moulton 1926）。ミルン（1926）も同様のアイデアを持っていた。アダムズはニュートン法を用いて陰方程式を解いた（Hairer, Nørsett & Wanner 1993, §III.1）。

後方微分公式（BDF）

BDF法は、とその他の係数を、 次数s （可能な限り最大の次数）に達するように選択した暗黙法です。これらの法は、特にスティッフな微分方程式を解くために使用されます。 ${\displaystyle b_{s-1}=\cdots =b_{0}=0}$

分析

線形多段階法、そして実際には微分方程式のあらゆる数値解法の解析における中心概念は、収束、順序、および安定性です。

一貫性と秩序

最初の疑問は、この方法に一貫性があるかどうか、つまり差分方程式が微分方程式の良い近似になっているかどうかです。より正確には、h がゼロになるときに、ローカル打ち切り誤差がステップ サイズhよりも速くゼロになる場合、マルチステップ法は一貫性があります。ここで、ローカル打ち切り誤差は、以前の値がすべて正確であると仮定した方法の結果と、時刻 における方程式の正確な解との差として定義されます。テイラー級数を使用した計算により、線形マルチステップ法が一貫性を持つのは、次の場合のみであることが示されています。上記の方法はすべて一貫性があります (Hairer、Nørsett、Wanner 1993、§III.2)。 $${\displaystyle {\begin{aligned}&a_{s}y_{n+s}+a_{s-1}y_{n+s-1}+a_{s-2}y_{n+s-2}+\cdots +a_{0}y_{n}\\&\qquad {}=h{\bigl (}b_{s}f(t_{n+s},y_{n+s})+b_{s-1}f(t_{n+s-1},y_{n+s-1})+\cdots +b_{0}f(t_{n},y_{n}){\bigr )},\end{aligned}}}$$ ${\displaystyle y'=f(t,y)}$ ${\displaystyle y_{n+s}}$ ${\displaystyle y_{n+s-1},\ldots ,y_{n}}$ ${\displaystyle t_{n+s}}$ $${\displaystyle \sum _{k=0}^{s-1}a_{k}=-1\quad {\text{and}}\quad \sum _{k=0}^{s}b_{k}=s+\sum _{k=0}^{s-1}ka_{k}.}$$

手法が整合的であれば、次の問題は、数値解法を定義する差分方程式が微分方程式をどの程度近似するかです。多段階法は、hがゼロに近づくにつれて局所誤差がp次となる場合、p次 を有すると言われます。これは、手法の係数に関する以下の条件と等価です。sステップAdams –Bashforth法はs次ですが、sステップAdams–Moulton法はp次です(Hairer, Nørsett & Wanner 1993, §III.2)。 ${\displaystyle O(h^{p+1})}$ $${\displaystyle \sum _{k=0}^{s-1}a_{k}=-1\quad {\text{and}}\quad q\sum _{k=0}^{s}k^{q-1}b_{k}=s^{q}+\sum _{k=0}^{s-1}k^{q}a_{k}{\text{ for }}q=1,\ldots ,p.}$$ ${\displaystyle s+1}$

これらの条件は、多くの場合、特性多項式を使用して定式化されます 。これらの多項式を使用すると、メソッドがp次を持つための上記の条件は次のようになります 。特に、メソッドが少なくとも 1 次を持つ場合、つまり およびの場合に、メソッドは一貫しています。 $${\displaystyle \rho (z)=z^{s}+\sum _{k=0}^{s-1}a_{k}z^{k}\quad {\text{and}}\quad \sigma (z)=\sum _{k=0}^{s}b_{k}z^{k}.}$$ $${\displaystyle \rho (e^{h})-h\sigma (e^{h})=O(h^{p+1})\quad {\text{as }}h\to 0.}$$ ${\displaystyle \rho (1)=0}$ ${\displaystyle \rho '(1)=\sigma (1)}$

安定性と収束

1ステップ法の数値解は初期条件に依存しますが、 sステップ法の数値解はs個の開始値すべてに依存します。したがって、数値解が開始値の摂動に対して安定しているかどうかが重要です。線形多段階法は、特定の時間間隔における特定の微分方程式に対して零安定であるとは、サイズ ε の開始値の摂動によって、その時間間隔における数値解が、ステップサイズhに依存しないKの値に対してK ε 以下しか変化しない場合を指します。これは、微分方程式の条件を確認すれば十分であるため、「零安定性」と呼ばれます(Süli & Mayers 2003, p. 332)。 ${\displaystyle y_{0}}$ ${\displaystyle y_{0},y_{1},\ldots ,y_{s-1}}$ ${\displaystyle y'=0}$

特性多項式ρの根がすべて 1 以下の絶対値を持ち、かつ絶対値 1 の根の重複度が 1 である場合、根条件が満たされるという。線形多段階法が零安定であるのは、根条件が満たされる場合のみである(Süli & Mayers 2003, p. 335)。

ここで、十分に滑らかな微分方程式に、矛盾のない線形多段階法を適用し、すべての初期値がとして初期値に収束すると仮定します。すると、数値解は、この方法が零安定であるかのように、かつその場合に限り、厳密解に収束します。この結果は、ゲルムント・ダールキストにちなんで名付けられたダールキスト同値定理として知られています。この定理は、有限差分法のラックス同値定理と精神的に類似しています。さらに、この方法がp次を持つ場合、グローバル誤差（固定時刻における数値解と厳密解の差）は です（Süli & Mayers 2003、p. 340）。 ${\displaystyle y_{1},\ldots ,y_{s-1}}$ ${\displaystyle y_{0}}$ ${\displaystyle h\to 0}$ ${\displaystyle h\to 0}$ ${\displaystyle O(h^{p})}$

さらに、この方法が収束する場合、が法1の唯一の根であるとき、この方法は強安定であると言われます。この方法が収束し、法1のすべての根が重複していないが、そのような根が複数存在するとき、この方法は相対安定であると言われます。この方法が収束するためには、1が根でなければならないことに注意してください。したがって、収束法は常にこれら2つのいずれかです。 ${\displaystyle z=1}$

スティッフな方程式に対する線形マルチステップ法のパフォーマンスを評価するために、線形テスト方程式y' = λ yについて考えます。ステップ サイズhでこの微分方程式にマルチステップ法を適用すると、特性多項式を持つ線形回帰関係が得られます。この多項式は、マルチステップ法の安定性多項式と呼ばれます。その多項式のすべての根の絶対値が 1 未満の場合、マルチステップ法の数値解は 0 に収束し、マルチステップ法はそのh λの値に対して絶対安定していると言えます。実部が負のすべてのh λに対して絶対安定している場合、その方法はA 安定であると言えます。絶対安定の領域は、マルチステップ法が絶対安定であるすべてのh λの集合です(Süli & Mayers 2003、347 & 348 ページ)。詳細については、スティッフな方程式とマルチステップ法のセクションを参照してください。 $${\displaystyle \pi (z;h\lambda )=(1-h\lambda \beta _{s})z^{s}+\sum _{k=0}^{s-1}(\alpha _{k}-h\lambda \beta _{k})z^{k}=\rho (z)-h\lambda \sigma (z).}$$

例

アダムズ・バッシュフォースの3段階法を考える。したがって、一つの特性多項式は であり 、その根は であり、上記の条件を満たす。は法1の唯一の根であるため、この法は強安定である。 $${\displaystyle y_{n+3}=y_{n+2}+h\left({23 \over 12}f(t_{n+2},y_{n+2})-{4 \over 3}f(t_{n+1},y_{n+1})+{5 \over 12}f(t_{n},y_{n})\right).}$$ $${\displaystyle \rho (z)=z^{3}-z^{2}=z^{2}(z-1)}$$ ${\displaystyle z=0,1}$ ${\displaystyle z=1}$

もう一つの特性多項式は $${\displaystyle \sigma (z)={\frac {23}{12}}z^{2}-{\frac {4}{3}}z+{\frac {5}{12}}}$$

第1および第2ダルキスト障壁

これら2つの結果はゲルムント・ダールキストによって証明され、線形多段階法の収束順序とA安定性に関する重要な上界を示している。最初のダールキスト障壁はダールキスト（1956）で、2番目のダールキスト障壁はダールキスト（1963）で証明された。

最初のダルキスト障壁

第一のダールキスト障壁は、零安定かつ線形なqステップ多段階法は、 qが奇数の場合はq + 1を超える収束次数を達成できず、 qが偶数の場合はq + 2を超える収束次数を達成できないことを述べています。この方法が明示的である場合、 qを超える収束次数を達成することはできません(Hairer, Nørsett & Wanner 1993, Thm III.3.5)。

第二ダルキスト障壁

第二のダルキスト障壁は、明示的な線形多段階法はどれもA安定ではないというものです。さらに、（暗黙的な）A安定線形多段階法の最大次数は2です。次数2のA安定線形多段階法の中で、台形則は誤差定数が最も小さいです（Dahlquist 1963、Thm 2.1および2.2）。

参照

• デジタルエネルギーゲイン

参考文献

• バッシュフォース、フランシス（1883）「毛細管現象の理論検証の試み：理論上の液滴形状と実測値の比較による。液滴形状の理論値を示す表を作成する際に用いられた積分法の説明付き」JCアダムズ著、ケンブリッジ{{citation}}: CS1 maint: location missing publisher (link)。
• ブッチャー、ジョン・C.（2003）、常微分方程式の数値解析法、ジョン・ワイリー、ISBN 978-0-471-96758-3。
• Dahlquist, Germund (1956)、「常微分方程式の数値積分における収束と安定性」、Mathematica Scandinavica、4 : 33– 53、doi : 10.7146/math.scand.a-10454。
• Dahlquist, Germund (1963)、「線形多段階法の特殊な安定性問題」、BIT、3 : 27–43、doi :10.1007/BF01963532、ISSN  0006-3835、S2CID  120241743。
• ゴールドスタイン、ハーマン H. (1977)、「16世紀から19世紀までの数値解析の歴史」、ニューヨーク：シュプリンガー・フェアラーク、ISBN 978-0-387-90277-7。
• ハイラー、エルンスト。ノーセット、シベール・ポール。 Wanner、Gerhard (1993)、常微分方程式の解法 I: Nonstiff 問題(第 2 版)、ベルリン: Springer Verlag、ISBN 978-3-540-56670-0。
• ヘアラー、エルンスト、ワーナー、ゲルハルト（1996年）、常微分方程式の解法II：スティフ問題と微分代数問題（第2版）、ベルリン、ニューヨーク：シュプリンガー・フェアラーク、ISBN 978-3-540-60452-5。
• イザーレス、アリエ（1996年）、微分方程式の数値解析入門、ケンブリッジ大学出版局、Bibcode：1996fcna.book.....I、ISBN 978-0-521-55655-2。
• ミルン, WE (1926)、「常微分方程式の数値積分」、アメリカ数学月刊誌、33 (9)、アメリカ数学会: 455– 460、doi :10.2307/2299609、JSTOR  2299609。
• モールトン、フォレスト・R.（1926）、『外部弾道学の新手法』、シカゴ大学出版局。
• クアルテローニ, アルフィオ;サッコ、リッカルド。サレリ、ファウスト (2000)、Matematica Numerica、Springer Verlag、ISBN 978-88-470-0077-3。
• Süli, Endre; Mayers, David (2003), 『数値解析入門』 , Cambridge University Press , ISBN 0-521-00794-1。

外部リンク

• ワイスタイン、エリック・W.「アダムス法」。MathWorld。

Retrieved from "https://en.wikipedia.org/w/index.php?title=Linear_multistep_method&oldid=1321447195"