ode.py now contains the euler-richardson function