ode: numeric and power series



Sam Steingold wrote:
> how do I solve an 1st order ODE IVP on maxima numerically?

load(dynamics);

There's an rk function for solving odes.  I've never used it.  "?? rk"
gives some info.  Looking in dynamics.mac gives an example:


results: rk(t-x^2,x,1,[t,0,8,0.1])$
plot2d([discrete, results]));

Don't know if that's right or not.

> how do I solve an 1st order ODE IVP on maxima using power series?

Don't know about that one.

Ray