From 7a9c4bf00aaf9db43bd63ba19d6968621e748114 Mon Sep 17 00:00:00 2001 From: Abhimanyu Susobhanan Date: Fri, 15 Aug 2025 14:05:18 +0200 Subject: [PATCH] scipy.integrate.simps -> simpson --- pygedm/yt2020.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/pygedm/yt2020.py b/pygedm/yt2020.py index bca4aa8..45fd828 100644 --- a/pygedm/yt2020.py +++ b/pygedm/yt2020.py @@ -48,7 +48,7 @@ # integrated function excluding n_0^{sphe} (left-hand side of eq.[4]) y_sphe = 4.0*np.pi*r_array**2*mu_e*m_p/c_scaling*np.exp(-Upsilon*(1-np.log(1+r_array/r_s)*r_s/r_array)) -I_sphe = integrate.simps(y_sphe,r_array) # integration +I_sphe = integrate.simpson(y_sphe,r_array) # integration n_0_sphe = M_b/Z_halo/I_sphe # central electron number density for spherical comp. [cm^{-3}] #print('Upsilon = %f' % Upsilon) #print('n_0_sphe = %f [cm^{-3}]' % n_0_sphe)