The Moho discontinuity is the dividing line between the crust and the mantle. Several methods, including the Parker-Oldenburg, Bot, and Euler methods, are used to determine the Moho depth. The Parker-Oldenburg method’s results are based on inverse iteration, which yields the Fourier transform. This algorithm’s great speed, which is due to its based on the fast Fourier transform, is one of its important features. The Parker-Oldenburg method has greatly and largely resolved the issue of the decrease in sensitivity to depth that exist in inversion of geophysical data and that prevents us from having adequate resolution and precision in depth for investigations. The MATLAB program was used to apply this approach, which is based on actual data from the Mahalat area in the Zagros region. After modifications and re-gridding, the algorithm stopped after 10 iterations, and The results were in accordance with the previous investigations in this area .