
[t1, y1] = eul('fun',0.0,1.0,1.0,0.1);
[t2, y2] = eul('fun',0.0,1.0,1.0,0.05);
[t3, y3] = eul('fun',0.0,1.0,1.0,0.025);
[t4, y4] = eul('fun',0.0,1.0,1.0,0.001);
t = 0.0:0.01:1.0;
y = func(t);
plot(t,y,'w-.',t1,y1,'m:',t2,y2,'g+',t3,y3,'yo',t4,y4,'c:');
xlabel('t');
ylabel('y');
title('Comparison: Exact Solution vs. Solutions Computed by Euler Method');
