function [time,y]=euler(time,y,dt);
% implements a standard forward Euler step
% f is the rhs of the equations
y=y+dt*f(time,y);
time=time+dt;
