Abstract
A numerical analysis of the well-known linear partial differential equation describing the relativistic wave is presented in this work. Three different operators of fractional differentiation with power law, exponential decay law and Mittag-Leffler law are employed to extend the Klein-Gordon equation with mass parameter to the concept of fractional differentiation. The three models are solved numerically. The stability and the convergence of the numerical schemes are investigated in detail.