OFFSET
0,1
LINKS
FORMULA
G.f.: (3x^5+2x^4+x^3+4x^2-x+6)/(-x^6-x^3+x^2-2x+1) (conjectured). - Ralf Stephan, May 12 2004
MATHEMATICA
RecurrenceTable[{a[n] == Ceiling[a[n - 1]^2/a[n - 2] - 1/2], a[0] == 6, a[1] == 11}, a, {n, 0, 31}] (* or *)
First@ Transpose[NestList[{#2, Round[#2^2/#1]} & @@ # &, {6, 11}, 31]] (* Michael De Vlieger, Aug 08 2016, after Harvey P. Dale at A021008 *)
PROG
(PARI) pisotP(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]-1/2));
a
}
pisotP(50, 6, 11) \\ Colin Barker, Aug 08 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved