function [mseq]=m_seq(init);
n=length(init);
N=2^n-1;
for i=1:N
    for k=1:n             %加上去的
    cinit(k)=init(k);    %加上去的
    end                   %加上去的
    init(1)=madd(init(3),init(10));
    for j=2:n
        init(j)=cinit(j-1);
    end
    mseq(i)=init(10);
end