function C = mk_DCT_matrix_C(N)
% 
% Written by:
% -- 
% John L. Weatherwax                2005-08-04
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
% Epage 
% 
%-----

C = zeros(N,N); 
C(:,1) = ones(N,1); 

n_values = 0:(N-1);
k_values = 1:(N-1); 

n_times_k  = ( 2 * n_values + 1 )' * k_values;
C(:,2:end) = sqrt(2) * cos( pi * n_times_k / (2 * N) );
C          = C/sqrt(N);