function randmtx=randCS(row,col,seed)
A=4;
Xb=2.5;
randmtx(1,1)=seed;
for i=2:row*col
    randmtx(l,i)=A*sin(randmtx(l,i-1)-Xb)*sin(randmtx(l,i-1)-Xb);
end
randmtx=reshape(randmtx,row,col);