Force multipeak fitting

Hi there, I am new to this program, here we have written a program to calculate peak area of a series spectra in one file ,for example, wave0 to wave 8000, to construct one intensity image, , while we found that there are two peaks so close that general fitting can not recognize these two peaks as one, we can force fit the specific peak range into two peaks, can any one can tell me how to do batch fitting so that we needn't to do so one by one ?
Thanks
Thanks to John's help, we made this script to do force fitting with two Gaussian functions, well, so far, the baseline is a little problem, over all, it works, only some times the fitting is not that reasonable, check the code, and hope it helps.
multipeakfit.ipf