Replace scipy's romberg with quad.

Romberg will be removed from scipy in v1.15.0

sbpy/activity/gas/core.py

@@ -24,7 +24,7 @@
 try:
     import scipy
     from scipy import special
-    from scipy.integrate import quad, dblquad, romberg
+    from scipy.integrate import quad, dblquad
     from scipy.interpolate import CubicSpline, PPoly
 except ImportError:
     scipy = None
@@ -1635,8 +1635,8 @@ def _column_density_at_rho(self, rho: np.float64) -> np.float64:
         def column_density_integrand(z):
             return self._volume_density(np.sqrt(z**2 + rhosq))
 
-        c_dens = 2 * romberg(column_density_integrand, 0,
-                             z_max, rtol=0.0001, divmax=50)
+        c_dens = 2 * quad(column_density_integrand, 0,
+                          z_max, rtol=0.0001)
 
         # result is in 1/m^2
         return c_dens
@@ -1745,13 +1745,12 @@ def _calc_num_fragments_grid(self) -> np.float64:
         def vol_integrand(r, r_func):
             return r_func(r) * r**2
 
-        r_int = romberg(
+        r_int = quad(
             vol_integrand,
             0,
             max_r,
             args=(self.vmr.volume_density_interpolation,),
             rtol=0.0001,
-            divmax=20,
         )
         return 4 * np.pi * r_int