I am trying to plot a Gaussian normal probability density function with 500 points, mean=1.5 standard deviation 0.5. The Below is my code. My code doesn't have a error but the plot is not accurate.

I am trying to plot a Gaussian normal probability density function with 500 points, mean=1.5 standard deviation 0.5. The Below is my codecode doesnt have a error but the plot is not accurate.

function

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

code:

x = [xmin:xmax];
m = 1.5
s = 0.6
N=500;
xmax= m+ N/2*s;
xmin = m- N/2*s;
f = gauss_distribution(x, m, s);
plot(x,f,'.')
grid on
xlabel('Randomly produced numbers')
ylabel('Gauss Distribution')

Plot: is attached

Answer 1

[Jose Marques]

Hello Sam. I used:

Y = normpdf(x,m,s)

and give the same result.

Seems that is correct.

[Jose Marques]

The way to generate the x values is correct?

[Sam17]

i was asked to generate x values in similar manner in my question.

[Sam17]

I am not getting the exact plot.