The curve fitting method for analyzing temperature-programmed desorption spectra has been developed. The method was applied to butene-1 desorption from the NaY zeolite and found to be useful for determining the surface coverage dependence of the apparent activation energy of desorption as well as the pre-exponential factor.