function E=sinm1(A)
%SINM1 函數用于由 Taylor 展開方法計算矩陣的正弦。

%Designed by Prof D Xue (c) 2000
E = zeros(size(A)); F = A; k = 1;
while norm(E+F-E,1) > 0
   E = E + F;  F = -A^2*F/((k+2)*(k+1)); k = k+2;
end