流过平板的超声速流动的C++求解求助

[xls]

最近在做安德森的计算流体力学第十章的流过平板的超声速流动，用C++写的，但是迭代到一百多步会报nan的计算错误，一直没有找到代码的问题，求助哪位大神能指出大概是什么的问题
附上求解过程的输出的密度变化图：
1.zip
书上的描述：

[xls]

还有监测的一个点的密度变化