close all;


fs = 1000000;
t = -5:1/fs:5;
g = @(t, sigma) 1./sqrt(2.*pi.*sigma.^2).*exp(-t.^2./(2.*sigma.^2));

for s = [0.1 0.001 0.00001 0.0000001]
    figure;
    subplot(2,1,1);
    plot(t, g(t, s));
    title(strcat('Gaussian s=', num2str(s)));

    subplot(2,1,2);
    PlotModSpectrum(g(t, s), fs);
    title(strcat('Centered CDFT s=', num2str(s)));
end