I am using the lsode function in Octave to solve a set of ODEs.
According to the documentation, there is a function, lsode_options,
which can be used to set tolerances, etc. for lsode to use.
Unfortunately, I haven't been able to find any information on the
usage of lsode_options, for example, what options are available
and what are their default values.
I appreciate any guidance you may be able to provide.