m=1021;
a=2:100;
for i=1:length(a)
disp(['a=', num2str(a(i)), ', c=0,m=1021   period=', num2str(period(1,a(i),0,m))])
end