By integrating the classical Boussinesq expression we derive analytically the vertical stress distribution induced by pressures distributed with arbitrary laws, up to the third order, over polygonal domains. Thus, one can evaluate in closed form either the vertical stress produced by shell elements, modelling raft foundations by finite elements, acting over a Winkler soil or those induced by a linear pressure distribution simulating axial force and biaxial bending moments over a pad foundation. To this end we include charts and tables, both for rectangular and circular domains, which allow the designer to evaluate the vertical stresses induced by linear load distributions by hand calculations. The effectiveness of the proposed approach is witnessed by the comparison between the analytical results obtained with the proposed formulas and the numerical ones of a FEM discretization of the soil associated with the loading distribution induced by a foundation modeled by plate elements resting on a Winkler soil.The classical approach to finding stress and displacements in an elastic, homogeneous and isotropic half-space due to surface tractions has been first developed by Boussinesq [4] who provided the solution for a point load by making use of the potential theory [3].The case of uniform pressure applied within a circular domain was addressed by Lamb [19] while Love [20] considered the case of isotropic uniform pressure applied within a rectangular region.However, even in many natural soils deposited through a geologic process of sedimentation over a period of time stiffness becomes greater with depth. Analytical continuum models reflecting this property were first studied by Gibson [16] by assuming incompressibility and a linear increase of soil stiffness with depth. In later papers Gibson further refined his model and presented on account of his results in a Rankine lecture [17].Subsequently Vrettos [30] solved the Boussinesq problem for an half-space exhibiting an exponential behavior of stiffness yet with a bounded value at infinite depth. Half-spaces of this kind were later studied by Selvadurai in the framework of contact static problems and a survey account of his result were presented in [24].Since the pioneering work by Gibson an abundant literature has been published on the evaluation of stresses and displacements in inhomogeneous and/or anisotropic half-spaces, see e.g. [31] for a fairly complete account.Nevertheless, Boussinesq solution is still widely used mainly for its elegance and conciseness, for instance to evaluate the order of magnitude of the depth of influence of a foundation load [7] and the related settlement. Actually, it is usually assumed that the variation of vertical stress with depth below a loaded region on the surface can be predicted with sufficient accuracy using linear, homogeneous, isotropic elastic theory.Further applications of this theory range, to quote only a few, from space geodesy to contact mechanics, rock mechanics and geomechanics.For instance, in the former case position...