function ex0801()
[x1,y1]=MyEuler('myfun01',0,2*pi,1,16);%歐拉法所的的解
h1=2*pi/15 %計算取16點的步長
[x11,y11]=MyEuler('myfun01',0,2*pi,1,32);%歐拉法所的的解
h2=pi/15 %計算取32點的步長
y=dsolve('Dy=y+sin(t)','y(0)=1');%該常微分方程的解析解
for k=1:33
t(k)=x11(k);
y2(k)=subs(y,t(k));%求其對應點的離散解
end
plot(x1,y1,'+b',x11,y11,'og',x11,y2,'*r')
legend('h=0.4的歐拉法解','h=0.2的歐拉解','符號解')