% C = dct4mtx(M)
% Provide the type IV discrete cosine transform matrix

function C = dct4mtx(M);
for I = 1:M
   for L = 1:M
     C(I,L) = sqrt(2/M)*cos((pi/M)*((I-1)+0.5)*((L-1)+0.5));
   end
end