AbstractUsing the basic Boussinesq's equation, the expression for the vertical stress distribution (σz) underneath any point on the ground surface due to a general triangular loaded region in a preferred orientation with a linearly varied loading has been successfully derived. When the triangle is not in a preferred orientation, a simple axis transformation is required and the expression will be equally applicable. Based on this expression, σzdue to an arbitrarily shaped loaded foundation can simply be determined by first triangulating the loaded area and summing up the contributions from each generated triangular region. The procedures for triangulating and calculating the stress distribution can be simply automated through computer program