cp's OEIS Frontend

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.

A095003 a(n) = 9*a(n-1) - 9*a(n-2) + a(n-3).

This page as a plain text file.
%I A095003 #17 Nov 12 2022 18:24:59
%S A095003 1,6,45,352,2769,21798,171613,1351104,10637217,83746630,659335821,
%T A095003 5190939936,40868183665,321754529382,2533168051389,19943589881728,
%U A095003 157015551002433,1236180818137734,9732430994099437,76623267134657760,603253706083162641,4749406381530643366
%N A095003 a(n) = 9*a(n-1) - 9*a(n-2) + a(n-3).
%C A095003 a(n)/a(n-1) tends to 7.87298... = 4 + sqrt(15) = C (having the property that C + 1/C = 8). Eigenvalues of M are C, 1/C, 1; being roots of x^3 - 9x^2 + 9x - 1.
%H A095003 <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (9,-9,1).
%F A095003 a(n+3) = 9*a(n+2) - 9*a(n+1) + a(n); given a(1) = 1, a(2) = 6, a(3) = 45.
%F A095003 Let M be the 3 X 3 matrix [1 1 1 / 1 2 3 / 1 3 6]. M^n * [1 0 0] = [A095002(n) a(n) A095004(n)].
%e A095003 a(4) = 352 since M^4 * [1 0 0] = [145, 352, 640].
%p A095003 a:= n-> (<<1|1|1>, <1|2|3>, <1|3|6>>^n)[1, 2]:
%p A095003 seq(a(n), n=1..23);  # _Alois P. Heinz_, Jun 06 2021
%t A095003 a[n_] := (MatrixPower[{{1, 1, 1}, {1, 2, 3}, {1, 3, 6}}, n].{{1}, {0}, {0}})[[2, 1]]; Table[ a[n], {n, 20}]; (* _Robert G. Wilson v_, May 29 2004 *)
%t A095003 LinearRecurrence[{9,-9,1},{1,6,45},30] (* _Harvey P. Dale_, Nov 12 2022 *)
%Y A095003 Cf. A095002, A095004, A076765.
%K A095003 nonn
%O A095003 1,2
%A A095003 _Gary W. Adamson_, May 27 2004
%E A095003 Edited and extended by _Robert G. Wilson v_, May 29 2004
%E A095003 Definition corrected and edited by _Georg Fischer_, Jun 06 2021