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