有限差分法求解二维椭圆型方程的Matlab程序

有限差分法求解二维椭圆型方程的Matlab程序 % 用有限差分法求解矩形域上的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=h^2/l^2;     s=2*(1+r);     k=1;     % 应用Gauss-Seidel法求解差分方程