Re: [Help-gsl] solving ODE in a loop

[Frank Reininghaus]

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