This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A094432 #21 Jun 13 2015 00:51:20 %S A094432 0,3,24,165,1104,7347,48840,324597,2157216,14336355,95275896, %T A094432 633179973,4207956720,27965034003,185848661544,1235103986325, %U A094432 8208193936704,54549615616707,362523179503320,2409238895476197,16011202548279696 %N 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]. %C A094432 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. %D A094432 Carl D. Meyer, "Matrix Analysis and Applied Linear Algebra", SIAM, 2000, pp. 86.-87. %H A094432 Vincenzo Librandi, <a href="/A094432/b094432.txt">Table of n, a(n) for n = 1..1000</a> %H A094432 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (8,-9). %F A094432 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 %F A094432 G.f.: 3*x^2/(1-8*x+9*x^2). - _Bruno Berselli_, Aug 03 2011 %e A094432 a(4) = 165 since M^4 * [1 0 0] = [38 -203 165]. %t A094432 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 *) %Y A094432 Cf. A094431, A154245. %K A094432 nonn,easy %O A094432 1,2 %A A094432 _Gary W. Adamson_, May 02 2004 %E A094432 More terms from _Robert G. Wilson v_, May 08 2004