※小説ではない※専門書 要約資料集 為替(換算)3.9万円でもらう 紐解集生成 専門 初入門 資料 作：｛作者名｝
> 学問の宇宙・形式科学の恒星群 数学冠・派生第1巻。ガイド役: Fable 5 監修 / Sonnet 5 執筆(脚本班)
> トーン規約: GAKUMON_UNIVERSE.md準拠。専門用語は初出で必ず説明する。
> 重複回避の方針: →BOOK-0063『微分積分』第1巻は「瞬間の変化(微分)」と「積み重ねの全体(積分)」という二つの道具そのものを扱う総論書だった。本冊はその道具を使って「何が起きるかを予言する式」を組み立てる技術、微分方程式(びぶんほうていしき)を扱う。微分の計算そのものの復習は最小限に留め、詳しい導出は→BOOK-0063『微分積分』第1巻に委ねる。
> 接続先: →BOOK-0063『微分積分』第1巻(微分・積分という道具そのもの)、→BOOK-0058『ばねと弾性』第1冊(単振動の力学的な起源)、→BOOK-0017『速度』第1冊(変化率という考え方の入口)、→BOOK-0064『電磁気』第1巻(場を記述する式の親戚)、→BOOK-0067『代数』第1巻(方程式を解くという営みそのもの)、→BOOK-0068『解析幾何』第1巻(解を座標平面上の曲線として見る視点)。
> 水準: 一〜四(基礎と応用の入口)。数値解法とカオスの入口までを扱い、多変数・偏微分方程式などの高度な話題は後半(水準五以降)の別巻に譲る。
# BOOK-0078 微分方程式 — 変化の法則を書く言語(数学派生 第1巻)
> 学問の宇宙・形式科学の恒星群 数学冠・派生第1巻。ガイド役: Fable 5 監修 / Sonnet 5 執筆(脚本班)
> トーン規約: GAKUMON_UNIVERSE.md準拠。専門用語は初出で必ず説明する。
> 重複回避の方針: →BOOK-0063『微分積分』第1巻は「瞬間の変化(微分)」と「積み重ねの全体(積分)」という二つの道具そのものを扱う総論書だった。本冊はその道具を使って「何が起きるかを予言する式」を組み立てる技術、微分方程式(びぶんほうていしき)を扱う。微分の計算そのものの復習は最小限に留め、詳しい導出は→BOOK-0063『微分積分』第1巻に委ねる。
> 接続先: →BOOK-0063『微分積分』第1巻(微分・積分という道具そのもの)、→BOOK-0058『ばねと弾性』第1冊(単振動の力学的な起源)、→BOOK-0017『速度』第1冊(変化率という考え方の入口)、→BOOK-0064『電磁気』第1巻(場を記述する式の親戚)、→BOOK-0067『代数』第1巻(方程式を解くという営みそのもの)、→BOOK-0068『解析幾何』第1巻(解を座標平面上の曲線として見る視点)。
> 水準: 一〜四(基礎と応用の入口)。数値解法とカオスの入口までを扱い、多変数・偏微分方程式などの高度な話題は後半(水準五以降)の別巻に譲る。
---
## 入口の物語 — 未来を予言する式
コーヒーカップから立ちのぼる湯気を見たことがあるだろう。淹れたてのコーヒーは、放っておけば必ず室温まで冷めていく。人口が増える国もあれば減る国もある。惑星は軌道を描いて太陽の周りを回り続け、振り子は規則正しく左右に揺れる。バネにおもりを吊るして引っ張ると、手を離した瞬間から上下に振動を始める。
これらはまったく異なる現象に見える。しかし、もし「今この瞬間の状態から、次の瞬間にどれだけ変化するか」という規則さえ分かっていれば、その先どうなるかを丸ごと予言できるとしたらどうだろうか。
これこそが**微分方程式(びぶんほうていしき、水準一: ある量そのものと、その量の変化率(微分)との関係を表した方程式のこと。この式を「解く」と、時間とともにその量がどう変わっていくかを表す関数が得られる)**という道具の正体である。
→BOOK-0067『代数』第1巻で学んだ普通の方程式、たとえば `2x + 3 = 11` は、「x という**数**」を求める式だった。ところが微分方程式が求めるのは数ではない。微分方程式が求めるのは「**関数そのもの**」である。つまり、時間とともに変化していく量の、変化の全履歴を一気に言い当てる式なのだ。
この道具を使うと、次のようなとても異なる問いに、驚くほど似た数式で答えられる。
- バクテリアの数は、なぜ最初はゆっくり増え、途中から爆発的に増えるのか。
- 放射性物質は、なぜ「半分に減るまでの時間」が常に一定なのか。
- 熱いコーヒーは、なぜ最初は急激に冷め、室温に近づくほどゆっくり冷めるのか。
- バネにつるしたおもりは、なぜ行ったり来たりを規則正しく繰り返すのか。
- 太陽と地球と月という、たった3つの天体の未来の動きが、なぜ厳密には予言できないのか。
本冊では、この「変化の法則を書く言語」を、実際に手を動かして検算しながら学んでいく。
```
[現在の状態] + [変化の規則(微分方程式)] ──解く──> [未来の全履歴(関数)]
```
---
## 第一階段(水準一) — 最も簡単な微分方程式 dy/dx = ky
### 微分の復習は最小限に
→BOOK-0063『微分積分』第1巻ですでに学んだ通り、**微分(びぶん、水準一: ある量が別の量に対してどれくらいの割合で変化しているかを表す量)**は `dy/dx` という記号で書かれる。これは「x が変化したときに、y がどれくらいの割合で変化するか」を表している。本冊ではこの記号の意味そのものの詳しい説明は→BOOK-0063に譲り、ここでは「dy/dx は y の変化率である」という事実だけを使って先に進む。
### 「増えるほど速く増える」という式
さて、次のような式を考えてみよう。
```
dy/dx = k y
```
これは日本語に訳すと、「y の変化率(dy/dx)は、その時点の y の値に比例する(k倍である)」という意味である。k は比例定数(ひれいていすう、水準一: 二つの量の比が常に一定であるときの、その一定の比の値)と呼ばれる、決まった数だ。
この式が言っていることを具体的に想像してみよう。たとえばバクテリアの数を y としよう。dy/dx = ky (ここでは x の代わりに時間 t を使い、dy/dt = ky と書くことが多い)は、「今いるバクテリアの数が多ければ多いほど、次の瞬間に増える数も多い」ということを意味する。10匹のバクテリアより100匹のバクテリアの方が、同じ時間で生まれる子供の数が多いのは直感的にも納得できるだろう。**「増えれば増えるほど、増える速さも速くなる」**という、雪だるま式の変化を表しているのがこの式である。
### 答えは指数関数
この微分方程式 `dy/dx = ky` を満たす関数は、次の形をしていることが知られている。
```
y = C・e^(kx)
```
ここで `e` は**ネイピア数(水準二: およそ2.71828...という値を持つ、自然対数の底と呼ばれる特別な定数)**、C は最初の量(x=0のときの y の値)を表す定数である。この `e^(kx)` という形の関数を**指数関数(しすうかんすう、水準一: xが変数として指数の位置にある関数)**と呼ぶ。
### 検算その1 — y = e^x が dy/dx = y を満たすことの確認
なぜ `y = e^(kx)` がこの微分方程式の答えなのか、まずは最も単純な k=1 の場合、つまり `y = e^x` について確かめてみよう。→BOOK-0063『微分積分』第1巻の手順に従えば、微分は「x をわずかに Δx だけ動かしたときの変化率の極限」として計算できる。
理屈だけでなく、実際に小さな数を代入して数値で確かめる。x = 1 の付近で、Δx = 0.0001 という非常に小さな変化を考える。
```
y(1) = e^1 ≈ 2.718281828
y(1.0001) = e^1.0001 ≈ 2.718553674
差 ≈ 0.000271846
変化率 = 差 / 0.0001 ≈ 2.71846
```
真の値 e^1 ≈ 2.718281828 と、計算した変化率 2.71846 を比べると、ごくわずかな誤差(Δxを有限の小さな値にしたことによる近似誤差)を除いてほぼ一致している。つまり `y = e^x` という関数は、「その値そのもの」と「その変化率」が常にほとんど等しいという、非常に特別な性質を持っている。これが `dy/dx = y` を満たすということの数値的な裏付けである。`k` が1以外の値でも、`y = e^(kx)` の変化率は `k・e^(kx) = ky` となることが同様の手順で確かめられる(詳細な導出は→BOOK-0063『微分積分』第1巻を参照)。
---
## 休憩所①(ここまでのまとめ箱)
- 微分方程式とは、ある量とその変化率の関係を表す式であり、これを解くと未来の全履歴を表す関数が得られる。
- 最も基本的な微分方程式は `dy/dx = ky`(変化率が量そのものに比例する)であり、答えは指数関数 `y = C・e^(kx)` になる。
- `y = e^x` の変化率が y 自身にほぼ等しいことは、Δx=0.0001での数値差分計算で検算済み。
次の階段では、この指数関数がマルサスの人口論と放射性崩壊という、まったく違う二つの現象で同じ数式として使われることを見る。
---
## 第二階段(水準二) — 人口爆発と半減期、同じ式の二つの顔
### マルサスの人口論 — 増えるほど速く増える人口
k が正の数、つまり `dy/dt = ky`(kは正)の場合、y はどんどん加速しながら増えていく。この考え方を歴史上はじめて経済・社会の議論に持ち込んだのが、イギリスの経済学者トマス・ロバート・マルサス(1766年〜1834年)である。マルサスは1798年に『人口論』を著し、「人口は制約がなければ等比数列的(幾何級数的)に増加するが、食料生産は等差数列的(算術級数的)にしか増加しない」という主張を展開した。これは、人口の増加を `dy/dt = ky` の形で捉える発想の先駆けとして知られている。
もっとも、マルサスの主張がそのまま現実の人口動態を正確に予言できたかどうかについては、後の歴史における農業技術の急速な進歩などもあり、経済学・人口学の分野で議論が続いている。この点は**諸説あり**、本冊では「1798年に『人口論』が発表され、指数関数的な増加という考え方を広めた」という事実のみを断定し、その予言がどこまで的中したかという評価には立ち入らない。
### 放射性崩壊 — 減るほどゆっくり減る
k が負の数の場合、つまり `dy/dt = -ky`(kは正の定数として、マイナス符号を式の外に出す書き方をする)の場合はどうなるか。この場合、y はどんどん減っていくが、減るペースも y 自身が小さくなるにつれて遅くなっていく。
これは放射性崩壊(ほうしゃせいほうかい、水準二: 不安定な原子核が、時間とともに一定の割合で別の原子核に変わっていく現象)という自然現象そのものである。放射性物質の量を N(t) とすると、次の式が成り立つ。
```
dN/dt = -kN
```
この式を解くと `N(t) = N₀・e^(-kt)` となる(N₀は最初の量)。この現象で特に重要なのが**半減期(はんげんき、水準二: 放射性物質の量がちょうど半分になるまでにかかる時間)**という概念である。
### 検算その2 — 半減期を2回経過すると量は1/4になる
半減期を T とすると、定義から `N(T)/N₀ = 1/2` が成り立つ。ここで「半減期を2回経過した後」、つまり時刻 2T での量を計算してみよう。
```
N(2T) = N₀・e^(-k・2T) = N₀・(e^(-kT))² = N₀・(1/2)² = N₀・(1/4)
```
具体的な数値でも確かめる。k = 0.1(1秒あたりの崩壊の割合と仮定する)としたとき、半減期は `T = ln2 / k` という式で求まる(ln2は2の自然対数、およそ0.6931)。
```
T = 0.6931 / 0.1 ≈ 6.931(秒)
N(2T) / N₀ = e^(-0.1 × 2 × 6.931) = e^(-1.3863) ≈ 0.2500
```
計算結果は 0.2500、つまりぴったり1/4になっている。半減期という考え方の核心は、「半分になるまでの時間」がどの時点から測っても常に同じ T であることだ。100個が50個になるのにかかる時間と、50個が25個になるのにかかる時間は、まったく同じ T である。これは`(1/2)`を何回掛けても、その「1回分の時間」が変わらないという、指数関数ならではの性質による。
同じ `dy/dx = ky` という一本の式が、k の符号が違うだけで「爆発的に増える人口」にも「規則正しく半分になっていく放射性物質」にもなる。これが微分方程式という言語の経済性であり、本冊冒頭で述べた「一つの言語で全く違う現象を語れる」ことの最初の実例である。
```
k > 0: y = C・e^(kt) 時間とともに急上昇 ─────╮
├─ 同じ形の式
k < 0: y = C・e^(kt) 時間とともに0へ漸近 ──────╯
グラフの概形(方向場のイメージ):
y
|＼ ／
| ＼ ／
| ＼ (k<0: 崩壊) (k>0: 増殖)／
| ＼＿＿＿＿ ＿＿＿／
+-------------------------------- t
```
---
## 休憩所②(ここまでのまとめ箱)
- k > 0 の `dy/dt = ky` は爆発的な増加を表し、マルサスが1798年『人口論』で人口の増え方をこの発想で論じた(その予言の的中度は諸説あり)。
- k < 0 の `dy/dt = -ky` は放射性崩壊を表し、半減期という「半分になるまでの一定時間」が生まれる。
- 半減期を2回経過すると量は (1/2)² = 1/4 になることを、k=0.1・T≈6.931秒の具体的な数値で検算済み。
次の階段では、コーヒーが冷める現象を微分方程式で表し、さらに「行ったり来たり」を繰り返すバネと振り子の運動を見ていく。
---
## 第三階段(水準二〜三) — 冷めるコーヒーと、揺れるバネ
### ニュートンの冷却法則
淹れたてのコーヒーはなぜ、最初は急激に冷めて、ぬるくなってからはなかなか冷めきらないのだろうか。この現象を説明するのが**ニュートンの冷却法則(れいきゃくほうそく、水準二: 物体の温度が周囲の温度に近づく速さは、物体と周囲の温度差に比例するという法則)**である。アイザック・ニュートン(1642年〜1727年、イングランドの自然哲学者)がこの法則を論じたことにちなんでこう呼ばれている。
コーヒーの温度を T(t)、周囲の室温を `T_env`(一定とする)とすると、この法則は次の微分方程式で書ける。
```
dT/dt = -k (T - T_env)
```
これは「コーヒーの温度の変化率は、コーヒーと室温の差に比例する」という意味だ。コーヒーが熱々(室温との差が大きい)なときほど急激に冷め、室温に近づく(差が小さくなる)ほど冷め方はゆっくりになる。この式を解くと次のようになる。
```
T(t) = T_env + (T₀ - T_env)・e^(-kt)
```
T₀ は最初のコーヒーの温度である。試しに T₀ = 90℃(熱々のコーヒー)、`T_env` = 20℃(室温)、k = 0.05(1分あたりの冷却の割合と仮定)として、時間ごとの温度を計算してみよう。
```
t=0分: T = 20 + 70×e^0 = 20 + 70.00 = 90.00℃
t=10分: T = 20 + 70×e^-0.5 = 20 + 42.46 = 62.46℃
t=20分: T = 20 + 70×e^-1.0 = 20 + 25.75 = 45.75℃
t=30分: T = 20 + 70×e^-1.5 = 20 + 15.62 = 35.62℃
```
最初の10分で90℃から62.46℃まで、27.5℃以上も下がっているのに対し、次の10分(10分〜20分)では62.46℃から45.75℃と、下がり幅は16.7℃程度に減っている。さらに次の10分(20分〜30分)では下がり幅は10.1℃程度とさらに小さくなる。温度差が小さくなるほど冷める速さが遅くなる、という法則通りの振る舞いが数値で確認できる。
### バネとおもり — 単振動という「行ったり来たり」
→BOOK-0058『ばねと弾性』第1冊で学んだように、バネにおもりをつるして引っ張ると、手を離した瞬間からおもりは上下に振動を始める。この運動は、ニュートンの運動方程式(質量×加速度=力)と、バネの伸びに比例して働く復元力(フックの法則)を組み合わせることで、次の微分方程式として書ける。
```
m (d²x/dt²) = -k x
```
ここで x はおもりの、つり合いの位置からのずれ、m はおもりの質量、k はバネ定数(バネの硬さを表す定数)である。`d²x/dt²` は**二階微分(にかいびぶん、水準三: ある量を2回続けて微分したもの。位置を2回微分すると加速度になる)**と呼ばれ、位置 x を時間で2回微分したもの、つまり加速度を表す。
この式が言っているのは、「おもりがつり合いの位置から離れれば離れるほど、それを引き戻そうとする力(加速度)が強くなる」ということだ。この微分方程式の答えは、驚くべきことに三角関数になる。
```
x(t) = A cos(ωt), ω = √(k/m)
```
A は最初にどれだけ引っ張ったか(振幅)、ω(オメガ)は**角振動数(かくしんどうすう、水準三: 振動が1秒間にどれだけの角度に相当する速さで揺れるかを表す量)**と呼ばれる量で、バネ定数と質量から `√(k/m)` として決まる。この運動を**単振動(たんしんどう、水準二: 一定の周期でくり返し行ったり来たりする、最も基本的な振動運動)**と呼ぶ。→BOOK-0058『ばねと弾性』第1冊ではこの運動を力学的な側面(復元力・エネルギー)から扱ったが、本冊ではこれを微分方程式の解として、数式の側面から捉え直す。
### 検算その3 — なぜ解が三角関数になるのか、数値で確かめる
「なぜ解が三角関数になるのか」を、具体的な数値を使って検算してみよう。質量 m=2kg、バネ定数 k=8 N/m とすると、`ω = √(8/2) = √4 = 2` である。このとき解は `x(t) = cos(2t)` となる(A=1とする)。
この関数が本当に元の微分方程式 `m(d²x/dt²) = -kx` を満たすか、t=0.3 の時点で数値的に確かめる。二階微分は、微小な幅 h=0.0001 を使って次の近似式で計算できる。
```
d²x/dt² ≈ (x(t+h) - 2x(t) + x(t-h)) / h²
```
実際に計算すると、
```
x(0.2999) = cos(2×0.2999) ≈ 0.573059
x(0.3000) = cos(2×0.3000) ≈ 0.572484
x(0.3001) = cos(2×0.3001) ≈ 0.571909
数値二階微分 ≈ -3.301342
理論値(-ω²・x(0.3) = -4×0.572484) ≈ -3.301342
```
両者はぴったり一致している。さらに `m・(d²x/dt²) + k・x` を計算すると、
```
2 × (-3.301342) + 8 × 0.572484 ≈ -6.602684 + 6.602684 ≈ 0(誤差の範囲内)
```
見事に0になった。これは `m(d²x/dt²) = -kx`、すなわち `m(d²x/dt²) + kx = 0` が成り立っていることの数値的な検算である。
### 検算その4 — 周期の次元(単位)を確かめる
単振動の周期は `T = 2π√(m/k)` という式で表される。この式が本当に「時間」の単位を持つか、次元(単位)だけを取り出して確かめてみよう。
質量 m の単位は kg、バネ定数 k の単位は N/m(ニュートン毎メートル、力÷距離)である。まず `m/k` の単位を計算する。
```
[m/k] の単位 = kg / (N/m) = kg・m / N
```
ここで、力の単位 N(ニュートン)は「質量×加速度」であり `N = kg・m/s²` と書ける。これを代入すると、
```
kg・m / (kg・m/s²) = kg・m × s²/(kg・m) = s²
```
`m/k` の単位は `s²`(秒の2乗)になった。これに平方根を取ると `√(m/k)` の単位はちょうど `s`(秒)になる。2π は単なる数(単位を持たない)なので、`T = 2π√(m/k)` 全体の単位も秒になる。周期という「時間の長さ」を表す量が、正しく時間の単位を持つ式で表されていることが、この次元だけの検算からも確認できる。
具体的な数値でも確かめておくと、先ほどの m=2kg, k=8N/m の場合、
```
T = 2π√(2/8) = 2π√(0.25) = 2π × 0.5 = π ≈ 3.14159(秒)
```
これは `ω = 2` から `T = 2π/ω = 2π/2 = π` として求めた値ともぴったり一致する。
### 振り子との接続
単振動は、振り子(おもりを糸でつるして揺らす運動)にも現れる。振れ幅が小さいとき、振り子の周期は `T = 2π√(L/g)` という、バネの式と非常によく似た形の式で表される(Lは糸の長さ、gは重力加速度)。長さ L=1メートルの振り子の周期を計算すると、`T = 2π√(1/9.8) ≈ 2.007`秒となる。バネの `k/m` が、振り子では `g/L` に置き換わっているだけで、数式の構造そのものは同一である。バネの振動も振り子の揺れも、「つり合いの位置から離れるほど強く引き戻される」という同じ論理構造から生まれる、同じ数式の異なる衣装なのである。
```
単振動の一般形: d²x/dt² = -ω²x
│
├─ バネ: ω² = k/m
└─ 振り子: ω² = g/L (振れ幅が小さいときの近似)
```
---
## 休憩所③(ここまでのまとめ箱)
- ニュートンの冷却法則 `dT/dt = -k(T-T_env)` により、コーヒーは室温との差が大きいほど急激に冷え、差が縮まるほどゆっくり冷える(90℃→10分後62.46℃→20分後45.75℃→30分後35.62℃と検算済み)。
- バネの運動 `m(d²x/dt²) = -kx` の解 `x(t)=A cos(ωt)` は、t=0.3での数値二階微分が理論値-3.301342と一致することで検算済み。
- 周期の式 `T=2π√(m/k)` は次元(単位)を計算すると正しく秒になることを検算済み(m=2kg,k=8N/mでT=π秒と数値一致)。
- 振り子の周期 `T=2π√(L/g)` はバネの式と同じ数学的構造を持つ。
次の階段では、すべての微分方程式がきれいな数式で解けるわけではないという事実と、その代わりに使われる数値解法、そして「3つの天体」ですら未来を完全には予言できないという驚くべき事実を見ていく。
---
## 第四階段(水準三〜四) — 解けない方程式と、一歩ずつ進む数値解法
### いつでもきれいに解けるわけではない
ここまで見てきた `dy/dx = ky` や `m(d²x/dt²) = -kx` は、指数関数や三角関数というすでに知られた関数を使って、ぴったりと式で解くことができた。これを**解析的に解ける(かいせきてきにとける、水準三: 既知の関数を組み合わせた式で、微分方程式の答えをぴったりと書き表せること)**という。
しかし、世の中の微分方程式の多くは、このように綺麗な式で解くことができない。たとえば、空気抵抗が速度の2乗に比例するような、より現実に近い運動の式や、複数の変数が絡み合う式になると、既知の関数を組み合わせただけでは答えを書き表せないことが非常に多い。この場合でも、微分方程式そのものは「その時点でどう変化するか」という規則を教えてくれる。そこで登場するのが**数値解法(すうちかいほう、水準三: 式でぴったりと解を求める代わりに、ごく短い時間ごとに変化を積み上げて、近似的に未来の値を計算していく方法)**である。
### オイラー法 — 一歩ずつ未来を予測する
最も基本的な数値解法が**オイラー法(水準三: 現在の値とその変化率を使って、ごく短い時間だけ先の値を近似的に予測し、それを繰り返して未来を追跡する方法)**である。レオンハルト・オイラー(スイスの数学者)にちなんで名付けられたこの方法は、次のような単純な考え方に基づいている。
「現在の位置」と「現在の変化率(傾き)」が分かっていれば、ごく短い時間 `h` だけ先の値は、次の式でおおよそ予測できる。
```
y(次) ≈ y(今) + h × (今の変化率)
```
これは、崖の上から一定の傾きで飛び出したボールが、ごく短い時間なら、ほぼ真っ直ぐ進むという近似に似ている。刻み幅 h を十分小さく取れば、この近似の積み重ねによって、解析的に解けない複雑な微分方程式でも、近似的な答えの推移を追いかけることができる。
### 検算その5(補足) — オイラー法を手計算で2歩進めてみる
オイラー法がどのように働くか、すでに答えを知っている `dy/dx = y`(つまり k=1 の場合。答えは `y = e^x`)を例にして、実際に手計算で確かめてみよう。初期条件を `y(0) = 1`、刻み幅を `h = 0.5` として、2歩分計算する。
**1歩目(x=0 → x=0.5):**
```
今の値: y(0) = 1
今の変化率: dy/dx = y = 1 (dy/dx=yなので、変化率は今のyの値そのもの)
次の値: y(0.5) ≈ y(0) + h × 1 = 1 + 0.5 × 1 = 1.5
```
**2歩目(x=0.5 → x=1.0):**
```
今の値: y(0.5) = 1.5
今の変化率: dy/dx = y = 1.5
次の値: y(1.0) ≈ y(0.5) + h × 1.5 = 1.5 + 0.5 × 1.5 = 1.5 + 0.75 = 2.25
```
2歩の計算の結果、x=1 における近似値は `y ≈ 2.25` となった。真の答え(解析的な解 `y=e^x` から) `e^1 ≈ 2.71828` と比べると、2.25という近似値はかなり粗く、誤差は約0.468(相対誤差にしておよそ17%)にもなる。これは刻み幅 h=0.5 が大きすぎるためで、刻み幅をもっと細かく(たとえば h=0.01 など)すれば、近似の精度は大きく向上することが知られている。この検算からわかる大事な教訓は、「オイラー法は正しい方向には進むが、刻み幅が粗いと誤差が積み重なる」ということだ。この誤差と刻み幅の関係を扱う理論は、数値解析という分野でさらに詳しく研究されている。
```
オイラー法のイメージ(方向場に沿って一歩ずつ進む):
y
| ・(真の解 y=e^x のカーブ)
| ・
| ・
| o (2.25: オイラー法2歩目)
| o (1.5: オイラー法1歩目)
| o (1.0: 出発点)
+------------------------------ x
0 0.5 1.0
```
### 三体問題 — 3つだけでも解けない
指数関数やオイラー法の話から、「複雑でもコンピュータで数値的に追いかければ何でも予言できる」と思うかもしれない。しかし、ここに数学史上もっとも驚くべき発見の一つがある。それが**三体問題(さんたいもんだい、水準四: 3つの天体が互いの重力で運動するとき、その未来の軌道を、指数関数や三角関数のような既知の関数の組み合わせだけでぴったりと表す一般解が存在しないという問題)**である。
2つの天体(たとえば太陽と地球だけ)が互いの重力だけで運動する場合、その軌道は解析的に、つまりきれいな式でぴったりと解くことができる(これを二体問題と呼ぶ)。ところが、そこにもう一つ、たとえば月を加えて3つの天体になった途端、事情が一変する。
フランスの数学者アンリ・ポアンカレ(1854年〜1912年)は、1880年代末から1890年代にかけて三体問題を研究し、一般には解析的な一般解が存在しないことを明らかにした。この研究のなかでポアンカレは、初期条件(最初の位置や速度)のごくわずかな違いが、時間が経つにつれて、まったく異なる軌道へと指数関数的に広がっていく現象があることを見出した。
この「わずかな初期条件の違いが、やがて予測不能なほど大きな違いに育っていく」性質は、後に**カオス(水準四: 決定論的な(偶然のいたずらが入り込む余地のない)方程式に従っているにもかかわらず、初期条件のわずかな違いによって、長期的な予測が事実上不可能になる現象)**と呼ばれる考え方の入口として位置づけられている。三体問題やカオスの厳密な数学的取り扱いは非常に高度であり、その全容は後半(水準五以降)の別巻に譲るが、ここで押さえておきたいのは次の一点である。「決まった規則(微分方程式)に従っているからといって、必ずしも遠い未来まで正確に予言できるとは限らない」。
なお、ポアンカレの三体問題研究の詳しい経緯(スウェーデン国王オスカル2世の懸賞論文をめぐる逸話など)には版によって細部の伝わり方に違いがあり、本冊では「ポアンカレが1880年代末〜1890年代に三体問題を研究し、一般解の非存在とカオス的な鋭敏性を明らかにした」という大枠のみを断定し、逸話の細部には立ち入らない。
---
## 休憩所④(ここまでのまとめ箱)
- すべての微分方程式が解析的に(きれいな式で)解けるわけではなく、その場合は数値解法を使う。
- オイラー法は「現在の値+刻み幅×変化率」を繰り返す最も基本的な数値解法である。`dy/dx=y, y(0)=1, h=0.5`で2歩計算すると `y(1.0)≈2.25`(真値e≈2.71828との誤差は約17%)と検算済み。
- 2つの天体の運動は解析的に解けるが、3つになると一般解は存在しない(ポアンカレ、1880年代末〜1890年代)。これがカオスという考え方の入口になっている。
---
## なぜ面白いのか — 一つの言語で語られる、すべての変化
微分方程式の面白さは、あまりにも異なって見える現象──人口の爆発的増加、放射性物質の静かな崩壊、コーヒーが冷めていく穏やかな変化、バネと振り子の規則正しい往復運動、そして3つの天体だけでも予言できなくなる宇宙の複雑さ──これらすべてが、同じ「変化の規則を書く」という一つの文法で語られるという点にある。
`dy/dx = ky` というたった一行の式が、バクテリアの増殖にも放射性崩壊にも当てはまる。`m(d²x/dt²) = -kx` という式が、バネにも振り子にも当てはまる。そして、その延長線上には、もはや式では書ききれない複雑さ、つまりカオスという壁が立ちはだかっている。
→BOOK-0063『微分積分』第1巻で学んだ微分・積分という二つの道具は、単独では「今この瞬間どうなっているか」を教えてくれるだけだった。しかし、その道具を「方程式」という形に組み上げた瞬間、私たちは未来を語る言語を手にすることになる。天気予報も、人口予測も、薬が体内で分解される速さの計算も、根っこにあるのはすべてこの微分方程式という同じ文法なのである。
---
## 次への一歩 — この先に待っている風景
本冊(第1巻)では、水準一〜四にあたる基礎と応用の入口──最も単純な `dy/dx=ky` の指数関数解、人口増加と放射性崩壊という対照的な二つの応用、コーヒーの冷却とバネ・振り子の単振動、そして解析的に解けない方程式に対する数値解法(オイラー法)と、三体問題・カオスの入口──を扱った。
次に待っているのは、水準五以降の風景である。複数の変数が絡み合う連立微分方程式、空間全体に広がる量を扱う偏微分方程式(熱の伝わり方や波の伝播を記述する)、そしてカオス理論のより厳密な数学的取り扱い。これらはCOVERAGE_MAPにおいて後続の巻で扱われる予定である。
コーヒーが冷めていく、その素朴な観察から出発した本冊の旅は、ここでいったん休憩する。だが、この言語を使いこなす先には、天気のような複雑系から宇宙の運行まで、あらゆる変化を記述しようとする、さらに広大な風景が待っている。
```
[微分積分]──変化と積み重ねを扱う道具そのもの(→BOOK-0063)
│
[微分方程式]──道具を使って未来を予言する式を組む ← 本冊はここ
│
[偏微分方程式・カオス理論(厳密版)]──空間に広がる変化、予測不能性の数理(後半)
```
---
## 参照リンク一覧(整合チェック用)
- →BOOK-0063『微分積分』第1巻(微分そのものの定義・導出、本冊冒頭・検算1で参照)
- →BOOK-0058『ばねと弾性』第1冊(単振動の力学的起源、第三階段で参照)
- →BOOK-0017『速度』第1冊(変化率という考え方の入口、入口の物語で参照)
- →BOOK-0064『電磁気』第1巻(場を記述する方程式の親戚として、なぜ面白いのかで言及)
- →BOOK-0067『代数』第1巻(方程式を解くという営みそのものとの対比、入口の物語で参照)
- →BOOK-0068『解析幾何』第1巻(解を座標平面上の曲線として見る視点、第三階段の図で参照)
- →DIC-0001 数学用語辞典(専門用語の相互参照)
- →DIC-0002 物理用語辞典(物理的な応用語彙の相互参照)
# BOOK-0078 微分方程式 — 変化の法則を書く言語(数学派生 第1巻)