Higher‐order algorithms for stable solutions of fractional time‐dependent nonlinear telegraph equations in space.

In this article, three novel algorithms are developed and successfully applied to investigate the stable solutions of time‐fractional nonlinear telegraph equations. Firstly, we presented the shifted Gegenbauer polynomials through appropriate transformations. The approximation of a function u(||x||, t) is defined via shifted Gegenbauer polynomials (SGPs) and then developed the operational matrices of a positive integer and non‐integer derivatives assisted by multi‐dimensional shifted Gegenbauer polynomials vectors. Linearized fully spectral method (LFSM) is proposed by means of function approximations and developed operational matrices whilst Picard iterative scheme is used to linearize the discussed problem. To introduce a semi‐discrete method (SDM), the temporal derivative is approximated through a central difference scheme whereas spatial derivatives are approximated assisted by novel operational matrices. A fully numerical method (FNM) is proposed by introducing the central difference for temporal and spatial variables. Proposed computational methods transform the time‐fractional nonlinear telegraph equation into a system of linear algebraic equations that are easier to tackle. An inclusive comparative study is witnesses that suggested schemes are effective, accurate and well‐matched to investigate the approximate solutions of the problems udnerstudy. Convergence and error bound of the suggested methods are investigated theoretically whereas stability analysis is performed numerically. Moreover, the developed methods can be used conveniently to examine the numerical solution of other multi‐dimensional highly nonlinear fractional or variable order problems of physical nature. [ABSTRACT FROM AUTHOR]