level 3
SUPER逍遥游
楼主
感觉做的很迷呀，求指点！a1 = 1; a2 = 10; a3 = 30; lam = 0.5; r = 10;
s = NDSolve[{Q'[t] == lam r l[t] - (1 - (1 - 1/(a2 - a1))/E^t) Q[t],
l'[t] == (1 - (1 - 1/(a2 - a1))/E^t) Q[
t] - (1 - (1 - 1/(a2 - a1))/E^t) l[t],
T'[t] == (1 - (1 - 1/(a2 - a1))/E^t) l[t],
Q[0] == T[0] == l[0] == 1}, {Q, T, l}, {t, 0, 2}]
Plot[Evaluate[l[t] /. s], {t, 0, 2}]
