
[t3, y3] = eul('fun',0.0,1.0,1.0,0.025);
hold on;
plot(t3,y3,'yo');

