A Maxima function for solving initial value problems with adaptive step size and error control.
Subject: A Maxima function for solving initial value problems with adaptive step size and error control.
From: Jaime Villate
Date: Wed, 26 Oct 2011 09:18:54 +0100
On 10/25/2011 06:09 PM, Panagiotis Papasotiriou wrote:
> Actually, rkf45 is not based on rk; I tried to keep it compatible with
> rk's syntax to avoid confusion for new users, and of course both rk
> and rkf45 are implementing Runge-Kutta methods, but otherwise their
> implementation has little in common.
> If rkf45 is included in share, rk could still be useful for
> educational purposes (teaching students about the original Runge-Kutta
> method of fourth order, and using rk for examples of use.) Otherwise,
> rkf45 should be used for real work. I believe we should keep them
> separated, as it is the case in most Numerical Analysis packages I am
> aware of.
OK, I do not have any objections then.
I should also point out at that in future releases rk and plotdf will be
more related to each other, once I figure out an efficient way to pass
data among themselves. plotdf already uses an adaptive step-size and it
can also use Adams-Moulton instead of Runge-Kutta (even though it is
currently disabled due to some upstanding bugs).
Regards,
Jaime