【发布时间】：2018-03-02 12:29:18
【问题描述】：
我正在尝试使用 np.trapz 打印积分的值。
我只得到一个值列表/数组而不是积分。
l0=0.5
h=0.3
k=500
time = 10
dt = 0.01
n=int(time/dt)
def F(x2):
return (-k*x2*(1-(l0/(sqrt(x2**2+h**2)))))-(k*h*(1-(l0/(sqrt(x2**2+h**2))))-g*m)
force=F(x2)*array([1,0])
x2=linspace(0.4,0.75,n+1)
W=trapz(force)
print(W)
【问题讨论】：
-
不完整的代码，缺乏预期的输出（与不是最小的可行示例有关），缺乏
F的定义... -
请添加
g、x2和m -
仍然不是 MVCE 并且缺少导入和缩进问题。
-
您对
F的调用很奇怪，因为您没有定义x2。反转force=...和x2=...行，看看会发生什么。 -
请检查您包含在问题中的代码是否可运行。目前看起来有些行的顺序不正确（例如，
x2在定义之前使用）。查看minimal, complete and verifiable example 的建议。