next up previous
: この文書について... : 計算物理学: No.4 : Euler 法の誤差

課題

課題 1.


\begin{displaymath}
\frac{d y}{d x} = 1- y^2, y(0) = 0
\end{displaymath} (10)

の数値解をオイラー法によって求めよ。この問題の解析解 $\tanh x$ と比較せ よ。刻み幅 $h$ に対する誤差の振舞いを確認せよ。


課題 2.

単振動の微分方程式

\begin{displaymath}
\frac{d^2 y}{d x^2} + y = 0
\end{displaymath} (11)

は2元1階常微分方程式系
$\displaystyle \frac{d y}{d x}$ $\textstyle =$ $\displaystyle z,$  
$\displaystyle \frac{d z}{d x}$ $\textstyle =$ $\displaystyle -y$ (12)

に置き換えられる。

この数値解をオイラー法によって求めよ(初期条件は適当に選べ)。結果を解析解 と比較せよ。刻み幅 $h$ に対する誤差の振舞いを確認せよ。



Kiyohide Nomura 平成17年6月6日