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 A161007 #25 Jan 31 2023 06:00:16
%S A161007 0,1,2,20,72,464,2080,11584,56448,298240,1499648,7771136,39536640,
%T A161007 203411456,1039409152,5333401600,27297349632,139929124864,
%U A161007 716615843840,3672097685504,18810048872448,96373660712960,493708103385088,2529394778177536,12958119210516480
%N A161007 a(n+1) = 2*a(n) + 16*a(n-1), a(0)=0, a(1)=1.
%H A161007 Colin Barker, <a href="/A161007/b161007.txt">Table of n, a(n) for n = 0..1000</a>
%H A161007 J. Borowska, L. Lacinska, <a href="https://doi.org/10.17512/jamcm.2014.4.03">Recurrence form of determinant of a heptadiagonal symmetric Toeplitz matrix</a>, J. Appl. Math. Comp. Mech. 13 (2014) 19-16, remark 2 for permanent of tridiagonal Toeplitz matrices a=2, b=4.
%H A161007 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (2,16).
%F A161007 a(n) = ((1+sqrt(17))^n - (1-sqrt(17))^n)/(2*sqrt(17)).
%F A161007 Limit_{n -> oo} a(n+1)/a(n) = 1 + sqrt(17).
%F A161007 G.f.: x / (1 - 2*x - 16*x^2). - _Colin Barker_, Jul 01 2015
%F A161007 a(n) = 2^(n-1)*A006131(n-1). - _R. J. Mathar_, Mar 08 2021
%F A161007 a(n) = (4*i)^n*ChebyshevU(n, -i/4). - _G. C. Greubel_, Oct 15 2022
%t A161007 LinearRecurrence[{2, 16}, {0, 1}, 50] (* _T. D. Noe_, Nov 07 2011 *)
%o A161007 (PARI) concat(0, Vec(-x/(16*x^2+2*x-1) + O(x^40))) \\ _Colin Barker_, Jul 01 2015
%o A161007 (Magma) [n le 2 select n-1 else 2*(Self(n-1) +8*Self(n-2)): n in [1..41]]; // _G. C. Greubel_, Oct 15 2022
%o A161007 (SageMath)
%o A161007 A161007=BinaryRecurrenceSequence(2,16,0,1)
%o A161007 [A161007(n) for n in range(41)] # _G. C. Greubel_, Oct 15 2022
%Y A161007 Cf. A006131, A010473 (sqrt(17)).
%K A161007 nonn,easy
%O A161007 0,3
%A A161007 _Sture Sjöstedt_, Jun 02 2009