cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A048605 Numerators of coefficients in function a(x) such that a(a(x)) = arctan(x).

This page as a plain text file.
%I A048605 #31 Jun 02 2024 11:06:31
%S A048605 1,-1,7,-43,4489,-49897,20130311,-319053131,329796121169,
%T A048605 -62717244921977,14635852695795623,-33233512260583073,
%U A048605 149490010959849868177,-3562767949848393597053
%N A048605 Numerators of coefficients in function a(x) such that a(a(x)) = arctan(x).
%C A048605 A recursion exists for coefficients, but is too complicated to use without a computer algebra system.
%D A048605 W. C. Yang, Polynomials are essentially integer partitions, preprint, 1999
%D A048605 W. C. Yang, Composition equations, preprint, 1999
%H A048605 Dmitry Kruchinin, Vladimir Kruchinin, <a href="http://arxiv.org/abs/1302.1986">Method for solving an iterative functional equation A^{2^n}(x) = F(x)</a>, arXiv:1302.1986
%H A048605 W. C. Yang, <a href="http://dx.doi.org/10.1016/S0012-365X(99)00412-4">Derivatives are essentially integer partitions</a>, Discrete Math., 222 (2000), 235-245.
%F A048605 a(n) = numerator(T(2*n-1,1)), T(n,m)=1/2*(2^(-m-1)*m!*((-1)^(n+m)+1)*(-1)^((3*n+m)/2)*sum(i=m..n, (2^i*stirling1(i,m)*binomial(n-1,i-1))/i!)-sum(i=m+1..n-1, T(n,i)*T(i,m))), n>m, T(n,n)=1. - _Vladimir Kruchinin_, Mar 12 2012
%e A048605 x - x^3/6 + x^5 * 7/120 + ...
%t A048605 n = 28; a[x_] = Sum[c[k] k! x^k, {k, 1, n, 2}];
%t A048605 sa = Series[a[x], {x, 0, n}];
%t A048605 coes = CoefficientList[ComposeSeries[sa, sa] - Series[ArcTan[x], {x, 0, n}], x] // Rest;
%t A048605 eq = Reduce[((# == 0) & /@ coes)]; Table[c[k] k!, {k, 1, n, 2}] /. First[Solve[eq]] // Numerator
%t A048605 (* _Jean-François Alcover_, Apr 26 2011 *)
%o A048605 (Maxima)
%o A048605 T(n, m):=if n=m then 1 else 1/2*(2^(-m-1)*m!*((-1)^(n+m)+1)*(-1)^((3*n+m)/2)*sum((2^i*stirling1(i, m)*binomial(n-1, i-1))/i!, i, m, n)-sum(T(n, i)*T(i, m), i, m+1, n-1));
%o A048605 makelist(num(T(2*n-1, 1), n, 1, 5)); /* _Vladimir Kruchinin_, Mar 12 2012 */
%Y A048605 Cf. A048604, A095885.
%K A048605 frac,sign,nice
%O A048605 0,3
%A A048605 Winston C. Yang (yang(AT)math.wisc.edu)