Is there a numerical ODE solver facility in Maxima



Using a computer algebra system there is another approach, which is to 
solve for the taylor series coefficients in f(x).
There are various packages to do this, though I don't know of one that 
is set up in Maxima.

A trivial setup (but expensive to run) is to define f, using deftaylor,  
with symbolic coefficients a0, a1, ...
and then expand f to some order, substitute in the DE, and solve for a0, 
a1, etc.

RJF