[Help-gsl] ODE fail

[Richard Harke]

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