from mpmath import mp
mp.dps = 120
f = lambda x: mp.diff(mp.scorergi, x)
roots = [mp.findroot(f, (0, 1))]
x = mp.mpf(0)
y = f(x)
while len(roots) < 13:
t = x - mp.mpf(1)/16
z = f(t)
if y*z < 0:
roots.append(mp.findroot(f, (t, x)))
x, y = t, z
for n, root in enumerate(roots):
print(n, mp.nstr(root, 100))
The values are written to 100 significant digits. Arb real balls at 700 bits enclose a 900-term Maclaurin expansion of the solution fixed by (1) and (2). With $\operatorname{Gi}(x)=\sum_{k\geq0}c_kx^k$, the coefficients satisfy $c_2=-1/(2\pi)$ and $c_{k+3}=c_k/((k+3)(k+2))$. Separate geometric tails for the three residue classes bound the function and its first two derivatives on $|x|\leq16$. Their bounds are respectively less than $7\cdot10^{-432}$, $4\cdot10^{-430}$ and $3\cdot10^{-428}$ and are added to the ball radii.
The generator covers $[-465/32,0]$ and $[0,1]$ by 994 adjacent rational intervals. Local enclosures follow from the differential equation and a Gronwall bound, with its derivation in the attached generator. On 976 intervals the enclosure of $\operatorname{Gi}^{\prime}$ excludes zero. On the other eighteen, $\operatorname{Gi}^{\prime\prime}$ has a fixed nonzero sign; five have equal endpoint signs and thirteen have opposite endpoint signs. This certifies twelve negative roots, one positive root and no omitted roots in the covered intervals. The uniqueness of the positive root on the entire nonnegative axis is established in [2].
Each root is further enclosed by exact rational endpoints separated by $2\cdot10^{-120}$, with opposite signs of $\operatorname{Gi}^{\prime}$ and a nonzero sign of $\operatorname{Gi}^{\prime\prime}$ throughout the bracket. The bracket is strictly inside the interval denoted by the stored decimal. Independent certified quadrature of (3) confirms the signs at both endpoints for all thirteen roots. Integrals of $t^2$ and $\operatorname{Ai}^{\prime}$ with known nonzero answers serve as controls, and each integral evaluation is finite with absolute radius less than $10^{-145}$.
Independent mpmath root computations at 150 and 200 decimal digits [7] agree to 140 significant digits and lie inside every certified bracket. Checks of the initial values, differential equation and differentiated sum identity use mpmath computations independent of the generator. These comparisons supplement the rigorous enclosures. The attached checks repeat the bracket and integral tests on all intervals read from the database.