Matrix methods in synthetic seismograms

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.