function [seq]=mrand(connection)
m=length(connection);
L=2^m-1;
registers=[zeros(1,m-1) 1];
seq(1)=registers(m);
for i=2:L
    new_reg_cont(1)=connection(1)*seq(i-1);
    for j=2:m
        new_reg_cont(j)=bitxor(registers(j-1),connection(j)*seq(i-1));
    end
    registers=new_reg_cont;
    seq(i)=registers(m);
end