function script3b
% investigates theperiod of lcg generators with m=97  ( a prime)
for a=1:2:35
    v=period(1,a,1,97);
    disp(['(ax+1)(mod 97) where a= ' num2str(a) '  period= ' num2str(v(1)) ])
end
