偏微分方程离散计算不收敛求助

vvt

太长了，实在没精力看。
我觉得应该是楼主的代码写得不好，而不是算法问题。

你确定 dx2 就只算最后一个？而不是累加？

vvt

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

do i=M,N+M
do j=-N,N
  dx2=(u1(i,j)-u0(i,j))**2.d0+(v1(i,j)-v0(i,j))**2.d0+(zeta1(i,j)-zeta0(i,j))**2.d0
  end do
end do

这就是只算了最后一个。不是累加。累加应该是 dx2 = dx2 + ...