[Help-gsl] ODE solver questions

[Thanate Dhirasakdanon]

Hello, 

I would like to know if there is any way to get
an interpolation of solution between each integration step.
In LSODE, we can do that with a call to DINTDY. 

I guess I am looking for something like the function
gsl_odeiv_evolve_intp() below:

int status = gsl_odeiv_evolve_apply (e, c, s,
                                     &sys, 
                                     &t, t1,
                                     &h, y);
double tout = (t+t1)/2.0;
double yout[2];
status = gsl_odeiv_evolve_intp (e, tout, yout);
// gsl_odeiv_evolve_intp() returns error if tout<t or tout>t1.

I am new to gsl, so any suggestion is welcomed.

Thanate


__________________________________________________
Do You Yahoo!?
Tired of spam?  Yahoo! Mail has the best spam protection around 
http://mail.yahoo.com