Intégration numérique Python sans SciPy?

Question

J'ai besoin d'intégrer la fonction: y (x) = e ^-ax cos (x) sans utiliser SciPy/NumPy etc. Je me bats un peu. Je sais que cela:

def integrand(x): 
       return 

Pour le bit de retour, je veux utiliser cette équation: la somme i=1-N de y (xi) (xi - x (i-1))

Comment pourrais-je écrire cela?

Answer 1

définie simplement comme:

from math import cos, e 

def y(x): 
    return (e**(-a*x))*cos(x) 

def integrand(x): 
    return sum(y(x[i])*(x[i]-x[i-1]) for i in range(1,len(x)))

Les range(1,len(x)) est compris entre 1 (inclus) à la longueur de x (exclusif).

Ou utilisez:

from math import exp, cos 

def y(x): 
    return exp(-a*x)*cos(x)