二维拉普拉斯方程的数值解法

        拉普拉斯方程（Laplace's equation）又称调和方程、位势方程，是一种偏微分方程，因由法国数学家拉普拉斯首先提出而得名。

        两个自变量的拉普拉斯方程具有以下形式：

        对于这样一个偏微分方程，我们想要通过计算机求它的数值解：对于一个函数f(x)在x处的一阶导数我们可以用它的中心差分来表示：f'(x)=(f(x+dx/2)-f(x-dx/2))/dx.令dx=1可得f(x)步长为1的一阶导数f'(x)=f(x+0.5)-f(x-0.5).同理f(x)的二阶导数f"(x)可以看作f(x+0.5)-f(x-0.5)的一阶导数，再使用一次中心差分可得f"(x)=f(x+1)-2f(x)+f(x-1)。

        对于二元函数F(x,y),F(x,y)对自变量x的二阶偏导为F(x+1,y)-2F(x,y)+F(x-1,y)，对自变量y的二阶偏导为F(x,y+1)-2F(x,y)+F(x,y-1)。对拉普拉斯方程可以得到F(x+1,y)-2F(x,y)+F(x-1,y)+F(x,y+1)-2F(x,y)+F(x,y-1)=0，即F(x,y)=(F(x+1,y)+F(x-1,y)+F(x,y+1)+F(x,y-1))/4。在已知边界条件和初始条件的情况下可以迭代求和得到满足精度要求的数值解F(x,y)。