In this article, the multi-domain bivariate spectral collocation method is applied, for the first time, in solving hyperbolic partial differential equations (PDEs) that are defined over large time domains. The proposed ...