|
|
A094432
|
|
a(n) = rightmost term in M^n * [1 0 0]. M = the 3 X 3 stiffness matrix [1 -1 0 / -1 4 -3 / 0 -3 3].
|
|
4
|
|
|
0, 3, 24, 165, 1104, 7347, 48840, 324597, 2157216, 14336355, 95275896, 633179973, 4207956720, 27965034003, 185848661544, 1235103986325, 8208193936704, 54549615616707, 362523179503320, 2409238895476197, 16011202548279696
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
A094431(n) = left term in M^n * [1 0 0]. A stiffness matrix in Hooke's Law governs the force on nodes of stretched or compressed springs (refer to A094431). a(n)/a(n-1) tends to 4 + sqrt(7) = 6.6457513...; a(n)/A094431(n) tends to 2 + sqrt(7). A stiffness matrix is symmetric.
|
|
REFERENCES
|
Carl D. Meyer, "Matrix Analysis and Applied Linear Algebra", SIAM, 2000, pp. 86.-87.
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (3/(2*sqrt(7)))*((4+sqrt(7))^(n-1)-(4-sqrt(7))^(n-1)). For n>1, a(n) = 3*A154245(n-1). - Francesco Daddi, Aug 02 2011
|
|
EXAMPLE
|
a(4) = 165 since M^4 * [1 0 0] = [38 -203 165].
|
|
MATHEMATICA
|
Table[(MatrixPower[{{1, -1, 0}, {-1, 4, -3}, {0, -3, 3}}, n].{1, 0, 0})[[3]], {n, 21}] (* Robert G. Wilson v, May 08 2004 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|