fff.m

关于一的微分方程数值解,一维波动方程初边值问题Crank-Nicholson格式,及追赶法

clear;clc;
i=1;
for T=0.1:0.1:1              %产生时间等于0.1，0.2 ,0.3, ...,0.9, 1的结果。 
    A=fun1(0.00005,0.1,T);
    A1(i,:)=A(end,:);
    A=fun2(0.00005,0.1,T);
    A2(i,:)=A(end,:);
    A=fun3(0.00005,0.1,T);
    A3(i,:)=A(end,:);
    A=fun4(0.00005,0.1,T);
    A4(i,:)=A(end,:);
    i=i+1;
end