a=[9:2:55];
m=2^10;

for i=1:length(a)
figure
x=lcg(1,a(i),0,m,m);
n=length(x);
plot(x(1:n-1),x(2:n),'.');
title(['a=', num2str(a(i)), ', c=0,m=1024   period=', num2str(period(1,a(i),0,2^10))])
xlabel('x(n)')
ylabel('x(n+1)')
pause(1)
end