尝试此代码
tspan = [0 500];
ic = [1;0;0;0;0];
[t，p] = ode45（@odefun，tspan，ic）;
情节（t，p）
传奇（{'P'，，，，'p \ _dot1'，，，，'p \ _dot2'，，，，'p \ _dot3'，，，，'p \ _dot4'}）
功能dpdt = odefun（t，p）
a = [-1 0 0 0 0;
1 -1 0 0 0;
0 1 -1 0 0;
0 0 1 -1 0;
0 0 0 1 -1];
b = [0 1 0 0 0;
0 -1 1 0 0;
0 0 -1 1 0;
0 0 0 -1 1;
0 0 0 0 -1];
lambda = 0.2;
t = 100;
g = t/2;
S0 = 0.5;
如果（t floor（t/t）*t）
S = S0;
别的
s = 0;
结尾
dpdt =（lambda*a + s*b）*p;
结尾
