日志
一道DAE的maple解答和mmtc解答
||
Clear["Global`*"]
SetOptions[Plot, PlotRange -> Full, Frame -> True];
sys = {a'[t] == -(a[t]/10) + 10000 b[t] c[t],
b'[t] == a[t]/10 - 10000 b[t] c[t] - 10000000 b[t]^2,
a[t] + b[t] + c[t] == 1, a[0] == 1, b[0] == 0};
sol = First[NDSolve[sys, {a, b, c}, {t, 0, 40000}]];
Plot[a[t] /. sol, {t, 0, 40000}]
Plot[Evaluate[{a[t], b[t], c[t]} /. sol], {t, 0, 40000}]
SetOptions[Plot, PlotRange -> Full, Frame -> True];
sys = {a'[t] == -(a[t]/10) + 10000 b[t] c[t],
b'[t] == a[t]/10 - 10000 b[t] c[t] - 10000000 b[t]^2,
a[t] + b[t] + c[t] == 1, a[0] == 1, b[0] == 0};
sol = First[NDSolve[sys, {a, b, c}, {t, 0, 40000}]];
Plot[a[t] /. sol, {t, 0, 40000}]
Plot[Evaluate[{a[t], b[t], c[t]} /. sol], {t, 0, 40000}]
