Stability of pseudospectral and finite-difference methods for variable coefficient problems

Abstract

It is shown that pseudospectral approximation to a special class of variable coefficient one-dimensional wave equations is stable and convergent even though the wave speed changes sign within the domain. Computer experiments indicate similar results are valid for more general problems. Similarly, computer results indicate that the leapfrog finite-difference scheme is stable even though the wave speed changes sign within the domain. However, both schemes can be asymptotically unstable in time when a fixed spatial mesh is used.