pith. sign in

arxiv: 2606.02744 · v1 · pith:UVJ5MZK3new · submitted 2026-06-01 · ✦ hep-ph · hep-th

IterInt: Evaluating iterated integrals via differential equations

Pith reviewed 2026-06-28 13:17 UTC · model grok-4.3

classification ✦ hep-ph hep-th
keywords iterated integralsdifferential equationsnumerical evaluationmultiple polylogarithmsshuffle regularisationbanana integrals
0
0 comments X

The pith

IterInt converts iterated integrals into systems of first-order linear differential equations that standard solvers can evaluate numerically with high precision.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

IterInt is a package in Mathematica and C++ that evaluates iterated integrals with arbitrary kernels. It transforms these integrals into a system of first-order linear differential equations. These equations are then solved using established numerical libraries for high precision. The package also automatically handles shuffle-regularisation to manage cases with poles at the integration start.

Core claim

The central claim is that by transforming iterated integrals into a system of first-order linear differential equations, IterInt enables efficient and high-precision numerical evaluation of integrals involving arbitrary integration kernels, including those requiring shuffle-regularisation.

What carries the argument

The transformation of the iterated integrals, defined by user-specified kernels, into a closed system of first-order linear differential equations.

If this is right

  • Comparison with GiNaC shows agreement for ordinary and elliptic multiple polylogarithms.
  • Validation against known results for banana integrals up to four loops.
  • Supports evaluation where the integrand has a pole at the starting point via automatic regularisation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This method could allow users to define custom kernels for integrals not covered by existing special functions.
  • Integration with other computational physics tools might enable broader applications in multi-loop calculations.

Load-bearing premise

Any user-defined integration kernels can be transformed into a finite closed system of first-order linear differential equations that remains numerically stable.

What would settle it

Finding an iterated integral with arbitrary kernels for which no finite system of first-order linear differential equations can be constructed or where the numerical solution becomes unstable.

Figures

Figures reproduced from arXiv: 2606.02744 by Claude Duhr, Gideon Baur.

Figure 1
Figure 1. Figure 1: The resulting tree graph representing the relations between the iterated in￾tegrals I(ω1, ω3, ω2;t), I(ω1, ω4, ω4;t), I(ω2, ω1;t) and I(ω2, ω3, ω4;t). Each arrow represents taking a derivative. The edges passing through the dashed line are removed for the computation (see the discussion in the main text). by combining the two systems into a single larger system, d dt   I(ω1, ω2;t) I(ω1, ω3;t) I(ω1;t) I… view at source ↗
Figure 1
Figure 1. Figure 1: figure 1. Afterwards the systems of differential equations are as small as possible without [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Numerical evaluation of the unequal-mass two-loop sunrise integral I111 (blue line). The red dots correspond to an independent evaluation of the integral by numerically integrating the Feynman parameter representation. Note that the range on the x-axis corresponds to the region below threshold, where the integral is real. The equal-mass sunrise integral. If all three masses are equal, there are only three … view at source ↗
Figure 3
Figure 3. Figure 3: Numerical evaluation using IterInt of the master integral I111 of the equal￾mass sunrise family (blue dots). The red line correspond to an independent evaluation of the integral using a series representation of the result. Note that the range on the x-axis corresponds to the region below threshold, where the integral is real. proceed in a different fashion. Indeed, we may exploit the fact that Mathematica … view at source ↗
Figure 4
Figure 4. Figure 4: Numerical evaluation using IterInt of the master integral I1111 of the equal￾mass banana family (blue dots and triangles). The red lines correspond to an independent evaluation of the integral using a truncated series expansion. Note that the range on the x-axis corresponds to the region above threshold, and both the real (solid lines) and imaginary (dashed line) parts are shown. Unlike for the three-loop … view at source ↗
Figure 5
Figure 5. Figure 5: Numerical evaluation using IterInt of the master integral I11111 of the equal￾mass 4-loop equal-mass banana family (blue markers). The red lines correspond to an independent evaluation of the integral. Note that the range on the x-axis corresponds to the region above threshold. Both the real (solid lines/circles) and imaginary (dashed line/triangles) parts are shown. ematica and in C++, and we have describ… view at source ↗
read the original abstract

We introduce IterInt, a novel package implemented in both Mathematica and C++ for the numerical evaluation of iterated integrals involving arbitrary integration kernels. After the user has defined the integration kernels, IterInt transforms the iterated integrals into a system of first-order linear differential equations which can be solved efficiently and with high precision using well established libraries. IterInt is also able to automatically perform shuffle-regularisation. This makes it possible to evaluate also integrals where the integrand has a pole at the starting point of the integration path. As an illustration of our code, and also to validate it and gauge its performance, we compare the output of IterInt to the results obtained by GiNaC for ordinary and elliptic multiple polylogarithms, and also to existing results for the first few orders for banana integrals with up to four loops.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper introduces IterInt, a package implemented in both Mathematica and C++ for the numerical evaluation of iterated integrals with arbitrary user-defined kernels. The method transforms the integrals into a system of first-order linear differential equations solved using established libraries, with automatic shuffle regularisation to handle cases with poles at the integration starting point. Validation consists of direct numerical comparisons to GiNaC results for ordinary and elliptic multiple polylogarithms, as well as to published results for banana integrals up to four loops.

Significance. If the implementation is correct, IterInt supplies a general-purpose numerical tool for iterated integrals appearing in multi-loop calculations, extending beyond libraries limited to specific classes such as MPLs. The dual-language implementation and the reported agreement with independent codes constitute concrete strengths that support potential utility in the field.

major comments (1)
  1. [method description (paragraph following abstract)] The manuscript states that the transformation to a closed first-order linear DE system is always possible for arbitrary kernels, but provides no explicit algorithm, basis construction, or closure proof in the method description; this is load-bearing for the generality claim even though the finite-dimensionality argument is standard.
minor comments (2)
  1. [validation section] The validation paragraphs report agreement with GiNaC and banana results but do not tabulate achieved relative precision, CPU times, or the specific integration paths and weight values used; adding such a table would improve reproducibility.
  2. [implementation and solver paragraph] No discussion of numerical stability, step-size control, or error estimation in the DE solver is included, which would clarify the 'high precision' claim.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of IterInt and the recommendation for minor revision. We address the single major comment below.

read point-by-point responses
  1. Referee: [method description (paragraph following abstract)] The manuscript states that the transformation to a closed first-order linear DE system is always possible for arbitrary kernels, but provides no explicit algorithm, basis construction, or closure proof in the method description; this is load-bearing for the generality claim even though the finite-dimensionality argument is standard.

    Authors: We agree that an expanded description of the basis construction and closure procedure would strengthen the presentation of the generality claim. In the revised manuscript we will add a concise subsection (or short appendix) outlining the algorithmic steps: (i) the finite-dimensional vector space spanned by the iterated integrals with the user-supplied kernels, (ii) the explicit construction of a basis via repeated differentiation and reduction using the shuffle product, and (iii) the resulting closed linear system of first-order DEs. This will make the standard finite-dimensionality argument fully explicit without altering the existing code or validation results. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper describes a software package that converts user-specified iterated integrals into a first-order linear DE system via the finite-dimensional vector space of weight-n integrals under differentiation. This is a direct algebraic construction, not a fitted or self-referential prediction. Validation is performed against independent external codes (GiNaC) and published banana-integral results, providing external benchmarks. No self-citations, ansatze, or renamings of known results serve as load-bearing steps in the central claim. The method is self-contained against external checks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard existence and uniqueness theorems for linear ODEs and the algebraic properties of shuffle regularization; no new free parameters or invented entities are introduced in the abstract.

axioms (2)
  • standard math Linear first-order ODE systems derived from iterated integrals admit unique solutions given appropriate boundary conditions.
    Invoked implicitly when stating that the integrals are transformed into solvable DE systems.
  • domain assumption Shuffle regularization preserves the value of the iterated integral when a pole is present at the lower limit.
    Stated as an automatic feature without further justification in the abstract.

pith-pipeline@v0.9.1-grok · 5658 in / 1228 out tokens · 23563 ms · 2026-06-28T13:17:54.881737+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. CHESS: CHEbyshev pSeudo-Spectral transport for Feynman integral differential equations

    hep-ph 2026-06 unverdicted novelty 5.0

    CHESS package implements Chebyshev-Lobatto spectral collocation for transporting epsilon-factorized differential equations of Feynman master integrals with benchmarks showing rapid convergence and shorter wall times t...

Reference graph

Works this paper leans on

118 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    Chen,Iterated path integrals,Bull

    K.-T. Chen,Iterated path integrals,Bull. Am. Math. Soc.83(1977) 831

  2. [2]

    ’t Hooft and M.J.G

    G. ’t Hooft and M.J.G. Veltman,regularisation and Renormalization of Gauge Fields,Nucl. Phys. B44(1972) 189

  3. [3]

    Bollini and J.J

    C.G. Bollini and J.J. Giambiagi,Dimensional Renormalization: The Number of Dimensions as a Regularizing Parameter,Nuovo Cim. B12(1972) 20

  4. [4]

    Kotikov,Differential equations method: New technique for massive Feynman diagrams calculation,Phys

    A.V. Kotikov,Differential equations method: New technique for massive Feynman diagrams calculation,Phys. Lett. B254(1991) 158

  5. [5]

    Kotikov,Differential equations method: The Calculation of vertex type Feynman diagrams,Phys

    A.V. Kotikov,Differential equations method: The Calculation of vertex type Feynman diagrams,Phys. Lett. B259(1991) 314. 41

  6. [6]

    Kotikov,Differential equation method: The Calculation of N point Feynman diagrams,Phys

    A.V. Kotikov,Differential equation method: The Calculation of N point Feynman diagrams,Phys. Lett. B267(1991) 123

  7. [7]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Differential equations for two-loop four-point functions,Nucl. Phys. B580(2000) 485 [hep-ph/9912329]

  8. [8]

    Henn,Multiloop integrals in dimensional regularisation made simple,Phys

    J. Henn,Multiloop integrals in dimensional regularisation made simple,Phys. Rev. Lett.110(2013) 251601 [1304.1806]

  9. [9]

    Bourjaily et al.,Functions Beyond Multiple Polylogarithms for Precision Collider Physics, inSnowmass 2021, 3, 2022 [2203.07088]

    J.L. Bourjaily et al.,Functions Beyond Multiple Polylogarithms for Precision Collider Physics, inSnowmass 2021, 3, 2022 [2203.07088]

  10. [10]

    Goncharov,Multiple polylogarithms, cyclotomy and modular complexes, Math

    A.B. Goncharov,Multiple polylogarithms, cyclotomy and modular complexes, Math. Res. Lett.5(1998) 497 [1105.2076]

  11. [11]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Two loop master integrals for gamma* —>3 jets: The Planar topologies,Nucl. Phys. B601(2001) 248 [hep-ph/0008287]

  12. [12]

    Remiddi and J.A.M

    E. Remiddi and J.A.M. Vermaseren,Harmonic polylogarithms,Int. J. Mod. Phys. A15(2000) 725 [hep-ph/9905237]

  13. [13]

    Goncharov, M

    A.B. Goncharov, M. Spradlin, C. Vergu and A. Volovich,Classical Polylogarithms for Amplitudes and Wilson Loops,Phys. Rev. Lett.105(2010) 151605 [1006.5703]

  14. [14]

    Brown,Multiple zeta values and periods of moduli spacesM0,n,Annales Sci.Ecole Norm.Sup.42(2009) 371 [math/0606419]

    F.C. Brown,Multiple zeta values and periods of moduli spacesM0,n,Annales Sci.Ecole Norm.Sup.42(2009) 371 [math/0606419]

  15. [15]

    C. Duhr, H. Gangl and J.R. Rhodes,From polygons and symbols to polylogarithmic functions,JHEP10(2012) 075 [1110.0458]

  16. [16]

    Duhr,Hopf algebras, coproducts and symbols: an application to Higgs boson amplitudes,JHEP08(2012) 043 [1203.0454]

    C. Duhr,Hopf algebras, coproducts and symbols: an application to Higgs boson amplitudes,JHEP08(2012) 043 [1203.0454]

  17. [17]

    C. Duhr,Mathematical aspects of scattering amplitudes, inTheoretical Advanced Study Institute in Elementary Particle Physics: Journeys Through the Precision Frontier: Amplitudes for Colliders, pp. 419–476, 2015, DOI [1411.7538]

  18. [18]

    Weinzierl,Feynman Integrals

    S. Weinzierl,Feynman Integrals. A Comprehensive Treatment for Students and Researchers, UNITEXT for Physics, Springer (2022), 10.1007/978-3-030-99558-4, [2201.03593]

  19. [19]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Numerical evaluation of harmonic polylogarithms, Comput. Phys. Commun.141(2001) 296 [hep-ph/0107173]

  20. [20]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Numerical evaluation of two-dimensional harmonic polylogarithms,Comput. Phys. Commun.144(2002) 200 [hep-ph/0111255]

  21. [21]

    Vollinga and S

    J. Vollinga and S. Weinzierl,Numerical evaluation of multiple polylogarithms, Comput. Phys. Commun.167(2005) 177 [hep-ph/0410259]. 42

  22. [22]

    Buehler and C

    S. Buehler and C. Duhr,CHAPLIN - Complex Harmonic Polylogarithms in Fortran,Comput. Phys. Commun.185(2014) 2703 [1106.5739]

  23. [23]

    Frellesvig, D

    H. Frellesvig, D. Tommasini and C. Wever,On the Reduction and Evaluation of Generalized Polylogarithms,PoSLL2016(2016) 040 [1609.00148]

  24. [24]

    Naterop, A

    L. Naterop, A. Signer and Y. Ulrich,handyG —Rapid numerical evaluation of generalised polylogarithms in Fortran,Comput. Phys. Commun.253(2020) 107165 [1909.01656]

  25. [25]

    Brown and A

    F. Brown and A. Levin,Multiple Elliptic Polylogarithms,1110.6917

  26. [26]

    Broedel, C.R

    J. Broedel, C.R. Mafra, N. Matthes and O. Schlotterer,Elliptic multiple zeta values and one-loop superstring amplitudes,JHEP07(2015) 112 [1412.5535]

  27. [27]

    Broedel, C

    J. Broedel, C. Duhr, F. Dulat and L. Tancredi,Elliptic polylogarithms and iterated integrals on elliptic curves. Part I: general formalism,JHEP05(2018) 093 [1712.07089]

  28. [28]

    Enriquez and F

    B. Enriquez and F. Zerbini,Elliptic hyperlogarithms,Canadian Journal of Mathematics(2023)

  29. [29]

    Walden and S

    M. Walden and S. Weinzierl,Numerical evaluation of iterated integrals related to elliptic Feynman integrals,Comput. Phys. Commun.265(2021) 108020 [2010.05271]

  30. [30]

    C. Duhr, F. Lorkowski, R. Marzucca, S. Mauc and S. Weinzierl,Elliptic Multiple Polylogarithms with Arbitrary Arguments inGiNaC,2602.09956

  31. [31]

    Duhr and F

    C. Duhr and F. Brown,A double integral of dlog forms which is not polylogarithmic,PoSMA2019(2022) 005 [2006.09413]

  32. [32]

    Manin,Iterated integrals of modular forms and noncommutative modular symbols, inAlgebraic geometry and number theory, vol

    Y.I. Manin,Iterated integrals of modular forms and noncommutative modular symbols, inAlgebraic geometry and number theory, vol. 253 ofProgr. Math., (Boston), pp. 565–597, Birkhäuser Boston, 2006 [math/0502576]

  33. [33]

    Brown,Multiple Modular Values and the relative completion of the fundamental group ofM1,1,1407.5167

    F. Brown,Multiple Modular Values and the relative completion of the fundamental group ofM1,1,1407.5167

  34. [34]

    Matthes,Iterated primitives of meromorphic quasimodular forms for sl2(Z), Transactions of the American Mathematical Society375(2022) 1443

    N. Matthes,Iterated primitives of meromorphic quasimodular forms for sl2(Z), Transactions of the American Mathematical Society375(2022) 1443

  35. [35]

    Broedel, C

    J. Broedel, C. Duhr and N. Matthes,Meromorphic modular forms and the three-loop equal-mass banana integral,JHEP02(2022) 184 [2109.15251]

  36. [36]

    Bönisch, C

    K. Bönisch, C. Duhr, F. Fischbach, A. Klemm and C. Nega,Feynman integrals in dimensional regularisation and extensions of Calabi-Yau motives,JHEP09 (2022) 156 [2108.05310]. 43

  37. [37]

    Pögel, X

    S. Pögel, X. Wang and S. Weinzierl,Taming Calabi-Yau Feynman Integrals: The Four-Loop Equal-Mass Banana Integral,Phys. Rev. Lett.130(2023) 101601 [2211.04292]

  38. [38]

    Pögel, X

    S. Pögel, X. Wang and S. Weinzierl,Bananas of equal mass: any loop, any order in the dimensional regularisation parameter,JHEP04(2023) 117 [2212.08908]

  39. [39]

    Pögel, X

    S. Pögel, X. Wang and S. Weinzierl,The three-loop equal-mass banana integral in ε-factorised form with meromorphic modular forms,JHEP09(2022) 062 [2207.12893]

  40. [40]

    C. Duhr, S. Maggio, C. Nega, B. Sauer, L. Tancredi and F.J. Wagner,Aspects of canonical differential equations for Calabi-Yau geometries and beyond,JHEP06 (2025) 128 [2503.20655]

  41. [41]

    Klemm, C

    A. Klemm, C. Nega, B. Sauer and J. Plefka,Calabi-Yau periods for black hole scattering in classical general relativity,Phys. Rev. D109(2024) 124046 [2401.07899]

  42. [42]

    Driesse, G.U

    M. Driesse, G.U. Jakobsen, A. Klemm, G. Mogull, C. Nega, J. Plefka et al., Emergence of Calabi–Yau manifolds in high-precision black-hole scattering,Nature 641(2025) 603 [2411.11846]

  43. [43]

    Dlapa, G

    C. Dlapa, G. Kälin, Z. Liu and R.A. Porto,Local-in-Time Conservative Binary Dynamics at Fifth Post-Minkowskian and First Self-Force Orders,Phys. Rev. Lett.135(2025) 251401 [2506.20665]

  44. [44]

    Huang and Y

    R. Huang and Y. Zhang,On Genera of Curves from High-loop Generalized Unitarity Cuts,JHEP04(2013) 080 [1302.1023]

  45. [45]

    Hauenstein, R

    J.D. Hauenstein, R. Huang, D. Mehta and Y. Zhang,Global Structure of Curves from Generalized Unitarity Cut of Three-loop Diagrams,JHEP02(2015) 136 [1408.3355]

  46. [46]

    Marzucca, A.J

    R. Marzucca, A.J. McLeod, B. Page, S. Pögel and S. Weinzierl,Genus drop in hyperelliptic Feynman integrals,Phys. Rev. D109(2024) L031901 [2307.11497]

  47. [47]

    C. Duhr, F. Porkert and S.F. Stawinski,Canonical differential equations beyond genus one,JHEP02(2025) 014 [2412.02300]

  48. [48]

    Bargiela, H

    P. Bargiela, H. Frellesvig, R. Marzucca, R. Morales, F. Seefeld, M. Wilhelm et al., The spectrum of Feynman-integral geometries at two loops,JHEP05(2026) 057 [2512.13794]

  49. [49]

    F. Yang, J. Gong and Y. Zhang,Genus drop involving non-hyperelliptic curves in Feynman integrals,2605.07729. 44

  50. [50]

    D’Hoker, M

    E. D’Hoker, M. Hidding and O. Schlotterer,Constructing polylogarithms on higher-genus Riemann surfaces,Commun. Num. Theor. Phys.19(2025) 355 [2306.08644]

  51. [51]

    D’Hoker and O

    E. D’Hoker and O. Schlotterer,Fay identities for polylogarithms on higher-genus Riemann surfaces,2407.11476

  52. [52]

    Baune, J

    K. Baune, J. Broedel, E. Im, A. Lisitsyn and F. Zerbini,Schottky–Kronecker forms and hyperelliptic polylogarithms,J. Phys. A57(2024) 445202 [2406.10051]

  53. [53]

    Baune, J

    K. Baune, J. Broedel, E. Im, A. Lisitsyn and Y. Moeckli,Higher-genus Fay-like identities from meromorphic generating functions,SciPost Phys.18(2025) 093 [2409.08208]

  54. [54]

    D’Hoker, B

    E. D’Hoker, B. Enriquez, O. Schlotterer and F. Zerbini,Relating Flat Connections and Polylogarithms on Higher Genus Riemann Surfaces,Commun. Math. Phys.407(2026) 43 [2501.07640]

  55. [55]

    D’Hoker and O

    E. D’Hoker and O. Schlotterer,Meromorphic higher-genus integration kernels via convolution over homology cycles,J. Phys. A58(2025) 33LT01 [2502.14769]

  56. [56]

    Baune, J

    K. Baune, J. Broedel, E. Im, Z. Ji and Y. Moeckli,Higher-genus multiple zeta values,J. Phys. A58(2025) 465401 [2507.21765]

  57. [57]

    D’Hoker and O

    E. D’Hoker and O. Schlotterer,Single-valued flat connections in several variables on arbitrary Riemann surfaces,2602.01461

  58. [58]

    D’Hoker and O

    E. D’Hoker and O. Schlotterer,Equivalence of flat connections and Fay identities on arbitrary Riemann surfaces,2602.09108

  59. [59]

    Berger and J

    M. Berger and J. Broedel,Translating auxiliary symmetries between Schottky uniformization and Jacobi parametrization,2602.04715

  60. [60]

    Prausa and J

    M. Prausa and J. Usovitsch,The analytic leading color contribution to the Higgs-gluon form factor in QCD at NNLO,JHEP03(2021) 127 [2008.11641]

  61. [61]

    Görges, C

    L. Görges, C. Nega, L. Tancredi and F.J. Wagner,On a procedure to derive ϵ-factorised differential equations beyond polylogarithms,JHEP07(2023) 206 [2305.14090]

  62. [62]

    Forner, C.C

    F. Forner, C.C. Mella, C. Nega, L. Tancredi and F.J. Wagner,Integrand Analysis, Leading Singularities and Canonical Bases beyond Polylogarithms,2604.25270

  63. [63]

    Bauer, A

    C.W. Bauer, A. Frink and R. Kreckel,Introduction to the GiNaC framework for symbolic computation within the C++ programming language,J. Symb. Comput. 33(2002) 1 [cs/0004015]

  64. [64]

    Gaiotto, J

    D. Gaiotto, J. Maldacena, A. Sever and P. Vieira,Pulling the straps of polygons, JHEP12(2011) 011 [1102.0062]. 45

  65. [65]

    Deligne,Le groupe fondamental de la droite projective moins trois points, in Galois Groups overQ, Y

    P. Deligne,Le groupe fondamental de la droite projective moins trois points, in Galois Groups overQ, Y. Ihara, K. Ribet and J.-P. Serre, eds., vol. 16 of Mathematical Sciences Research Institute Publications, (New York), pp. 79–297, Springer (1989), DOI

  66. [66]

    Brown,The Massless higher-loop two-point function,Commun

    F. Brown,The Massless higher-loop two-point function,Commun. Math. Phys. 287(2009) 925 [0804.1660]

  67. [67]

    Runge,Ueber die numerische Auflösung von Differentialgleichungen,Math

    C. Runge,Ueber die numerische Auflösung von Differentialgleichungen,Math. Ann.46(1895) 167

  68. [68]

    Kutta,Beitrag zur näherungsweisen integration totaler differentialgleichungen, Zeitschrift für Mathematik und Physik(1901) 435

    W. Kutta,Beitrag zur näherungsweisen integration totaler differentialgleichungen, Zeitschrift für Mathematik und Physik(1901) 435

  69. [69]

    The C++ Alliance,The Boost C++ Libraries,https://www.boost.org/

  70. [70]

    The GNU Project,GSL - GNU Scientific Library, https://www.gnu.org/software/gsl/

  71. [71]

    Fousse, G

    L. Fousse, G. Hanrot, V. Lefèvre, P. Pélissier and P. Zimmermann,Mpfr: A multiple-precision binary floating-point library with correct rounding,ACM Trans. Math. Softw.33(2007) 13

  72. [72]

    A. Enge, M. Gastineau, P. Théveny and P. Zimmermann,mpc — A library for multiprecision complex arithmetic with exact rounding. INRIA, 1.3.0 ed., Dec., 2022

  73. [73]

    Maddock and C

    J. Maddock and C. Kormanyos,Boost.Multiprecision documentation, https://www.boost.org/doc/libs/latest/libs/multiprecision/doc/html/index.html (2020)

  74. [74]

    Ahnert and M

    K. Ahnert and M. Mulansky,Boost.Numeric.Odeint documentation, https://www.boost.org/doc/libs/latest/libs/numeric/odeint/doc/html/index.html (2015)

  75. [75]

    Bartels,Numerical Mathematics 3 x 9, Springer (2025)

    S. Bartels,Numerical Mathematics 3 x 9, Springer (2025)

  76. [76]

    Hanke-Bourgeois,Grundlagen der numerischen Mathematik und des wissenschaftlichen Rechnens, Springer (2006)

    M. Hanke-Bourgeois,Grundlagen der numerischen Mathematik und des wissenschaftlichen Rechnens, Springer (2006)

  77. [77]

    Galassi et al.,GNU Scientific Library Reference Manual (3rd Ed.)(2009)

    M. Galassi et al.,GNU Scientific Library Reference Manual (3rd Ed.)(2009)

  78. [78]

    Lappo-Danilevskij,Mémoires sur la théorie des systèmes des équations différentielles linéaires

    J.A. Lappo-Danilevskij,Mémoires sur la théorie des systèmes des équations différentielles linéaires. Vol. II,Travaux de l’Institut Physico-Mathématique Stekloff7(1935) 5

  79. [79]

    Broedel, N

    J. Broedel, N. Matthes, G. Richter and O. Schlotterer,Twisted elliptic multiple zeta values and non-planar one-loop open-string amplitudes,J. Phys. A51(2018) 285401 [1704.03449]. 46

  80. [80]

    Tarasov,Connection between Feynman integrals having different values of the space-time dimension,Phys

    O.V. Tarasov,Connection between Feynman integrals having different values of the space-time dimension,Phys. Rev. D54(1996) 6479 [hep-th/9606018]

Showing first 80 references.