A simple, numerically accurate and efficient method is given for the computation of Green's functions and their spatial derivatives for SH motion in a layered elastic medium. It is shown that with appropriately refined analysis the usual approximate method for the computation of Love waves can be modified to yield the complete solution containing both Love and body waves for SH motion, The complete solution is expressed as a sum of modes of both Love and leaking mode type. Contrary to commonly held notions, no branch line integrals are required. The method applies for both near-and far-field solutions. The application of the method is illustrated with a numerical example.