现有Allen的病毒扩散模型SIRS的马尔科夫过程程序，求指导
 

由Allen编写的SIRS模型的马尔可夫过程如下：
for k=1:4 % 解马尔科夫链的样本轨迹
clear t s i r
t(1)=0; i(1)=i0; s(1)=N-i0;r(1)=0;
j=1;
while i(j)>0 && t(j)<time % Stop hits zero or at time
nn=s(j)+i(j)+r(j);
u1=rand;u2=rand;
den=((beta/nn)*i(j)*s(j)+gam*i(j)+alpha*i(j)+delta*r(j));%所有patch的概率之和
t(j+1)=-log(u1)/den+t(j); % Time to next event
e1=(beta/nn)*s(j)*i(j)/den;%以下是每个patch概率与总概率之比
e2=e1+gam*i(j)/den;
e3=e2+alpha*i(j)/den;
e4=e3+delta*r(j)/den;
if (u2<=e1) %每个patch情况下的个体变化情况
s(j+1)=s(j)-1;
i(j+1)=i(j)+1;
r(j+1)=r(j);
elseif (u2>e1 && u2<=e2)
s(j+1)=s(j);
i(j+1)=i(j)-1;
r(j+1)=r(j)+1;
elseif (u2>e2 && u2<=e3)
s(j+1)=s(j);
i(j+1)=i(j)-1;
r(j+1)=r(j);
else
s(j+1)=s(j)+1;
i(j+1)=i(j);
r(j+1)=r(j)-1;
end
j=j+1;
end
（省略出图部分）
end
本人菜鸟，对上述标黄部分的分支概率不太理解，变量“den”是各个转换概率之和，e1-e4分别是四种转换概率与“den”的比值，请问各位大神，e1-e4的顺序有区别么？另附本程序的转换图，见“SIR模型”图，若想在本程序的基础上，改写成“SIRP模型”图所应用的程序，应该怎么写呢？
目前存在两个问题：1.S→I之间有两种转换概率的情况该怎么解决呢？
2. P作为独立变量，对S→I有影响，但其增加或减少的变动情况与S,I,R相对独立，P的每一步数量应该单独统计么？
感谢各位！论文卡在这里很久了，请大家帮忙。
[ATTACH]3388[/ATTACH]

[ATTACH]3389[/ATTACH]