login
This site is supported by donations to The OEIS Foundation.
Logo

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A154350 a(n) = ((9+2*sqrt(2))^n-(9-2*sqrt(2))^n)/(4*sqrt(2)). 0
1, 18, 251, 3204, 39349, 474390, 5666543, 67367304, 798953833, 9463355802, 112016774627, 1325476969740, 15681360907549, 185504677544862, 2194344849556439, 25956365831240976, 307027410944717521 (list; graph; refs; listen; history; internal format)
OFFSET

1,2

COMMENTS

lim_{n -> infinity} a(n)/a(n-1) = 9+2*sqrt(2) = 11.8284271247....

FORMULA

G.f.: x/(1-18*x+73*x^2). [From Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Jan 12 2009, corrected Oct 08 2009]

a(n) = 18*a(n-1)-73*a(n-2) for n>1; a(0)=0, a(1)=1. [From Philippe DELEHAM (kolotoko(AT)wanadoo.fr), Jan 12 2009]

MATHEMATICA

Join[{a=1, b=18}, Table[c=18*b-73*a; a=b; b=c, {n, 40}]] (*From Vladimir Joseph Stephan Orlovsky, Feb 09 2011*)

PROG

(MAGMA) Z<x>:=PolynomialRing(Integers()); N<r>:=NumberField(x^2-2); S:=[ ((9+2*r)^n-(9-2*r)^n)/(4*r): n in [1..17] ]; [ Integers()!S[j]: j in [1..#S] ]; [From Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Jan 12 2009]

CROSSREFS

Cf. A002193 (decimal expansion of sqrt(2)).

Sequence in context: A153886 A154241 A154250 * A001722 A060788 A144708

Adjacent sequences:  A154347 A154348 A154349 * A154351 A154352 A154353

KEYWORD

nonn

AUTHOR

Al Hakanson (hawkuu(AT)gmail.com), Jan 07 2009

EXTENSIONS

Extended beyond a(7) by Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Jan 12 2009

Edited by Klaus Brockhaus (klaus-brockhaus(AT)t-online.de), Oct 08 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 15 20:24 EST 2012. Contains 205852 sequences.