
clf;
t = 0.0:0.01:1.0;
y = func(t);
[t1, y1] = eul('fun',0.0,1.0,1.0,0.1);
plot(t,y,'w-.',t1,y1,'c:');
xlabel('t');
ylabel('y');
title('Comparison: Exact Solution vs. Solutions Computed by Euler Method');
