level 4
lwy1993414
楼主
请吧内大神帮忙运行一下这段求微分方程的代码,我查破头也不懂到底哪里出了问题,万分感谢~
s = NDSolve[{x'[
t] == -I*(t*x[t] - (3/(2*(t + 44.132*10^6)))*5.25*10^34*x[t]*
Abs[x[t]]^2 +
1.5*1.25*10^33*Conjugate[x[t]]*Conjugate[x[t]]*
y[t]/(t + 44.132*10^6)) - 0.5*993*x[t] +
0.25*100/(t + 44.132*10^6),
y'[t] == -I*((3*(t + 44.132*10^6) - 132.25*10^6) -
0.5*(-8.3*10^32)*y[t]*Abs[y[t]]*
Abs[y[t]]/(t + 44.132*10^6) +
0.16*1.25*10^33*x[t]^3/(t + 44.132*10^6)) - 1900*0.5*y[t]},
x[0] == 10^(-8), y[0] == 10^(-8), {x, y}, {t, 0, 1000}];
ParametricPlot[Evaluate[{x[t], y[t]} /. s], {t, 0, 1000}]
2020年06月30日 08点06分
1
s = NDSolve[{x'[
t] == -I*(t*x[t] - (3/(2*(t + 44.132*10^6)))*5.25*10^34*x[t]*
Abs[x[t]]^2 +
1.5*1.25*10^33*Conjugate[x[t]]*Conjugate[x[t]]*
y[t]/(t + 44.132*10^6)) - 0.5*993*x[t] +
0.25*100/(t + 44.132*10^6),
y'[t] == -I*((3*(t + 44.132*10^6) - 132.25*10^6) -
0.5*(-8.3*10^32)*y[t]*Abs[y[t]]*
Abs[y[t]]/(t + 44.132*10^6) +
0.16*1.25*10^33*x[t]^3/(t + 44.132*10^6)) - 1900*0.5*y[t]},
x[0] == 10^(-8), y[0] == 10^(-8), {x, y}, {t, 0, 1000}];
ParametricPlot[Evaluate[{x[t], y[t]} /. s], {t, 0, 1000}]


