function val = morlet_wavelet(t)
% evaluate Morlet wavelet
    w0 = 5;  % scale factor for the frequency of the sinusoidal component
    val = (pi^-0.25) * exp(-0.5*t^2) * (exp(1i*w0*t) - exp(-0.5*w0^2));
end