Rigorous numerical integration of algebraic functions

We present an algorithm which uses Fujiwara’s inequality to bound algebraic functions over ellipses of a certain type, allowing us to concretely implement a rigorous Gauss–Legendre integration method for algebraic functions over a line segment. We consider path splitting strategies to improve convergence of the method and show that these yield significant practical and asymptotic benefits. We implemented these methods to compute period matrices of algebraic Riemann surfaces and these are available in SageMath.