使用odeint求解具有快速振荡的不同方程

Question

我已经编写了此代码来求解方程，我知道此功能的行为具有非常快的振荡，当我运行时，它会为某些“ m [x]”和某些“ t” s的虚假值，带有此错误：

c：\ users \ dani \ anaconda3 \ lib \ lib \ site-packages \ scipy \ intemate \ odepack.py：247：odeintwarning：在此呼叫上完成的多余工作（也许是错误的dfun类型）。使用full_output = 1运行以获取定量信息。 WARNINGS.WARN（WARNNING_MSG，ODEINTWARNING）

我不知道有什么问题。 如何获得正确的结果？还是至少尽可能准确？还是我应该以另一种形式重写代码？ 谢谢。

import scipy as sio
import numpy as np
import mpmath as mp
import scipy.integrate as spi
import matplotlib.pyplot as plt
import time
from scipy.integrate import quad
initial_value=np.logspace(24,27,100)
t=np.logspace(-20,6,100)
m=np.logspace(0,6,100)
start_time=time.perf_counter()
phi_m={}
phi_m_prime={}
phi=[]
phi_prime=[]
j=0
i=np.pi*2.435*initial_value[0]
while i<(np.pi*(2.435*10**(27))):
    i=np.pi*2.435*initial_value[j]
    phi=[]
    phi_prime=[]
    for x in range (len(m)):
        def dzdt(z,T):
            return [z[1], -3*1.4441*(10**(-6))*m[x]*np.sqrt(0.69)*(mp.coth(1.5*np.sqrt(0.69)*(10**(-6))*1.4441*m[x]*T))*z[1] - z[0]]
        z0 = [i,0]
        ts = t/m[x]
        zs = spi.odeint(dzdt, z0, ts)
        phi.append(zs[99,0])
        phi_prime.append(zs[99,1])
    phi_m[j]=phi
    phi_m_prime[j]=phi_prime
    j+=1
end_time=time.perf_counter()
print(end_time-start_time,"seconds")