音速の壁

規格化した値での計算はうまくいくのですが、実際の値（密度がかなり小さい）を入れるとおかしくなってしまいます。
亜音速から超音速に変わる付近で密度が大きく揺れる・・・。
うまく安定化させないとダメなようです。