近似計算するプログラム

質問

近似計算法である４次のRunge-Kutta法を用い、プロット点を一斉にdos画面に表示したい場合、全体的にどのようなプログラムにすれば良いのでしょうか？最終的にはその多数のプロット点をエクセルに出力し、グラフを作りたいと思っています。
何を近似するのかは気にしていないので、色々探したところ、あるサイトから４次のRunge-Kutta法の一部分のプログラムを見つけ、以下のを使用することにしたのです。プログラムは全くの無知なので、どのくらいの量になるのか想像出来ませんが、もし短い容易なプログラムになるのであれば載せて頂けると有難いです。コンパイラはVBC++です。

反復 {
v1dt = v(x[j] ) * dt
v2dt = v(x[j] + v1dt/2 ) * dt
v3dt = v(x[j] + v2dt/2 ) * dt
v4dt = v(x[j] + v3dt ) * dt
x[j+1] = x[j] + (v1dt + 2*v2dt + 2*v3dt + v4dt)/6
}

回答No.1

neKo_deux · ベストアンサー

こちらのような書籍を購入しておくと、後々お役に立つと思います。

C言語による最新アルゴリズム事典
http://www.matsusaka-u.ac.jp/~okumura/algo/

--
> あるサイトから

こういうサイトのURLなども提示してくれると、いろいろ参考になるんですが…。

> もし短い容易なプログラムになるのであれば

比較的短いです。

Sinéad page - harm_osc_rk.c
http://www.maths.tcd.ie/~ryan/teaching/harm_osc_rk.c.html

参考URL：

http://www.matsusaka-u.ac.jp/~okumura/algo/

補足

早くのご回答、本当に有難う御座いました。この
http://www.maths.tcd.ie/~ryan/teaching/harm_osc_rk.c.html
にあるプログラムが、neKo deuxさんが作ってくださったプログラムですよね？
更にお勧めの書籍まで紹介して頂き、大変感謝しています。

>こういうサイトのURLなども提示してくれると、いろいろ参考になるんですが…。

確かに提示すべきでした。申し訳ありません。一応、そのサイトのURLをお教えしますので、そのサイトをneKo deuxさんに見て頂いた上でこのスレッドを締め切らせてもらいますね。

http://www.damp.tottori-u.ac.jp/~ooshida/edu/ode/euler2rk.html