%用乘同余法產生隨機數
A=6;  N=100;      %初始化；
x0=1;  M=255;
for k=1:N          %乘同余法遞推100次；
   x2=A*x0;        %x2和x0分別表示xi和xi-1；
  x1=mod (x2,M);    %將x2存儲器的數除以M，取余數放x1（xi）中；
   v1=x1/256;       %將x1存儲器的數除以256得到小于1的隨機數放v1中；
   v(:,k)=v1;        % 將v1中的數（ ）存放在矩陣存儲器v的第k列中，v(:,k)
%表示行不變、列隨遞推循環次數變化；     
   x0=x1;           %xi-1= xi；
   v0=v1;
end                 %遞推100次結束；
v2=v                %該語句末無‘；’，實現矩陣存儲器v中隨機數放在v2中，%且可直接顯示在MATLAB的window中；
k1=k;
%grapher             %以下是繪圖程序；
k=1:k1;
plot(k,v,k,v,'r');
xlabel('k'), ylabel('v');tktle('(0-1)均勻分布的隨機序列')