At present, the calculation of nuclear magnetic resonance sounding simulation is mainly based on homogeneous half space model or equivalent half space model. As the equivalent model makes the calculation result vary greatly from the actual response, in order to avoid the error brought by the equivalent model, the paper first implement the oscillating integral calculation by using gauss numerical integration and gets the scattering field of the underground space based on layered model harmonic field electromagnetic response derived by the recursive method, then obtains nuclear magnetic resonance (NMR) sounding response under one-dimensional medium case. In order to verify the effectiveness of the proposed algorithm, the layered model of calculation method is used to calculate the homogeneous half space of NMR response, being consistent with the actual uniform half-space model results, indicating that the calculation method given in the text reliable. The paper calculated the different orientation of the cutting angle under toppled layered model, and compared the differences in cross-sectional direction of cutting angle. Analysis is showed that NMR response is calculated more accurately by layered model recursive method than by the method of approximate half-space model.