这有个尾巴该如何解决啊

fcode

迭代过程如果误差积累了，试试用乘法代替加法。
比如：

[Fortran] 纯文本查看 复制代码

y = 0
zb=-(4.76_QP-4.9_QP)*0.14_QP            
Do i = 1 , 100
  y = y + zb
End Do

改写为：

[Fortran] 纯文本查看 复制代码

Do i = 1 , 100
  y = 0 + (i-1)*zb
End Do