login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A020742
Pisot sequence T(7,9).
4
7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91, 93, 95, 97, 99, 101, 103, 105, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133, 135, 137, 139
OFFSET
0,1
FORMULA
a(n) = 2*n + 7.
a(n) = 2*a(n-1) - a(n-2).
From Elmo R. Oliveira, Oct 30 2024: (Start)
G.f.: (7 - 5*x)/(1 - x)^2.
E.g.f.: (7 + 2*x)*exp(x).
a(n) = A016825(n+3)/2 = A028560(n+1) - A028560(n). (End)
MATHEMATICA
T[x_, y_, z_] := Block[{a}, a[0] = x; a[1] = y; a[n_] := a[n] = Floor[a[n - 1]^2/a[n - 2]]; Table[a[n], {n, 0, z}]]; T[7, 9, 66] (* Michael De Vlieger, Aug 08 2016 *)
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]));
a
}
pisotT(50, 7, 9) \\ Colin Barker, Aug 08 2016
CROSSREFS
Subsequence of A005408, A020735. See A008776 for definitions of Pisot sequences.
Sequence in context: A295299 A066669 A059808 * A193071 A329112 A347886
KEYWORD
nonn,easy
STATUS
approved