A020747 - OEIS

A020747

Pisot sequence T(4,6), a(n) = floor(a(n-1)^2/a(n-2)).

4, 6, 9, 13, 18, 24, 32, 42, 55, 72, 94, 122, 158, 204, 263, 339, 436, 560, 719, 923, 1184, 1518, 1946, 2494, 3196, 4095, 5246, 6720, 8608, 11026, 14123, 18089, 23168, 29673, 38004, 48674, 62339, 79840, 102254, 130960, 167724, 214808, 275109, 352337, 451244, 577915

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

OFFSET

0,1

LINKS

Colin Barker, Table of n, a(n) for n = 0..1000

PROG

(PARI) pisotT(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]));

}

pisotT(50, 4, 6) \\ Colin Barker, Jul 29 2016

CROSSREFS