Code:

x = -40:0.01:40;
K = 10;
mu = 0;
sigma = 5;

y = K * exp(-(((x - mu).^ 2)/(2*sigma.^2)));

plot (x, y);