MATLAB: Using Linear congruential to generate 10,000 uniform random variables

lcgrng

I want to generate uniform random variables between 0 and 1 by using X(n+1)=(1664525*X(n)+1013904223) mod^32, but i couldn't.
Please help…
this is my code:
clear all; clc;
a=1664525; c=1013904223; m=2^32;
for n=1:10000
X(n)=mod((a*X(n)+c), m);
end
disp(X);

Best Answer

You need an intial value for X(1); and you need to set
for n=2:10000
X(n+1)=mod((a*X(n)+c), m);
end