A052134
Numerators of coefficients in function a(x) such that a(a(a(x))) = sinh x.
Original entry on oeis.org
1, 1, -7, 643, -13583, 29957, -24277937, 6382646731, 2027394133729, -10948179003324221, 177623182156029053, -126604967848904128751, -2640658729595838040517543, 423778395125199663867841, 134802774190189008299452419971
Offset: 0
- W. C. Yang, Composition equations, preprint, 1999.
-
n = 15; m (* = maximal degree *)= 2n - 1; a[x_] = Sum[c[k]*x^k, {k, 1, m, 2}]; coes = DeleteCases[ CoefficientList[ Series[ a@a@a@x - Sinh[x], {x, 0, m}], x] // Rest, 0]; Do[s[k] = Solve[coes[[1]] == 0] // First; coes = coes /. s[k] // Rest, {k, 1, n}]; (CoefficientList[a[x] /. Flatten@Array[s, n], x] // Numerator // Partition[#, 2] &)[[All, 2]] (* Jean-François Alcover, Jan 04 2013 *)
A052132
Numerators of coefficients in function a(x) such that a(a(a(x))) = sin x.
Original entry on oeis.org
1, -1, -7, -643, -13583, -29957, -24277937, -6382646731, 2027394133729, 10948179003324221, 177623182156029053, 126604967848904128751, -2640658729595838040517543, -423778395125199663867841
Offset: 0
- W. C. Yang, Composition equations, preprint, 1999.
-
n = 14; m = 2 n - 1 (* m = maximal degree *); a[x_] = Sum[c[k] x^k, {k, 1, m, 2}] ; coes = DeleteCases[ CoefficientList[Series[a @ a @ a @ x - Sin[x], {x, 0, m}], x] // Rest , 0]; Do[s[k] = Solve[coes[[1]] == 0] // First; coes = coes /. s[k] // Rest, {k, 1, n}]
(CoefficientList[a[x] /. Flatten @ Array[s, n], x] // Numerator // Partition[#, 2] &)[[All, 2]]
(* Jean-François Alcover, May 04 2011 *)
-
T(n,m):=if n=m then 1 else 1/3*((((-1)^(n-m)+1)*sum((2*i-m)^n*binomial(m,i)*(-1)^((n+m)/2-i),i,0,m/2))/(2^m*n!)-sum(T(k,m)*sum(T(n,i)*T(i,k),i,k,n),k,m+1,n-1)-T(m,m)*sum(T(n,i)*T(i,m),i,m+1,n-1));
makelist(num(T(2*n-1,1)),n,1,7); /* Vladimir Kruchinin, Mar 10 2012 */
A072350
E.g.f. A(x) satisfies A(A(x)) = tan(x), where A(x) = Sum_{n>=1} a(n)*x^(2n-1)/(2n-1)!.
Original entry on oeis.org
1, 1, 3, 17, 225, 3613, -42997, 8725357, 2116966081, -549193907111, -114757574954509, 117893333517545097, 14433599120070484321, -65568697910890921624715, 2968238619232726100394235, 86999609037195113208781248165
Offset: 1
a(x) = x/1!+x^3/3!+3*x^5/5!+17*x^7/7!+225*x^9/9!+3613*x^11/11!-42997*x^13/13!+...
-
a[n_] := Module[{A, B, F}, F = Tan[x + O[x]^(2n+1)]; A = F; For[i = 0, i <= 2n-1, i++, B = InverseSeries[A, x]; A = (A + (B /. x -> F))/2]; If[n<1, 0, (2n-1)!*SeriesCoefficient[A, {x, 0, 2n-1}]]]; Table[a[n], {n, 1, 16}] (* Jean-François Alcover, Oct 29 2015, adapted from PARI *)
-
{a(n)=local(A,B,F);F=tan(x+O(x^(2*n+1)));A=F; for(i=0,2*n-1,B=serreverse(A);A=(A+subst(B,x,F))/2); if(n<1,0,(2*n-1)!*polcoeff(A,2*n-1,x))} \\ Paul D. Hanna, Dec 09 2004
Showing 1-3 of 3 results.
Comments