Updates 
14 Apr 2009 
1.1 
Bug fixes. No major features added. 
27 Apr 2009 
1.2 
Expanded description 
13 Apr 2011 
1.5 
Added ability to access the xy data for the individual fitted peaks, via the optional output parameters xi and yi. 
25 Apr 2011 
1.8 
Added two new peakshapes: 6=Equalwidth Gaussians, and 7=Equalwidth Lorentzians. 
04 May 2011 
1.9 
Expanded the description 
25 Jul 2011 
1.10 
Several bug fixes. Reformatted peak table on graph. 
10 Aug 2011 
1.15 
Version 1.8: Aug, 2011. Takes AUTOZERO setting as 9th input argument; improved data input argument flexibility.

24 Aug 2011 
1.16 
Bug fix 
28 Sep 2011 
1.17 
Version 2.1: Sept, 2011. Accepts AUTOZERO 0 (none), 1 (linear), or 2 (quadratic). 
20 Oct 2011 
1.19 
Version 2.2: October, 2011. Adds exponential pulse and sigmoid models 
19 Jan 2012 
1.20 
Version 2.3: January, 2012. Bug fixes in background subtraction modes and
in handlng very small data sets. 
19 Jan 2012 
1.21 
Version 2.3: January, 2012. Bug fixes in background subtraction modes and
in handlng very small data sets. 
18 May 2012 
1.23 
Version 2.4: May, 2012, Exponential broadening uses normal rather than circular convolution. 
18 May 2012 
1.24 
Version 2.4: May, 2012, Exponential broadening uses normal rather than circular convolution. 
08 Jun 2012 
1.25 
Version 2.5: June, 2012, Allows zeros as placeholders for unspecified input arguments. 
18 Jun 2012 
1.26 
Version 2.6: June, 2012. Added fixedwidth Gaussian and Lorentzian peak shapes (shape numbers 11 and 12). 
04 Sep 2012 
1.27 
Version 3: September, 2012, adds the ability to estimate the uncertainty of peak parameters using the bootstrap sampling method. 
12 Sep 2012 
1.29 
Version 3.1: September, 2012. Unlimited Number of peaks. Bug fixes. 
20 Sep 2012 
1.30 
Version 3.3: September, 2012. Added 11th input argument ('plots') to turn plotting OFF (plots=0) or on (plots=1); added Gaussian/Lorentzian blend, bifurcated (asymmetrical) Gaussian and Lorentzian. 
15 Oct 2012 
1.31 
Version 3.4: October, 2012. Works in Matlab or Octave 3.6.1 
02 Nov 2012 
1.33 
Version 3.4.2; Slight improvement in speed of exponentiallybroadened shapes. Works in Matlab or Octave 3.6.1 
15 Jan 2013 
1.36 
Version 3.51: January, 2013. Improved accuracy of linear autozero calculation. Improved calculation of default "start" guess when not specified in the input arguments. 
21 Feb 2013 
1.37 
Version 3.6: February, 2013. Addition of fixedposition Gaussian shape (16) and fixedposition Lorentzian shape (17). 
13 Sep 2013 
1.39 
Version 3.9 adds exponentially broadened Lorentzian(peak shape 18); and alpha function (peak shape 19). 
11 Nov 2013 
1.40 
Version 4.2 ,corrects some bugs, adds an additional autozero mode that subtracts a flat baseline without requiring that the signal return to the baseline at both ends of the signal segment, and adds a Voigt profile peak shape. 
03 Feb 2014 
1.41 
Version 4.31: Jan, 2014. Adds 12th input argument, for + or +/ peak mode. 
11 Feb 2014 
1.42 
Version 5: Feb, 2014. Adds multipleshape models, designated by using a vector as the 5th input argument. See examples 17 and 18. 
20 Jun 2014 
1.43 
Version 5.4: June, 2014. Replaces bifurcated Lorentzian with BreitWignerFano resonance peak (Shape=15). 
29 Jul 2014 
1.44 
Version 5.5: July, 2014. Adds shape 21 (triangular) and 25 (Lognormal distribution). Moves peak table to lower panel and plots residuals as red dots. 
02 Sep 2014 
1.45 
Version 5.7: August, 2014. Adds minimum width constraint as 13th input argument (See example 19); Can be a vector for multiple peak shapes. The default if not specified is the independent variable (x) interval. 
24 Jan 2015 
1.46 
Adds an additional optional input argument, new peak shapes, and new output arguments. 
24 Jan 2015 
1.47 
Adds an additional optional input argument, new peak shapes, and new output arguments. 
24 Apr 2015 
1.48 
Version 7.1: April, 2015, adds peak shapes with three unconstrained iterated variables: 30=Voigt (variable alpha), 31=ExpGaussian (variable time constant), 32=Pearson (variable shape factor), 34=Gaussian/Lorentzian blend (variable percent). 
24 Apr 2015 
1.49 
Version 7.1: March, 2015, adds peak shapes with three iterated variables: 30=voigt (variable alpha), 31=ExpGaussian (variable time constant), 32=Pearson (variable shape factor), 34=Gaussian/Lorentzian blend (variable percent). See Examples 2528. 
22 Jun 2015 
7.2 
Version 7.2: June, 2015, fixes bug when the bootstrap is used with multiple shape models. 
02 Aug 2015 
7.4 
August, 2015, Allows fixed width models (shapes 11 ands 12)
to have different fixed widths for each peak (specified by a vector);
also adds a fixedheight Gaussian model (shape 34) with heights specified
by a vector (input argument 10) 
24 Oct 2015 
7.45 
This version corrects the Voigt profile model so that it reports the FWHM of the peak. 
18 Apr 2016 
7.9 
Version 7.9: April, 2016, Added shapes 38=variable time constant ExpLorentzian; 40=sine wave; 41=rectangle;
42=flattened Gaussian; 43=Gompertz function (3 variable logistic); 44=1exp(k*x). "extra" defaults to 1 if not specified in input arguments. 
18 Apr 2016 
7.9 
Edited discription 
15 Jul 2016 
8.0 
Version 8: July 2016, adds shape 46 (quadslope), for fitting curved baselines. See example 33 in the builtin help. 
29 Jul 2016 
8.1 
Minor bug fix 