fun1.m

来自「关于一的微分方程数值解,一维波动方程初边值问题Crank-Nicholson格式」· M 代码 · 共 18 行

function A=fun1(t,h,T)
%上机实习一显式格式;
%T是时间，t迭代的时间步长。h空间步长。
r=t^2/h^2;
D=100;E=5000;F=3000;L=1;
N=T/t+2;
J=L/h+2;
A(2,J)=0;                         %另第一,二行为零
for i=3:N
    A(3,1)=0;
    A(3,end)=0;
    for k=2:J-1
        A(3,k)=r*D*A(3-1,k+1)+(2-2*r*D)*A(3-1,k)+r*D*A(3-1,k-1)-A(3-2,k)+t^2*(E*(k-1)*h*(L-(k-1)*h)+F);
    end
   A(1,:)=A(2,:);
   A(2,:)=A(3,:);
end