# This script shows the distribution of points produced by repeated calls to randcong.m
clear dat
a=23098209;
npts = 10000;
for i=1:npts
	a=randcong(a);
	dat(2*i-1)=a;
	x(i)=a;
	a=randcong(a);
	dat(2*i)=a;
	y(i)=a;
endfor

figure(1)
hist(dat)
title("Random Congruential Histogram")
figure(2)
plot(dat,'.')
title("Random Congruential Successive Points")
figure(3)
plot(x,y,'.')
title("Random Congruential Successive x-y Points")