% Demo2

% The vectors 
%               cVec = cos((0:n-1)*(2pi/n))' 
%               sVec = sin((0:n-1)*(2pi/n))'
%
% are very important. So let's learn something about them...
   clc
   disp('    n     cTc     sTs      cTs ')
   disp('----------------------------------')
   for n = 5:5:50
       [cVec,sVec] =  csVecs(n);
       cTc = cVec'*cVec;
       sTs = sVec'*sVec;
       cTs = cVec'*sVec;
       fprintf('%5d  %6.2f  %6.2f  %10.2e\n',n,cTc,sTs,cTs)
   end
   
% Conclusion: cVec'cVec = sVec'sVec = n/2  and cVec'sVec = 0
       