直列多重振り子シミュレーション(改良版)

http://www.aihara.co.jp/~taiji/pendula-equations/present.html さんリスペクトです。
逆行列計算をLU分解にしました。コードはhttp://code.google.com/p/face-recognition-library-as3/source/browse/trunk/FaceRecognitionLib/src/com/karthik/math/lingalg/LUDecomposition.as?r=3を参考にさせて頂きました。
数が多くなっても軽い。