ルンゲクッタ法で Rutherford 散乱
謝辞
大学での Fortran77 での計算機実習でのコードを Java 化しただけです。
青山学院大学理工学部物理学科の先生方に厚く感謝したいと思います。
Rutherford 散乱とは
中心にある点電荷によって散乱される荷電粒子の軌道。
点電荷 Q は原点に固定されているとして、電荷 q、質量 M を持ち、x 軸の負の領域から、x 軸に沿って、速度 v、y 座標 b で入射すると、
極座標 (r,φ)
1/r = sin(φ)/b - D(1 + cos(φ))/2b
2
D=qQ/2πε
0
Mv
基本方程式は、
中央の円が、半径 1.0 の円で、y の値(衝突係数)を入力。
赤線が、理論値
青線が、ルンゲクッタによる近似値
緑線が、ルンゲクッタによる近似値(あなたが入力した衝突係数)
(PARAM の dB は理論解の衝突係数の幅)
ソースコードなどはここ
(ruth.lzh is 8,070byte)
ルンゲクッタ法の JavaClass
mail to active@window.goukaku.com