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. 