On Wed, Jan 24, 2007 at 05:09:16PM -0500, Dan Solomon wrote:
...HTML email Elided...
Your email was a bit difficult to read from a non HTML email reader,
but I gather that the problem seems to be something about how you're
defining your expressions for the components since heavily nested
lists are being formed .
rk expects as its first argument, a list of expressions. Each
expression is the derivative of one of the unknown dependent
variables.
here's an example:
solve the system dx/dt = .1*x + .1*y, dy/dt = -.1*y^2, x0=1,y0=1, for
x and y when t in the interval 0,1 with step size .01
rk([.1*x+.1*y,-.1*y^2],[x,y],[1,1],[t,0,1,.01]);
Now perhaps you can convert your equations into a similar form given
the template?
Hope that helps.
--
Daniel Lakeland
dlakelan at street-artists.org
http://www.street-artists.org/~dlakelan