gnuplot，MATLABを使ったグラフの作成

プログラミングの方法として，可能な限り単純な仕様 を満たすものをまず書くべきである．今回は数値シミュレーションの正当性をま ずチェックすることにした．文献[1]の結果を参考にした．特に回転 エネルギーが時間に関して不変であることを確かめた．一般論として，実時間 アニメーションを行なうための複雑なGUI(Grahical user interface)は最初は避 けるのがよい．とりあえずシミュレーションの計算結果をちょっと見て みよう，というとき，グラフ作成ツールgnuplotを使うのが良い考えである．我々 はMATALAB が使える環境にあるのでMATLAB を使った．次の図 2.1$\sim$2.4は計算結果をMATLABでグラフ化し PostScript(EPS)ファイルに変換したものである．

図 2.1: 角運動量ベクトルの成分$La$の符号は変化しない

図 2.2: 角運動量ベクトルの成分$Lb$の符号は変わる．

図 2.3: 角運動量ベクトルの成分$Lc$の時間変化　　　　　　　　　 　　　　　　　　　　　　　　　　

図 2.4: c軸がz軸となす角の余弦の時間変化．　-1から1までで変化してい るからペンチが上下反転している．