next up previous contents
Next: 5.10.2.2 Runge-Kutta 法の絶対安定領域 Up: 5.10.2 スキームの数値的安定性 Previous: 5.10.2 スキームの数値的安定性

5.10.2.1 中点公式の不安定性

中心差分近似に基づく中点公式 $x_{n+1}=x_{n-1}+h f(t_n,x_n)$$x'=-\lambda x$, $x(0)=1$ を解く( $\lambda = -4.35$)。左側は $x_1=e^{-\lambda h}$ としたもの、右側は $x_1=x_0+h f(t_0,x_0)$ とした もの。

\includegraphics[height=7cm]{figure/test-4.35.eps}     \includegraphics[height=7cm]{figure/test-4.35-by-midpoint2.eps}


next up previous contents
Next: 5.10.2.2 Runge-Kutta 法の絶対安定領域 Up: 5.10.2 スキームの数値的安定性 Previous: 5.10.2 スキームの数値的安定性
Masashi Katsurada
平成17年6月2日