オイラー法とルンゲ・クッタ法を使って微分方程式を解こう

＞＞前回 今回は、微分方程式を数値的に解く方法として、オイラー法とルンゲ・クッタ法について説明します。 オイラー法 微分方程式の数値解法では、初期条件からわずかに値を増加させた点の値を計算し、さらに計算された点を元にわずかな値を増加させた点の値を計算します。その操作を繰り返します。 この際、オイラー法では微分の定義に基づいた（1）式を用います。 数値解法のプログラムでは、（1）式のf（x, y）の部分にsin（x）＋yなどの具体的な式が入りますが、f（x, y）が複雑な式だと（1）式自体が複雑になります。そこで、「feval」コマンドを用いて、あらかじめ定義したf（x, y）を（1）式では呼び出して計算する方法を使ってみます。fevalコマンドはfeval（関数, 変数1，変数2）という形で用いて、関数に変数1と変数2を代入した値を返します。コマンドウィンドウ上で、試してみましょう。