[go: up one dir, main page]

login
A020734
Pisot sequence L(4,10).
2
4, 10, 25, 63, 159, 402, 1017, 2573, 6510, 16472, 41679, 105461, 266850, 675216, 1708513, 4323086, 10938795, 27678662, 70035898, 177213300, 448406526, 1134612428, 2870933600, 7264383443, 18381221638, 46510390257, 117686215012, 297783895756, 753488830979
OFFSET
0,1
LINKS
PROG
(PARI) pisotL(nmax, a1, a2) = {
a=vector(nmax); a[1]=a1; a[2]=a2;
for(n=3, nmax, a[n] = ceil(a[n-1]^2/a[n-2]));
a
}
pisotL(50, 4, 10) \\ Colin Barker, Aug 07 2016
CROSSREFS
See A008776 for definitions of Pisot sequences.
Sequence in context: A020748 A021004 A020709 * A225301 A298412 A289245
KEYWORD
nonn
STATUS
approved