Quite some time ago, I imported colnew into maxima (with permission from
the authors). colnew is a solver for boundary value ODEs.
I never got around to implementing a maxima interface to the routines.
Until now. There is a very rudimentary and very raw interface to colnew
in place. There's no documentation except what is in the Fortran code.
However, there are two sample programs, prob1.mac and prob2.mac, that
demonstrate how to use the raw interface.
Hope this is still useful to someone. Sorry for taking so long.
I haven't figured out what a good simple interface would be, but perhaps
a user can help with suggestions on this.
Ray