% call: paired_matrix
% calculate the matrix NxN of the paired transform 
% when N is a power of two 

    function PM=paired_matrix(N)

    PM=zeros(N);    
    for m=1:N
        y=zeros(1,N);
        y(m)=1;
        PM(:,m)=fst_1d(y);
        % PM(:,m)=fastpaired_1d(y);
    end;