本程序利用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
5星
