diff --git a/spaceKLIP/utils.py b/spaceKLIP/utils.py
index 561da6cf..8f8f9a74 100644
--- a/spaceKLIP/utils.py
+++ b/spaceKLIP/utils.py
@@ -16,7 +16,7 @@
 import importlib
 import scipy.ndimage.interpolation as sinterp
 
-from scipy.integrate import simps
+from scipy.integrate import simpson
 from scipy.ndimage import fourier_shift, gaussian_filter
 from scipy.ndimage import shift as spline_shift
 
@@ -654,8 +654,8 @@ def _get_tp_comsubst(instrume,
             # Compute COM substrate transmission averaged over the respective
             # filter profile.
             bandpass_throughput = np.interp(comsubst_wave, bandpass_wave, bandpass_throughput)
-            int_tp_bandpass = simps(bandpass_throughput, comsubst_wave)
-            int_tp_bandpass_comsubst = simps(bandpass_throughput * comsubst_throughput, comsubst_wave)
+            int_tp_bandpass = simpson(bandpass_throughput, comsubst_wave)
+            int_tp_bandpass_comsubst = simpson(bandpass_throughput * comsubst_throughput, comsubst_wave)
             tp_comsubst = int_tp_bandpass_comsubst / int_tp_bandpass
     
     # Return.