Beschreibung

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

Lizenz

Ich, der Urheberrechtsinhaber dieses Werkes, veröffentliche es als gemeinfrei. Dies gilt weltweit. In manchen Staaten könnte dies rechtlich nicht möglich sein. Sofern dies der Fall ist: Ich gewähre jedem das bedingungslose Recht, dieses Werk für jedweden Zweck zu nutzen, es sei denn, Bedingungen sind gesetzlich erforderlich.