有限差分法的Matlab程序(椭圆型方程)3页.doc

有限差分法的Matlab程序（椭圆型方程）function FD_PDE(fun,gun,a,b,c,d) % 用有限差分法求解矩形域上的Poisson方程 tol=10(-6); % 误差界 N=1000; % 最大迭代次数 n=20; % x轴方向的网格数 m=20; % y轴方向的网格数 h=(b-a)/n; % x轴方向的步长 l=(d-c)/m; % y轴方向的步长 for i=1:n-1 x(i)=a+i*h; end % 定义网格点坐标 for j=1:m-1 y(j)=c+j*l; end % 定义网格点坐标 u=zeros(n-1,m-1); %对u赋初值 % 下面定义几个参数 r=h2/l2; s=2*(1+r); k=1; % 应用Gauss-Seidel法求解差分方程 while k=N % 对靠近上边界的网格点进行处理 % 对左