|
|
A014003
|
|
Pisot sequence E(9,15), a(n) = floor( a(n-1)^2/a(n-2) + 1/2 ).
|
|
1
|
|
|
9, 15, 25, 42, 71, 120, 203, 343, 580, 981, 1659, 2806, 4746, 8027, 13576, 22961, 38834, 65680, 111085, 187879, 317761, 537431, 908960, 1537329, 2600093, 4397552, 7437605, 12579264, 21275381, 35983173, 60858545, 102930403, 174086776, 294433954, 497977820
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
LINKS
|
|
|
FORMULA
|
Known not to satisfy any linear recurrence.
|
|
MATHEMATICA
|
nxt[{a_, b_}]:={b, Floor[b^2/a+1/2]}; NestList[nxt, {9, 15}, 40][[All, 1]] (* Harvey P. Dale, Jan 31 2023 *)
|
|
PROG
|
(PARI) pisotE(nmax, a1, a2) = {
a=vector(nmax); a[1]=a1; a[2]=a2;
for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]+1/2));
a
}
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|