cp's OEIS Frontend

A012494 Expansion of e.g.f. arctan(sin(x)) (odd powers only).

%I A012494 #51 Aug 17 2018 02:41:53
%S A012494 1,-3,45,-1743,125625,-14554683,2473184805,-579439207623,
%T A012494 179018972217585,-70518070842040563,34495620120141463965,
%U A012494 -20515677772241956573503,14578232896601243652363945,-12198268199871431840616166443,11871344562637111570703016357525
%N A012494 Expansion of e.g.f. arctan(sin(x)) (odd powers only).
%C A012494 arctan(sin(x)) = x - 3*x^3/3! + 45*x^5/5! - 1743*x^7/7! + 125625*x^9/9! + ....
%C A012494 Absolute values are coefficients in expansion of
%C A012494 arctanh(arcsinh(x)) = x + 3*x^3/3! + 45*x^5/5! + 1743*x^7/7! + ....
%C A012494 arccot(sin(x)) = Pi/2 - x + 3*x^3/3! - 45*x^5/5! + 1743*x^7/7! - ....
%H A012494 Vincenzo Librandi, <a href="/A012494/b012494.txt">Table of n, a(n) for n = 0..127</a>
%F A012494 a(n) = n!*sum(k=1..ceiling(n/2), (1+(-1)^(n-2*k+1))*2^(1-2*k)*sum(i=0..(2*k-1)/2, (-1)^((n+1)/2-i)*binomial(2*k-1,i)*(2*i-2*k+1)^n/n!)/(2*k-1)), n>0. _Vladimir Kruchinin_, Feb 25 2011
%F A012494 G.f.: cos(x) /(1 + sin^2(x)) = 1 - 3*x^2/2! + 45*x^4/4! - ... . - _Peter Bala_, Feb 06 2017
%F A012494 a(n) ~ (-1)^n * (2*n)! / (log(1+sqrt(2)))^(2*n+1). - _Vaclav Kotesovec_, Aug 17 2018
%p A012494 a:= n-> (t-> t!*coeff(series(arctan(sin(x)), x, t+1), x, t))(2*n+1):
%p A012494 seq(a(n), n=0..20);  # _Alois P. Heinz_, Aug 16 2018
%t A012494 Drop[ Range[0, 25]! CoefficientList[ Series[ ArcTan[ Sin[x]], {x, 0, 25}], x], {1, 25, 2}] (* Or *)
%t A012494 f[n_] := n!Sum[(1 + (-1)^(n - 2k + 1))2^(1 - 2k)Sum[(-1)^((n + 1)/2 - j)Binomial[2k - 1, j]((2j - 2k + 1)^n/n!)/(2k - 1), {j, 0, (2k - 1)/2}], {k, Ceiling[n/2]}]; Table[ f[n], {n, 1, 25, 2}] (* _Robert G. Wilson v_ *)
%o A012494 (Maxima) a(n):=n!*sum((1+(-1)^(n-2*k+1))*2^(1-2*k)*sum((-1)^((n+1)/2-i)*binomial(2*k-1,i)*(2*i-2*k+1)^n/n!,i,0,(2*k-1)/2)/(2*k-1),k,1,ceiling((n)/2)); /* _Vladimir Kruchinin_, Feb 25 2011 */
%o A012494 (Maxima) a(n):=sum(sum((2*i-2*k-1)^(2*n+1)*binomial(2*k+1,i)*(-1)^(n-i+1),i,0,k)/(4^k*(2*k+1)),k,0,n); /* _Vladimir Kruchinin_, Feb 04 2012 */
%Y A012494 Bisection of A003704, A013208.
%Y A012494 Cf. A101923, A001209, A000364, A000281, A156134, A002437.
%Y A012494 Cf. other sequences with a g.f. of the form cos(x)/(1 - k*sin^2(x)): A000364 (k=1), A001209 (k=1/2), A000281 (k=2), A156134 (k=3), A002437 (k=4).
%K A012494 sign,easy
%O A012494 0,2
%A A012494 Patrick Demichel (patrick.demichel(AT)hp.com)