The leading Lyapunov vector is computed as follows:
1) Start with an arbitrary perturbation of arbitrary size
2) Evolve it from  t  to               using the TLM
3) Repeat 2) for the succeeding time intervals
After a sufficiently long time , the perturbation converges to the leading Lyapunov
vector. The direction of this vector is independent of the initial perturbation, the
length of the time interval or the choice of norm of the perturbation,
properties not shared by singular vectors. If during the repeated application of
the TLM the LV becomes too large it may be scaled down to avoid computational
blow up.
Additional LVs can be obtained by the same procedure, except that after
each time step the perturbation has to be orthogonalized with respect to
the subspace of the previous LVs, since otherwise all the LVs would
converge to the leading LV.