#DFT Implementation
%X=0.0;
%A(k)=abs(X(k))/N;
% X1(k+1,n+1)=summe;
%n=0:N-1;
xn=input('Please put the vector where the data is going to be evaluated: ');
N=0:length(xn);
% N could maybe be N=input('Please put the "N" factor where the data is going to be evaluated: ');
L=length(xn);
%need zeros function maybe?
for k=0:N-1
  for n=0:N-1
    summe=xn(k+1)+xn(n+1)*exp(-j*2*pi*k*n)/N;
end
end
disp('x(k)=');
disp(summe);
k=0:N-1;
n=0:N-1;
xk=zeros(
subplot(2,1,1), stem (L,abs(summe)), xlabel('k'), ylabel('|Xk|'), title('Magnitude Plot')
subplot(2,1,2), stem(L,angle(summe)), xlabel('k'), ylabel('Angle (Xk)'), title('Phase Plot')