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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A117648 Simplest bounching ball system at: f[n]=(4/3)*f[n-1]-f[n-2]. 0
2, 3, 2, 1, 3, 3, 2, 0, 2, 2, 0, 2, 3, 3, 1, 2, 2, 1, 1, 3, 3, 2, 1, 2, 2, 0, 2, 3, 3, 1, 2, 2, 1, 1, 3, 3, 2, 1, 2, 2, 0, 2, 3, 3, 0, 2, 2, 1, 1, 3, 3 (list; graph; refs; listen; history; internal format)
OFFSET

0,1

COMMENTS

There is a large family of this type of linear harmonic equation. General equation form is: (x - (a + I*Sqrt[2*a + 1])/(a + 1))*(x - (a - I*Sqrt[2*a + 1])/(a + 1)) Rational numbers are Coefficients: Table[Coefficient[ExpandAll[(x - (a + I*Sqrt[2*a + 1])/( a + 1))*(x - (a - I*Sqrt[2*a + 1])/(a + 1))], x], {a, 1, 285}] Binet solution for this one is: f[n]= 2*Cos[n*ArcTan[Sqrt[5]/2]]+Sqrt[5]*Sin[n*ArcTan[Sqrt[5]/2]]

FORMULA

f[0]=2;f[1]=3; f[n]=(4/3)*f[n-1]-f[n-2] a(n) = Abs[Floor[f[n]]]

MATHEMATICA

Clear[f, M, v] f[0] = 2; f[1] = 3; f[n_] := f[n] = (4/3)*f[n - 1] - f[n - 2] Table[Abs[Floor[f[n]]], {n, 0, 50}] ListPlot[%, PlotJoined -> True] M = {{0, 1}, {-1, (4/3)}}; v[0] = {2, 3}; v[n_] := v[n] = M.v[n - 1] Table[Abs[Floor[v[n][[1]]]], {n, 0, 80}] ListPlot[%]

CROSSREFS

Sequence in context: A118105 A125211 A139367 * A037222 A102288 A107357

Adjacent sequences:  A117645 A117646 A117647 * A117649 A117650 A117651

KEYWORD

nonn,uned

AUTHOR

Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Apr 10 2006

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 17 12:38 EST 2012. Contains 206021 sequences.