pub fn vectors<A, TEO>(
    teo: TEO,
    x: Array1<A>,
    alpha: A::Real,
    duration: usize
) -> Vec<(Array1<A>, Array2<A>, Array1<A::Real>)>where
    A: Scalar + Lapack,
    TEO: TimeEvolution<Scalar = A, Dim = Ix1> + Clone,Expand description
Calculate the Covariant Lyapunov Vectors at once
This function saves the time series of QR-decomposition, and consumes many memories.