%ex0806.m 用ode113(),MyHamming(),解微分方程y'=-y+1
clear, clf
x0 = 0; %積分區(qū)域開(kāi)始值
xt = 10; %積分區(qū)域終止值
y0 = 0; %y初值
N = 50; %離散點(diǎn)數(shù)
fun66 = inline('-y+1','x','y'); %需要求解的微分方程
f66 = inline('1-exp(-t)','t'); %方程的精確解析解
tic
[x113,y113] = ode113(fun66,[x0,xt],y0); 
t113 = toc
tic
[x1,yH] = MyHamming(fun66,x0,xt,y0,N); 
tH = toc
yt1 = f66(x1); %在精確解上取離散點(diǎn)
subplot(1,2,1) %將窗口分為一行兩列，在第一個(gè)圖中作出各種方法的解的曲線
plot(x1,yt1,'*r',x1,yH,'og')
legend('精確解','Hamming解')
title('Hamming法所得的解')
yt113=f66(x113);
subplot(122 ) 
plot(x113,yt113,'*r', x113, y113,'og')
legend('精確解','ABM解')
title('ode113即ABM算法所得的解')