function y=gaussian(x,a,m,s);
% x is time point vector
% m is peak max time point
% s is stdev (peak width)

y=a/(s*sqrt(2*pi))*exp(-1*(x-m).^2/(2*s^2));
