newseq=zeros(1,31);
sim('msequence');
mseque=zeros(31,1);
for i=1:31
    if mseq(i)==0
        mseque(i)=1;
    else 
        mseque(i)=-1;
    end
end
for i=0:100
    for j=0:30
    newseq(j+1)=mseque(mod(i+j,31)+1);
    end
autocorrelate(i+1)=newseq*mseque/31;
end
k=1:101
plot(k-1,autocorrelate);