val, abserr
=> 0.7853981633971.63822994214e-13
注意到我们为y积分的边界传参的方式,这样写是因为y可能是关于x的函数。...在这个例子的实现中,我们会加上额外的参数到 RHS 方程中:
def dy(y, t, zeta,w0):
"""
The right-hand side of the dampedoscillator...t, args=(0.0, w0)) # undamped
y2 = odeint(dy, y0, t, args=(0.2, w0)) # under damped
y3 = odeint(dy,...y0, t, args=(1.0, w0)) # critialdamping
y4 = odeint(dy, y0, t, args=(5.0, w0)) # over damped
fig, ax...= plt.subplots()
ax.plot(t, y1[:,0], 'k', label="undamped", linewidth=0.25)
ax.plot(t, y2[:,0], 'r',