擅长:python、mysql、java
<p>当然,你可以一起破解一些可能有用的东西。
您可以记录<code>t</code>,但必须注意
<em>可能不会持续增加。这取决于ODE算法及其工作方式(正向、反向和中心有限差分)</p>
<p>但它会让你知道你在哪里</p>
<pre><code>logger = [] # visible in odefunc
def odefunc(P,t):
x = P[0]
y = P[1]
z = P[2]
print(t)
logger.append(t)
if logger: # if the list is not empty
if logger[-1] > 2.5: # then read the last value
print('hua!')
if A * x - B * x * y < 0.6:
dxdt = A/5 * x
dydt = -C * y + D * x * y
dzdt = - B * z * y
else:
dxdt = A * x - B * x * y
dydt = -C * y + D * x * y
dzdt = 0
dPdt = np.ravel([dxdt, dydt, dzdt])
return dPdt
print(logger)
</code></pre>