数値計算の手順

1. u0 を仮定し分母 = 分子 = 0 とおいて Critical point
   と,そこでの速度勾配を求める(de l'Hopital の定理)

2. 運動方程式を Critical point から内側(太陽側)へ積分し
   u0 を再計算する

3. 再計算した u0 を用いて再び Critical point を求める.

4. u0 が 収束するまで 2,3 を繰り返す

結果 
反復計算の途中で負の平方根がでてきて計算がストップ 
→ Alfven 波が入ると方程式が複雑になるので Gas pressure
   のみの計算からはじめることにした