[seqfan] Re: A063747 sign of power series expansion of the Gamma function

[Allan Wechsler]

One of you said "e-38" and the other said "e-58". Do you really have 20
orders of magnitude of disagreement?

On Thu, May 18, 2023 at 11:07 PM Brendan McKay via SeqFan <
seqfan at list.seqfan.eu> wrote:

> I also get -0.444..e-58  * x^46 in Maple using 60 digits.  The calculation
> is numerically unstable though.
>
> B/
>
> On 19/5/2023 12:31 pm, Sean A. Irvine wrote:
> > Hi,
> >
> > I believe the terms in A063747 are incorrect for n>=46.
> >
> > https://oeis.org/A063747
> >
> > My attempts to reproduce the calculation in Maple and using constructible
> > reals in Java both yield a coefficient, like -0.44458...e-38 * x^46, from
> > which a(46) should be -1 rather than 1.
> >
> > Is someone able to make an independent check?
> >
> > I was unable to run the existing Pari program beyond n=15, but I assume
> > there is some way to tell Pari to try harder?
> >
> > (I previously emailed the author of the sequence, but received no
> response.)
> >
> > Sean.
> >
> > --
> > Seqfan Mailing list - http://list.seqfan.eu/
>
>
> --
> Seqfan Mailing list - http://list.seqfan.eu/
>