ただのダジャレですが、
bash と forth を使って Adams-Bashforth 法を実装してみるというのはどうだろうか。
追記: BashForth! そういうのもあるのか。
追記: む。BashForth には実数がないのか?
追記: とりあえず固定小数点数
\ fixedpoint.fs \ fixed point real variable fpprec 10000 fpprec ! : i>fr fpprec @ * ; : fr. fpprec @ /mod . ." ." . ; : fr* fpprec @ */ ; : fr/ fpprec @ swap */ ; \ test \ print 1.0 1 i>fr fr. cr \ 2.0*2.0 2 i>fr dup fr* fr. cr \ 1.0/2.0 1 i>fr 2 i>fr fr/ fr. cr bye
% ./bashforth-0.58a fixedpoint.fs 0 .1409 4 .1409 0 .1409
あれ? Gforth なら上手くいくのに。
% gforth fixedpoint.fs 1 .0 4 .0 0 .5000
追記: emit 使えばいいか。
: fr. fpprec @ /mod . 46 emit . ;
追記: いやこれじゃ負の数とか小数部が0.1未満の場合とかダメじゃないか。もうちょっとまじめにやる。
: fr. dup 0< if 45 emit negate then fpprec @ /mod . 46 emit 5 0 do 10 * fpprec @ /mod 48 + emit loop ;