Abstract
The mathematical description of various processes such as the nonlinear Klein-Gordon equation occurring in mathematical physics leads to a nonlinear partial differential equation. The mathematical model is only the first step, however, towards finding the solution of the problem under consideration. It has become possible to develop realistic mathematical models with the currently available computing power for complicated problems in science and engineering. To the best of our knowledge, systematically using the collocation method to acquire the numerical solution has not been previously used for the Klein-Gordon equation. The main aim of this paper is to systematically use the collocation method to acquire the numerical solution of the two coupled nonlinear non-homogeneous Klein-Gordon partial differential equations. We examine and analyze their stability, in detail. To this aim, we use the Von Neumann stability method to show that the proposed method is conditionally stable. A numerical example is introduced to demonstrate the performance and the efficiency of the proposed method for solving the coupled nonlinear non-homogeneous Klein-Gordon partial differential equations. The numerical results demonstrated that the proposed algorithm is efficient, accurate, and compares favorably with the analytical solutions.