wave equations. By employing a fourth-order accurate scheme to approximate the solution
at first time level, it is shown that the discrete solution is conditionally convergent in the
maximum norm with the convergence order of two. Since the asymptotic expansion of the
difference solution consists of odd powers of the mesh parameters (time step and spacings),
an unusual Richardson extrapolation formula is needed in promoting the second-order …