物理学コミュの計算ソフト Python sf を作っています

Python sf と名付けた下の計算ソフトを作っています。
概要;;;http://www.nasuinfo.or.jp/FreeSpace/kenji/sf/sfShotConcise/sfShotConcise.htm
マニュアルhttp://www.nasuinfo.or.jp/FreeSpace/kenji/sf/pysf/manual/readMeOfPythonSf.htm

傲慢に聞こえるでしょうが、物理・工学の方には Matlab/Mathematica より使えるソフトにできたと自負してます。多くを scipy, sympy, vpython などの凄いソフトがあるからこそなのですが、自分の付け加えた部分もそれなりにあると自負しています。

Python sf の応用例として「量子論の行列による力学」を下のブログに書いてみました。

http://loboskobayashi.wordpress.com/category/quantum-mechanics/mechanics-with-matrix/
http://loboskobayashi.wordpress.com/2010/12/27/量子力学での行列力学-2/
http://loboskobayashi.wordpress.com/2010/12/27/位置・運動量-演算子行列-5/
http://loboskobayashi.wordpress.com/2010/12/27/調和振動子演算子：行列と固有値-4/
http://loboskobayashi.wordpress.com/2010/12/27/調和振動子での波動関数の変化-3/

位置演算子の行列 X,運動量演算子の行列 Px を実際に構成し、 Hamiltonian H = Px^2 + X^2 の固有値・固有ベクトルを検討しています。また波動関数の初期値を与えて、その時間変化 expm( i t H) ψ0 を計算させています。

添付図のような波動関数の時間変化を得られます。RGB の混ざり具合で位相回転も表しているので、複素数値波動関数の時間変化を表せています。位相の時間変化の様子をイメージできている方は少ないと思っています。

皆様の御意見を伺わせてください。