divergence function

Returns the rate of divergence of close trajectories needed for the maximum Lyapunov exponent estimation.