fft/dst.m

   1 function Y = dst(X)\r
   2 % Y = dst(X)\r
   3 % apply sine transform to columns of X\r
   4 [n,m] = size(X);\r
   5 X1 = [ zeros(1,m); X ; zeros(n+1,m) ];\r
   6 X2 = imag( fft(X1) );\r
   7 Y = X2 (2:n+1,:);\r