n = 1001;
R = 1:n;
R(1) = 0;

for i = 2:1:n
    R(i) = (11 * R(i - 1) + pi) - fix(11 * R(i - 1) + pi);
end

x = (R / n);

subplot(2,2,1); hist(x);
subplot(2,2,2); plot(sort(x), [1:n] / n);
T = 1000; G = rand(1,T);
subplot(2,2,3); hist(G);
subplot(2,2,4); plot(sort(G), [1:T] / T);

dn = 1:n;
for i = 1:1:n
    dn(i) = max(abs(R(i) - (i - 1) / n), abs(R(i) - i / n));
end


P = 1:n;
for i = 1:1:n
    P(i) = 1 - 2 * exp(-2 * n * dn(i) * dn(i));
end

u = unique(R);
F = foo(u,R);
j = 0:0.001:1;
f = unifcdf(j, 0, 1);

for i = 1:1:n
    ks(i) = abs((f(i) - R(i)));
    kv(i) = f(i) - R(i);
end

kss = max(ks);
ksv = max(kv);


%plot(j, f)
%hold on
%plot(R, [1:n] / n);
%xlabel('zn CB')
%ylabel('ver')
%hold on
%plot(R, [1:n] / n - kss, 'k');