概要

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

ライセンス

この著作物の著作権者である私は、この著作物における権利を放棄しパブリックドメインとします。これは全世界で適用されます。 一部の国では、これが法的に可能ではない場合があります。その場合は、次のように宣言します。 私は、あらゆる人に対して、法により必要とされている条件を除き、如何なる条件も課すことなく、あらゆる目的のためにこの著作物を使用する権利を与えます。