Parameter r,a;
Variable I,R,,E,t;
InitialValue I=346,D=0,R=0,E=852,S=10110000-346-852;
ODEFunction S'=-r*0.03*I*S/10110000-r*0.03*E*S/10110000
E'=r*0.03*I*S/10110000-a*E+r*0.03*E*S/10110000
I'=a*E-0.0677*I
R'=0.0677*I
Data;
//确诊 治愈 潜伏 天数
346 0 852 1
464 0 844 2
635 0 817 3
687 0 965 4
833 1 941 5
895 1 1039 6
1030 1 1014 7
1197 1 1103 8
1293 15 1231 9
1435 17 1383 10
1636 16 1423 11
1759 16 1701 12
1919 16 1794 13
2029 16 2012 14
2285 16 2163 15
2523 16 2437 16
2819 16 2762 17
3059 16 3085 18
3460 16 3049 19
3713 16 3172 20
4041 16 3331 21
4448 16 3650 22
4960 16 3864 23
5494 103 3965 24
5955 2