Jacobian Free Computation of Lyapunov Exponents

The purpose of this paper is to present new algorithms to approximate Lyapunov exponents of nonlinear differential equations, without using Jacobian matrices. We first derive first order methods for both continuous and discrete QR approaches, and then second order methods. Numerical testing is given, showing considerable savings with respect to existing implementations.