“…In this study, we consider a numerical solution of the general nonlinear wave problem αu tt + βu t − γ∆u t − λ∆u + ρ∇u = g(u, ∇u, ∆u) + f , x ∈ Ω, t ≥ 0 u(x, 0) = g 1 (x), u t (x, 0) = g 2 (x), u(x, t) = h(x, t) for x ∈ Γ (1) in which α, β, γ, and ρ are constants; Γ is the boundary of the problem domain Ω = [0, T] 2 ; u is the unknown function of time and space; g can be an arbitrary continuous function; f is the source term; g 1 (x), g 2 (x), and h(x, t) describe the initial and boundary conditions. The general form (1) covers several types of well-known wave problems arising in many scientific and engineering fields, such as nonlinear optics [1][2][3], solid state physics [4][5][6][7][8], and quantum field theory [9][10][11]. For instance, when α = 1, γ = 0, β= 0, λ= 1, ρ = 0, and g(u) = au + bu c with constants a, b, and c, Equation (1) degrades to the so-called nonlinear Klein-Gordon equation:…”