Resumen

sirsys.m program

beta=0.001;v=0.1;%;%before modified, beta=0.005，v=0.05，
                       the figure produced by running this program in matlab is different with Fig.Sirsys-p9.png
                     so let beta=0.001,v=0.1, and try again, then the resoult is right.
options = odeset ('RelTol',1e-4,'AbsTol',1e-4);
[T,Y] = ode45(@sirsys1,[0 60],[500 1 0],options,beta,v);
plot(T,Y(:,1),'.',T,Y(:,2),'.',T,Y(:,3),'.')

sirsys1.m function

function dy = sirsys1(t,y,beta,v)
dy = zeros (3,1);
dy(1)=-beta*y(1)*y(2);            %(force of infection) lambda = beta*y(2)
dy(2)=beta*y(1)*y(2)-v*y(2);
dy(3)=v*y(2);

Licencia

Yo, el titular de los derechos de autor de esta obra, lo libero al dominio público. Esto aplica en todo el mundo. En algunos países esto puede no ser legalmente factible; si ello ocurriese: Concedo a cualquier persona el derecho de usar este trabajo para cualquier propósito, sin ningún tipo de condición al menos que éstas sean requeridas por la ley.