J'essaie d'écrire un programme qui me dessine la courbe d'une fonction définie par une intégrale à paramètre :
$ F(t)=\int_{0}^{\infty }exp(-xt)(1+4x)dx $
J'ai écrit ce code mais apparemment il fonctionne pas et j'ignore la raison ... :
Code : Tout sélectionner
import numpy as np
import scipy.integrate as integrate
import matplotlib.pyplot as plt
def F(t):
def g(x):
return np.exp(-t*x)*(1+4*x)
return integr.quad(F,0,np.inf)[0]
T=np.linspace(1,10)
Y=F(T)
plt.plot(T,Y)
plt.show()