In this article, a numerical method based on the shifted Chebyshev functions for the numerical approximation of the coupled nonlinear variable-order fractional sine-Gordon equations is shown. The variable-order fractional derivative is considered in the sense of Caputo-Prabhakar. To solve the problem, first, we obtain the operational matrix of the Caputo-Prabhakar fractional derivative of shifted Chebyshev polynomials. Then, this matrix and collocation method are used to reduce the solution of the nonlinear coupled variable-order fractional sine-Gordon equations to a system of algebraic equations which is technically simpler for handling. Convergence and error analysis are examined. Finally, some examples are given to test the proposed numerical method to illustrate the accuracy and efficiency of the proposed method.
On the sub–diffusion fractional initial value problem with time variable order
