Issue #2, although seemingly fixed revealed that it is not a good idea to hard code the parameters of the strategy that decides whether or not the orbit has closed. We have to pass it as a user options.
For starters, perhaps a parameter tolerance_coefficient, with default value ==1000
so that closing_orbit_tolerance = tolerance_coeff * ode_coeff will do the work for us.
Issue #2, although seemingly fixed revealed that it is not a good idea to hard code the parameters of the strategy that decides whether or not the orbit has closed. We have to pass it as a user options.
For starters, perhaps a parameter tolerance_coefficient, with default value ==1000
so that closing_orbit_tolerance = tolerance_coeff * ode_coeff will do the work for us.