##
## This computes the time pi (period of the potential term) map for the
## Hill's differential equation (see hill.m)
##
function ypi = hillstrob(y)
  ypi = lsode("hill", y, [0, 2 * pi]);
endfunction

## ;;; Local Variables: ***
## ;;; mode: octave ***
## ;;; End: ***
