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代码
最新推荐文章于 2021-03-18 05:32:00 发布