function y=gen5(x,N)
a=16807;
rm=2147483647;
q=127773;
r=2836;
y=zeros(N,1);
for i=1:N
A=floor(x/q);
test=a*x - A*(a*q+r);
if test>0
x=test;
else
x=test+rm;
end
y(i)=x/rm;
end