[Fix] trapz

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