A symmetrical isolated rectangular footing with centered biaxial overturning develops soil pressure that shifts to counter balance the loads. The highest soil pressure is at a corner. The objective of this paper is to extend the uniaxial soil pressure solution to include biaxial loads and to document a simple and understandable way to directly calculate the shape of the soil pressure distribution. Another objective is to make the solution suitable for automation. In uniaxial overturning there are two transition shapes, trapezoidal and triangular. In biaxial overturning there are three transition shapes and they form 4, 5 & 6 sided polyhedrons. This analysis calculates those volumes and compares them to the design vertical load to determine the characteristic shape of the soil pressure distribution. The calculation then proceeds to converge on the precise shape and calculate its centroid and moment capacity. Assemblies of tetrahedrons are used to model all of the soil pressure shapes. The advantage of this methodology is that matrix algebra can be used to organize the calculations and make them computationally efficient. The assumed soil pressure and footing dimensions can be adjusted until the calculated moment capacity matches the overturning moment.