Hello,
I am trying to solve a system of stiff ODEs using the GSL. The part of
my code that solve the equations looks like:
for (i = 0; i < time_steps; i++)
{
while (t < tim[i])
gsl_odeiv_evolve_apply (evolve, control, step, &syst,
&t, tim[i], &h, y);
}
I would like to check what is the error for y[] reached at the end of each
gsl_odeiv_evolve_apply() call.
How can I do that? An example would be appreciated.
Sébastien
_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl
