
[t2, y2] = eul('fun',0.0,1.0,1.0,0.05);
hold on;
plot(t2,y2,'g+');

