% Simple function to trace a Gaussian curve
mu = 0;
sigma = sqrt(1);

for i = 1:20001
    x(i) = i/1000-10;
    PG(i) = (1/(sigma*sqrt(2*pi)))*exp(-((x(i)-mu)^2)/(2*sigma^2));
end

plot(x,PG);
axis([-5 5 -1 1]);
axis 'auto y';
xlabel('x');
ylabel('P(x)');
title('Distribution Gaussienne / normale standard');

text(2,.35,'\mu = 0');
text(2,.33,'\sigma^{2} = 1');

FWHM = sigma * sqrt(8*log(2));
PGFWHM = (1/(sigma*sqrt(2*pi)))*exp(-((FWHM-mu)^2)/(2*sigma^2));

text(2,.31, ['FWHM = ' num2str(FWHM)]);
text(2,.29, ['P(FWHM) = ' num2str(PGFWHM)]);