A020734 Pisot sequence L(4,10).

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

Author

David W. Wilson

Keywords

• nonn

Links

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

Crossrefs

See A008776 for definitions of Pisot sequences.

Programs

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