A117655 a(n) = floor(f(n)), where f(n) = (5 - 4.37/(f(n-1) + f(n-2))^(1/3))^3, f(0) = 0, and f(1) = 1.

0, 1, 0, 0, 0, 0, 1, 3, 12, 34, 54, 64, 69, 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72

Offset

0,8

Links

G. C. Greubel, Table of n, a(n) for n = 0..1000

Formula

a(n) = floor(f(n)), where f(n) = (5 - 4.37/(f(n-1) + f(n-2))^(1/3))^3, f(0) = 0, and f(1) = 1.

Mathematica

f[n_]:= f[n]= If[n<2, n, (5 -4.37/(f[n-1]+f[n-2])^(1/3))^3];
a[n_]:= Floor[f[n]];
Table[a[n], {n, 0, 100}]

Prog

(Magma)
f:=[n le 2 select n-1 else (5 - 4.37/(Self(n-1) + Self(n-2))^(1/3))^3 : n in [1..101]];
A117655:= func< n | Floor(f[n+1]) >;
[A117655(n): n in [0..100]]; // G. C. Greubel, Mar 30 2024
(SageMath)
@CachedFunction
def f(n): return n if n<2 else (5 - 4.37/(f(n-1)+f(n-2))^(1/3))^3
def A117655(n): return floor(f(n))
[A117655(n) for n in range(101)] # G. C. Greubel, Mar 30 2024

Crossrefs

Sequence in context: A054610 A183468 A196234 * A081423 A184705 A257890

Adjacent sequences: A117652 A117653 A117654 * A117656 A117657 A117658

Keyword

nonn,less

Author

Roger L. Bagula, Apr 11 2006

Extensions

Edited by G. C. Greubel, Mar 30 2024

Status

approved