いろいろと手抜きですが、BashForth を使って、 を時間刻み δt=0.01、2次の Adams-Bashforth 法で解いてみました。 図は解析解と比較したグラフです。 プログラム adams.fs \ adams.fs \ fixed point real variable fpprec 10000 fpprec ! : i>fr fpprec @ *…