[Maxima] ode: numeric and power series
raymond.toy at ericsson.com
Wed Sep 17 16:03:56 CDT 2008
Sam Steingold wrote:
> Raymond Toy wrote:
>> Sam Steingold wrote:
>>> how do I solve an 1st order ODE IVP on maxima numerically?
>> 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.
> runge-kutta sounds good but it has been running on my amd64 ws for 20+
> minutes now... (maxima/sbcl):
I've never actually used rk before, but I think you should have written
This runs pretty fast and seems close to the Lisp code you gave.
More information about the Maxima