updates for scipy 1.14

CHANGELOG-unreleased.md

@@ -22,6 +22,7 @@ the released changes.
 - The following observatories no longer have a default of `include_bipm=False`: magic, lst, virgo, lho, llo, geo600, kagra, hess, hawc
 - New algorithm for TCB <-> TDB conversion
 - Reordered plotting axes in `pintk`
+- Changed `scipy.integrate.simps` to `scipy.integrate.simpson` to work with scipy 1.14
 ### Added
 - `bayesian_information_criterion()` function 
 - `dmx_setup` function

src/pint/templates/lcprimitives.py

@@ -16,7 +16,7 @@
 from math import atan, cosh, tan
 
 import numpy as np
-from scipy.integrate import quad, simps
+from scipy.integrate import quad, simpson
 from scipy.interpolate import interp1d
 from scipy.special import erf, i0, i1, owens_t
 from scipy.stats import cauchy, norm, vonmises, skewnorm
@@ -1593,7 +1593,7 @@ def integrate(self, x1=0, x2=1):
         x = self.interpolator.x
         y = self.interpolator.y
         mask = (x >= x1) & (x <= x2)
-        return simps(y[mask], x=x[mask])
+        return simpson(y[mask], x=x[mask])
         # return self.interpolator.y[mask].sum()/len(mask)