微分方程式の数値的解法 (1) Excel のみの方法

前回、VBA から Excel のシートにアクセスする方法を学んだ。 この機能を用いると、様々な数学的な操作を Excel 上で行うことができるようになる。 ここでは「微分方程式の数値解法」を取り上げ、VBA (および Excel) 学習の総仕上げとする。 微分方程式の数値的に解くのは、「Excel のみを使う方法」と「VBA と Excel を組み合わせて使う方法」のどちらでもできる。 ここではその両方を取り上げ、その違いを体感してもらう。 準備 Euler 法で数値解を求めてみる 数値解と解析解を比較してみる 準備 本日取り扱う微分方程式は以下のものである。 この微分方程式は、ばね系に周期的な外力が加わっている系の運動方程式 において、質量 m、ばね定数 k、外力の大きさ F、外力の振動数 ω に対し、 以下の制限を設けたものである。 この微分方程式 (運動方程式) は 2 階微分を含