こんにちは、
下記ＨＰのプログラムは、コンプトン散乱を計算したものです。
通常は、
sl[q] = (?u[0]*q0 + ?u[1]*(-q1) + ?u[2]*(-q2) + ?u[3]*(-q3) + ms);
で計算して
{p1 -> 0, p2 -> 0, k0 -> p3, k1 -> 0, k2 -> 0, k3 -> -p3, q0 -> p0, q1 -> p3*Sqrt[1 - z^2],で、纏めます。
すると、正常の解が得られます。

ここで、
sl[q] = (?u[0]*q0^0.5 + ?u[1]*(-q1^0.5) + ?u[2]*(-q2^0.5) + ?u[3]*(-q3^0.5) + ms);
{p1^0.5 -> 0, p2^0.5 -> 0, k0^0.5 -> p3, k1^0.5 -> 0, k2^0.5 -> 0, k3^0.5 -> -
のように変更します。
すると、全く違う答えになってしまいます。どこが間違っているのでしょうか？
計算時間は５分程度です。

http://www.geocities.jp/dirac_equation/index.htm