Numerical simulations of time-dependent behaviour of advances structures need the analysis of systems of partial differential equations of hyperbolic type, whose semi-discretization, using the Fourier multiplicative decomposition together with the finite element or similar techniques, leads to large sparse systems of ordinary differential equations. Effective and robust methods for numerical evaluation of their solutions in particular time steps are required; thus still new computational schemes...