|
Search: id:A091141
|
|
|
| A091141 |
|
a(n) = 2a(n-1) + 4a(n-2) - 2a(n-3). |
|
+0 4
|
|
| 1, 4, 13, 40, 124, 382, 1180, 3640, 11236, 34672, 107008, 330232, 1019152, 3145216, 9706576, 29955712, 92447296, 285304288, 880486336, 2717295232, 8385927232, 25880062720, 79869243904, 246486884224, 760690618624, 2347590286336
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
COMMENT
|
One of 3 related sequences generated from finite difference operations. Let r(1)=s(1)=t(1)=1. Given r(n), s(n), and t(n), let f(x) = r(n) x^2 + s(n) x + t(n), and let r(n+1), s(n+1), and t(n+1) be the 0-th, 1-st, and 2-nd differences of f(x) at x=1. I.e. r(n+1) = f(1) = r(n)+s(n)+t(n), s(n+1) = f(2)-f(1) = 3r(n)+s(n), and t(n+1) = f(3)-2f(2)+f(1) = 2r(n). This sequence gives s(n).
|
|
FORMULA
|
Let v(n) be the column vector with elements r(n), s(n), t(n); then v(n) = [1 1 1 / 3 1 0 / 2 0 0] v(n-1).
The limit as n->infinity of a(n+1)/a(n) is the largest root of x^3 - 2x^2 - 4x + 2 = 0, which is about 3.086130197651494.
|
|
MATHEMATICA
|
a[n_] := (MatrixPower[{{1, 1, 1}, {3, 1, 0}, {2, 0, 0}}, n-1].{{1}, {1}, {1}})[[2, 1]]
|
|
CROSSREFS
|
Cf. r(n) = A091140(n), t(n) = A091142(n).
Sequence in context: A025567 A076040 A003462 this_sequence A098183 A094628 A034742
Adjacent sequences: A091138 A091139 A091140 this_sequence A091142 A091143 A091144
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Gary W. Adamson (qntmpkt(AT)yahoo.com), Dec 21 2003
|
|
|
Search completed in 0.002 seconds
|