Static Green's functions in multilayered half-spaces (Q1374171)

We present a method for calculation of three-dimensional Green's functions in multilayered transversely isotropic or isotropic half-spaces. The formulation is based on the cylindrical system of vector functions combined with the propagator matrix method. The Green's solutions are expressed in terms of infinite integrals involving Bessel functions and are numerically integrated by an adaptive Gauss quadrature with continuous fraction expansions. Numerical examples show that the method is both accurate and efficient.