The line‐integral method of Hubbert (1948) is extended to obtain the gravity anomalies of two‐dimensional bodies of arbitrary cross‐sections with density contrast varying linearly with depth. The cross‐section is replaced by an N‐sided polygon. The coordinates of two vertices of any given side are used to determine the associated contribution to the gravity anomaly. The gravity contribution of each side is then summed to yield the total gravity effect. The case where density contrast varies exponentially with depth is also considered. This technique is used to obtain the structure of the San Jacinto Graben, California, where sediments filling the graben have an exponential increase in density with depth.