function seq=rand(n,a,c,m,seed)
seq=seed;
x=seed;
for i=1:
x=rem(a*x+c,m)
seq=[seq x];
end