// An highlighted blockimport numpy as np##定义未知量和待定参数amplitude和detuningamplitude=20detuning=20#郎之万方程经典形式def dfun(t,x,y,q,p):dvec=np.zeros((4))dvec[0]=-0.5*x-(detuning+0.02*q)*y+amplitudedvec[1]=-0.5*y+(detuning+0.02*q)*xdvec[2]=20*pdvec[3]=0.02*(x**2+y**2)-20*q-0.01*preturn dvec#四阶龙格库塔:dfun=微分方程；t为自变量，x,y,q,p为因变量；求解区间(a,b)；x0,y0,q0,p0=初值,h=步长def main(dfun,a,b,t,x,y,q,p,x0,y0,q0,p0,h):K1=dfun(t,x,y,q,p)[0]L1=dfun(t,x,y,q,p)[1]
