clear; clc;

a = 8121;
b= 28411;
M = 134456;
x0 = mod(tic, M);
x = zeros(1, 1e3);
x(1) = x0;
for i = 2 : 1 : length(x)
      x(i) = mod(a * x(i - 1) + b, M);
end
x = x / M;