在我们使用Python包scipy.integrate中的simps函数计算积分时,如果积分上下限一致,在数学上我们知道,此时积分值应该为0。然而当我们写程序时,如果采用一样的上下限来计算积分,会出问题,代码如下:
import numpy as np from scipy.integrate import simps num_ = 100 y_min = 1.0 y_max = 1.0 y_ = np.linspace(y_min, y_max, num_

在我们使用Python包scipy.integrate中的simps函数计算积分时,如果积分上下限一致,在数学上我们知道,此时积分值应该为0。然而当我们写程序时,如果采用一样的上下限来计算积分,会出问题,代码如下:
import numpy as np from scipy.integrate import simps num_ = 100 y_min = 1.0 y_max = 1.0 y_ = np.linspace(y_min, y_max, num_