A precise integration algorithm is proposed to solve the differential equat
ion of surface wave propagation for stratified anisotropic material. The me
thod can easily be extended to 3D wave problems. The numerical result is ex
act in the sense that it depends only on the precision of the host computer
.