I need ziggurat matlab code to draw

Hello, i need a ziggurat matlab code to generate somehing like this

i'm working with (gauss) random number generator ziggurat rejection method

Answer 1

[Chunru]

where is a normalizing constant that we can ignore.

x = [0:.1:2]

x = 1×21

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000 1.1000 1.2000 1.3000 1.4000 1.5000 1.6000 1.7000 1.8000 1.9000 2.0000

y = 1/(2*pi)*exp(-x.^2/2);
stairs(x, y, 'b');
hold on
xl = xlim;
plot([xl(1)*ones(1, length(x)-1); x(2:end)], [y(2:end); y(2:end)], 'b')
% longer tais
x = [0:.1:4]

x = 1×41

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000 1.1000 1.2000 1.3000 1.4000 1.5000 1.6000 1.7000 1.8000 1.9000 2.0000 2.1000 2.2000 2.3000 2.4000 2.5000 2.6000 2.7000 2.8000 2.9000

y = 1/(2*pi)*exp(-x.^2/2);
plot(x, y, 'r')

[Mar One]

do you have a code with the tail ? basically the tail when i=0

i'm working with (gauss) random number generator ziggurat rejection method

[Chunru]

I don't have the equation for the so called ziggurat method. I am only helping the plotting part here. Do you have further info on the ziggurat rejection method?

[Mar One]

Yes this is a simple resume of what ziggurat method is: https://blogs.mathworks.com/cleve/2015/05/18/the-ziggurat-random-normal-generator/

[Chunru]

See the update above.

[Mar One]

I need a Matlab code from this principal

[Chunru]

The matlab code is given above. What else do you need?

[Mar One]

[Mar One]

from this above