main.m

用Admas四阶预估-校正算法的PECE模式求下列初值问题的数值解

clear;clc;
[ex41,t41,u41,e41]=adams4(0,1,1,0.1);
subplot(1,2,1);
x=0:0.01:1;
plot(x,exat(x),'r-','linewidth',1.5);
hold on
plot(t41,u41,'b*');
hold off
legend('精确解','步长0.1数值解');
xlabel('自变量');
ylabel('函数值');
subplot(1,2,2);
plot(t41,e41,'b','linewidth',2.0);
legend('步长0.1时误差');