有限差分法的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轴方向的步长 fori=1:n-1 x(i)=a+i*h;