lua-users home
lua-l archive

[Date Prev][Date Next][Thread Prev][Thread Next] [Date Index] [Thread Index]

On 16 February 2011 08:47, Francesco Abbate <> wrote:
> I'm working now to implement the system to control the error and
> adjust the step of automatically. This is something required in real
> cases where a fixed step size would be a too much coarse approach.

If you are serious about writing RK (explicit) ODE integrator with
auto time step control, then the Dormand–Prince method with embedded
RK 4/5 pair is the one to go for:–Prince_method

See also:

"A survey of the MATLAB ODE suite"

Double time step error control is not exactly the state-of-the-art any more.