Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)

Revision History for A109084

(Underlined text is an addition; strikethrough text is a deletion.)

Showing all changes.
A109084 G.f. A(x) satisfies: A(x) = 1/G000041(x/A(x)) where G000041(x) is the g.f. of the partition numbers A000041.
(history; published version)
#9 by Vaclav Kotesovec at Mon Oct 02 13:28:29 EDT 2023
STATUS

editing

approved

#8 by Vaclav Kotesovec at Mon Oct 02 13:28:12 EDT 2023
MATHEMATICA

(* Calculation of constant c: *) val = Sqrt[r*s^5*(-1 + s/r)*(Log[r/s]^2 / (2*Pi*(2*s^3*(-s*Log[1 - r/s] + ArcTanh[1 - 2*r/s] * (2*r - (r - s)*(Log[1 - r/s] - 2*Log[r/s]))) + (r - s)*(s^3*(2 - 2*Log[1 - r/s] + 3*Log[r/s]) * QPolyGamma[0, 1, r/s] - s^3*QPolyGamma[0, 1, r/s]^2 + s^3*QPolyGamma[1, 1, r/s] + r*Log[r/s]*(r*Log[r/s] * Derivative[0, 2][QPochhammer][r/s, r/s] - 2*s^2*Derivative[0, 0, 1][QPolyGamma][0, 1, r/s])))))] /. FindRoot[{QPochhammer[r/s] == s, (Log[1 - r/s] + QPolyGamma[0, 1, r/s])/Log[r/s] == 1 + (r*Derivative[0, 1][QPochhammer][r/s, r/s])/s^2}, {r, 1/5}, {s, 1/2}, WorkingPrecision -> 1000]; N[Chop[val], -Floor[Log[10, Abs[Im[val]]]] - 3] (* Vaclav Kotesovec, Oct 02 2023 *)

STATUS

approved

editing

#7 by Vaclav Kotesovec at Sun May 13 07:45:08 EDT 2018
STATUS

editing

approved

#6 by Vaclav Kotesovec at Sun May 13 07:44:58 EDT 2018
FORMULA

a(n) ~ -c * d^n / n^(3/2), where d = A270915 = 5.35270133348664268777241581416... and c = 0.146705445870000769931272287955221766131167... - Vaclav Kotesovec, May 13 2018

#5 by Vaclav Kotesovec at Sun May 13 07:38:46 EDT 2018
LINKS

Vaclav Kotesovec, <a href="/A109084/b109084.txt">Table of n, a(n) for n = 0..1000</a>

STATUS

approved

editing

#4 by Russ Cox at Fri Mar 30 18:36:49 EDT 2012
AUTHOR

_Paul D. Hanna (pauldhanna(AT)juno.com), _, Jun 18 2005

Discussion
Fri Mar 30 18:36
OEIS Server: https://oeis.org/edit/global/213
#3 by N. J. A. Sloane at Sun Jun 29 03:00:00 EDT 2008
EXAMPLE

The initial terms [x^0] thruthrough [x^n] of n-th self-convolution

KEYWORD

sign,new

sign

#2 by N. J. A. Sloane at Sat Nov 10 03:00:00 EST 2007
KEYWORD

sign,new

sign

AUTHOR

Paul D . Hanna (pauldhanna(AT)juno.com), Jun 18 2005

#1 by N. J. A. Sloane at Tue Jul 19 03:00:00 EDT 2005
NAME

G.f. A(x) satisfies: A(x) = 1/G000041(x/A(x)) where G000041(x) is the g.f. of the partition numbers A000041.

DATA

1, -1, -2, -5, -17, -63, -253, -1062, -4615, -20570, -93538, -432211, -2023567, -9578815, -45767162, -220431025, -1069079067, -5216655257, -25592441875, -126157044454, -624560659184, -3103962569509, -15480272621533, -77450458331100, -388627340240958, -1955249529839424

OFFSET

0,3

COMMENTS

Note: coefficient [x^n] A(x)^n = -A000203(n) (sum of divisors of n) for n>0.

FORMULA

G.f.: A(x) = x/series_reversion(x*eta(x)). G.f.: A(x) = 1/G109085(x) where G109085(x) is g.f. of A109085.

EXAMPLE

The initial terms [x^0] thru [x^n] of n-th self-convolution

are persistently small:

A^0: 1;

A^1: 1,-1;

A^2: 1,-2,-3;

A^3: 1,-3,-3,-4;

A^4: 1,-4,-2,0,-7;

A^5: 1,-5,0,5,0,-6;

A^6: 1,-6,3,10,3,6,-12;

A^7: 1,-7,7,14,0,7,0,-8;

A^8: 1,-8,12,16,-10,0,-8,8,-15;

A^9: 1,-9,18,15,-27,-9,-21,0,0,-13;

PROG

(PARI) a(n)=polcoeff(x/serreverse(x*eta(x+x*O(x^n))), n)

CROSSREFS

Cf. A109085, A000041, A000203.

KEYWORD

sign

AUTHOR

Paul D Hanna (pauldhanna(AT)juno.com), Jun 18 2005

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 18 20:50 EDT 2024. Contains 375284 sequences. (Running on oeis4.)