2017-03-10 19 views
0

我需要整合功能:Y(X)= E -ax COS(x)的不使用SciPy的/ NumPy的等我掙扎了不少。我知道這一點:沒有SciPy的Python數值整合?

def integrand(x): 
       return 

對於返回位我想用這個公式:從i=1Y(XI)(XI - X(I-1))N總和

我該怎麼寫這個?

+0

總和大概是從'1'到**'的N- 1' **? –

回答

2

只需將其定義爲:

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

從1(含)的範圍range(1,len(x))到的x(不包括)的長度。

或使用:

from math import exp, cos 

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

@ Kasramvd:謝謝,修改。 –

+0

謝謝。但是,我現在失去了如何繼續。我很習慣能夠使用SciPy。當'a'= 0時,我希望它打印整數的值。 – dtfd1998

+0

@webberwizard:你首先將'a'設置爲'a = 0' ... –