今天是个好日子。
由于代码太大，我附上repl的链接：
https://replit.com/@PabloJunior1/Explicit-scheme-ENG
我正在尝试为二维方程实现显式差分方案。我处理了所有公式并且程序运行良好并且在 256 次迭代后收敛，如果我使用时间步长 τ = 0.25 * h ^ 2。但是如果我使用步长 τ = h ^ 2，那么在大约第 370 次迭代之后这个错误出去，这个过程进入一个永恒的循环。不幸的是，这是需要使用的值。.
我会注意到，在当前参数 N1 = N2 = 10 和误差精度 eps = 10 ^ (- 5) 的情况下，步长 τ = h ^ 2 程序应该在大约 75-85 步内收敛。我用隐式方案和交替方向隐式方案检查了它。如有必要，我可以附上原始公式。
此外，我已经尝试将 float64 类型分配给所有数组，正如在类似问题中所建议的那样。但是我的代码仍然不起作用。我不知道如何解决它。