A097719 - OEIS

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

A097719

a(n)=(a^n-b^n)/(a-b), where a=1.3802775690976141157... and b=-0.8191725133961644397... are the real roots of x^4-x^3-1=0.

0, 1, 2, 1, 2, 3, 6, 7, 9, 13, 20, 27, 37, 51, 71, 98, 136, 187, 258, 357, 494, 681, 940, 1298, 1792, 2474, 3415, 4714, 6506, 8981, 12396, 17110, 23617, 32599, 44996, 62107, 85725, 118324, 163320, 225427, 311153, 429477, 592798, 818226

(list; graph; refs; listen; history; text; internal format)

OFFSET

0,3

LINKS

Table of n, a(n) for n=0..43.

MATHEMATICA

NSolve[x^4 - x^3 - 1 == 0, x] b01 = 1.380277569097614 b11 = -0.8191725133961645 b21 = 0.2194474721492751 - 0.9144736629677264*I b31 = 0.2194474721492751 + 0.9144736629677264*I f[n_] := N[(b01^n - b11^n - b21^n - b31^n)/(b01 - b11 - b21 - b31)]; digits = 200 c = Table[Floor[Re[f[n]]], {n, 1, digits}]

CROSSREFS

Cf. A001644.

Sequence in context: A301336 A128474 A108618 * A249050 A341649 A056493

Adjacent sequences: A097716 A097717 A097718 * A097720 A097721 A097722

KEYWORD

nonn

AUTHOR

Roger L. Bagula, Sep 21 2004; corrected Oct 06 2004

STATUS

approved