%Peter Parente
%make a non-orthnormal DFT kernel matrix
function M = dft(N)
for n = 0:N-1
    for k = 0:N-1
        M(n+1,k+1) = exp(-j*n*k*2*pi/N);
    end
end