Abstract: The half space boundary value problem in classical elasticity can be transformed into the half space harmonic equation problem with Neumann boundary conditions. Based on the Dirac \delta function in generalized functions and its related properties, a simple proof for the solution of this problem can be given along with its related engineering applications. This proof does not need the use of very complicated mathematical tools and can help the study of the related problems in elasticity.