We show constructively that, under certain regularity assumptions, any system of coupled linear differential equations with variable coefficients can be tridiagonalized by a time-dependent Lanczos-like method. The proof we present formally establishes the convergence so-called ?-Lanczos algorithm and yields full characterization algorithmic breakdowns. From there, solution original is available...