function [y]=congruential1(a,imax)
b=0;
m=8192;
x0=1;
y(1)=mod(a*(x0+b),m)
for i=2: imax
    
    y(i)=mod(a*(y(i-1)+b),m);
end
