function DFT=slowFT(f,N)

DFT=zeros(N,1);
for k=0:N-1
    for n=0:N-1
        DFT(k+1)=DFT(k+1)+f(n+1)*exp(-i*k*n*2*pi/N);
    end
    DFT(k+1)=DFT(k+1)/N;
end
