The problem of tracing seismic rays between specified source and receiver is discussed for Earth models consisting of layers, in which velocity varies linearly, that are separated by material interfaces of arbitrary shape. The calculation of travel times, amplitudes, and phase shifts is considered. Fast and efficient numerical algorithms are developed. Computed examples are presented.