|
%I #3 Jun 15 2015 16:46:57
%S 1,-24,-169344,-25255286784,-23089632627769344,
%T -79051067969864491597824,-766667475511149432871084621824,
%U -17578325209217134578862801556544159744,-839197248407269659950832532302025663168118784
%N E.g.f.: S(x) = Series_Reversion( Integral 1/(1-x^5)^(1/5) dx ), where the constant of integration is zero.
%F Let e.g.f. C(x) = Sum_{n>=0} a(n)*x^(5*n)/(5*n)! and e.g.f. S(x) = Sum_{n>=0} a(n)*x^(5*n+1)/(5*n+1)!, then C(x) and S(x) satisfy:
%F (1) C(x)^5 + S(x)^5 = 1,
%F (2) S'(x) = C(x),
%F (3) C'(x) = -S(x)^4/C(x)^3,
%F (4) C(x)^4 * C'(x) + S(x)^4 * S'(x) = 0,
%F (5) S(x)/C(x) = Integral 1/C(x)^5 dx,
%F (6) S(x)/C(x) = Series_Reversion( Integral 1/(1+x^5) dx ) = Series_Reversion( Sum_{n>=0} (-1)^n * x^(5*n+1)/(5*n+1) ).
%e E.g.f. with offset 0 is C(x) and e.g.f. with offset 1 is S(x) where:
%e C(x) = 1 - 24*x^5/5! - 169344*x^10/10! - 25255286784*x^15/15! - 23089632627769344*x^20/20! +...
%e S(x) = x - 24*x^6/6! - 169344*x^11/11! - 25255286784*x^16/16! - 23089632627769344*x^21/21! +...
%e such that C(x)^5 + S(x)^5 = 1:
%e C(x)^5 = 1 - 120*x^5/5! + 604800*x^10/10! + 13208832000*x^15/15! +...
%e S(x)^5 = 120*x^5/5! - 604800*x^10/10! - 13208832000*x^15/15! -...
%e Related Expansions.
%e (1) The series reversion of S(x) is Integral 1/(1-x^5)^(1/5) dx:
%e Series_Reversion(S(x)) = x + 24*x^6/6! + 435456*x^11/11! + 115075344384*x^16/16! +...
%e 1/(1-x^5)^(1/5) = 1 + 24*x^5/5! + 435456*x^10/10! + 115075344384*x^15/15! +...
%e (2) d/dx S(x)/C(x) = 1/C(x)^5:
%e 1/C(x)^5 = 1 + 120*x^5/5! + 3024000*x^10/10! + 858574080000*x^15/15! +...
%e S(x)/C(x) = x + 120*x^6/6! + 3024000*x^11/11! + 858574080000*x^16/16! + 1226178516326400000*x^21/21! +...+ A258925(n)*x^(5*n+1)/(5*n+1)! +...
%e where
%e Series_Reversion(S(x)/C(x)) = x - 1/6*x^6 + 1/11*x^11 - 1/16*x^16 + 1/21*x^21 +...
%o (PARI) /* E.g.f. Series_Reversion(Integral 1/(1-x^5)^(1/5) dx): */
%o {a(n)=local(S=x); S = serreverse( intformal( 1/(1-x^5 +x*O(x^(5*n)))^(1/5) )); (5*n+1)!*polcoeff(S, 5*n+1)}
%o for(n=0, 15, print1(a(n), ", "))
%o (PARI) /* E.g.f. C(x) with offset 0: */
%o {a(n)=local(S=x, C=1+x); for(i=1, n, S=intformal(C +x*O(x^(5*n))); C=1-intformal(S^4/C^3 +x*O(x^(5*n))); ); (5*n)!*polcoeff(C, 5*n)}
%o for(n=0, 15, print1(a(n), ", "))
%o (PARI) /* E.g.f. S(x) with offset 1: */
%o {a(n)=local(S=x, C=1+x); for(i=1, n+1, S=intformal(C +x*O(x^(5*n+1))); C=1-intformal(S^4/C^3 +x*O(x^(5*n+1))); ); (5*n+1)!*polcoeff(S, 5*n+1)}
%o for(n=0, 15, print1(a(n), ", "))
%Y Cf. A258925 (S/C), A258878, A258900, A258927.
%K sign
%O 0,2
%A _Paul D. Hanna_, Jun 15 2015
| 