求特征值给出以下结果，说明这个极限不可能存在
ans =
   1.000008538586829
  -0.043025817257730
  -0.072282721329098
  -0.058200000000000
方法也不正确，矩阵幂的极限一般使用特征值来求，具体可参考常见线性代数书籍中关于矩阵特征值的部分，我编了以下函数实现

function M = matlim(A)
%compute limit of matrix exponential A^n by eig, make sure A be square
norma = norm(A);
if norma>1
   disp('norm larger than 1, no limit exist')
   return
end

[U,eiga] = eig(A);

eiglim = round(eiga);

Uinv = U\eye(size(A));

M = U*diag(eiglim)*Uinv;
end