Description

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);

Conditions d’utilisation

Moi, propriétaire des droits d’auteur sur cette œuvre, la place dans le domaine public. Ceci s'applique dans le monde entier. Dans certains pays, ceci peut ne pas être possible ; dans ce cas : J’accorde à toute personne le droit d’utiliser cette œuvre dans n’importe quel but, sans aucune condition, sauf celles requises par la loi.