Remove["`*"];
zeroPolynomial[A_] := Module[{n = Length[A], base, m, s},
base = NestList[Dot[A, #] &, IdentityMatrix[n], n^2];
m = Transpose[Flatten /@ base];
s = NullSpace[m][[1]];
Total[Table[x^i*s[[i]], {i, Length[s]}]]
]
求解n阶方阵零化多项式的mathematica代码
最新推荐文章于 2024-08-03 16:33:59 发布