If you have the Statistics Toolbox, then this is extremely easy, just using the normpdf function which encodes the formula for you.
MU = 10;
SIGMA = 1;
x = 5 : 0.01 : 15;
y = normpdf(x,MU,SIGMA)
figure
plot(x,y)
If you don't have the toolbox, then just use
y =1./sqrt(2*pi*SIGMA) * exp(-(x-MU).^2 /(2 * SIGMA));
instead.
