clear all
s=20;
a=13;
b=1;
m=16;
for I=1:m
       n=a*s+b;
u=mod(n,m);
v=u/m
s=u;
end
