% C = dst4mtx(M)
% Provide the type IV discrete sine transform matrix

function S = dst4mtx(M);
for I = 1:M
   for L = 1:M
     S(I,L) = sqrt(2/M)*sin((pi/M)*((I-1)+0.5)*((L-1)+0.5));
   end
end