sinm - matrix sine function
sinm(x) is matrix sine of x matrix.
A=[1,2;2,4]; sinm(A)+0.5*%i*(expm(%i*A)-expm(-%i*A))