The Parker-Oldenburg inversion algorithm is defined as follows:
$$h(x)=F^{-1}(\frac{-F[\Delta g(x)]e^{-kz_0}}{2\pi G \rho}-\displaystyle\sum_{n=2}^{\infty} \frac{k^{n-1}}{n!}F[h^{n}(x)])$$
Where F is the Fourier transform operator; ∆g(x) is gravity anomaly; h(x) is the depth to the interface; k is the wavenumber; G is the gravitational constant; ρ is the density contrast across the interface; and z0 is the average depth of the horizontal interface. According to this algorithm, the topography of the interface density can be estimated by an iterative inversion procedure.