产生0-1均匀分布伪随机数的乘同余法（matlab）

%%乘同余法计算随机数列
%%公式：Sk+1=A*Sk*mod(M)
M=64;%
A=11;
x0=1;
Sk=x0
%%输出n位随机数列
n=10;
for i=1:n
    Sk=mod(A*Sk,M)
end

%%混合同余法计算随机数列
%%迭代公式：Sk+1=(A*Sk+C)mod(M),当C为0时，混合同余法退化为乘同余法
clc;clear
M=64;
A=13;
C=7;
x0=1;
Sk=x0
%%输出n位随机数列
n=10;
for i=1:n
    Sk=mod(A*Sk+C,M)
end