Shan Chen R-P 验证

各位大佬好;
请问:
我想验证R-P 方程,但是只有100x100勉强可以通过,400x400 就差的有点多,尤其是增长的情况。另外psi是用哪一个?我用的代码里的第一个(也是很多文献中的公式)发散了,只有第二个可以出结果。Guo的力模型在c-s eos时也不适用,不知道是不是我代码错了。
代码如下:
https://github.com/xiongxin9000/shan-chen/blob/master/s-c.cpp

谢谢