Storing values in a loop

plot of x does not give a normal distribution as expected, seems to be an error in my code

n = 20000;

x = zeros(n,1)

x(1) = Mint(1)

for j = 1:n-1

x_c = normrnd(x(j),25);

if rand > min(1,normpdf(x_c)/normpdf(x(j)))

x(j+1) = x_c;

else

x(j+1) = x(j);

end

end

Answer 1

[Bhaskar R]

Already your 20000 iterations of the x is stored in the variable in x

and for the histogram use command as

hist(x)