function [] = lcg2d(polynomial, seed)
modulus = polynomial(2);
upper = modulus-1;
r = seed;
x = zeros(upper,1);
for i =1:modulus
       r = lcg(r,polynomial);
       x(i) = r;
end;
x1 = x(1:(modulus-1));
x2 = x(2:modulus);
plot(x1,x2,'.');
