Hi,
On Tuesday 09 September 2008 19:04:34 Lado Samushia wrote:
> I am trying to solve am ode system in a double loop (for different values
> of parameters).
> ...
> It works fine for the first time, but brakes down on second (and
> consequent) loop executions (gives "nan" for function values and jumps
> trough interval in one step).
I think you should call
gsl_odeiv_step_reset (s);
gsl_odeiv_evolve_reset (e);
before each loop execution, see
http://www.gnu.org/software/gsl/manual/html_node/Stepping-Functions.html
http://www.gnu.org/software/gsl/manual/html_node/Evolution.html
I hope that helps!
Regards
Frank
