Abstract: Levenberg-Marquardt method (LM) was proposed to
be applied as a non-linear least-square fitting in the analysis of a
natural gamma-ray spectrum that was taken by the Hp (Ge) detector.
The Gaussian function that composed of three components, main
Gaussian, a step background function and tailing function in the lowenergy
side, has been suggested to describe each of the y-ray lines
mathematically in the spectrum. The whole spectrum has been
analyzed by determining the energy and relative intensity for the
strong y-ray lines.