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

figure(4)
hist(dat,[-10:.5:5])
title("Iterative Log Map Histogram")
figure(5)
plot(dat,'.')
title("Iterative Log Map Successive Points")
axis([1,2*npts,-10,5])
figure(6)
plot(x,y,'.')
axis([-10,5,-10,5])
title("Iterative Log Map Successive x-y Points")