Hello,

I'm in the process of solving a single ODE using odeint.h and steppersie.h. The solution to the ODE is actually imbedded in a fitting algorithm such that the ODE is solved many times using different parameters for an input parameter set that defines the exact values of the ODE to solve. As part of the fitting algorithm, the ODE solver successfully finds a solution the majority of the time. However, there are a few cases when I get an out of bounds error. The error occurs at line 275 in steppersie.h; the call is

fsave[ipt][i]=del[i];

The value of ipt is found to be out of bounds.

I realize that without my specific code it may be hard to track down this problem, but I was wondering if anyone might have an answer to why this is happening before I go delving into the code to track down the error (whether due to my own code or a true bug in the NR code).

Sincerely,
Kyle

PS - To the authors, I find your book invaluable in my astronomy research. This book and Binney & Tremaine's Galactic Dynamics are always within arms reach. Thank you for all your work here and in the book itself.