MaximaでMathematica組み込みシンボルを実装していく
自作関数一覧PowerMod(a,b,p):=map(lambda([u],power_mod(a,u,p)),b)$
MemberQ(a,b):=map(lambda([u],member(u,b)),a)$
Primitiveroot(n):=block(
[a:[],i,j:makelist(i,i,1,n-1)],for i:2 thru n-1 do (
if is(sort(PowerMod(i,j,n))=j) then a:endcons(i,a)),a)