login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A015583 a(0) = 0, a(1) = 1; for n >= 2, a(n) = 9*a(n-1) + 7*a(n-2). 3
0, 1, 9, 88, 855, 8311, 80784, 785233, 7632585, 74189896, 721137159, 7009563703, 68134033440, 662273246881, 6437397456009, 62572489832248, 608214190682295, 5911935144966391, 57464915639473584, 558567786770026993, 5429364490406558025, 52774254921049211176 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,3
LINKS
FORMULA
G.f.: 1/(1 - 9*x - 7*x^2). - Zerinvary Lajos, Dec 20 2008
MAPLE
a:=proc(n) option remember: if n=0 then 0 elif n=1 then 1 else 9*procname(n-1)+7*procname(n-2) fi: end; seq(a(n), n=0..22); # Muniru A Asiru, Jul 15 2018
MATHEMATICA
Join[{a=0, b=1}, Table[c=9*b+7*a; a=b; b=c, {n, 60}]] (* Vladimir Joseph Stephan Orlovsky, Jan 27 2011 *)
LinearRecurrence[{9, 7}, {0, 1}, 30] (* Vincenzo Librandi, Nov 15 2012 *)
PROG
(Sage) [lucas_number1(n, 9, -7) for n in range(0, 19)] # Zerinvary Lajos, Apr 26 2009
(Magma) [n le 2 select n-1 else 9*Self(n-1) + 7*Self(n-2): n in [1..30]]; // Vincenzo Librandi, Nov 15 2012
(PARI) x='x+O('x^30); concat([0], Vec(1/(1-9*x-7*x^2))) \\ G. C. Greubel, Jan 06 2018
(GAP) a:=[0, 1];; for n in [3..22] do a[n]:=9*a[n-1]+7*a[n-2]; od; a; # Muniru A Asiru, Jul 15 2018
CROSSREFS
Sequence in context: A245491 A160466 A362509 * A152266 A260041 A084022
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Extended by T. D. Noe, May 23 2011
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 03:51 EDT 2024. Contains 371264 sequences. (Running on oeis4.)