function [time,y]=euler(time,y,dt);
y=y+dt*f(time,y);
time=time+dt;
