diff --git a/mbtrack2/utilities/beamloading.py b/mbtrack2/utilities/beamloading.py
index a2771cd83eabe06b6b5515ee40f9a7c697c94b73..d4f0cda386b1d154ff31b11122b508e96e2f0aab 100644
--- a/mbtrack2/utilities/beamloading.py
+++ b/mbtrack2/utilities/beamloading.py
@@ -414,7 +414,7 @@ class BeamLoadingEquilibrium():
 
         """
         rho0 = gaussian_bunch(self.tau0, self.ring.sigma_0) / c
-        R = trapezoid(rho0**2, self.tau0) / trapz(self.rho0**2, self.tau0)
+        R = trapezoid(rho0**2, self.tau0) / trapezoid(self.rho0**2, self.tau0)
         return R
 
     @property