(19) |
について、初期条件
(20) |
なお、牧野が書いたプログラムが
http://grape.astron.s.u-tokyo.ac.jp/~makino/pcphysics/programs/index.htmlの hermonic5a.C にあるので、今回はこれをコピーしてコンパイル・実行する だけでも一応いいことにする。但し、これは古い C++ の規格にそって書かれ ているので色々変更が必要である。
但し、シュテルマー公式とは、2階の微分方程式用の公式である。時間の2階微 分に、偏微分方程式の時に使ったのと同じ3点を使う中心差分を使う。