微分の数値計算です。減衰振動をルンゲ・クッタ法で数値計算するプログラムです。
振幅が最初の1/2になるときの時間を求めたいのですが、どうもうまく行きません。
できれば早急に指摘お願いします…
http://kansai2channeler.hp.infoseek.co.jp/cgi-bin/joyful/img/1224.txt