同SIR模型一样，假设网路中的节点总数为1000000，最初被感染的计算机台数 ， =0。感染率以及排除率不变。在该模型中，由于采取了了反病毒措施，所以假设采用反病毒措施的概率 。
（2）模型编程
   function y=ill(t,x)
N=10^6;
a=0.8/N;b=0.05;k=0.005;
y=[-a*x(1)*x(2)-k*x(1),a*x(1)*x(2)-(b+k)*x(2),k*x(1),(b+k)*x(2)]';
>>  ts=0:40;
>> x0=[10^6-1,1,0,0];
>> [t,x]=ode45('ill',ts,x0);[t,x]
>> plot(t,x(:,2)+x(:,4)),grid,pause
（3）程序结果
ans =

  1.0e+005 *

         0   10.0000    0.0000         0         0
    0.0000    9.9501    0.0000    0.0499    0.0000
    0.0000    9.9005    0.0000    0.0995    0.0000
    0.0000    9.8510    0.0001    0.1489    0.0000
    0.0000    9.8018    0.0002    0.1980    0.0000
    0.0001    9.7527    0.0004    0.2469    0.0000
    0.0001    9.7036    0.0008    0.2955    0.0001
    0.0001    9.6543    0.0017    0.3439    0.0001
    0.0001    9.6042    0.0034    0.3921    0.0003
    0.0001    9.5525    0.0070    0.4400    0.0005
    0.0001    9.4972    0.0141    0.4876    0.0011
    0.0001    9.4344    0.0285    0.5349    0.0022
    0.0001    9.3565    0.0571    0.5819    0.0045
    0.0001    9.2487    0.1139    0.6284    0.0090
   。。。
    0.0002    1.3907    6.4180    0.9113    1.2800
    0.0002    0.8194    6.6239    0.9167    1.6400
    0.0002    0.4797    6.5960    0.9199    2.0043
    0.0003    0.2833    6.4319    0.9218    2.3631
    0.0003    0.1700    6.1966    0.9229    2.7106
    0.0003    0.1041    5.9282    0.9235    3.0441
    0.0003    0.0652    5.6483    0.9239    3.3625
    0.0003    0.0418    5.3686    0.9242    3.6654
    0.0003    0.0273    5.0951    0.9244    3.9531
    0.0003    0.0183    4.8311    0.9245    4.2261
    0.0003    0.0125    4.5782    0.9246    4.4848 传染病模型在计算机病毒传播中的应用(10):http://www.youerw.com/shuxue/lunwen_3552.html