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

%I #10 Jul 29 2016 22:03:19

%S 4,6,9,13,18,24,32,42,55,72,94,122,158,204,263,339,436,560,719,923,

%T 1184,1518,1946,2494,3196,4095,5246,6720,8608,11026,14123,18089,23168,

%U 29673,38004,48674,62339,79840,102254,130960,167724,214808,275109,352337,451244,577915

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

%H Colin Barker, <a href="/A020747/b020747.txt">Table of n, a(n) for n = 0..1000</a>

%o (PARI) pisotT(nmax, a1, a2) = {

%o a=vector(nmax); a[1]=a1; a[2]=a2;

%o for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]));

%o a

%o }

%o pisotT(50, 4, 6) \\ _Colin Barker_, Jul 29 2016

%Y See A008776 for definitions of Pisot sequences.

%K nonn

%O 0,1

%A _David W. Wilson_