numerical solutions to differential equations


> Thanks for your reply. Two more questions. Am I correct in assuming
> that maxima as of yet does not have a mechanism to return an
> interpolating function as the numerical solution to an ode?

I have just commited (ten minutes ago) some changes in the interpolation
package (don't use older versions: you can get overflow errors):

http://maxima.cvs.sourceforge.net/maxima/maxima/share/numeric/interpol.mac

Documentation is here:

http://maxima.cvs.sourceforge.net/maxima/maxima/doc/info/interpol.texi

Download interpol.mac and execute something like this:

load("your_path/interpol.mac")$
load("your_path/dynamics.mac")$
points: rk(t-x^2,x,1,[t,0,4,0.1]);
lagrange(points); /* Lagrange polynomial */
cspline(points);  /* Cubic splines interpolation */
linearinterpol(points); /* Linear interpolation */

In the examples in interpol.texi you'll see how to handle these three
latest results to make interpolations and plotting.

-- 
Mario Rodriguez Riotorto
www.biomates.net