LA10 -- Sage

#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

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]

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)

LA10

166 days ago by Tina