Add a forward euler function that approximates ode45