Thanks for your reply. I have found it. My mistake, hstart was way too large. Now I'm only puzzled that the failure mode was rather obscure.
On Sat, Sep 19, 2015 at 3:55 AM, Juan Pablo Amorocho D. <[email protected]> wrote: > Richard- > > I agree with Matthias. You'll need to provide us with your source code, > and the GSL version, if you want us to help you. > - Juan > > > > 2015-09-18 19:38 GMT+02:00 Richard Harke <[email protected]>: > >> I have a set of ODE's I'm trying to solve. After I installed gsl, I built >> and ran >> the example, vanderpol.c. No problem. I then changed func and jac for my >> problem, changed the number of eqs. to 4. The output was simply the >> initial values repeated over and over. I put beak points in func and jac. >> They are never called!? Is there possibly some simple explanation or am >> I doomed to slog through the library code to find this? >> >> The eqs are: >> dx/dt =xdot >> dy/dt = ydot >> d(xdot)/dt = a function of x, y,xdot, ydot >> d(ydot)/dt = a function of x, y, xdot, ydot >> >> Thanks for any suggestion >> > >
