初期条件をうまく設定して、一定時間後の未来で収束するように計算しています。今回は固いばねと質点による振り子。これによって逆行列を解かなくてよくなり、シンプレクティック積分法が使える(エネルギーがほぼ厳密に保存する)という利点も。あとちょっとだけキモさを削減。
http://www.nicovideo.jp/watch/sm31293675