Fixed a bug in HaloFeedback

bradkav · Feb 6, 2020 · bf10d27 · bf10d27
1 parent 2768f65
commit bf10d27
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/HaloFeedback.py b/HaloFeedback.py
@@ -304,7 +304,7 @@ def dfdt_minus(self, r0, v_orb, v_cut=-1, n_kick = 1):
            step = np.append(0, step)
 
            #Make sure that the integral is normalised correctly
-           renorm = simps(self.P_delta_eps(v_orb, delta_eps_list), delta_eps_list)
+           renorm = np.trapz(self.P_delta_eps(v_orb, delta_eps_list), delta_eps_list)
            frac_list = 0.5*(step[:-1] + step[1:])/renorm
 
         #Sum over the kicks