MatrixFunction(A) コマンドは、重複を含み、A の固有値 lambdaそれぞれに対して、[lambda, F( lambda )] を補間して得られる行列を返します。 ここで、行列多項式は r(lambda) = F(lambda) - p(lambda)*q(lambda) であり、p(x) は特性多項式、q(lambda) は商、r(lambda) は 剰余です。

output オプション outopts の詳細は、 Matrix を参照してください。

with(LinearAlgebra):
A := Matrix([[-13, -10], [21, 16]]);

MatrixFunction(A, v^2, v);

MatrixFunction(A, sin(x), x);

MatrixFunction(A, exp(x), x);

MatrixExponential(A, t);

MatrixFunction(A, exp(x*t), x);

MatrixFunction(A, ln(x), x, outputoptions=[readonly=true]);

B := Matrix([[y]]);

MatrixExponential(B);

MatrixFunction(B, exp(y), y);

MatrixFunction(B, exp(y));