In this paper we develop a finite-difference scheme to approximate radially symmetric solutions of the initial-value problem with smooth initial conditions (1)in an open sphere D around the origin, where the internal and external damping coefficients-β and γ, respectively-are constant, and the nonlinear term has the form G ′ (w) = w p , with p > 1 an odd number. The functions φ and ψ are radially symmetric in D, and φ, ψ, rφ and rψ are assumed to be small at infinity. We prove that our scheme is consistent order O(∆t 2 ) + O(∆r 2 ) for G ′ identically equal to zero, and provide a necessary condition for it to be stable order n. Part of our study will be devoted to compare the physical effects of β and γ.