%is it a generator?
%question: is g a generator for Z*p?
clc
clear

g1 = 6;
p1 = 13;
is_generator(g1,p1)
find_generator(29)

function glist = find_generator(p)
glist = [1:p-1;1:p-1];
for i=1:length(glist)
    glist(1,i) = is_generator(i,p);
end
end

function val = is_generator(g,p)
modlist = zeros(1,p-1);
for i=1:p-1
    modlist(i)=mod(g^(i-1),p);
end
val = length(unique(modlist)) == p-1;
end