clear all; close all;

gen = @(x,z) mod(z*x*pi, 1);

Ta = -1; Tb = 1; Tc = 0;
c = 1;
N = 1024;
tooths_U = 17;
tooths_G = 31;

U = zeros(1, N*5);    % rozklad rownomierny [0  1]
G = zeros(1, N*5);    % rozklad rownomierny [-1 1]
T = zeros(1, N);      % rozklad trojkatny [-1 1]

U(1) = 0.218731; G(1) = 0.218731;
for k=2:N*5
    U(k) = gen(U(k-1), tooths_U);
end

for k=2:N*5
    G(k) = gen(G(k-1), tooths_G);
end
G = G*(Tb-Ta) - (Tb-Tc); % [0,1] -> [-1,1]

iter = 1; k = 1;
while(k<5*N)&&(iter < N)
    if(G(k) < Tc)
        temp = 2*(G(k)-Ta) / (Tb-Ta)*(Tc-Ta);
    else
        temp = 2*(Tb-G(k)) / (Tb-Ta)*(Tb-Tc);
    end

    if(c*U(k) <= temp)
        T(iter) = c*G(k);
        iter = iter + 1;
    end
    k = k+1;
end
pkt_T = (Tc-Ta)/(Tb-Ta);
figure(1);
histogram(T, 100);