Avec quad il est possible d'intégrer une fonction qui prend plusieurs arguments, exemple avec une fonction f qui prend deux arguments arg1 et arg2:
quad( f, x_min, x_max, args=(arg1,arg2,))
Exemple de code:
from scipy.integrate import quad
import matplotlib.pyplot as plt
import scipy.stats
import numpy as np
def normal_distribution_function(x,mean,std):
value = scipy.stats.norm.pdf(x,mean,std)
return value
x_min = 0.0
x_max = 30.0
mean = 15.0
std = 4.0
ptx = np.linspace(x_min, x_max, 100)
pty = scipy.stats.norm.pdf(ptx,mean,std)
plt.plot(ptx,pty, color='gray')
plt.fill_between(ptx, pty, color='#e1b1b4', alpha='1.0')
plt.grid()
plt.title('How to integrate a function that takes parameteres in python ?', fontsize=10)
plt.xlabel('x', fontsize=8)
plt.ylabel('Probability Density Function', fontsize=8)
res, err = quad(normal_distribution_function, x_min, x_max, args=(mean,std,))
print(res)
plt.savefig("integrate_function_takes_parameters.png")
plt.show()
Références
Links | Site |
---|---|
How to evaluate single integrals of multivariate functions with Python's scipy.integrate.quad? | stackoverflow |
quad | doc scipy |