function retval = sinm(A)
  retval = -.5 * i * (expm(i * A) - expm(-i * A));
endfunction

% ;;; Local Variables: ***
% ;;; mode: octave ***
% ;;; End: ***
