Python Numerische Integration ohne SciPy?

Question

Ich brauche die Funktion zu integrieren: y (x) = e ^ax cos (x) ohne SciPy/NumPy etc. Ich bin kämpfen, um eine ganze Menge. Ich weiß, das:

def integrand(x): 
       return 

Für die Rücklieferbit Ich mag diese Gleichung verwenden: die Summe i=1-N von y (xi) (xi - x (i-1))

Wie würde ich das schreiben?

Answer 1

einfach definieren als:

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)))

Die range(1,len(x)) reicht von 1 (einschließlich) zur Länge des x (exklusiv).

oder -nutzung:

from math import exp, cos 

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