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 A087798 #33 Sep 08 2022 08:45:11 %S A087798 2,9,83,756,6887,62739,571538,5206581,47430767,432083484,3936182123, %T A087798 35857722591,326655685442,2975758891569,27108485709563, %U A087798 246952130277636,2249677658208287,20494051054152219,186696137145578258 %N A087798 a(n) = 9*a(n-1) + a(n-2), starting with a(0) = 2 and a(1) = 9. %C A087798 a(n+1)/a(n) converges to (9 + sqrt(85))/2. %C A087798 For more information about this type of recurrence follow the Khovanova link and see A054413 and A086902. - _Johannes W. Meijer_, Jun 12 2010 %H A087798 Vincenzo Librandi, <a href="/A087798/b087798.txt">Table of n, a(n) for n = 0..1000</a> %H A087798 Tanya Khovanova, <a href="http://www.tanyakhovanova.com/RecursiveSequences/RecursiveSequences.html">Recursive Sequences</a> %H A087798 <a href="/index/Rea#recur1">Index entries for recurrences a(n) = k*a(n - 1) +/- a(n - 2)</a> %H A087798 <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (9,1). %F A087798 a(n) = ((9 + sqrt(85))/2)^n + ((9 - sqrt(85))/2)^n. %F A087798 G.f.: (2 - 9*x)/(1 - 9*x - x^2). - _Philippe Deléham_, Nov 02 2008 %F A087798 From _Johannes W. Meijer_, Jun 12 2010: (Start) %F A087798 a(2n+1) = 9*A097840(n), a(2n) = A099373(n). %F A087798 a(3n+1) = A041150(5n), a(3n+2) = A041150(5n+3), a(3n+3) = 2*A041150(5n+4). %F A087798 Lim_{k->infinity} a(n+k)/a(k) = (A087798(n) + A099371(n)*sqrt(85))/2. %F A087798 Lim_{n->infinity} A087798(n)/A099371(n) = sqrt(85). (End) %e A087798 a(4) = 9*a(3) + a(2) = 9*756 + 83 = 6887. %t A087798 RecurrenceTable[{a[0] == 2, a[1] == 9, a[n] == 9 a[n-1] + a[n-2]}, a, {n, 30}] (* _Vincenzo Librandi_, Sep 19 2016 *) %t A087798 LinearRecurrence[{9,1}, {2,9}, 30] (* _G. C. Greubel_, Nov 07 2018 *) %o A087798 (Magma) I:=[2,9]; [n le 2 select I[n] else 9*Self(n-1)+Self(n-2): n in [1..30]]; // _Vincenzo Librandi_, Sep 19 2016 %o A087798 (PARI) x='x+O('x^30); Vec((2-9*x)/(1-9*x-x^2)) \\ _G. C. Greubel_, Nov 07 2018 %Y A087798 Cf. A014511. %K A087798 easy,nonn %O A087798 0,1 %A A087798 Nikolay V. Kosinov, Dmitry V. Poljakov (kosinov(AT)unitron.com.ua), Oct 10 2003 %E A087798 More terms from _Ray Chandler_, Nov 06 2003