OFFSET
1,2
COMMENTS
The matrix is derived from the polynomial 8x^3 + 4x^2 - 4x - 1 shown on page 204 of "Advanced Trigonometry"; cos(2*Pi/7), cos(4*Pi/7), and cos(6*Pi/7) are roots of this polynomial.
REFERENCES
C. V. Durell & A. Robson, "Advanced Trigonometry", Dover 2003, p. 204.
LINKS
Index entries for linear recurrences with constant coefficients, signature (4, 4, -8).
FORMULA
a(n)/a(n-1) tends to 1/cos(3*Pi/7) = 4.4939592074...
From Colin Barker, Feb 01 2012: (Start)
a(n) = 4*a(n-1) + 4*a(n-2) - 8*a(n-3), n > 4.
G.f.: x*(1 + 4*x - 16*x^2 + 16*x^3)/(1 - 4*x - 4*x^2 + 8*x^3). (End)
EXAMPLE
a(3) = 20 since M^3 * [1 1 1] = [120 -20 8]. Take the absolute value of the center term.
MATHEMATICA
Table[ Abs[ MatrixPower[{{-4, 4, 8}, {1, 0, 0}, {0, 1, 0}}, n].{1, 1, 1}][[2]], {n, 23}] (* Robert G. Wilson v, Apr 28 2004 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Gary W. Adamson, Apr 25 2004
EXTENSIONS
More terms from Robert G. Wilson v, Apr 28 2004
STATUS
approved