我有一个微分方程系统，我试图绘制，它具有后勤增长，应该围绕各自的承载能力进行平衡。然而，由于某些不可思议的原因，他们会无限制地走向无限。
我的代码是：
eqn1 =（ - （S [t] / K2）T [t] + r2（1 - S [t] / K2）* S [t]）
eqn2 =（（S [t] / K2）T [t] + r1（1 - T [t] / K1）* T [t]）
K2 = 2.5 * 10 ^ 10
K1 = 10 ^ 11
r1 = .001
r2 = .001
solns = NDSolve [{S＆＃39; [t] == eqn1，T＆＃39; [t] == eqn2，S [0] == 10000，T [0] == 10000}，{S， T}，{t，0,5000}]
绘图[评估[{S [t]，T [t]} /。 solns]，{t，0,100000}，PlotStyle - ＆gt; {红色，蓝色}]
非常感谢任何见解。