A vector–matrix differential equation is formulated using normal mode analysis from the governing equations of a three-dimensional anisotropic half space in presence of heat source and gravity. The corresponding solution is obtained with the help of eigenvalue approach. Numerical computations for displacement, thermal strain and stress component, temperature distribution are evaluated and presented graphically.