“…Among these equations, the hyperbolic telegraph equation is an important type of hyperbolic partial differential equation commonly used in modeling of mathematical problems such as signal analysis, random walking theory and wave propagation [1]. Many methods have been applied to solve numerically the second order linear hyperbolic telegraph equation, such as the implicit three-level difference scheme [2], the collocation methods based on the thin plate splines, the quartic B-spline, the septic B-spline, the cubic B-spline, the modified cubic B-spline, the extended cubic B-spline and the cubic trigonometric B-spline [3][4][5][6][7][8][9][10][11], the Rothe-Wavelet method [12], the Legendre multiwavelet Galerkin method [13], the Chebyshev tau method [14], the finite difference methods based on the quartic spline functions, the non-polynomial splines and the cubic Hermite interpolation polynomial function [15][16][17], the unconditionally stable parallel difference scheme [18], the boundary integral equation method and the dual reciprocity technique [19], the differential quadrature method [20,21], the Rothe-wavelet-Galerkin method [22], the reduced differential transform method [23], the reproducing Kernel Hilbert space method [24], the sinc-collocation method [25], the Chebyshev wavelets method [26], the meshfree method based on the radial basis functions [1], the polynomial scaling functions method [27], the high-order shifted Gegenbauer pseudospectral method [28], the uniformly convergent Euler matrix method [29], the Euler method [30], the Bessel collocation method [31], the Galerkin method [32] and the adaptive Monte Carlo method …”