我正在尝试求解二维拉普拉斯方程，

$\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2} = 0; \qquad 0 \lt x \lt 1, \quad0 \lt y \lt 1$

受限于边界条件，

${\frac{d T}{d x}}|_{x=0,y}=0;$

${\frac{d T}{d y}}|_{x,y=0}=0;$

${\frac{d T}{d y}}|_{x,y=1}=0;$

$T|_{x=1,y} = 1;$

我正在使用有限差分方法来离散前两个 B.Cs 中的微分方程和前向差异以及第三个 BC 的后向差异的解在 Gauss Seidal 和 Jacobi 中完美出现迭代方法，但使用 BicgStab 算法，即使与 Jacobi 或 Gauss Seidal 矩阵相同，它也不会收敛。任何想法都可能有所帮助。谢谢你。$T=1$