function p=period(x0,a,c,m)
x=lcg(x0,a,c,m,m+2);
p=find(x(2)==x(3:m+2));
