球状分布電荷と核力による散乱

ルンゲクッタ法で球状分布電荷と核力による散乱

謝辞大学での Fortran77 での計算機実習でのコードを Java 化しただけです。

青山学院大学理工学部物理学科の先生方に厚く感謝したいと思います。

球状分布電荷と核力による散乱 球状に一様に分布している電荷によって散乱される荷電粒子の軌道については、ここを参照。

核力としては、次の力を仮定する。

無次元化すると、

内側の円が、半径 1.0 の円で、y の値(衝突係数)を入力。
外の円は、半径が 3.0 で、その中に一様に電荷が分布している。

(PARAM の dB は理論解(ラザフォード散乱)の衝突係数の幅)
(衝突係数 = 3.2422384 辺りが引力と斥力がつりあう辺り)
ソースコードなどはここ(snuc.lzh is 8,745byte)