OFFSET
0,2
COMMENTS
Generating function: AQM(1,1-8x) where AQM(u,v) (arithmetic-quadratic mean of u and v) is the fixed point obtained by iterating ((u+v)/2, sqrt((u^2+v^2)/2)) (we choose 1-8x in order to avoid denominators, as in A060691).
EXAMPLE
First steps of iteration of ((u+v)/2, sqrt((u^2+v^2)/2)) are (1, 1-8x), (1 - 4*x, 1 - 4*x + 8*x^2 + 32*x^3 + 96*x^4 + O(x^5)), then (1 - 4*x + 4*x^2 + 16*x^3 + 48*x^4 + O(x^5), 1 - 4*x + 4*x^2 + 16*x^3 + 56*x^4 + O(x^5)) and (1 - 4*x + 4*x^2 + 16*x^3 + 52*x^4 + O(x^5), 1 - 4*x + 4*x^2 + 16*x^3 + 52*x^4 + O(x^5)), so the first terms of this sequence are 1, -4, 4, 16, 52.
PROG
(Sage)
R.<x> = PowerSeriesRing(QQ, default_prec=50)
(a, b) = (1, 1-8*x)
for i in range(50):
(a, b) = ((a+b)/2, sqrt((a^2+b^2)/2))
a.coefficients()
(PARI) seq(n)={my(p=1, q=1-8*x+O(x*x^n)); while(p!=q, my(t=p+q); q = sqrt((p^2 + q^2)/2); p=t/2); Vec(p)} \\ Andrew Howroyd, Mar 22 2021
CROSSREFS
KEYWORD
sign
AUTHOR
David A. Madore, Mar 22 2021
STATUS
approved