|
MATHEMATICA
|
Clear[m11, p, f]; m11 = {{2, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0}, {-1, 2, -1, 0, 0, 0, 0, 0, 0, 0, 0}, {0, -1, 2, -1, 0, 0, 0, 0, 0, 0, -1}, {0, 0, -1, 2, -1, 0, 0, 0, 0, 0, 0}, {0, 0, 0, -1, 2, -1, 0, 0, 0, 0, 0}, {0, 0, 0, 0, -1, 2, -1, 0, 0, 0, 0}, {0, 0, 0, 0, 0, -1, 2, -1, 0, 0, 0}, {0, 0, 0, 0, 0, 0, -1, 2, -1, 0, 0}, {0, 0, 0, 0, 0, 0, 0, -1, 2, -1, 0}, {0, 0, 0, 0, 0, 0, 0, 0, -1, 2, 0}, {0, 0, -1, 0, 0, 0, 0, 0, 0, 0, 2}}; p[x_] = ExpandAll[(x + Sum[m11[[n, n]], {n, 1, Length[m11]}])*CharacteristicPolynomial[m11, x]]; f[x_] = ExpandAll[1/(x^12*p[1/x])]; a = Table[SeriesCoefficient[Series[f[t], {t, 0, 35}], n], {n, 0, 35}]
|