
[t4, y4] = eul('fun',0.0,1.0,1.0,0.00001);
hold on;
plot(t4,y4,'c:');

