We present an explicit numerical scheme to solve the variable coefficient wave equation in one space dimension with minimal restrictions on the coefficient and initial data.
The algorithmA range of physical and biological applications involve the one dimensional wave equation 1where: there exist x − < x + such that ζ takes constant values ζ − and ζ + on the respective intervals (−∞, x − ) and (x + , ∞); and 0 < c < ζ < C for some c < C ∈ R. Applications include imaging of layered media such as seismic imaging and the acoustic imaging of laminated structures, microwave imaging of skin tissue, and modelling the human vocal tract and cochlea; see