OFFSET
1,2
FORMULA
G.f. A(x) also satisfies:
(1) A(x) = 1+x - cos(A(x) + x).
(2) A(x) = x + 2 * sin( (A(x) + x)/2 )^2.
(3) A(x) = -x + 2 * Series_Reversion(x - sin(x)^2).
a(n) = 2 * A143134(n) for n>1.
EXAMPLE
E.g.f.: A(x) = x + 4*x^2/2! + 24*x^3/3! + 224*x^4/4! + 2880*x^5/5! + 47104*x^6/6! + 935424*x^7/7! + 21853184*x^8/8! + 587089920*x^9/9! + 17829167104*x^10/10! +...
such that A(x - sin(x)^2) = x + sin(x)^2.
RELATED SERIES.
A(x - sin(x)^2) = x + 2*x^2/2! - 8*x^4/4! + 32*x^6/6! - 128*x^8/8! + 512*x^10/10! - 2048*x^12/12! +...
which equals x + sin(x)^2.
cos(A(x) + x) = 1 - 4*x^2/2! - 24*x^3/3! - 224*x^4/4! - 2880*x^5/5! - 47104*x^6/6! +...
which equals 1+x - A(x).
PROG
(PARI) {a(n) = my(A=x); for(i=1, 21, A = subst(x + sin(x +x*O(x^n) )^2, x, serreverse(x - sin(x +x*O(x^n) )^2))); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {a(n) = my(A=x); for(i=0, n, A = 1+x - cos(A + x +x*O(x^n))); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 06 2016
STATUS
approved