有限差分法解一阶微分方程

% Copyright(C), 2024

%%% Test#12:    Compute the Poisson equation definite solution problem 
%%%             used by Finite-Difference-Method

% Author:       YANG Yk    #8211221221
% Class:        2201, Geophysics
% Advisor:      TONG Xz
% Date:         2024/04/05 05:18
%

% clear
clear;

% definition
a = 1;
n = 100;
dx = a/n;
x = 0:dx:a;
A = sparse(n+1,n+1);
b = sparse(n+1,1);

% boundary conditions
A(1,1) = 1;
A(n+1,n+1) = 1;
b(1) = 0;
b(n+1) = 0;

% operator
dx2 = dx^2;
k1 = 1/dx2;
k2 = -2*k1;
for i = 2:n
    A(i,i-1:i+1) = [k1 k2 k1];
    b(i) = exp(x(i));
end

% solute
u = A\b;

% figure
figure;
plot(x,u,'ro');grid on
xlabel('x');
ylabel('u');
legend('Solution of FDM')

方程是，边界上都是0.

效果如下