3810.sagews

#AUFGABE 10c)

b = 1/11 
c = 6.8 
a = b*c

S0 = 10000 
I0 = 1000 
R0 = 19000 
t = 40 

S = [S0] 
I = [I0] 
R = [R0] 

n = S0+I0+R0 

for i in range(1,t+1): 
    S.append(S[i-1]-a*I[i-1]*S[i-1]/n) 
    I.append(I[i-1]+a*I[i-1]*S[i-1]/n -b*I[i-1]) 
    R.append(R[i-1]+b*I[i-1]) 
    print 't=',i,'S=',S[i],'I=',I[i],'R=',R[i] 

print 'Liste S=',S 
print 'Liste I=',I 
print 'Liste R=',R 

#sinnvoller wäre vielleicht die Ausgabe in aufgerundeten Zahlen

t= 1 S= 9793.93939393939 I= 1115.15151515152 R= 210000/11 t= 2 S= 9568.88563875671 I= 1238.82785986588 R= 19192.2865013774 t= 3 S= 9324.61723149398 I= 1370.47555259535 R= 19304.9072159107 t= 4 S= 9061.28908700643 I= 1509.21501048332 R= 19429.4959025103 t= 5 S= 8779.49227539982 I= 1653.81045750054 R= 19566.6972670996 t= 6 S= 8480.30018531227 I= 1802.65614236076 R= 19717.0436723270 t= 7 S= 8165.29399293897 I= 1953.78450361036 R= 19880.9215034507 t= 8 S= 7836.56087438859 I= 2104.90084910525 R= 20058.5382765062 t= 9 S= 7496.66012062421 I= 2253.44698022370 R= 20249.8928991521 t= 10 S= 7148.55521895157 I= 2396.69306551237 R= 20454.7517155361 t= 11 S= 6795.51379317122 I= 2531.85330351887 R= 20672.6329033099 t= 12 S= 6440.98149161010 I= 2656.21712294191 R= 20902.8013854480 t= 13 S= 6088.43970912327 I= 2767.28462152493 R= 21144.2756693518 t= 14 S= 5741.25961846345 I= 2862.89338295521 R= 21395.8469985814 t= 15 S= 5402.56575068653 I= 2941.32421591802 R= 21656.1100333955 t= 16 S= 5075.12107553511 I= 3001.37578053144 R= 21923.5031439335 t= 17 S= 4761.24244141464 I= 3042.40207096723 R= 22196.3554876181 t= 18 S= 4462.75100421424 I= 3064.31150171607 R= 22472.9374940697 t= 19 S= 4180.95778368618 I= 3067.53094936084 R= 22751.5112669530 t= 20 S= 3916.68057666766 I= 3052.94170643748 R= 23030.3777168949 t= 21 S= 3670.28571943597 I= 3021.79640853848 R= 23307.9178720256 t= 22 S= 3441.74686429642 I= 2975.62649926544 R= 23582.6266364382 t= 23 S= 3230.71292012213 I= 2916.14894350651 R= 23853.1381363714 t= 24 S= 3036.57827628046 I= 2845.17913793850 R= 24118.2425857811 t= 25 S= 2858.54996626757 I= 2764.55479904789 R= 24376.8952346846 t= 26 S= 2695.70814024505 I= 2676.07346152060 R= 24628.2183982344 t= 27 S= 2547.05781146935 I= 2581.44438470352 R= 24871.4978038271 t= 28 S= 2411.57114790104 I= 2482.25428602605 R= 25106.1745660729 t= 29 S= 2288.22053225876 I= 2379.94542112051 R= 25331.8340466207 t= 30 S= 2176.00322361576 I= 2275.80405511619 R= 25548.1927212681 t= 31 S= 2073.95877716211 I= 2170.95722383200 R= 25755.0839990059 t= 32 S= 1981.18048817233 I= 2066.37576520069 R= 25952.4437466270 t= 33 S= 1896.82209192602 I= 1962.88181915604 R= 26140.2960889180 t= 34 S= 1820.10083232231 I= 1861.15927701829 R= 26318.7398906594 t= 35 S= 1750.29785252091 I= 1761.76595890893 R= 26487.9361885702 t= 36 S= 1686.75669134893 I= 1665.14657836192 R= 26648.0967302892 t= 37 S= 1628.88050799961 I= 1571.64580004197 R= 26799.4736919584 t= 38 S= 1576.12851459858 I= 1481.52090253010 R= 26942.3505828713 t= 39 S= 1528.01197548261 I= 1394.95372323424 R= 27077.0343012832 t= 40 S= 1484.09003378093 I= 1312.06169009644 R= 27203.8482761226 Liste S= [10000, 9793.93939393939, 9568.88563875671, 9324.61723149398, 9061.28908700643, 8779.49227539982, 8480.30018531227, 8165.29399293897, 7836.56087438859, 7496.66012062421, 7148.55521895157, 6795.51379317122, 6440.98149161010, 6088.43970912327, 5741.25961846345, 5402.56575068653, 5075.12107553511, 4761.24244141464, 4462.75100421424, 4180.95778368618, 3916.68057666766, 3670.28571943597, 3441.74686429642, 3230.71292012213, 3036.57827628046, 2858.54996626757, 2695.70814024505, 2547.05781146935, 2411.57114790104, 2288.22053225876, 2176.00322361576, 2073.95877716211, 1981.18048817233, 1896.82209192602, 1820.10083232231, 1750.29785252091, 1686.75669134893, 1628.88050799961, 1576.12851459858, 1528.01197548261, 1484.09003378093] Liste I= [1000, 1115.15151515152, 1238.82785986588, 1370.47555259535, 1509.21501048332, 1653.81045750054, 1802.65614236076, 1953.78450361036, 2104.90084910525, 2253.44698022370, 2396.69306551237, 2531.85330351887, 2656.21712294191, 2767.28462152493, 2862.89338295521, 2941.32421591802, 3001.37578053144, 3042.40207096723, 3064.31150171607, 3067.53094936084, 3052.94170643748, 3021.79640853848, 2975.62649926544, 2916.14894350651, 2845.17913793850, 2764.55479904789, 2676.07346152060, 2581.44438470352, 2482.25428602605, 2379.94542112051, 2275.80405511619, 2170.95722383200, 2066.37576520069, 1962.88181915604, 1861.15927701829, 1761.76595890893, 1665.14657836192, 1571.64580004197, 1481.52090253010, 1394.95372323424, 1312.06169009644] Liste R= [19000, 210000/11, 19192.2865013774, 19304.9072159107, 19429.4959025103, 19566.6972670996, 19717.0436723270, 19880.9215034507, 20058.5382765062, 20249.8928991521, 20454.7517155361, 20672.6329033099, 20902.8013854480, 21144.2756693518, 21395.8469985814, 21656.1100333955, 21923.5031439335, 22196.3554876181, 22472.9374940697, 22751.5112669530, 23030.3777168949, 23307.9178720256, 23582.6266364382, 23853.1381363714, 24118.2425857811, 24376.8952346846, 24628.2183982344, 24871.4978038271, 25106.1745660729, 25331.8340466207, 25548.1927212681, 25755.0839990059, 25952.4437466270, 26140.2960889180, 26318.7398906594, 26487.9361885702, 26648.0967302892, 26799.4736919584, 26942.3505828713, 27077.0343012832, 27203.8482761226]

#AUFGABE 10d)

s=list_plot(S,plotjoined=True,color='#ffa500',linestyle='--',legend_label='$S - Anfaellige Personen$')
i=list_plot(I,plotjoined=True,color='#ee4000',linestyle=':',legend_label='$I - Infizierte Personen$')
r=list_plot(R,plotjoined=True,color='#7cfc00',legend_label='$R - Genesene Personen$')

p=s+i+r
p.set_legend_options(loc = 0)
p.axes_labels(['$Zeit. in. Tagen$','$Anzahl. der. Personen$'])
show(p)