%function c=gauss(dt,f0)
% GAUSS generate a Gaussian wavelet
% input (dt,f0)

dt=.01;
f0=5;

T=1/f0;
nt=4*T/dt;
c=1:nt;
c=c*0;

        t0=2*T/dt;
        b=1/(T*sqrt(2*pi));
        a=T*sqrt(2);

for it=1:nt
            t=(it-t0)*dt;c(it) = b*exp(-(a*t)^2);
end



