[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Help-gsl] [ODE] GSL and Octave comparison
|
From: |
Brian Gough |
|
Subject: |
Re: [Help-gsl] [ODE] GSL and Octave comparison |
|
Date: |
Sat, 26 Aug 2006 10:56:11 +0100 |
|
User-agent: |
Wanderlust/2.14.0 (Africa) Emacs/21.3 Mule/5.0 (SAKAKI) |
At Thu, 24 Aug 2006 21:24:58 +0200,
Marco Maggi wrote:
> and then recomputing with GSL step function rk2
> control function 'y' eps_abs = eps_rel = 0.001,
Those are large tolerances -- they are applied to each step, so the
global error is O(Nsteps * tolerance). Also rk2 is the least accurate
method, try one of the others, e.g. rk4 at least.
--
Brian Gough
Network Theory Ltd,
Publishing the GSL Manual - http://www.network-theory.co.uk/gsl/manual/