Python:Simpson的规则,曲线下正面积的否定答案?

时间:2016-04-22 21:41:32

标签: python scipy integrate simpsons-rule

我有

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]

所有y值均为正值,因此曲线下面积应为正值。但是当我尝试使用辛普森的规则来整合

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

  print b

我得到-2.45630795891e-05作为答案。我究竟做错了什么?

1 个答案:

答案 0 :(得分:4)

x正在减少,因此您从右到左进行整合。这翻转了结果的标志。