Python的解决与四功能二阶ODE

Question

我学习阻尼的动态，用二阶ODE驱动摆定义像so，具体我预设电台：

d^2Y/dt的^ 2 + C * DY/DT +罪（Y）= A * cos（重量）

import numpy as np 
import matplotlib.pyplot as plt 
from scipy import integrate 


def pendeq(t,y,a,c,w): 
    y0 = -c*y[1] - np.sin(y[0]) + a*np.cos(w*t) 
    y1 = y[1] 
    z = np.array([y0, y1]) 
    return z 

a = 2.1 
c = 1e-4 
w = 0.666667  # driving angular frequency 

t = np.array([0, 5000]) # interval of integration 
y = np.array([0, 0])  # initial conditions 

yp = integrate.quad(pendeq, t[0], t[1], args=(y,a,c,w)) 

这个问题看起来很相似Need help solving a second order non-linear ODE in python，但我得到我是什么错误

Supplied function does not return a valid float. 

做错了？

Answer 1

integrate.quad要求提供的功能（pendeq，在你的情况下）只返回一个浮点数。你的函数正在返回一个数组。