% ch8example9prog1.m
clear;
r(1)=13;    % 隨機種子（為奇數）
N=100;      % 樣本數
for n=2:N
    r(n)=mod(7^5*r(n-1), 2^31);% 乘同余法產生隨機數
end
r=r./2^31; % 隨機數歸一化
[f,x] = ecdf(r);   % 計算隨機數樣本r的經驗分布
P = unifcdf(x,0,1);% 均勻分布的理論概率分布函數
Kn=sqrt(N)*max(abs(f-P)) % 計算統計量
u_alpha=1.3581 % 顯著性水平為0.05時候的分位點
H=(Kn>u_alpha)
stairs(x,f);hold on;plot(x,P);% 經驗分布和理論分布曲線對比
[H,P,KSSTAT,CV] = kstest(r,[x,P],0.05)