a b s t r a c tThe dynamic system of anisotropic elasticity from three second order partial differential equations is written in the form of the time-dependent first order symmetric hyperbolic system with respect to displacement velocity and stress components. A new method of deriving the time-dependent fundamental solution of the obtained system is suggested in this paper. This method consists of the following. The Fourier transform image of the fundamental solution with respect to a space variable is presented as a power series expansion relative to the Fourier parameters. Then explicit formulae for the coefficients of these power series are derived successively. Using these formulae the computer calculation of fundamental solution components (displacement velocity and stress components arising from pulse point forces) has been made for general anisotropic media (orthorhombic and monoclinic) and the simulation of elastic waves has been obtained. These computational examples confirm the robustness of the suggested method.