%I #14 Aug 28 2018 12:45:18
%S 1,2,6,40,360,4592,70896,1279360,26497920,619457792,16166151936,
%T 466022394880,14708199367680,504453778491392,18681868054910976,
%U 742996971891097600,31583887537425776640,1429076863804079931392,68575244394079858262016,3478457209493103235563520,185971933231431479545036800
%N E.g.f. A(x) satisfies: cos(A(x)) + sin(A(x)) = 1/( cos(A(-x)) + sin(A(-x)) ).
%H Paul D. Hanna, <a href="/A318006/b318006.txt">Table of n, a(n) for n = 1..300</a>
%F E.g.f. A(x) satisfies:
%F (1) A(-A(-x)) = x.
%F (2a) 1 = Sum_{n>=0} (-1)^floor(n/2) * ( A(x) - (-1)^n*A(-x) )^n/n!.
%F (2b) 1 = Sum_{n>=0} (-1)^floor(n/2) * ( A(A(x)) + (-1)^n*x )^n/n!.
%F (3a) 1 = cos( A(x) - A(-x) ) + sin( A(x) + A(-x) ).
%F (3b) 1 = ( cos(A(x)) + sin(A(x)) ) * ( cos(A(-x)) + sin(A(-x)) ).
%F (4a) 1 = cos(A(A(x)) + x) + sin(A(A(x)) - x).
%F (4b) 1 = ( cos(A(A(x))) + sin(A(A(x))) ) * (cos(x) - sin(x)).
%F (5a) A(x) = arcsin( 2*sin(2*x)/(2 - sin(2*x)) )/2.
%F (5b) A(A(x)) = arcsin( sin(2*x)/(1 - sin(2*x)) )/2, which is the e.g.f. of A318005.
%F (6) cos(A(x)) + sin(A(x)) = sqrt( (2 + sin(2*x))/(2 - sin(2*x)) ).
%F a(n) ~ sqrt(3) * 5^(1/4) * 2^(n-2) * n^(n-1) / (exp(n) * (arcsin(2/3))^(n - 1/2)). - _Vaclav Kotesovec_, Aug 28 2018
%e E.g.f.: A(x) = x + 2*x^2/2! + 6*x^3/3! + 40*x^4/4! + 360*x^5/5! + 4592*x^6/6! + 70896*x^7/7! + 1279360*x^8/8! + 26497920*x^9/9! + 619457792*x^10/10! + ...
%e such that
%e cos(A(x)) + sin(A(x)) = 1/( cos(A(-x)) + sin(A(-x)) )
%e and
%e sin(2*A(x)) = 2*sin(2*x)/(2 - sin(2*x)) = 2*sin(A(x))*cos(A(x)).
%e RELATED SERIES.
%e (a) A(A(x)) = x + 4*x^2/2! + 24*x^3/3! + 224*x^4/4! + 2880*x^5/5! + 48064*x^6/6! + 989184*x^7/7! + 24218624*x^8/8! + ... + A318005(n)*x^n/n! + ...
%e where A(A(x)) = arcsin( 2*sin(2*x)/(2 - sin(2*x)) ) /2.
%e (b) cos(A(x)) + sin(A(x)) = 1/(cos(A(-x)) - sin(A(-x))) = 1 + x + x^2/2! - x^3/3! - 7*x^4/4! - 59*x^5/5! - 239*x^6/6! - 421*x^7/7! + 4913*x^8/8! + 108361*x^9/9! + 1000321*x^10/10! + ...
%e where cos(A(x)) + sin(A(x)) = sqrt( (2 + sin(2*x))/(2 - sin(2*x)) ).
%t nmax = 25; Rest[CoefficientList[Series[ArcSin[4/(2 - Sin[2*x]) - 2]/2, {x, 0, nmax}], x] * Range[0, nmax]!] (* _Vaclav Kotesovec_, Aug 28 2018 *)
%o (PARI) /* A(x) = arcsin( 2*sin(2*x)/(2 - sin(2*x)) )/2 */
%o {a(n) = my(A = asin( 2*sin(2*x +x*O(x^n))/(2 - sin(2*x +x*O(x^n))) )/2 ); n!*polcoeff(A,n)}
%o for(n=1,25, print1(a(n),", "))
%o (PARI) /* From 1 = cos(A(A(x)) + x) + sin(A(A(x)) - x) */
%o {a(n) = my(A=[1,1]); for(i=1, n, A = concat(A,0);
%o A[#A] = -Vec(cos(subst(x*Ser(A),x,x*Ser(A)) + x) + sin(subst(x*Ser(A),x,x*Ser(A)) - x))[#A+1]/2; ); n!*A[n]}
%o for(n=1, 25, print1(a(n), ", "))
%Y Cf. A318005, A318001, A200560.
%K nonn
%O 1,2
%A _Paul D. Hanna_, Aug 27 2018