|
%I #16 Oct 27 2023 21:05:23
%S 1,1,1,0,-1,-2,-2,0,3,6,6,2,-6,-14,-16,-8,11,32,42,26,-18,-74,-108,
%T -82,18,162,268,238,16,-344,-656,-664,-189,694,1570,1792,826,-1294,
%U -3668,-4698,-2866,2110,8364,12034,8960,-2432,-18508,-30134,-26254,-910,39492,73862,73560,19120
%N G.f. A(x) satisfies A(x) = 1 + x * A(x) / A(x^2).
%H Alois P. Heinz, <a href="/A218033/b218033.txt">Table of n, a(n) for n = 0..1000</a> (terms up to n=200 from Vincenzo Librandi)
%t m = 60; A[_] = 1;
%t Do[A[x_] = 1 + x A[x]/A[x^2] + O[x]^m // Normal, {m}];
%t CoefficientList[A[x], x] (* _Jean-François Alcover_, Oct 04 2019 *)
%o (PARI)
%o N=66; R=O('x^N); x='x+R;
%o F = 1 + x;
%o for (k=1,N+1, F = 1 + x * F / subst(F,'x,'x^2) + R; );
%o Vec(F)
%K sign
%O 0,6
%A _Joerg Arndt_, Oct 19 2012
|
