Moved internal gaussian functions to top of method.
Added check of the length of the peak vector and case for a length less than 2