Matrix methods in synthetic seismograms

Abstract

Some methods are examined for the matrix problems arising in the computation of seismograms in a stratified medium. In particular, the Thomson‐Haskell and Kennett algorithms are placed in a common setting, and their connections with standard methods in numerical analysis are given. Finally, a new formulation is presented which has advantages both for physical insight and for numerical accuracy. Copyright © 1984, Wiley Blackwell. All rights reserved