Fortran修正了原代码的错误，与原代码的数值差异不大。

Matlab代码运行正常，两者趋势一致，但存在精度差异。