The additional pressure at the bottom of a building’s foundation produces an additional stress in the foundation soils under the building’s foundation. In order to overcome the limitations of traditional elastic theory methods and the finite element method when calculating the additional stress in foundation soils, we use the Gauss-Simpson formula to derive the Gauss-Simpson Quadrature Algorithm based on the elasticity theory. The Gauss-Simpson Quadrature Algorithm is a method designed to calculate the additional stress in foundation soils under an irregularly shaped foundation and an irregular load distribution. This new method is based on the fact that the Gaussian quadrature formula and the Simpson formula are independent of the specific type of integrand. The finite element method with n interpolation points can only achieve an algebraic accuracy of n. The interpolation points of the Gaussian quadrature formula are n zeros of orthogonal polynomials, which can achieve an algebraic accuracy of 2n + 1. Moreover, the weights of the nodes in the quadrature formula are all positive, and thus, it has a high numerical stability. In the proposed method, the Simpson formula is necessary. The Simpson formula is used to transform the implicit additional stress formula with the integral sign into an explicit cumulative integral, which can be considered similar to the rectangular domain case to obtain an explicit analytical algebraic formula for solving the additional stress approximation. In engineering applications, we only need to provide the field engineers with the locations of the interpolation points of the Gauss-Legendre formula, the interpolated weight coefficients, and the specific type of Simpson's formula, and then, the results of the additional stress can be calculated manually, which is nearly impossible using the traditional methods and finite element methods. From the point of view of academic rigor and theoretical completeness, it is possible to use the compound Gauss-Simpson Quadrature Algorithm in conjunction with the looping function in computer programs. Under standard conditions, the proposed Gauss-Simpson Quadrature Algorithm is in good agreement with the results of the traditional elasticity theory.