Mathematica

■Mathematica
Mathematicaでいっぱい遊ぶ．スゴく楽しい．計算計算．

■調和振動子ポテンシャル中の波動関数とか
理論値は
n+0.5
結構良い．ずれるのはRunge-Kuttaで微分方程式を積分するときに初期値の取り方が悪いせい．eigenvalueに対して初期値を変えることが必要．
# eigenvalue[0] is 5.000000e-01
# eigenvalue[1] is 1.500000e+00
# eigenvalue[2] is 2.500000e+00
# eigenvalue[3] is 3.500001e+00
# eigenvalue[4] is 4.500012e+00
# eigenvalue[5] is 5.500094e+00
# eigenvalue[6] is 6.500576e+00
# eigenvalue[7] is 7.502852e+00
# eigenvalue[8] is 8.511172e+00
# eigenvalue[9] is 9.535938e+00