拉普拉斯方程五点差分格式算法

五点差分格式解二维拉普拉斯方程。程序思路更清晰，已经调试通过。

%****二维椭圆方程的差分解法****%

%二维椭圆方程为：uxx+uyy=0

% x方向的结点数 nx
% x方向的步长 hx
% 求解区间x的左端：minx
% 求解区间x的右端：maxx
% y方向的结点数 ny
% y方向的步长 hy
% 求解区间y的左端：miny
% 求解区间y的右端：maxy
% 求解区间上的数值解 u

%%%处理求解区域上边界的函数%%%
clear;
minx=100;
maxx=400;
miny=0;
maxy=300;
nx=8;
ny=8;
hx=（maxx-minx）/（nx-1）;
hy=（maxy-miny）/（ny-1）;
u0=zeros（nxny）; %u0[nxny]二维数组存放边界条件的值
x=minx:hx:maxx; %差分网格x向坐标
y=miny:hy:maxy; %差分网格y向坐标
Z=zeros（nxny）; %Z[nxny]二维数组存放所有差分网格上的数值解（包括已知边界网格点上的值）
[XY]=meshgrid（xy）;
for i=1:nx
    %u0（1i）=UpBorder（maxyminx+（i-1）*hx）; %处理求解区域的上边界
    %u0（nyi）=DownBorder（minyminx+（i-1）*hx）; %处理求解区域的下边界
    u0（1i）=8+（i-1）*0.1;
    u0（nyi）=6.5+（i-1）*0.1;
end
for j=1:ny
    %u0（j1）=LeftBorder（maxyminx+（i-1）*hx）; %处理求解区域的左边界
    %u0（jnx）=RightBorder（）; %处理求解区域的右边界
    u0（j

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3780  2011-08-02 22:09  difference.m

----------- ---------  ---------- -----  ----

                 3780                    1