On Aug 30, 2008, at 10:41 PM, Jason Merrill wrote:
* Please don't make me state the problem in terms of first order equations only. Turning this second order equation into a coupled set of first order equations is a simple (but boring) algebraic manipulation, and Sage should be able to do it. In cases where Sage can't pull off the necessary manipulation, fail with a descriptive error message.
I have MATLAB code to numerically solve a system of second order linear ODEs directly. It runs much faster than MATLABs ODE solvers on my finite element problems. Once I get the symbolic side coded, it should be straightforward to convert to SciPy. That said, if someone wants it, I'd be happy to supply it under GPL if someone else wants to convert it to Python. It's a particular implementation of the Newmark-Beta method that doesn't require iteration and instead reduces it to matrix manipulation. --- Tim Lahey PhD Candidate, Systems Design Engineering University of Waterloo
smime.p7s
Description: S/MIME cryptographic signature