Write programs with AMPL to solve odes