2016-04-22 1 views
2

J'aiPython: la règle de Simpson, réponse négative pour la zone positive sous la courbe?

y1=[ 9.49110000e-004 4.74145420e-004 1.41847155e-008 3.33228420e-028 
     3.76352289e-081 4.48206815e-193 0.00000000e+000 0.00000000e+000 
     0.00000000e+000 0.00000000e+000 0.00000000e+000] 

x=[ 112. 111.97667396 111.90666665 111.78989038 111.62619837 
    111.41538384 111.15717866 110.85125168 110.49720654 110.09457901 
    109.64283388] 

Toutes les valeurs y sont positives, de sorte que la surface sous la courbe devrait être positive. Mais quand j'essaie d'utiliser la règle de Simpson pour intégrer

from scipy.integrate import simps 
    b= simps(y1, x) 

    print b 

Je reçois -2.45630795891e-05 comme la réponse. Qu'est-ce que je fais mal?

Répondre

3

diminue, vous intégrez ainsi de droite à gauche. Cela retourne le signe du résultat.

+0

J'ai inversé à la fois «y1» et «x», j'ai toujours eu la même réponse. – HuShu

+0

@HuShu: Eh bien, [quand je les inverse] (http://ideone.com/o3tepD), la réponse est positive. – user2357112

+0

Wow! Merci. Je ne sais pas comment/pourquoi cela s'est passé, mais après avoir redémarré mon noyau et inversé la commande, cela m'a donné un nombre positif. – HuShu