mercredi 2 juillet 2014

Mathematica Not Graphing Diff Eqs Properly


Vote count:

0




I have a system of differential equations that I'm trying to graph, which have logistic growth and should be leveling out around their respective carrying capacities. However, they're going off to infinity for some indiscernible reason.


My code is:



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'[t] == eqn1, T'[t] == eqn2, S[0] == 10000 , T[0] == 10000}, {S, T}, {t, 0, 5000}]


Plot[Evaluate[{S[t], T[t]} /. solns], {t, 0, 100000}, PlotStyle -> {Red, Blue}]



Any insights would be much appreciated.



asked 1 min ago






Aucun commentaire:

Enregistrer un commentaire