On 03/25/2010 10:25 AM, dabu wrote:
Hi,
I am new in sage. I was wondering about Sage's capability to solve
odes numerically.
I was expecting to find something which is like ndsolve of
Mathematica.
For example it should not only as for the first order equations, nor
that one has to supply jacobians manually. It should also tackle
boundary conditions as equations.
I wonder any such things exist. Without such a construct. Sage would
not be useful (read compete with Mathematica :)) to a large part of
theoretical physics community.
I would be happy to write such a thing in case it does not exist.
I've always just used the scipy functions to do this sort of thing,
though I'm not an expert in the area, so I'm not sure what the scipy
functionality is missing. See:
http://www.sagenb.org/home/pub/258/ - basic example
http://www.sagenb.org/home/pub/106/ - calculating streamlines in Sage
http://www.sagenb.org/home/pub/879/ - epidemic modeling
http://www.sagenb.org/home/pub/42/ - spring-mass systems
Docs:
http://docs.scipy.org/doc/scipy/reference/tutorial/integrate.html#ordinary-differential-equations-odeint
http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.odeint.html#scipy.integrate.odeint
http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html#scipy.integrate.ode
It would be great to have something like this (or wrap the scipy stuff)
in Sage, so if you want to contribute, go for it! I would be glad to
see it.
Thanks,
Jason
--
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
To unsubscribe from this group, send email to sage-support+unsubscribegooglegroups.com or
reply to this email with the words "REMOVE ME" as the subject.
