%constructing linear sequence, feedback polynomial and initial vector given;
function d=lseq(f,din); %primitive polynomial and initial state are input variables;
n=length(f)-1; N=2^n-1; %length of sequence;
d=[din zeros(1,N-n)]; %template for a sequence;
for k=n+1:N
    d(k)=sum(f(1:n).*d(k-n:k-1)); d(k)=rem(d(k),2);
end; %linear sequence is ready;