update for scipy > 1.14

eazy/photoz.py

@@ -2587,8 +2587,11 @@ def show_fit(self, id, id_is_idx=False, zshow=None, show_fnu=0, get_spec=False,
         """
         import matplotlib.pyplot as plt
         from matplotlib.gridspec import GridSpec
-        from scipy.integrate import cumtrapz
-
+        try:
+            from scipy.integrate import cumtrapz
+        except ImportError:
+            from scipy.integrate import cumulative_trapezoid as cumtrapz
+
         import astropy.units as u
         from cycler import cycler
 
@@ -3824,7 +3827,10 @@ def pz_percentiles(self, percentiles=[2.5,16,50,84,97.5], oversample=5,
             
         """
         import scipy.interpolate 
-        from scipy.integrate import cumtrapz
+        try:
+            from scipy.integrate import cumtrapz
+        except ImportError:
+            from scipy.integrate import cumulative_trapezoid as cumtrapz
 
         interpolator = scipy.interpolate.Akima1DInterpolator