Abstract: The algorithms behind a toolbox for approximating Lyapunov exponents of nonlinear differential systems by QR methods are described. The basic solvers perform integration of the trajectory and ...