clear all;
a = 7 ^ 5;
b = 0;
m = 2 ^ 31 - 1;
x0 = 4294967291;
x = x0;
N = 100000;

val = [1 2 3 4];
prob = [0.15 0.45 0.05 0.35];
freq = zeros(1, 4);

for j = 1:N
    i = 0;
    F = 0;
    x = mod(a * x + b, m);
    U = x / m;
    while F < U
        i = i + 1;
        F = F + prob(i);
    end
    xg(i) = val(i);
    freq(i) = freq(i) + 1;
end
freq = freq / m
subplot(1, 2, 1);
bar(val, prob);
title("Variable aléatoire");
subplot(1, 2, 2);
bar(val, freq);
title("variable générées");