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: Raymond Toy
Date: Mon, 24 Oct 2011 22:15:28 -0700
On 10/24/11 9:27 PM, dlakelan wrote:
> On 10/24/2011 02:12 PM, Panagiotis Papasotiriou wrote:
>> Dear community,
>>
>> I have written a Maxima package for solving initial value problems. The
>> package implements a Maxima function called rkf45, which is an
>> implementation of the Runge-Kutta-Fehlberg method of 4th-5th order.
> ...
>
> Thank you very much for this contribution to Maxima, I hope one of the
> maxima developers will include it into the standard distribution as a
> shared package. However, rather than having a rkf45 specific share
> package, perhaps we should create a "numode" or other package named
> after the type of problem, where we can put several related numerical
> differential equation codes. For example perhaps the next task someone
> will take on will be some kind of implicit scheme for stiff systems, or
> someone else will produce boundary value problem solvers.
colnew will do boundary value problems. This is already included in maxima.
I have a translation via f2cl of odepack, which can do stiff odes. This
is not in maxima but could be.
>
> Organizing shared code by the task they solve rather than the method
> used helps the naive user to find the code.
>
It's unfortunate that ?? doesn't work with share packages. Then
appropriately named share packages could be found.
Ray