数値計算の手順

1. (nu)E を用いて,全てのグリッドにおける (nu) を計算する.

2. u(r) と T1を仮定し,エネルギーの式を外側から内側へ積分. 積分して得られる
   T0が境界条件のT0と合うように T1を決定する(はさみ打ち法)



3. 得られた T(r), dT(r)/dr を用いて運動方程式を Critical point  から両側境界へ積
  分する. 

4. 得られた u(r) を用いて 手順 2 と同様に T(r), dT(r)/dr を計算する.

5. 3,4 の積分を u0の変化率が定められた範囲内になるまで反復する.