有没有黄平热弹流润滑的正确代码

书中代码运行得不到书中的结果，初学fortran,希望交流学习一下,最好能有MATLAB的程序，希望师兄师姐指点一二