function output = dft(input)

%Direct computation of DFT
%usage: y=dft(x);

l=length(input);
output=zeros(1,l);				%construct empty output matrix

for i=1:l
   for k=1:l
      output(i)=output(i)+input(k)*exp(-j*2*pi/l*(k-1)*(i-1));	
   end								%compute from formula
end

%David Phillips 09507351
%Curtin University Linear Systems 202
%23 October 2000