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 205 3946 25

6320 205 4047 26

6696 205 4003 27

7183 205 4013 28

7909 205 3769 29

8635 205 3656 30

9357 205 3610 31

9725 381 3672 32

10191 381 3769 33

10523 381 3829 34

11020 381 3877 35

11502 381 4077 36

12115 381 4215 37

12622 550 4351 38

13228 550 4355 39

13791 550 4255 40

14183 550 4149 41

14728 550 4299 42

15410 550 4298 43

16161 550 4337 44

16973 550 3953 45

17128 1005 3905 46

17591 1005 3880 47

17877 1005 3998 48

18572 1005 3798 49

19253 1005 3819 50

20043 1005 3734 51

20471 1005 3948 52

21033 1005 4042 53

21471 1005 4005 54

21875 1005 3949 55

19301 4074 4056 56

20003 4074 3991 57

20708 4074 3959 58

20453 4971 3942 59

21109 4971 3756 60

21510 4971 3821 61

21858 4971 3707 62

22460 4971 3447 63

23097 4971 3232 64

23770 4971 3283 65

24395 4971 3307 66

24865 4971 3474 67

25331 4971 3387 68

25565 4971 3424 69

25907 4971 3466 70

26329 4971 3641 71

27053 4971 3565 72

27702 4971 3555 73

28339 4971 3667 74

28718 4971 3925 75

28989 4971 4083 76

29373 4971 3971 77

29970 4971 4149 78

30618 4971 4488 79

31257 4971 5014 80

32006 4971 5345 81

32643 4971 6712 82

33072 4971 7126 83

33344 4971 7257 84

34119 4971 7273 85

35106 4971 7176 86

36271 4971 7485 87

37351 4971 7801 88

39355 4971 7044 89

40198 4971 6884 90

40601 4971 7250 91

41392 4971 7399 92

42282 4971 7748 93

43756 4971 7755 94

45152 4971 7745 95

46399 4971 7793 96

47082 4971 8749 97

47851 4971 9733 98

48791 4971 10698 99

50030 4971 10946 100

51511 4971 11031 101

52897 4971 10892 102

54192 4971 10746 103

55831 4971 10181 104

57584 4971 9365 105

59489 4971 8286 106

60976 4971 7583 107

62542 4971 6756 108

63789 4971 6456 109

64938 4971 6087 110

66012 4971 5740 111

66949 4971 5487 112

67775 4971 5410 113

68559 4971 5162 114

69298 4971 4842 115

70245 4971 4220 116

71025 4971 3733 117

71752 4971 3370 118

72436 4971 3079 119

73185 4971 2677 120

73721 4971 2405 121

74140 4971 2213 122

74465 4971 2273 123

74758 4971 2384 124

75122 4971 2435 125

75515 4971 2388 126

75862 4971 2159 127

76126 4971 2233 128

76353 4971 2143 129

76738 4971 1880 130

77142 4971 1710 131

77557 4971 1496 132

77903 4971 1277 133

78021 4971 1248 134

78359 4971 1009 135

78496 4971 1160 136

78618 4971 1678 137

78852 4971 1582 138

79053 4971 1538 139

79180 4971 1593 140

79269 4971 1673 141

79368 4971 1933 142