ts = 0:50;i0 = 0.02;s0 = 0.98;x0=[i0;s0];[t,x]=ode45(@sirfun,ts,x0);figure; plot(x(:,2),x(:,1)); grid on; function dy = sirfun(t,x)% x(1)--> i, x(2)-->s% di/dt = λsi-μi% ds/dt = -λsi lambda = 1; mu = 0.3; dy = [lambda*x(1)*x(2)。
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024