function y=mca(x)

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

for i=2: 10
    
y(i)=a*mod(y[i-1]+c,m)

end