[seqfan] Lambert W function and converging factors: merging A001662 and A217538?
Paolo Bonzini
bonzini at gnu.org
Wed Jun 22 23:39:29 CEST 2016
While preparing a new series for OEIS, I stumbled upon some weirdness in
the existing series A001662 and A051711:
- A001662 ("Numerators in expansion of W(exp(x)) about x=1, where W is
the Lambert function") says in the comments section "Please note: a(17)
is not 10125320047141". A217538 is the same sequence as A001662 except
that a(17) is indeed 10125320047141.
- A051711 ("a(0) = 1; for n>0, a(n) = n!*4^n/2") says "For n <= 16,
denominators in expansion of W(exp(x)) about x=1, where W is the Lambert
function".
>From the history, the original definition of A001662 was "Coefficients
of Airey's converging factor" and, based on Airey's paper, a(17) was
10125320047141. Edit #9 changed it to the current form. However, I
believe that the edit was incorrect.
First of all, let's factor the old and new 17-th term:
10125320047141 = 13*15683*49663379
778870772857 = 15683*49663379
This suggests that the coefficient of the series is A001662/A051711, but
the 17th term of A001662/A051711 is not a reduced fraction. So as
things stand A001662 is *not* the coefficient of Airey's converging
factor, despite saying so in the comments.
In addition, in the comments, "(-1)^n times the polynomials with
coefficients in triangle A008517, evaluated at -1" is wrong. The values
in the sequence are simply the polynomials with coefficients in triangle
A008517, evaluated at -1:
1-2= -1
1-8+6= -1
1-22+58-24= 13
1-52+328-444+120= -47
1-114+1452-4400+3708-720= -73
without the need for another sign adjustment. This is also visible from
the paper "A Sequence of Series for The Lambert W Function" (Section 2.2),
and is referenced in the comments section of A008517.
Another "yellow flag" is the 1 + 2*x - 9346449274284*x^17 - x*Q(0)
generating function, where of course the huge x^17 term is simply
778870772857-10125320047141.
Luckily, A217638 already has good quality content for the formula,
maple, mathematica and programs sections. Therefore, my suggestion is:
1) change A001662 as follows
- restore the title to "Coefficients of Airey's converging factor"
- restore the 17th term to 10125320047141
- replace the comments section with:
n!*4^n/2 times the coefficient in expansion of W(exp(x)) about x=1,
where W is the Lambert function. - Paolo Bonzini, Jun 22 2016
The polynomials with coefficients in triangle A008517, evaluated at
-1.
- add to the links the following paper (cited in A217538)
R. M. Corless, D. J. Jeffrey and D. E. Knuth, A sequence of series
for the Lambert W Function (section 2.2).
http://www.apmaths.uwo.ca/~rcorless/frames/PAPERS/LambertW/CorlessJeffreyKnuth.ps
- replace the formula, maple, mathematica and prog sections of A001662
with those in A217538
2) delete A217538 which is now a duplicate of A001662 only with fewer
links and references.
3) create two new series:
- "Numerators in expansion of W(exp(x)) about x=1, where W is the
Lambert function, A001662/gcd(A001662,A051711)". This series can reuse
most of the content and programs of A001662.
- "Denominators in expansion of W(exp(x)) about x=1, where W is the
Lambert function, A051711/gcd(A001662,A051711)"
If this is okay, I can submit edits for (1) and (3), but I thought it'd
be clearer to post here first and get approval.
Thanks,
Paolo Bonzini
More information about the SeqFan
mailing list