%I #10 Aug 07 2016 09:28:29
%S 4,10,25,63,159,402,1017,2573,6510,16472,41679,105461,266850,675216,
%T 1708513,4323086,10938795,27678662,70035898,177213300,448406526,
%U 1134612428,2870933600,7264383443,18381221638,46510390257,117686215012,297783895756,753488830979
%N Pisot sequence L(4,10).
%H Colin Barker, <a href="/A020734/b020734.txt">Table of n, a(n) for n = 0..1000</a>
%o (PARI) pisotL(nmax, a1, a2) = {
%o a=vector(nmax); a[1]=a1; a[2]=a2;
%o for(n=3, nmax, a[n] = ceil(a[n-1]^2/a[n-2]));
%o a
%o }
%o pisotL(50, 4, 10) \\ _Colin Barker_, Aug 07 2016
%Y See A008776 for definitions of Pisot sequences.
%K nonn
%O 0,1
%A _David W. Wilson_