function res = congruence(p0,n)
p = zeros(n,1);
 
prem=1500007;   %nombre premier
A=1300;         %racine primitive de P
 
p(1)=max(1,rem(p0,prem)); %s'assure que 1<=p0<prem
 
for k=1:n-1
    p(k+1)=rem(A*p(k),prem);
end
 
res=p/prem;