在Python中集成数组

如何解决在Python中集成数组

我是第一次使用Simpson方法在python中进行集成。 我已经编写了代码，它可以成功运行，但是我不知道它能否给出正确的输出。

这是我的输入文件：我想在x的范围内执行y的积分（在这种情况下为57.3到101.8）。

 X        Y
57.3    794.158
59.4    816.297
61.5    748.34
63.8    834.679
66.1    810.91
68.5    882.33
71.0    837.228
73.7    898.849
76.4    882.926
79.1    928.2
82.0    1122.0341
85.1    1120.606
88.2    1213.764
91.4    1268.479
94.7    1303.058
98.2    1433.6331
101.8   1642.562

以下是我的python脚本：

import pandas as pd
import numpy as np
from scipy.integrate import simps

def print_integration(filename='infile'):

    df=pd.read_csv(filename,delim_whitespace=True,header=None)
    x=np.array(df[0])
    y=np.array(df[1])

    length=len(x)
    upper_bound=x[1]
    lower_bound=x[length-1]
    del_x=(lower_bound - upper_bound)/length

    intg=simps(y,dx=del_x)
    print('{}'.format(intg))

if __name__ == '__main__':
    print_integration()

这给了我一个答案：

40741.8545537255

这是正确的答案吗？ 预先谢谢你。