[SciPy-User] Maximum Lyapunov exponent of my previous system

[Barrett B]

I need some help setting up the calculation of the maximum Lyapunov
exponent of the system I was describing in my previous thread, "Heaviside
function in vector form." Note that I list the actual function being
integrated, f(X, t), in my first follow-up post on that thread.

Here is a brief discussion of the Lyapunov exponent:
https://en.wikipedia.org/wiki/Lyapunov_exponent

A simple Google search may come up with a better description.

I tried simply introducing a fourth differential variable and starting to
track it after some arbitrary time t_crit, but because of the nature of the
integrator, that doesn't work. So I need help trying to calculate the MLE.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20140709/beeb2d36/attachment.html>