Very recently, Atangana and Seda made use of the Newton polynomial interpolation to suggest new numerical scheme to solve linear and non-linear problems. In this paper, we applied such numerical method to solve the well-known telegraph equations with different differential operators. The stability analysis of each case is presented in detail and numerical simulations are presented for different values of fractional orders and fractal dimensions.