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.