login
E.g.f. A(x) satisfies: cos(A(x)) + sin(A(x)) = 1/( cos(A(-x)) + sin(A(-x)) ).
3

%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