Abstract
In this article, a novel variable order fractional nonlinear Klein Gordon model is presented where the variable-order fractional derivative is defined in the Caputo sense. The merit of nonstandard numerical techniques is extended here and we present the weighted average nonstandard finite difference method to study numerically the proposed model. Special attention is paid to study the convergence and to the stability analysis of the numerical technique. Moreover, the truncation error is analyzed. Three test examples are provided. Comparative studies are done between the used numerical technique and the weighted average finite difference method. It is found that the stability regions are larger by using the weighted average nonstandard finite difference method.