我正在尝试使用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
的第一个值需要由trapz
设定。在你的情况下force
不是自我计算力的功能,但已经评估点x2
。