function X=mdft(x)
     N=length(x);
     X=zeros(0,N);
     n=0:(N-1);
     for k=0:(N-1);
         e=exp(-%i*2*%pi*n*k/N);
         X(k+1)=sum(e.*x);
     end;
endfunction
fp=20;
t=[-4:0.05:4];
 x=sin(2*%pi*3*t)+2*sin(2*%pi*5*t);
  X=mdft(x);
   plot(abs(X));