[seqfan] Re: A281181 - Need Formula for Terms
Paul Hanna
pauldhanna.math at gmail.com
Fri Jan 20 05:00:51 CET 2017
Seqfans,
The e.g.f. of A281181 is one of a family of nice functions that I wish
to describe here.
A281181 begins (zero-valued coefficients omitted):
[1, 1, 13, 493, 37369, 4732249, 901188997, ...]
Notice that the e.g.f. of A281181 satisfies:
C(x) = cosh( Integral C(x)^3 dx )
as well as
C(x)^5 = d/dx Series_Reversion( Integral C(i*x)^5 dx ).
This generalizes to a nice family of interrelated functions,
characterized by the following observation.
OBSERVATION.
If
C(x) = cosh( Integral C(x)^m dx ),
then
C(x)^(2*m-1) = d/dx Series_Reversion( Integral C(i*x)^(2*m-1) dx ).
I give some examples below (given in quick PARI code, which I hope is
intelligible).
I hope to submit these sequences as a table sometime.
In case you are interested ...
Paul
EXAMPLES.
----------------------------------------------------------------------
\\ CASE m=4: C(x) = cosh( Integral C(x)^4 dx ) begins:
C=1; for(i=1,31, C = cosh( intformal( C^4 +O(x^31)))); Vec(serlaplace(C))
\\ [1, 0, 1, 0, 17, 0, 865, 0, 88865, 0, 15335425, 0, 3993275825, ...]
\\ ( d/dx Series_Reversion( Integral C(i*x)^k dx ) )^(1/k) for k=3..7
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^3 )))^(1/3)))
\\ [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^4 )))^(1/4)))
\\ [1, 0, 1, 0, 5, 0, 61, 0, 1385, 0, 50521, 0, 2702765, 0, 199360981, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^5 )))^(1/5)))
\\ [1, 0, 1, 0, 9, 0, 225, 0, 11025, 0, 893025, 0, 108056025, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^6 )))^(1/6)))
\\ [1, 0, 1, 0, 13, 0, 493, 0, 37369, 0, 4732249, 0, 901188997, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^7 )))^(1/7)))
\\ [1, 0, 1, 0, 17, 0, 865, 0, 88865, 0, 15335425, 0, 3993275825, ...]
----------------------------------------------------------------------
\\ CASE m=5: C(x) = cosh( Integral C(x)^5 dx ) begins:
C=1; for(i=1,31, C = cosh( intformal( C^5 +O(x^31)))); Vec(serlaplace(C))
[1, 0, 1, 0, 21, 0, 1341, 0, 173961, 0, 38032281, 0, 12572222301, ...]
\\ ( d/dx Series_Reversion( Integral C(i*x)^k dx ) )^(1/k) for k=4..10
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^4 )))^(1/4)))
\\ [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^5 )))^(1/5)))
\\ [1, 0, 1, 0, 5, 0, 61, 0, 1385, 0, 50521, 0, 2702765, 0, 199360981, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^6 )))^(1/6)))
\\ [1, 0, 1, 0, 9, 0, 225, 0, 11025, 0, 893025, 0, 108056025, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^7 )))^(1/7)))
\\ [1, 0, 1, 0, 13, 0, 493, 0, 37369, 0, 4732249, 0, 901188997, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^8 )))^(1/8)))
\\ [1, 0, 1, 0, 17, 0, 865, 0, 88865, 0, 15335425, 0, 3993275825, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^9 )))^(1/9)))
\\ [1, 0, 1, 0, 21, 0, 1341, 0, 173961, 0, 38032281, 0, 12572222301, ...]
Vec(serlaplace( deriv( serreverse( intformal( subst(C,x,I*x)^10 )))^(1/10)))
\\ [1, 0, 1, 0, 25, 0, 1921, 0, 301105, 0, 79715041, 0, 31953352585, ...]
----------------------------------------------------------------------
\\ CASE m=6: C(x) = cosh( Integral C(x)^6 dx ) begins:
C=1; for(i=1,31, C = cosh( intformal( C^6 +O(x^31)))); Vec(serlaplace(C))
\\ [1, 0, 1, 0, 25, 0, 1921, 0, 301105, 0, 79715041, 0, 31953352585, ...]
----------------------------------------------------------------------
[END]
More information about the SeqFan
mailing list