function [y]=mca(x)

a=13;
c=0;
m=31;
y(1)=mod(a*(x+c),m);

for i=2: 100
    
y(i)=mod(a*y(i-1)+c,m);

end
