findpeaks(M31,'MinPeakProminence',1,'Annotate','extents')
[pks,locs,widths,proms] = findpeaks(M31,'MinPeakProminence',1,'Annotate','extents');
widths
proms

W1 = widths*(1e-11)

findpeaks(M37,'MinPeakProminence',1,'Annotate','extents')
[pks,locs,widths,proms] = findpeaks(M37,'MinPeakProminence',1,'Annotate','extents');
widths
proms

W1 = widths*(1e-11)