A method for computing the deformation in a layered, self-gravitating half space due to internal point sources is formulated. The coupled elastic-gravitational field equations appropriate to the first-order perturbation problem are presented. Although the first-order perturbation problem is not strictly coupled to a physically reasonable zeroth-order problem, the perturbation equations together with boundary conditions present a well-posed problem. The general homogeneous solution for deformation in a uniform infinite medium is found in the form of eight independent solution vectors. It is easily shown that the solution to the coupled elastic gravitational equations reduces to that for the uncoupled equations in the limit of vanishing density. The solutions for various internal point sources in a layered half space are constructed out of the solution found for the infinite medium by the user of well-known techniques. The half space is divided into a series of horizontal layers, and a matrix method is used to propagate the solution from one layer to the next. The point source is represented by a discontinuity in several of the various propagated quantities across a plane z=constant. The solution is ultimately determined by the requirements of vanishing normal tractions on the free surface and vanishing of all perturbation quantities at infinite depth. Finally, expressions are given for the direct calculation of surface displacements from several types of point sources.