……まずい(-_-!)。いいかげんに書いたら
定常振動のはずが発振しやがった。

それはさておきこの問題、拡張性を無視すればそう難しくなく書けるが、
パッケージ化とか考えると(とっくに公開されてはいるのだが)
正規形に直した式を渡す部分はともかくも(C でいうなら
値をセットする領域の配列と、その配列を対象とした
関数のポインタを渡すような形になる)、
初期値をセットする部分でけっこうプログラム的に
嬉しくないことになる。
なんかこう、「こういう形で書け」みたいな情報は
ないのか?

ついでに>>407は
×とりあえず、f(x) + �冉(x)*�凅 で書いてもらって、
○とりあえず、f(t+�冲) = f(t) + f'(t)*�冲 の形で書いてもらって、
だ。