a=2;
b=a;
x=1;
k=1;
clear iterates

for m = 10^7 + 1 : 2 : 10^7 + 99
    a=2;
    b=a;
    x=1;
    n=m-1;
    while n > 0
        d=rem(n,2);
        if d==1
                x=rem(x*b,m);
                iterates(k)=x;
                k=k+1;
        end
        b=rem(b*b,m);
        n=(n-d)/2;
    end
end