function [] = tw2d(polynomial)
p = polynomial(1);
r = ones(1,p);
upper = power(2,p)-1;
x = zeros(upper,1);
for i =1:upper
        n = 0;
       for bits = 1:p
           [b,r] = tausworthe(r,polynomial);
           n = 2 * n + b;
       end
       x(i) = n/upper;
end;
x1 = x(1:(upper-1));
x2 = x(2:upper);
plot(x1,x2,'.');