Doha, E. H.Bhrawy, A. H.Hafez, R. M.Abdelkawy, M. A.2019-11-202019-11-202014Cited References in Web of Science Core Collection: 592325-0399https://doi.org/10.12785/amis/080211http://www.naturalspublishing.com/Article.asp?ArtcID=4681Accession Number: WOS:000331386900011A numerical approximation of the initial-boundary system of nonlinear hyperbolic equations based on spectral collocation method is presented in this article. A Chebyshev-Gauss-Radau collocation (C-GR-C) method in combination with the implicit Runge-Kutta scheme are employed to obtain highly accurate approximations to the mentioned problem. The collocation points are the Chebyshev interpolation nodes. This approach reduces this problem to solve system of nonlinear ordinary differential equations which are far easier to be solved. Indeed, by selecting a limited number of collocation nodes, we obtain an accurate results. The numerical examples demonstrate the accuracy, efficiency, and versatility of the method.enUniversity for System of nonlinear hyperbolic equationsCollocation methodChebyshev-Gauss-Radau quadratureImplicit Runge-Kutta methodSPECTRAL GALERKIN METHODBOUNDARY-VALUE-PROBLEMSFINITE-ELEMENT METHODSPARTIAL-DIFFERENTIAL-EQUATIONSINTEGRAL-EQUATIONSCOLLOCATION METHODINTEGRODIFFERENTIAL EQUATIONSPOLYNOMIAL SOLUTIONSVISCOELASTICITYDOMAINSArticlehttps://doi.org/10.12785/amis/080211