I am playing with a coupled system of differential equations. This
system has an initial transient behavior that involves a lot of fast
dynamics, and then a long slow oscillation after the initial behavior.
Using rk required a tiny time step and didn't allow me to get to the
long slow oscillation efficiently.
I installed the latest maxima sources thinking that the rkf45 code that
was posted to the list would be included, but it doesn't seem to be as
"load(rkf45)" doesn't work. I then went directly to
https://sites.google.com/site/pjpapasot/maxima/libraries/rkf45 and
downloaded the code myself and used batch to load it, and it was very
helpful (THANK YOU Panagiotis !)
Will this code be included in the next Maxima release? It is very
helpful for getting numerical simulations of equations that are derived
via Maxima.