In this work, we address the problem of solving nonlinear general Klein–Gordon equations (nlKGEs). Different fourth- and sixth-order, stable explicit and implicit, finite difference schemes are derived. These new methods can be considered to approximate all type of Klein–Gordon equations (KGEs) including phi-four, forms I, II, and III, sine-Gordon, Liouville, damped Klein–Gordon equations, and many others. These KGEs have a great importance in engineering and theoretical physics.
The higher-order methods proposed in this study allow a reduction in the number of nodes, which might also be very interesting when solving multi-dimensional KGEs. We have studied the stability and consistency of the proposed schemes when considering certain smoothness conditions of the solutions. Additionally, both the typical Dirichlet and some nonlocal integral boundary conditions have been studied. Finally, some numerical results are provided to support the theoretical aspects previously considered.