Gaussian distribution – how to plot it in Matlab

where 

The formula above can me coded in Matlab easily, like this: 

function f = gauss_distribution(x, mu, s)
p1 = -.5 * ((x - mu)/s) .^ 2;
p2 = (s * sqrt(2*pi));
f = exp(p1) ./ p2; 

Now, let’s use it in an example. 

We produce 500 random numbers between -100 and 100, with mean m = 0 and standard deviation s = 30. The code is: 

a = -100; b = 100;
x = a + (b-a) * rand(1, 500);
m = (a + b)/2;
s = 30; 

Then, we plot this information using our bell curve: 

f = gauss_distribution(x, m, s);
plot(x,f,'.')
grid on
title('Bell Curve')
xlabel('Randomly produced numbers')
ylabel('Gauss Distribution') 

The produced shape is:

An important property of this bell-shaped curve is that the values less than one standard deviation from the mean (between green lines below) represent approximately 68.2% of the area under the curve, while two standard deviations from the mean (between red lines below) take about 95.4%, and three standard deviations account for about 99.7% of the area.