授業資料/Q のバックアップ(No.2)

各種例題 †

いくつか、やってみると良い課題を書いておこう．レポート課題として使ってみるのも良い．

f(x) = x^2 - 2 の近似根を求めるプログラムを書いて動かしてみよう．なお、解は明らかに ±√2 であるが、求めるのはこのどちらかだけで良い．方法は多数あるので、以下、列挙しておこう．全部やる必要はないが、複数やっておいたほうがいいだろう．
- 二分法で…
- はさみうち法で…
- 割線法で…
- 縮小写像法で(ただし、問題を少し書き換えないといけないだろう)
- Newton 法で…
f(x) = e^{-x} - x の近似根を、同様に求めてみよう．
f(x) = 3x^2 + 1 + (log (Pi - x))^2/(Pi^4) について近似根を求められるだろうか、チャレンジしてみよう．たぶん、ちょっと難しい．
連立非線形方程式である、
2x^2 + y^2 - 1 = 0,
x - (√3)y = 0,
の近似解を求めてみよう．たぶん Newton 法がいいだろう．
f(x) = e^{-x} - x の近似根を、(簡易版)ホモトピー法で求めてみよう．
f(x) = x^3 - 3x + 3 の近似根を、ホモトピー法で求めてみよう．

LU 分解を使って連立一次方程式の求解を行うプログラムを書いてみよう．
いきなり全体のプログラムを書くよりは、各種の線形計算、たとえば、ベクトルとベクトルの内積を計算するルーチンを書いてから全体を書いたほうが楽だろう．
同様に CG 法を使うプログラムを書いてみよう．

du/dt = u(1-u) という単純な方程式の近似解を求めるプログラムを書いて動かしてみよう．方法としては、
- Euler 法,
- Runge-Kutta 法,
- 線形多段階法, などがあるが、まずは Runge-Kutta 法でプログラムが組めるといいだろう．
連立常微分方程式
du/dt = (2-v)u,
dv/dt = (2u - 3)v,
の近似根を求めるプログラムを書いて、動かしてみよう．この問題の解は時間発展すると (u,v) 平面を一周して同じ点へ戻る性質があるが、その性質が「どれくらい再現されるか」グラフを描いて確かめてみよう．