        function d=F_DM(m,N)

        j=[0:1:N-1];  x=[2*pi*j/N]; 

        for j=1:N
            for k=1:N
               A(j,k)=(-1)^(j-1)*exp(i*(k-1)*x(j))*(i*(k-1-N/2))^m;
               B(j,k)=(-1)^(k-1)*exp(-i*(j-1)*x(k));
            end
        end

        d=((1/N)*A)*B;
