On 9/4/10, Tuomo Keskitalo <[email protected]> wrote: > Thanks for the link. That is a special code meant for second order ODEs. > Apparently your case requires precise solution, and therefore you should > modify the Newton iteration code if you want to use GSL for your problem. > However, I think your case is special. At least I do not know of a generally > accepted all-purpose tolerance for the Newton iteration.
I just want to use, e.g. the implicit RK method of order 2 (midpoint rule) and rely on the proved properties of the stepper function (symmetric, symplectic). Either the method shouldn't be named "rk2imp". Hope, this opinion is clear enough. _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
