We adapt convergent look-ahead and backward finite difference formulas to compute future eigenvectors and eigenvalues of piecewise smooth time-varying matrix flows A(t). This is based on the Zhang Neural Network model for time-varying problems and uses the associated error function
E(t)=A(t)V(t)−V(t)D(t)
with the Zhang design stipulation
˙E(t)=−ηE(t).
Here E(t) decreased exponentially over time for η>0. It leads to a discrete-time differential equation of the form P(tk)˙z(tk)=q(tk) for the eigendata vector z(tk) of A(tk). Convergent high order look-ahead difference formulas then allow us to express z(tk+1) in terms of earlier discrete A and z data. Numerical tests, comparisons and open questions follow.