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

[Brendan McKay]

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/