-
Notifications
You must be signed in to change notification settings - Fork 0
/
SIR with dempgraphy.nb
15836 lines (15760 loc) · 854 KB
/
SIR with dempgraphy.nb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 8.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 157, 7]
NotebookDataLength[ 874361, 15827]
NotebookOptionsPosition[ 871429, 15747]
NotebookOutlinePosition[ 871889, 15764]
CellTagsIndexPosition[ 871846, 15761]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]",
RowBox[{"Needs", "[", "\"\<PlotLegends`\>\"", "]"}], "\n",
RowBox[{
RowBox[{"eq1", "=",
RowBox[{
RowBox[{"\[Mu]", " ", "n"}], "-",
RowBox[{"\[Mu]", " ",
RowBox[{"s", "[", "t", "]"}]}], "-",
RowBox[{"\[Beta]", " ", "*", " ",
RowBox[{"s", "[", "t", "]"}], "*",
FractionBox[
RowBox[{"i", "[", "t", "]"}], "n"]}]}]}], ";"}], "\n",
RowBox[{
RowBox[{"eq2", "=",
RowBox[{
RowBox[{"\[Beta]", "*",
RowBox[{"s", "[", "t", "]"}], "*",
FractionBox[
RowBox[{"i", "[", "t", "]"}], "n"]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{"\[Mu]", "+", "\[Gamma]"}], ")"}], "*",
RowBox[{"i", "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"eq3", "=",
RowBox[{
RowBox[{"\[Gamma]", "*",
RowBox[{"i", "[", "t", "]"}]}], "-",
RowBox[{"\[Mu]", " ",
RowBox[{"r", "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Mu]", "=", "0.6"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Beta]", "=", "0.5"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Gamma]", "=", "0.4"}], ";"}], "\[IndentingNewLine]",
RowBox[{"R0", "=",
RowBox[{"N", "[",
FractionBox["\[Beta]",
RowBox[{"(",
RowBox[{"\[Gamma]", "+", "\[Mu]"}], ")"}]],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"tf", "=", "100"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"szero", "=", "0.5"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"izero", "=", "0.4"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rzero", "=", "0.1"}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"sol", "=",
RowBox[{"NDSolve", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"s", "'"}], "[", "t", "]"}], "\[Equal]", "eq1"}], ",",
RowBox[{
RowBox[{
RowBox[{"i", "'"}], "[", "t", "]"}], "\[Equal]", "eq2"}], ",",
RowBox[{
RowBox[{
RowBox[{"r", "'"}], "[", "t", "]"}], "\[Equal]", "eq3"}], ",",
RowBox[{
RowBox[{"s", "[", "0", "]"}], "\[Equal]", "szero"}], ",",
RowBox[{
RowBox[{"i", "[", "0", "]"}], "\[Equal]", "izero"}], ",",
RowBox[{
RowBox[{"r", "[", "0", "]"}], "\[Equal]", "rzero"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"s", ",", "i", ",", "r"}], "}"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "tf"}], "}"}]}], "]"}]}], ";"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"plot1", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"a", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"s", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",",
RowBox[{"\"\<S\>\"", "[", "t", "]"}]}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Blue"}]}]}], "]"}]}], " ",
RowBox[{"(*", " ",
RowBox[{"For", " ",
RowBox[{"S", "[", "t", "]"}], " ", "versus", " ", "t"}], " ",
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"plot1a", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"a", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"s", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",", "\"\<Density\>\""}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Blue"}]}], ",",
RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",",
RowBox[{"GridLinesStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Gray", ",", "Dotted"}], "]"}]}]}], "]"}]}], ";"}], " ",
RowBox[{"(*", " ",
RowBox[{"For", " ",
RowBox[{"S", "[", "t", "]"}], " ", "versus", " ", "t"}], " ", "*)"}],
"\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"plot2", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"b", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"i", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",",
RowBox[{"\"\<I\>\"", "[", "t", "]"}]}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Red"}]}]}], "]"}]}], " ",
RowBox[{"(*", " ",
RowBox[{"For", " ",
RowBox[{"I", "[", "t", "]"}], " ", "versus", " ", "t"}], " ",
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"plot2a", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"b", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"i", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",", "\"\<Density\>\""}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Red"}]}], ",",
RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",",
RowBox[{"GridLinesStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Gray", ",", "Dotted"}], "]"}]}]}], "]"}]}], " ", ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]"}], "\n",
RowBox[{
RowBox[{"plot3", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"c", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"r", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",",
RowBox[{"\"\<R\>\"", "[", "t", "]"}]}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Green"}]}]}], "]"}]}], " ",
RowBox[{"(*", " ",
RowBox[{"For", " ",
RowBox[{"R", "[", "t", "]"}], " ", "versus", " ", "t"}], " ",
"*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"plot3a", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{
RowBox[{"c", "[", "t", "]"}], "=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"r", "[", "t", "]"}], "/.", "sol"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"t", ",", "\"\<Density\>\""}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Green"}], ",",
RowBox[{"Mesh", "\[Rule]", "10"}], ",",
RowBox[{"MeshStyle", "\[Rule]",
RowBox[{"Darker", "@", "Green"}]}], ",",
RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",",
RowBox[{"GridLinesStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Gray", ",", "Dotted"}], "]"}]}]}], "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"ShowLegend", "[",
RowBox[{
RowBox[{"Show", "[",
RowBox[{"plot1a", " ", ",", "plot2a", ",", "plot3a"}], "]"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Graphics", "[",
RowBox[{"{",
RowBox[{"Blue", ",",
RowBox[{"Line", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}]}], "}"}], "]"}]}], "}"}],
"]"}], ",", "\"\<S(t)\>\""}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Graphics", "[",
RowBox[{"{",
RowBox[{"Red", ",",
RowBox[{"Line", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}]}], "}"}], "]"}]}], "}"}],
"]"}], ",", "\"\<I(t)\>\""}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Graphics", "[",
RowBox[{"{",
RowBox[{"Green", ",",
RowBox[{"Line", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}]}], "}"}], "]"}]}], "}"}],
"]"}], ",", "\"\<R(t)\>\""}], "}"}]}], "}"}], ",",
RowBox[{"LegendShadow", "\[Rule]", "None"}], ",",
RowBox[{"LegendSpacing", "\[Rule]", "0"}], ",",
RowBox[{"LegendSize", "\[Rule]", " ",
RowBox[{"{",
RowBox[{"0.5", ",", "0.5"}], "}"}]}]}], "}"}]}], "]"}],
"\[IndentingNewLine]"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"plotall", "=",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"s", "[", "t", "]"}], "/.", "sol"}], "]"}], ",",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"i", "[", "t", "]"}], "/.", "sol"}], "]"}], ",",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"r", "[", "t", "]"}], "/.", "sol"}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "tf"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "0.05"}], ",", "1.1"}], "}"}]}], "}"}]}], ",",
RowBox[{"PlotLegend", "\[Rule]",
RowBox[{"{",
RowBox[{"\"\<S(t)\>\"", ",", "\"\<I(t)\>\"", ",", "\"\<R(t)\>\""}],
"}"}]}], ",", " ",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"Thick", ",", "Blue"}], "}"}], ",",
RowBox[{"{",
RowBox[{"Thick", ",", "Red"}], "}"}], ",",
RowBox[{"{",
RowBox[{"Thick", ",", "Green"}], "}"}]}], "}"}]}], ",",
RowBox[{"LegendShadow", "\[Rule]", "None"}], ",",
RowBox[{"LegendSpacing", "\[Rule]", "0"}], ",",
RowBox[{"LegendPosition", "\[Rule]",
RowBox[{"{",
RowBox[{"0.2", ",", "0"}], "}"}]}], ",",
RowBox[{"LegendSize", "\[Rule]", " ",
RowBox[{"{",
RowBox[{"0.2", ",", "0.3"}], "}"}]}], ",",
RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",",
RowBox[{"GridLinesStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Gray", ",", "Dotted"}], "]"}]}], ",",
RowBox[{"Frame", "\[Rule]",
RowBox[{"{",
RowBox[{"True", ",", "True", ",", "True", ",", "True"}], "}"}]}], ",",
RowBox[{"FrameLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"\"\<Density \>\"", ",", "\"\<\>\""}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"\"\<time \>\"", "[", "t", "]"}], ",", "\"\<\>\""}],
"}"}]}], "}"}]}]}], "]"}]}], "\[IndentingNewLine]"}], "\n",
RowBox[{
RowBox[{"(*",
RowBox[{"plotP", "=",
RowBox[{
RowBox[{
RowBox[{"ParametricPlot", "[",
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x1", "[", "t", "]"}], ",",
RowBox[{"x2", "[", "t", "]"}]}], "}"}], "/.", "sol"}], "]"}], ",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"S", "[", "t", "]"}], ",",
RowBox[{"\"\<E\>\"", "[", "t", "]"}]}], "}"}]}]}], "]"}],
"\[IndentingNewLine]", "plotQ"}], "=",
RowBox[{
RowBox[{
RowBox[{"ParametricPlot", "[",
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x2", "[", "t", "]"}], ",",
RowBox[{"x3", "[", "t", "]"}]}], "}"}], "/.", "sol"}], "]"}],
",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"\"\<E\>\"", "[", "t", "]"}], ",",
RowBox[{"\"\<I\>\"", "[", "t", "]"}]}], "}"}]}]}], "]"}], "\n",
"plotR"}], "=",
RowBox[{
RowBox[{
RowBox[{"ParametricPlot", "[",
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x1", "[", "t", "]"}], ",",
RowBox[{"x3", "[", "t", "]"}]}], "}"}], "/.", "sol"}], "]"}],
",",
RowBox[{"{",
RowBox[{"t", ",", "0", ",", "tf"}], "}"}], ",",
RowBox[{"PlotRange", "\[Rule]", "All"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"S", "[", "t", "]"}], ",",
RowBox[{"\"\<I\>\"", "[", "t", "]"}]}], "}"}]}]}], "]"}],
"\[IndentingNewLine]",
RowBox[{"h1", "[", "t_", "]"}]}], ":=",
RowBox[{
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"x1", "[", "t", "]"}], "/.", "sol"}], "]"}], "\n",
RowBox[{"h2", "[", "t_", "]"}]}], ":=",
RowBox[{
RowBox[{
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"x2", "[", "t", "]"}], "/.", "sol"}], "]"}],
"\[IndentingNewLine]",
RowBox[{"h3", "[", "t_", "]"}]}], ":=",
RowBox[{"Evaluate", "[",
RowBox[{
RowBox[{"x3", "[", "t", "]"}], "/.", "sol"}], "]"}]}]}]}]}]}]}],
"*)"}]}]}], "Input",
CellChangeTimes->{{3.564274624437256*^9, 3.564274722470863*^9},
3.5642747638122272`*^9, {3.5642747970531287`*^9, 3.5642748024324365`*^9}, {
3.5642749052103148`*^9, 3.564274962159572*^9}, {3.5642750086032286`*^9,
3.56427508275747*^9}, 3.564275136914568*^9, {3.564275234073125*^9,
3.564275265760937*^9}, 3.5642753049601793`*^9, {3.5642753485936747`*^9,
3.5642753488586903`*^9}, {3.5642755120030217`*^9, 3.564275543131802*^9}, {
3.5642757403530827`*^9, 3.5642757496856165`*^9}, 3.564277040585451*^9, {
3.5642780252857733`*^9, 3.564278058153653*^9}, {3.5642782372008944`*^9,
3.5642782974883423`*^9}, 3.564278465501952*^9, {3.5642786024247837`*^9,
3.564278605454957*^9}, {3.5642788326109495`*^9, 3.564278851760045*^9}, {
3.564957552972752*^9, 3.5649575726158757`*^9}, {3.564957828559515*^9,
3.5649578311986656`*^9}, 3.5649581785675344`*^9, {3.5649582463214097`*^9,
3.5649582480745096`*^9}, {3.5649582928020678`*^9,
3.5649583123031836`*^9}, {3.5649583613849907`*^9,
3.5649584176172075`*^9}, {3.5649584601716413`*^9, 3.5649585199930625`*^9}}],
Cell[BoxData["0.5`"], "Output",
CellChangeTimes->{{3.564274740722907*^9, 3.564274767446435*^9},
3.5642748028234587`*^9, 3.5642749163049498`*^9, {3.564274948295779*^9,
3.56427496316963*^9}, {3.564275021932991*^9, 3.564275051166663*^9},
3.56427508485859*^9, 3.564275139445712*^9, 3.564275235108184*^9,
3.564275273076356*^9, 3.5642753067792835`*^9, 3.564275350135763*^9, {
3.5642755243997307`*^9, 3.5642755441498604`*^9}, 3.56427575096869*^9,
3.564278420393372*^9, 3.5642784673480577`*^9, 3.5642786153265214`*^9,
3.5642786640853105`*^9, {3.564278835838134*^9, 3.564278852500087*^9},
3.5649575786322203`*^9, 3.56495779176241*^9, 3.56495783214172*^9,
3.564958181733715*^9, 3.56495824912757*^9, {3.5649582932170916`*^9,
3.564958312814213*^9}, {3.564958366562287*^9, 3.5649584302139273`*^9},
3.5649584611786985`*^9, 3.564958491854453*^9, 3.5649585231852455`*^9,
3.570935090474802*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwtmXc81l8bxymUUEJGIXvvvc+VUUokuySSUNkVGspKMgsVJaMlkYwQiqzE
T7IKKdnbfWdvnvPl+Yv36+a45ue6zsFr725ybgsFBUULJQUF8ZUvevV5Ene5
lq07bycFRTc68mmb0aDTOhIMEpHlrhlG3xiUaKlsGMF4LrQsZd8E0o7W+5x2
jB0Y9xQpiDSR0DnpJQZZBh5gVu3vWw4mI0PRc5OshwTgc+/dkzHNZJTb6nel
VkUE1CaafI6SyEidSewd1S0JWI0tImUvkpHTcszTl1YyQPuf6q7iVTJy/zIk
/85QDop3rSfzrJNR2EkawVUXBeBEhXo5mM/WXHeUdVYCPfl0m17M5ds4jl88
oQIVf5+lz2CW3e3/cOSyGmSLPKRYwWw3rSlh5KIB9ScO+q1hTk2q5pyP1IL6
OM+r65jzotQ6m00AbAWTmwi+XKxp/PToAQh+EU0iuKzu5dRda22w5ghdJPj8
IdVP5610oKcmYJXgxNey+Rcv6MKH3ZzrBN+f/hTg4aMHXElvN3jwzRVvZ6+D
8HXXVhztf6jqYEx8YOAhqA7l32DKY796o/31Iaeda4MpDhRQplw9DEd37d7g
ULZ5w+Y7RyCfg3WDPU6erZgPMoBIQZYN7glQcEAPjsKl4m0bzBKb8OtuqiEE
/qLe4FeMYldrnhhBuuYml3vt1xPOOAbM6ZQbnMSTwqlpZgzcpyk2+NONJQrn
RmPQs9tkXoNtK/yWx2FL5ya//y4a3tF6HOraNzniy5oMj40J2B3Y5NhLIxfs
fpuA+85NvrGHexffCVOQ2r3J6aLUdZe7TSFVcJNpZig5yk+bgSbfJqfJ5NxM
GjKDJNpNFrV/Y7Tqag4HKTe583BRrvmYOWRTbHLz82uS6JIF8P2fF/KOa/BN
WQDMrm/Ef/2DV36ImyWU/NtkA5mDejdmLYE8tck1NtLnL16yApf5Tc7K9sDf
WYHo0ibX3e/0HHYio9FXWdHxcgOIZ7UwcVVpO2RRPPPSchxDoQ2f+s0tWMDm
3kfhdwIkRDPPSaJS5ITphWTnAC0yOpH5bp1fjA94WJNp2fLJiOKWp/w/JARZ
hnRLj/vIyHvalGnIWgzOvsmIvjhLRsFUIrPC8VKgEKT3QXGFjMbsPjmLmcpC
159KeLRGRjcYfSQsbOSh8Gne8evYvjnGcr1HFxTh1ohJ1nfMMU511g12ysDx
QDZ9AvPViPpLYk6qMO7horuIWX+Lbc8TD3UYtNYNX8VMURvwgs9PE8aqefcR
/lKdfnq+IQEBy0J3NsFxAh/3fPgO4D3X1Ufw4wrPVObWA1Dy9PQswSx3yW2L
HdpAu7y4TPAWX8WiCz06UCEbtEYwa4PPPtI/Xfiq770RT43vHDfGZ/TAt/DX
Bv9YCLpzZekg8P23Wd+9Rwo071LqQ4yN0GY/6OtPsa3qw7F+5g3WzfwpXERx
BKa/MW3wiork9gFqAyj2YNtgEZJB7BP6o7BHgn6Dr+1wSnbYZwh2u2g2+E9T
T/prbiN4N7bZD+/LBWx9+I+Bp+xmf66PP7IdZjGGpv826ynL6FuUSYwxMOhu
MrnPf62a+Tj8eLrJFyLjK9ceHAe1uU3mMWDpWmE3galbm8zINWP4OdEE6rk3
ObRqq2b5XlMIpdvkmOT7Lv8lmkLjvk1ujwqafLzfDHaJbrJFTMd/0s/N4Arz
JjvXvh38JmYO81s3ecvPHP/0THMI+3/9/5zj8DgibQH9/+cB8cH32bkWoL28
Wc83lHp9zWQsYXpyk19ePuVHmWsJCv9n2hq1tksKViD0//5R22acnV9kBRcW
N7lWLaho3OkvovmPwj+ovRf1Tvw54zu7Fbzog5syVUeRGGlxLvMEE3BnPxAR
yZpAAjqylhbae0H6YNSXSgYy+jvRGvJmHy84Jy/bTCeTkZDZxZggHUGgKui6
wPQb11P/sI7VEVEYPFYUd2wK6zedqMnDcEmoZ7Ia2LdMRqJSIluHf8iAKq1k
mSLuB++L0tmBP+TgCB97sym273P+WcW6MQXwyOwyrsB8jY2/51a/EsjVRb4b
xsxV59V7pFsFrFQWdeYxf7xLLqmYUoNS59Y/xLyQ6bJ1/jCjAQupAT+JeWFy
U7aylxlB+LW38YT/L05c99W7D8AcwfKb4Dulbb5W9w6AZ5LtFMHsldd0JiLn
kZ6VnpCk8BCyoY9pC/Slh4OQUd5gNY4GsgKa10xY4Y2xyIJeMAkNv9l+O0+Y
GyLaU8+mniOjlgthvk9V+MH/8BD982oy6kj3aT2lIQw3fzwZjxkho/YtZd1X
XMXhZtD26ZZ5MkpfFYoveiENmYycLox4fr6NNWlwaZSFVxUX2PpxfILyS1+L
DsiD1vv7XxOwfRIu3BYOI4qwZZegyC/Mbwo9f3L+UYbbrQceTBL9aCREvYOk
Ctq5r78vYW691kBzhaQOvi6szEQ8GvOed0dQa4G0D4Mx4e+vJKrtf3YB3G8w
qCQ4mLySK7bzADgdmx0hOPMLXV0DmzakVr2dJ/hrRLP+rNN31Mgs6CZxswfp
DnpU5bdSgk8KXdf83RG0ssrmb2W9G4QNapNKnSeQ/PkSunV9Dkg9Nme+d5KE
nkU51bnG8cCjrwJN+rFkVKw1fMCwWgCumj7OjG/D+8GpvN/3ykWA8sm3l27/
yKghKt3fkFoSJDuD/nou4b8fZ9N48JEMTN1gvL+M4zUVTzcVf18O/vPo6lbH
9gluV3x7NlMBEnNHjxVjdnP//tzqlRLIpJALBjCTh6s5OhNVAMzj6ecwe39T
Fv2ZpwZ6pVu8ifr5Omx0JChbA6SnazOJeF2LWiaPN2rBtjH6MMJ/iZlfovRe
AC+75doIXtNTKP2XMo3uV4txO6kMolIqQbXPljvgF8uOv+yDY6igNNn6jdke
8NKJVmK0JKH5xhjVU1JcYO6QYSNoSkYP6bg+2uXwwfaFc/+6PpGRS9xMkFCp
EHA4f7heOUhGwtfvXs3qF4Oi3sveeXNYjwPSvnhwS8Nzg6j1dDxfji7defLf
fVnou+E2W4brhaLv8ShKlQefaIawcGKfKSV/W03H9RL5QuMH5qRd3AmzScoQ
zqfyiYw5zfWqxuM3qkCdN8NK1EtUKevl0HfqENV1qoWYL86lMRWW5ZowKiKv
TPj727BFo+o3Aq2B98UE04edcpgkAcS/5hgi+NEjIydy5BBSN137GcLcjwJm
j6yrsdKAW+FK2GjbKDrR1elKacUMjbnHbR+MT6A9XGnP3qrvg6TA9dc+wmTE
aaitM/SMF6gpHqmIZZBR7HlSb1G5IIj7OWXs7Caj6ovcGftbRGE4HZHmpnH+
77zm/7tbCrz4Kj/kYH0Jq5lSvMsiC0eOezicx/FY1DR1c+GQhyZ2nQ5nbJ/1
2ZhPMuKKENlo71SLOaTMmzuNTxka3PxMxoj9ImN/RjG/KgT8/Wa2QMzDg2lG
72TVIdCpSpWIR3yLR8iEiibEivMvE/Uh7ub0sN0CQc90wkvCf86yfTU3swHS
hXW7CZaF8q8j15ZRdmdRZxvFMFphXr+sl7ATTI5/lH1SNo5Yks8eNj3OBvl7
gq80ZZJQR/DNxh28+2HH0hMHSW8yYqybV/hVxg/rBUG1/vVkdLBsSSC2WBhO
cjzcyTyO/Q2bOaQ5Jw7JauKfmxbIiFuzRSFgjwyQXtUlXcD9sVOSi3LnLjmY
C3EUpcL20IT2OjKLKsBQ+O0jr4h65zDsqhNRgp2S9YFdmAtsq5l196nAHSdT
5WnMW2Vavl2VVQO7Hb88lzHvVQ5JOCOpATxvLhkQ/qukv18rP6wFaA+fI+Hv
l4fJk5ayAFMXuuoIXjF7SukgeQB6Rj3HNz6X8XSicS5DNB9+vh+u7EaX7Yuv
mFVRwLkqymY68RHUT52mK/WbEVwOnDgmpzOB1G5QOJrWswNrcP6v6B4Sovor
4/LdmAdGt/FXnIjA+877VxxyIQJQoX2dS+QHjkc8u6iknwh0Zl5BjmSsZ40j
bKIVEjDgcMeDBesHu/6C0pyPDJgYFH3txPHJfrfI9MJdDkbNw3KksH1WOVzn
syMU4KLhyEQ+Uf9PBKn17igBfYsQWz/xua+7Iv91FWjJ6babxWxG7vxrG6cG
IzPsmoR+sLgdLHGO0oChU00PiPhEPbSj8XqnBYNCVIGE/x7lYYHu9gCrurda
Cfbadv+615tJdNpBwDHwwwAi+7uNnT5ECysVNhP9aWMo+YsUj/8PFhjPL0qM
1CKhtb622OwcTnialcjDeZiMFqwHbSpc+eBpU5LP0yKcP9o8a7pbQpB8vXaL
+gAZGRWSqo+kiAFVZ28pC9aPbxVDNkc7pcBeWDDXE+tH4lu+Lde9ZMHGL10p
E/dL6WIoncEtedhFK5gThO2rcmq/eyRMEb6r8IY3E/Wvfnmp55Yy0FaylJAw
M67yWfWGqsI+9f05xH56sbtYjOke1o/trtlEv7zXCZPXTNYEiY+JooS//8xV
TA+XILD8N5tP8M2923unUvpQIz+b4furfSjXMy67g5oaAtbYtvWEjSLtwQty
be1MoIvo7eQbJhCnzPgKZdle4Lrm4D+/F//+U7oBfVteCKriDi1+ifWpK7M4
KVAQ0v9cMN/fRUbjaUHzJVGiYN0YRSGG9SK9dri4uEESxov/ip/HevHrWmSi
yqwM7JJ7ZGGM/f9gdiDaZ14O9ARJmaeJ/NFfv+pHpwgTUdHz1Zj9x+pvxW5V
BrK9AN0o4Z9a7ZctW1SBI6K/kdhH/joZUKsyqYMKizgD4X+/lfN8M6sm/Fjz
HCHqQcu7PtNGEYHg2ZPJG/sIJXfIxeBF9MWq5uLdh0NI4p6+UH8oA3gt3vVI
jhxHK6+mGESbWYFv79YzHAkk1P2oUsz2BTdomRpx87nh+tlb/lP9Bj8EKq8t
S9SS0T2FohWf68IQuTvZ+vcoGeV8ORHXkSMOUmXBoxZYH6ZKBim6hqThT9uO
aS1c/6e7C50+kWVBU2n7kxnsf1VS4x3ubQqwmnVdMwXbt1Sh+JiNVglss285
/sacHb7visyKMnBZpLtMYXYpuzFSvksN3Dss+Ql9mLcSMx5j0IBT+XwyhL/B
UpnuacJa4KFwxZrwt3Fv843HPAB1h2RrCAatJBH6dh5g2J1m9j0B3y94jO7N
UIuC1OL7SsdJMmIzPD/iqS8J5aatsbW4f317Ks4/KpABi+iI2xzYXoVArqbI
PDn4xF+moY/Pmzp9dJK1XgEqviS6lxLvATpXtV2qlWC2NFpriND7E5xrMaUq
wNx64S0x/xMyYr2OtqjBr/rbGUT/3jhr2fG0WQOaoo5VEPY7hjzscJrSgmur
D+8T9g4LiTe9DgCwXCjqIFihxHTv+gwfjG0Njf5bjucvEpBkEhUHg9uXzfTx
fmeHnKQfnJGG2LiRn12431YOTihw5sjC1uFf0c3Yfu36zFqFSnnginjRH4PP
O71nG1NOrSKks8RVtxH10VejTFuqDNIGp7L/EfXT9Z3HpF4VVPaaexPzOifR
69Nagzpk9ynMEPXGqfrkuFefJlxrXNbZ2Nf03ITQGgLpBIEygst9s16IOL9D
rjbPlhXPdqMIHtsGfg0KAIdt7R8Wh1H7XZJQQiIjUOXzVGVIT6BSav7Y1Uh2
GE3UtvftICGFO3s0KaR5YCHAjK0olIxek27ffuEoACuRQVS9LXi/UDI4tnxS
BKb+3UzoJGH/J+/2PUyRAOOBGJ++Rbzfp+qNGjrLwJv3+yX+w/XHWOLxPP20
HNgfXq4WxvZNUjPV/LyhAKm8P57lEfHPaVQs8VGCrKOFmn3EvgG6OocuqABK
/khD6K82hVTJi2A1SG8x3E/kr1im4pK1nwaw8Xy+S+Qvfm/s/fdJWmBTW++3
sb85o5Py7/A9cK4rwcR7AG1HzLpSZdth7u9CZuCdMRR3xHvL2AMWOB5GvVgn
S0KBrkufT17iBJJOkkePDhmJUP0nNW3ABxHJDHxnC/E+oRtT42crBPQf+and
+snI3px3nP2mGBTs3RlQPovj03dHtLhICpxd06TMcP51t3nkN9vLgq9xTV8q
zj8XNUWisZs8cMg2x9zC9tl8rB/yuKYIGc+O8zVhfnB8v7iVhzJYSQsrEHpb
VTV/jsZXFc5Rvb5O6O3JSXoaS391uPiHN4nI/0kamnHpKE0I0PLjI/xd6rjQ
Ff4GQU2Kdx7Bk50OyZff9CC+F2a+O8X6UJ20Km26HxXAcor9d/tRpHQhl0U7
gQlcoAA5lE6glKKdYpn+eyGj+Gf2ZRYyytL7l39Vkxf06Hao9D4jowCOCKHd
ZwWhvbMqi/MPrv+DBVxhLqIwntcSWorvf8GjoSEibyWBfipQQAfrbWDrr2P3
BmTAKDZJUwf7z2b3Yb/foBwEHWK2OIHtEz7j3vZ8WQHS9Q3EqzD7Ufpm084o
QYjjtO0IMY+3JLUKTarAa/mEcEJvG3yfdvttUYfq1I5ZIv9h+/rfVVNrwtsl
vm4i//pLa8opAgiGOhgTCf9P29c7SNxdQBSN/zz9bYeQrtHEzTwKBjC9+2Ip
9tI4+iem9i82hhUmWBQpW6NIyIF30JLWkRsWhTSb/50no2cVL9NcTvLDjlN7
Vt1qcP+7xzuPWguDOGUivRrW2yqHiOK0aHFIC8+6xor1tuF3SL3BF2mwT+Da
KYLrnUaNet+ebllgfFDvO479F3hziF1gTh52UroyPsX2xbKJLWguKkJ4WoFB
J+a3qV9evhpXhsajCVsIveW5c3x+YE0V9uhJLRH9P3j+hrjMGr7fhfzgJfwV
mxQWecyqBdFW7eaEv56DgYsaBTwwR7GjIv4hGRmX+hdU9IlAiNhU5ji+X025
0ub/FZKEsWdBF5OwvvKoM3KRX8lAhqFC1Q5sX2xk3JnYZ3Jwz9pORBufF5EH
xzQ/KkDa++9OH4nzTwhc6s1Xgo/PM/wGifcG3km6wrcqkChidprQ19CPA4b3
K9XAaUHyPpEfPqGW+/fKNUBOtKiQsNdjR63Dmx4tyGXviSTsPXdWZUy6hQ+O
pblyDZXhfeqPhWk3jTg8jqNiW8f7i7rdtuUtIA2K1mU/q3A/vUhRG76UKguk
sK2Otdje1Db+k2q58vCi8vXBaHzewvll97YiRah6uSPkJ+bURS/nW1nKACxJ
KoSeTr18IeVdogqJVK6HiHimvK2K1ylTh791tP1EP1Uss+4416wJzLXvNAn7
mmtTcgImsJ6qz30kuHC76H5X52R0aeWw24RsN5Jdt59l/beOfgqMTMZ1D6Nm
A6b4hhuMcFI5lapQcAL5H3mf3uzGDmL5h69Z/iCho513XrJz84BmXgxnVwju
1/ArJX8sBMDy29A5daynr2872PceFoFDruaP47CeNlCerP0WLQHBDaPn6rCe
xhmNV161k4GfiXNvK3B9fYbMixIWcmDx1VlNgJgXMvRzI5cVQGQ02jsXs4OG
kNGcuxIcELnTQLy/Dz6XOtJxRgXkCqPbiff34u1tIjU31GBmehcjkS8GJorG
1ssaQB15L4jIV+IpiTbzh1og9fn5dcL/0PiHtbGy/9Ahew/3NLMBxJJ07t+e
69vBJrJSYNV7DO9Pjy/xeLNAzJCOtb44CenmhrdJm3HCLTt5YaUDeP5baT6i
U+eDvoUwC1SA9aUnYLXgmBB8rf+jQ4311OWF2vq6ixgYuOp3P8R6yjpjuWvy
tRRQoJNMejj/7O9Zo9dOysKn4+6JiTj/9WHqDi3n5OFzzDtvP2zfm8EvrMFe
ikDq/HS+EXPc78s6sc7KoNovRSbeV5NEel3i3VWh6qGOHaGnTYv5l676qoN/
enUckf+3Qzr66bc1gVO6jZvwN0OD+abiu2404ctwO4SyD/2MCs+/p0oF3nUa
iXZGo0hG4udtn6tMMOeltG5aMIHs15Ih4MxeILvnMvHtxvefBHvuGAleeDrE
neaaSkatBddPHTEVhG+2+4Ljf5PRjgffZTytReGjLK//HayfU0tWLraJkhAl
/SBTGuunRUCA6Mc/MnhfOxiugf39MjA9tP5HDo5WNyVYEPtO+Hmawmncn1U9
qJLYZyK+fdUdVwJB//PbCf3Udgzt/DCsAh70e5wJ/dQLOGiUsqQG24cnRoh8
/1nhpUerGnA+6VQnke+SGP7WCLoF9Mh1f6i/zhAyfovEn3ykB8fL8T36juNo
dNzmTK0XKxh9qJF5eJeEhBx93LoOcgOXLtV4vROuJ1Z2xtuH+SHlxaG6sC94
/6IPTeI/KgzPpDvZh0ZwPee9/BZ1UxwkjX9cm8X7ksaERW5gvjTcFBE258T1
/GswYwbaZSG8nct4GPvbKh3U/WpCHjQ4h92fYPu6DwgIH59UhFeadCXEe1im
BGdpT78yFK9b1hPvYUuj6ms+c6rw4DvbANHfNl6CCSZz6qDhHLyX8M+v+Vww
804t+Dg9b0LkN0xt9oz5Sx6QGEinnInD85ThcW16swjQ/j4lmoX10u5mx2NW
VknwtEs3CcV6OTNUMLc3WQY+8zKobiX2OeZGUdbHcjCgy0mH8Hntf9RuGr9X
gAHykE0JoXf9V8O5svB9spNuH6GXLLc7mQ6/VIF1wXpFQi/Xq6JK8krUIFEi
L4jIx5Vtu0m0RRrAN3Q0l7D3ZuEbU/0KPvCu9aIeKMX3fQkPx5pZMThOaR7X
gfXxtUdbt6o0nje+7iGFuD98DXvKkxNkIVm6oKAK22dsY8v3KV0eXB21XCLx
ebpxs0mjOYqw7dajXuJ9KCX4q3/xK2WYiD2zTLwP1dNGK27PU4WWvLPSRPzY
ClxpHAvVQa6pspPojz1/zBy7ajUhpchHjYifpplspFTxQUjxZdh4vz601cXq
ctMxEHfc/H+VyImW3TaXzSCGa/N9+0U7+Ve4FB+UGjlqO+Rjvc1o8GrSMwDF
6j0bnzu+pUoOlzIBIZfNnx+lehYsVWwBEf///86Hldi0cCltoKY0XCDY8chE
UJOeOkhLtpsS9p16kvshXEoOSOvyV5kw/w8nEp5y
"], {{{}, {},
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwt02WUUFUUBtCRDkkJAYlh6O5u6QbpZigBEYvuFAODUEDBFjDBBINuqaER
FOwCCUUBA/dZyx97zvet9e6dmfvuS0we3fnumxISEg77EXOCH8XYK89jMEv0
ppyXX6MqD+nV+EJeylhW6l34R/6YVEzUi7Mv1jAknteb8av8OtV4WK/Ol/Iy
xrFK78q/8iekZpJegv2xhqHxvN6cC/IbVOcRvQZn5GcYz2q9GzfkT0nDZL0k
B2INw+J5vQUX5TepwaN6Tc7Kz8Y5xVno3ePQ2GCkZYpcioOxhg/kQWThzlhr
tuSS/BY12afPpxZf6cvZJMeLycfEOCezR7wjNurpmCqXJkWez/BYq7fisvw2
tXhMr83X8oo4wzgnvSep2KSnZ5pchkOxhg/lwWRlRKw1W/ObvIbaPK7X4Rv5
OTbLo8nP5Dg/sxep2axnYLpclsOxnnXyELIxMvYx2/C7vJY6PKHX5Vv5ebbI
91CAKXGWZm/SsEXPyAy5HEdiPevloWTnrtjHbMsV+R3q8qRej+/kF9gq38tt
TI0zNfuQlq3/fzeZzJlmeY7GHnwkDyMHo2Ivsx1/yO9SjwV6fb6XX2SbfB8F
mRZnbPYlHdv0zMySK3As1sd3FveKnHFfYh+zPX/K71GfhXoDfpBfYrt8P4WY
Hmds9iM92/WbmS1X5Hisj+9PHs4t8X5jH7MDV+X3acAivSE/yi+zQ36AwsyI
Mzb7k4EdehbmyJU4Eevju4y7Rq54v7GP2ZFr8S3RkMV6I36SX2GnPIYizIzz
NQeQkZ16VubKlTkZ69kQd43c8X5jH7MT1+Pecyh+h/mU2Zif5VfZJY8lkVlx
xubAuAPs0rPxoFyFz2M9G+O+kSfeb+xjduavuPs05mm9Cb/IK9ktj6Mos+OM
zWQys1vPzjy5KqdiPaPkvPFu5TvMv+PO04Tb9XPmKsbLScyJd23uMXNwmiXc
Gu8t7nCsY3WcH5+RkzE0jf+TlLgTZjNy0ZzctCAPLcnLHs+1ir1pTT7akJ+2
FKAd++NbM9tTkA4UoiOF6USR+FviDM3E+F8pSheS6EoxunHAc8XN7pSgByXp
SSl6UZq9nuttlqEPZelLOfpRnv4cjO/PHEBFBlKJZCoziCr8BwaK2jo=
"]]}}, {
{RGBColor[0, 0,
NCache[
Rational[2, 3], 0.6666666666666666]],
PointBox[{458, 459, 460, 461, 462, 463, 464, 465, 466, 467}]}, {}, {}}}],
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox[
RowBox[{"\"S\"", "(", "t", ")"}], TraditionalForm]},
AxesOrigin->{0, 0.5},
Method->{},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.564274740722907*^9, 3.564274767446435*^9},
3.5642748028234587`*^9, 3.5642749163049498`*^9, {3.564274948295779*^9,
3.56427496316963*^9}, {3.564275021932991*^9, 3.564275051166663*^9},
3.56427508485859*^9, 3.564275139445712*^9, 3.564275235108184*^9,
3.564275273076356*^9, 3.5642753067792835`*^9, 3.564275350135763*^9, {
3.5642755243997307`*^9, 3.5642755441498604`*^9}, 3.56427575096869*^9,
3.564278420393372*^9, 3.5642784673480577`*^9, 3.5642786153265214`*^9,
3.5642786640853105`*^9, {3.564278835838134*^9, 3.564278852500087*^9},
3.5649575786322203`*^9, 3.56495779176241*^9, 3.56495783214172*^9,
3.564958181733715*^9, 3.56495824912757*^9, {3.5649582932170916`*^9,
3.564958312814213*^9}, {3.564958366562287*^9, 3.5649584302139273`*^9},
3.5649584611786985`*^9, 3.564958491854453*^9, 3.5649585231852455`*^9,
3.5709350904818025`*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwUl/cjln8bhm2SSlZmyd4rKzz3fdmbx+MxshKytyQzIYRKKFJGUSTJSiUf
IdmbMpOEEkklJOP9vj9d/8B5ncdxHncJIJ2loqCgOElLQfH/K3B9pyj/aDPW
eUeJIfv2OG6E6M0WPPbwqeBJPjvyK7z3gPI+Gkdm+CoAcmPchbjWdd2mEnNO
oMbzA9+GJ+JnZbcOyB/gB5fXGr9nMgJxU/GzPzn0hcBa63eH0qA1Xj0SHdqp
KgbuvJIMbeuAq7NIPKW5KAUPzivL2uiJ4h7/MvIe2MqBa5NiXKExMx7Q9uXE
U1MFOOh7ri44ZhNLsaMT3vFVhB+xRhJPf89gru2R7vKeyjBx8nWZiHQn1kzP
ZeFzShX6YndpDhysweQPx95aPKcGF2+3nZW8n4s5/yZImflqwDsevoW7O/HY
vfy3vBtXMSgXTzbTe+iD1VxTmxwiAbhZvOW+8J2EnasnEPNMNOGgxq3fpl4a
2OuuB7+u2GsBmfPtAfVDwpiX/knkZasNr+vjb9PKMGN3S+Wf+XjrAGtQknQJ
+wbhxm90KTBMF75tMMvNPpohLJSFnvcM1gMDyTTRpVPFqFUvIycuTh+KzITo
r76fQpTmE7PXYw1g6lL2WJPMGKLQrKMsDDcEzUG5I4OZ3Sj5yIbpUJIRmPXb
T3z9OIgC7VxbNuKNIfB25F19kwH06ZKiG37TBDycuLm3Vl4gtszbE1fumUKE
g7Skzsdq9JBZIrz9jhks8sQPx3+rRM3Bx3RFH5sD/dJWLcH5JsrnL+QlkInA
kLUsf/BpMEJRWxSeA0SIyXBz5pQNQceN6bcFbSxg+UTe2Cv6RFTbL546PmIB
XRWvpSbeX0Zpbbty/I4k8Hd/W8Oo640yQxa9nadIYJYrFdqyqIei2I8eEjhl
CVwS234KvsbokTht17kZS4gCG8UQLSdEt0bJ1exEhslHr8rvX7dHJXJVMflf
yOBvx6HSOIMhcZcysx0/K4j6YsP+efkQmjR8WW21ZAW8XM4PKSgp0VBRhDQe
Yg3Oal4aFfyjDZs1FhoCv6zByj21/OdpecLei+Bnif42sFwRwnWzWJVgLKen
G/XHBg7p/QFfNiVCu6Osl0+ILbh4OH71r5AkVFQGbuzt2UK1WkiB+YgAoevG
ZNBXjx/4nS4RdvpbnTj/zvO7O8oMsJKSdr9W+Ame3IfmrKzZoPCL5GLrQgZO
t8G7QqPEC7tGVikCJyPxU+VP9wQlBEC/tkbmhK8rTnEx6MQqLgKXClNjP/Yb
4+d/W7J8sZeAyiUtidkPJ/AEGrE/ojkysFPl9Iw9hQdfckaeEpbyILLZt+2+
RY1HMYdJWTuegG0dTneRwWVsnblZN9tbCWLcWK8/KB3BMjy67PucVUClPCC6
7F8jFp7WEyLhcRISXxeI/b78EDOgOv3pTqA6EP3ybdlM0zGKzkvFAtEEqJBU
omfNCcNonPK8+m7joDB7UGRA3hnLEmpgf9EPYEz95Czjsh6W2xJ0j3VEE+y8
HkitWMhjbFd+jP4d1wKkVPOlSJQbo7qg9NL7kzZIXRwXxI2oMY6+MJ6VVR1g
erOhNf5piaDRzxW1vKYLQqOja7xH+gjvNuOTQrf0oFdVk505aQjNGtURrlAa
wC1c2rWD9yOiMDD4dWTHAPYvk1TZOQeQTvl70ZcURqAr+NO40L0PbatKM8zT
GkPOl3A+io0hJLZinHmHyQS8D3Oq6z1+gyIYPQrceEzh+00aK/rkGvRh8NOj
0qNmIGbrtmrPXYNqm4VOhwmaA8n1pEdsTjHaW84+/ZWNCPEZxxJZvRJQhVnv
NVIGEcwO0lSQ573Rj8+xu29ZLeCyzDdlm0cXkffVnDe7Ny2AtmOqjrcnGfEb
s01vc5Ig7UbspQfN4YiZb8206S4JHo+yyWU/sEXJrdSEZm5LwMZXr8TraKGM
ghu+3XctobYkTZHKwhqNXYv/mXuMDCP0CaEjl52RdcZ4t2wRGbhF6yzoTcyR
Z+eThV4JK/B0T7huTBZHVO+rYh+VW4EBJA4HKm41vF/nCjSStQaHrqbQOz50
aF5yobay2hpIlPmfhbqOEqKUZy+Q5WzA2ykq/9pTZcKDcw7RlNU2wFl7SOTw
JxXCvna10RBFW7jVlnSAGCdPUKMnVj57aQsypTX6lZgwoVMt/uWyx0ecZi23
hm1pEJ/9/uHMhT/UsHHoJsOZ8RpcYuXvevkpFpi0/1jwvSAHF9KWt7HW4gaN
x3m29BGX8I/fRxLLeI7D9EGCEZejFy5C9smI1xYGreTO5Oh+C5xt7qu2rZE4
+KR68K78OYk77xcn3UqVhkLNkImBheO4uIwY9dd3ciB711tKKWEfft5HtjLu
nQLAOnMVdP7Cmp65KnUtKULyo0/v8LAJLOKI4KeLc8ownR1I3fn6DcbXFTxr
NKMKdhcf3bI3fII1XPnxquWXGiidWzD/+C4Lk5s+7fliTQPy3hylF6qLwUgx
8m9mWXG43doud1vuLFZ8KvKC7g0A95kDM2cMTbGkxtELtumaEJP3Z2bhhDL2
7ZHRfHOuFhiklA7sPeDH6k6rqH+7qw2b72ruNP9mwE69In8tf6IDFJJKd5QW
fhLcPy3FyNTpAqksYVT55BhhenB7179ODywbvGzipNrQsIHQEm2LPry5fSPV
em8axRrf+DbcaAANN3NfKOoNo/vn27KkkCFQyiR+5iX2IC0xlseTbUYgb2Kw
T6pxCJkUj5YZdhhDsg/r3RfZnUimrZeSNGYCSLQzJPTsM3ShmqDFNWMKGsbn
Cx8eqUHSXqK+BR/NgGN0xK3Y6jFyfihKFfXNHIK07TiO7F1FD5/XvdMPIQLw
sq5w1/si60q0x/GLCD5DyCuOKhzpLfefvnHeArxU7idL/0lC3gcKFJfWLcDT
4MfsrNIl1BcUjq9GkYDa1yjaSMsFHby7Ppz+jwSJ3Y1KG95aKGliFK1HWIJ8
ZrF9Ub8FAofWD/QUZNhJpuMUUnNGMm/8gjdjydARfT5t4ZgVYoh+qqC8zwpW
Z+m9MWdFxMFb4WF51QrEdL0nDthTIdtCr8PEA9Zg4lT8JE2EFk15yP3cyLCG
9ZZH3j1kOoLan9t26odtwJrvqH3KRUXC0C0PLp10G0iPMN2ezVQllEXCwU02
W8h44ETwsz9BiK1+nZZ+yxZ4YhNnZrjECJxvIrS/X93AtaOo5I5mtOCOTBmj
cReYILyVLXHu1gN8vuLS0C6JAyAw4MAp3qv41zKGyzWiR2FEMlVhWS0UH/ZO
uZCnKggnTRJ7K9844uOPwkYcNERhTrJCXfCFLj5G9Xom1E8Skv8m8nZ6yuCP
dkRyXhbLws5DGktuWXb8SSapz3dAHoSTM3aqDfaw+GeNpeLzJ6CHcJXjOPcC
JuV71NptUQmm3nzt0+EfwMqeB73n/aACw3Rl6I7RC+ydmQgt48pJOAOSEkqk
e9hIRB9d6Io6tC+QZ1wTrmADNUUzabQY2NPfehzPG4xN5NMwfDgE0PhxKb2h
zxZL+LFdLXFQE2YCI5xerWli5W37u/qOaAHw90f78EhimGdWlzyXNhjRdcLW
bTbMPP2EypKwDlC6Bg3Ltu8QRmZyzWQVdeHhE3/W7px5wmGkH152Qg+W9R5E
Uom/IPQc288chekDs6/C9c60MfTLxCCdSdMArAunWA59n0KxOnt/7dUNYePq
5KLEix608GmhXlfHCAIHK/idrg+gOkODsys6xuDiVQzO+BCiUW2K/E40AYf5
f5o6LY2o3YyCEOZkChMlh/iEJ6qR1s/Lh1Ud/8u/TY++VVw1UuMWIJl7mAMx
0NPESyIf9dHxNmYoEqEI+EsOGEaiuvnP+o5VRMhJ43HeSPP/z9d0ds4pWcDf
nATlDOF4dF3mjKjvMwsgSVPsL1JMRu+TryQdVCeBEPj+ePwsGJE6/B3SG0jw
q4rlfGoIEVF9Z37486QluD5eszC+oYvy3Yx5ChstYRRRH6QutUMBzVq0MjgZ
mJvelNbfOI1izugKe7wlA9XNM/QJ/nqog70896yJFWTtyhRq2/GhlevPbu91
WwHne6L7xUP/GsimXX06FtZAPn7hqMA5SrTy484FryFruKDsQknAJQiNYRLZ
fuY2kKvEq/OdQZVwgWTxqXHYBkaIbqwbj5UJnixKhT8sbSGK1jf5i64sYc1+
ZnhiyhaeCf2TGnUTJHSkDRn88ejH6ZwFF37bv8N1FgJbn41Qwj0+rquXNZ/j
2ztHYm3tD8O1SCtUL3wXP+H1av+eARc8bOLjuzEQj9+/5tHll8UPYanYnwO7
vng99lXT9K0QfJs3Hk8WI+OuDjVT6c1ikEXDpXvpCQHvu/Yo1pRWGv4MlG+E
eAnhHVmOA3rZcsCxuu/ShfdM+K+c/b9ybigAaZWO62LIH0yYQemJa7kiJE2I
l9utfcD8A/qLbB8qQ7jI1LKtSRv24+tbrsm7qpCyur5/s6YSO9+rIv6+Rg2+
xC1tLV3Nxjq+mhnFV2rAtblfrimXL2ER1/79WB7AwP+u17cnyp6Y1NqEOFMw
gCdnhZBRMRGLaUiZi/DXBL+PkflSQicxstrmknCkFjw+wyKW1SCApVuqbj+K
0IbbEnceJpxmwrbk35gxpelAdnE80Zi0Roj6Ro7tv6ULVO/fGI/1TBFcsJf1
Pll60KNu7vPx1SvUoz+pqXNfH+ZjtJNwsWkEvc8WqgoMoFQl1L005h0atMeM
3e8YwuPUL5wir7rRfEHCtnWJEfg9bKVnOD2EnC7m1fE/NAbvJOLLJawXMQhR
1nTUmcBppzXj5eU61Hbf59pgkylYRhvxTfyrRt8aSoNrkBnIRAyoau5UoPS0
oqLBLnPAjCt1Tm/dQAK+kuEtrkQo46O+Hn8wEKVwOMTin4jwov3rwthqKKq9
s/Rkwt0CklTJd9p4ktB1/7nTs/MWEBRkuC++Kh5J82qoMviR4HjnwGLozlk0
QB9YsLlEgrl3NMcpDusgqs0IdyZvS5BZ5Tl06ZQ5WmuZs8RWLYFD49I1LOk0
yp18/nYqgAxs2v7cOYa2aMy60tZ7kwwzvtWKNXAS8eNS9qXRVlD8IMu04gc9
4g49PEvatgI5189F7GrUaPSy1XWeeGvwtTxclKK3qnF7JLS7gtIGyFUsMfih
E4QtqzfiarS2oGU5HXqNoEgQSt7Y2kmwBbb3caORg+KEXV3FxtXC3zjPe/nC
sea3eCONsFqTDSOQmHQZkw49wusaC+zLyOzwK6ZH62FWOr4xkHHSQYYP9rRF
i+L3XcBv7edrcK4SgHr9bQUeijO4b9ZavEijCFDKqbM+zTTARSOvhFfMScAP
g2lP7Iw8znGppC3wqCywtBMvmYly4iZbSXe6b8iD2QBDSPgTSpzic+43/N4J
uHDwR+pLjUXsRuOP3p1HSvC++u6hEhjC8g8dvf0nXwVWZSsrDANeYSV+4Rq5
ZSchosDyvGBTEXatkeNc8lN1EO5out39Pg3zbMxosWkmwN7A1t07T89hU6bD
Gq1TONh9Gs4TtXHAmFIc3H6uAJikxOiFGepgRa2NSj//aMJRXdJYbYUM9mhT
c/DmXy24Trs4FSh3BBNLUNwso9YBlS+Hj6WvUGChNvueXWPRhfInymeXU74S
RNoqI4zZ9CA2Vrqvk+otgZjW3n2B97+8h1tecWx4h74Z7y0YHTeAmwK1NhWk
//LfJrh+ktcQ8rxc7Mr5+xB7u1PQRQEjqAjeFKat70enmuvrvISNQaE/a2sF
DSEO1sSkVhkTuGerdoreuhkVlp9Y+a5hCtFb6u7fWWpQ+qcax1WCGfD32DNQ
fq1GG35Ftq665lAel1N0quoeWvg7fg0JEWGsO/qHRnUsUsi1XRq9RwROapoD
bUd9UUP0bHSekAU4Vun8UUi4hLY/RIU+e2AB9yTrzjyuSUb/FvuqaSVIUHP8
Ct9JxvOoKljeqrScBFypgtxpR8norW1K0CkJS9DyCW0eTNFGnHErb549sYS5
yKC0nym2SEjMNdJNlgwpTI/LhEyc0bfw6xeMa8kAo+/yl/KMUM6P8qbb6lag
bCUvxRIqiDarT1wlN/znS5UBxGcpfxvYBYIL9MEaVJeZouS+06BvB5o5rrZY
w3fLcca1Y8IEPhs6hk+4DcTUSb2ID1UhXPtsLTvTYgNNHNHl/WQVwovl56/k
df7by0ufJrpvyBGy5Tm1DLptIfP5fXEuaSFCdraZx4+rX/Dv97qd3JX78Et/
jPbUOOggkFukq3m7Ej81PelHacsK7189Ub6YfhNn5yu5/0SdB3YvHvH3GY3B
eU21tL/cPw4N7EqnL8m745leK7Mvm4WhSI59wT/cHH/rc/TxsWFx4OVi/PPm
mTL+K6lU8ONhGTCuqqKJcT+Gp7T/UrrCJg/ytNZPWLPo8L8ES39frhNQ7mx4
ZiRkFbN3zUBykkowtvdy32GXUSzx9fmjJQIq8OpqtF59dDNm/PjY43rBkzDQ
kKC7UPEIo9IrMXsqrw6FG/R/NfwzsJzhwMTvqgSIUbhUI70YgUn6e9was8ZB
bprmruVlF4z3NU97TCWAS1XBX9sPhlh6k9AB5WpNEDtP/Soq5wR2kStv2rFe
Czzn3w2LJfJhUz7KERKvtcE3/NBm8Tot9kJn4ZJpvw6UfK8qv5H3g2B3FrSD
x/TglX16ytb3XjRh6JfV/kkfSMaPi7MtP6KmTt+VPx8MIEjniYJD1CAKuWE2
/WvaEMKUJW7eU+9FBZUFp5wXjODWST6dtd9DSK3zcEHwN2Ng4qqo4PRqQ2TP
XOLkXxMw4b5EoDGsRRIth9ad90zhqbPxBL1IDQpkf3g2ldIc5N/vneTrK0Gc
M1ZyrLRE6BeTGtEfSUYS1FoKvXH/9X/D1dQdRh80plF72CDRAnJuR7/7nZyM
dMtP7VujJ8GUL2PGnZloFNBepzeWRoIvucJd00eckLFcNIP4fkuwfeEsSzOp
idL+FHkMHiJDCNczZ7toZ9RSJniMLYMMbfFfBsdxElItqmCi4raCiZRD3q4K
cojpr+y53btW0ERfVGP3bKchu7sZTI9aA7je8Zizo0MzhvstR+5bg92YnOgO
NRthkr1xsZ3fBj74OPCt1igRPkwlHRUstYUFaWPdvHQRgjw0dyxG/MOdf/vR
RPa/xrdZ987p3j4IwdVGvs337uNsBa6GlhZHoIFqatQs4wo+nhAzwHj8GNSP
WoXENAfjzF0bihOvBaEglNxy5JAdrvd6SyizXhSKbh7b6JXXxlNS1vQJ65JQ
436kJDRXAj9KGFa8xC4HxgK7rAFZLPhBaT7Kg4cU4LK+Bgqr+IfRJc+6s4or
gtgPw0A1+89YB5fpdJeYMjQ38jZ9ien5b2++ZdXhUQV+ytZDwvHPMGq54d5w
eTUIbE4Pdz6Tj3GrJN4+I60BJcaxf5hvJmKqj2p3mw0xuKXPyDn1IABru1Xw
00YeIDk4eb6i2grbJudRuklrwgF/TfsnkzhW12HrVH5SC3Lc/8xpLYhiI6ML
b3pVteHWjfcCkros2Pm7CuclDXTg9Id/r8bwLYKCt9TXfitd2P+mwDR0/2eC
QY0Z+aKzPjCm3GBTXZlAW6Y3TvY5G8AnN23hys8TqKllYifazhDKbzosZMx1
o+fb/jRM7kZwiOvY9Bg2iFiH3rGGuBmDoOxH8u+gQfTLYSuAI8QEuNuzyS+T
X6GV1Z5UilhTmNvxEObrq0ZrFLZvL0aZweVM4kjdVBWKaSs0Vr9iDsGH9E3U
Lt1GdaxO3z10iFD4HF/y8DqP1OQDimebiOBg5ivbThWEyFUPokz0LGChLfGk
1uMEJHBRd1zgrQXQK1d9odNJQlWiOg5HTEjgfWyOoiQmAP3iK3e90EMC8ww9
1zYJY4Rfyp64YmgJj9+fe26iboCE9EuDU/otYZW712X1pgOqtxO4vmFGhsDQ
pqvHFx0Ri+G/wbj3ZHD5kEa/VqeFNtUfn5N3tIL71BZa1G/Z0ZyvysOySSsY
k+KxEOHebci8ILCr6WINtgl/KKfFNxp+D56f0J21hoNZX64/aZchOOWYiWw5
2YBa3l+eRV9VgkuzcWbenA3Q3bMzivumRLBILD3rt2ILSi8KMhdXBAhtckEe
dJ6v8a2mldAVx1H8nEt9KLmVAizfPrexP/ESn6Mt0ZGZYoZ+kYwTLDz5uFoU
hbtlDyfsjDFz8nNfxmk+yvn2E/nBENomZvf8cZrah1wKiUJwMN8tWjfYCk/J
4RSXjhaD2Y5QOp1XOJ4wsHhEvEUKfl2/cutrtzDOabCpvB4mB7v6RhUypQfx
yqd/WYoDFEDiu0UC4/I6ZlvF51WZpghz/hgDNv0R87wjTKubpAw/OIIfJle1
Y7YXApQEI1XhVjBPLN9oFUb+MfnxdJYa5D7Osprcfxtj89d75XlNAwLHD60z
Po/Drt1ypgt+isH2peXx4zNeWGBzSlyAC0BbdWpSxWsLLFzdjTvNURO8Uzje
H29Qw1KozH/Y+GhBQtiD4yJTglhU74upK17akBKuZFNhfBAz9tui1IjUgfMz
ByOYH/0h0AXso+RP0oUIsTvLriwfCcXnsjm94/WgUOjCAV7mGmTQruznmqEP
ciPPDylf+4D4D3bM8acbgF6Lip9s+3vkaRYhEp9iCE77Bvmtc7pR3SThQVa2
EWQbzPR8EhhCX3xjNZhuGoNFsGT6wr5+tDaUu5/moQlgqXR4yYfn6HtuXK10
lSksN4UoTi9Xow73xVBUbgZ/TINtqrQqUenFM+e2682hZZ1DI7MgEzm3zXft
2RMh1fh9VCYehIoNHzcNjxEha9C9wCjhHPIO+Uvec7KAUMKfL2dzEpHAsaQ/
jh8sQHl5H9/mWAJyCMpTZDtLgi6b1IyThzxRNZK6LDVHgvS5hnncShe9TuxR
yz5jCc1ZkyxxVqbIbH3t2qmvlkCzRef6ecUJDXxx+GHmQQat1tnzmVunUNOb
Sh+aVTJsveEVW99QR8XEpZGhUCvIqWuZWb2zH0WcPms08tsKEq6O97rzU6ER
Nk8n6UhrcL4sdtV8d1DDoPVz/Ne/1tAYEj5YhisQBN4kqNjv2ADy10i6XaZI
WA9mrVaMsoX++CHixCsJQjD9jcjgsp/43yOc9fw32/Efsf5LTvr7wKrRaYSF
9jFe0CbDH/uODdZUGNiu1dzAdz+PZlZW8YK5trfSmaxwfNN+wbHFTwBm5465
IUkXvG5fjf3+iyLwzNHlzYqQEW72fOWtUaEE5GpxLqxkKeC9LV8cTSZloPGF
Ul9SFxd+94kAVWSwPHjPFL742kWFN/5N3m988QT8prvU0LbxDWv1GLtilKIE
7LFn8/l8hrEc9XNbny6qQLXcU4opUYQx7wjYziafhLvzflT2Rx5gPjP1Eizp
6nDLQu9E5KNrWK12yglCAQESUrxENo6cx1atVC0NX+EQkcmgEmHshBFtmBD5
I8BEhNwP2hu6WL52jg/VvCYsP+X6MmQqh0035Yy0L2pBywFBEf4CTox9t3Ys
5Zc2fF1IaieyUGFtGeDZSKkLxew5Lp3i3whVRmf/3aHRAybqw3OyVzoJGS0P
xMP268PE1pJ6lfQIcu+2Y/E+bAAxBg8dS7umET+98lopkyEkaWvmebP0o/KX
uSZ3DxtB40KZwnOmfiSxE3K6hN0YiqQ0lG7MDqECCq8DlXwmsJQbsn7kXAu6
at+rJy5lCly6j1k3dGpQ3XhXvImMGXzP6Qrvo6tB74/M2fspmkPgAEuyFU8R
Cv6tujHES4S7TeU2A2FxKHCNhcCeSwSFP+H1/bY+6H2lfKM8rwUobrGO2z+J
RTXPUngn8iwgv6frGH1PMrqWG2Rw6fh//rLPqdlgPAzNtGdY1BaT4LNCg1Ke
rjV6eFDf5qqAJejJkakme7XQO2bIgQeWsLj7SSko0AbF2IfUXRQlw/T+zocv
gpwRf6LiffSYDAbP7kg9mzNBsQ5Phk+csAK1AZKsCIig80xt9J0q/+U75r45
1Twt+np3eD6/3hqyXZDwEysBQo/euSZuVRsovlESwcSlQpg5Fd5J0rCFAuFJ
GSZGeYInr+LfmBZbqJe7y2D0RIgQw80w+6vwMy7H7mRlFjSAVwdlVY7T0kL+
fbKsuX01rrXgrTA6xgLHlhRfJTFm47xyy9uUr7nh5sVB6iDvWDwmb/+8wenj
4F8i21X40gOfmi6vz48ThlPrcaIfERFfLonfeHVNHKp5gH1ZXxV/1Pm1vr5P
GlQ/Ps4ceMOPT0Rcvav6Rw6aJ9+mFFIw4C/ImtfDNhSgNoptoCjrJ3aNKTI8
er8SLLNORiTMjGGxSz0XM6lVwDU4eGT9ZgtWq9bZRkV1Es4TuUJqfR9jHz2M
aU+yqMOpXTuKvpxMbM7Wc2OIgwDf8olBO41RGHa+p9xRCQdxg03urH5XzDmt
+9hsIcDz1giNzVhjLA3f1dK/pwnn6Rjtb00qYrWffn+6XqYFvH6eu+IOx7DO
+To7xgpt8Jz8HcDrQ48xxPIpRTfqgN33oEPzA6uEu3eO3gts14NrTNyawvxd
yLqz4/XCoD78UoyMslL7iAoZLy679RnAxa8+ZbTkIcSYkPLyUJ8hpD+5iWIa
elDgsi1V/nsjOK2g9Ghwegi9qVsTcho3BhmvBzbdP9uR5hHFRxyLJvBJqezk
u65aZPZ6wfPmqhnI/BnHfZMeIYMaL3mLLXOwzv2txWCRijwdt+MZI4ng/O4R
/YsKH3SdreWFb4wFJL/ru37OIBmNOWznB+z95zeTh5pkmWJRPi+n9lwCCfiQ
ctj2LWe0WddhkHiFDKHGMxdmH1oifuanx+4etoLmz5fDek8ooBuMP76L3bSC
WqHlT7mOFGh6jE9gIdcarEVS3te9OEhIoZW6T8llAyLb/ygZjisRKtqL3p0q
sAW6FNKJOTNRQjHl0USfhL/4tywtke/SzbhUuoHIXPIBCJQQqiAeKMa3H/46
ID7EASdElsKOGqTiM9lvJE4XHwVPyqU9wloIHszd/F49ShCWCvIS0sbt8XTF
l9thkaLAyxh19re3Dl7VdiprvEoSNA4qvbIJksJ/vVqgmP4iCw4FZq9m3rPi
TjPPPdAPeYhIPW1+zm8Ha80fSDpKrwjqzT89Dzydw7ZalHKP7FOGjZ81h7qJ
fVhlKk+o3LYKJGdsq3yarMN8X0ctNh9SgzpP24QnqoXYhq0EcemABrBNtO5y
UCZjCTLlASWiGIQsZm1TrARiA9xDUbn8AB7xYuvUSjYYCJjatfNpAqWB2LFr
XwCr0lfIQmJaQMUhzcUaKo7tKyS074lpg/KY+sUPHqyYC0eN9X4VHaC8txUy
v/6PEP/A9c6eji5UsDO8TtSdI6QsSlVHm+vDAuLQpuoZR+d7l92+Ew2AKuCo
qv3OJCqk/v5v2cgQHJ3pVmY1epDZj7KPYlZGYG8Zoyr3dQDNHbjeUEk2Bp29
0lyez4Mo0ujDH0oXE6Dj/3Ykhhkhnr/+S5O+ZhBmxRJnplqNHAI/meleMIfR
zWPrDq13kNYBZBWtQYTlC/r3xGsuoAup4uvHXhLBRUvviyl7IPKfW2uYIVhA
rr5YxovleDRmbRmrjyxA1/surtKdhHTGFn5maJMg0fx8O7VsEBIIOeCc0EqC
qIvsQ/mjpohfhdyENC2hBSdrPRjWQ5/e/dVjb7cEsRXPIqqH9ojtmkIXiz4Z
1jSvHM0bckIGVoX/1nrJcJC5aUomQQfNlD5O6iZbgdMRlUvXLLmQV8z6sXfD
VhD45lTvPf/thjWX8+5vTlmD2M73CMbKnYZV5yuPbCasoeViROTCTSnC9ytm
ajJTNoBpzV1mNlMm3IL1it9XR/EZ5wOEn2nDOKu3KOuReSqwzyzjUY94hrOu
J7+KmzgMRRMLbrOkXFzK+nxV2VsuMCVc5n10Jw4HLF+MaYwfrlxq65vl88ED
O8W9PWmEIdy9sv7xLxKuwW+WvkYrDu0iVfRdzer4EVOvxSADaeh63sp4nU8Q
v/CpxSu7Tg7G9gdrk/j344pxfINXa/7rT7c1BhbWNeyXk8lPjh5FGFc7CD9F
p7Bz2uFavm+VoZ42IU9luxVLPMW7m9GoChVpo6mM+k+x248zg02G1UB4SoW4
T/QWFuVqM543pAF9d/XtAj1iMffEW+MevzBQnbs8n1vijn0VkRwsvQTgT3Nq
7qayORZHqKrwuqgJMuptJv73VLCi8ec8XanaMJesevxaNiNWtp56Z+quDrig
+Uqbw78JxI6LHMLV+qA1N3HsftE0ck4r+qr11AAEs1vTGckjSM2LKtCw3BA2
6VXjnf52Iza3bkuu50aw/PxpmGzmEHp3o15O6bkx3N7bevHFvhvxhHrYXW03
gYXPzhyj3HWooIzsGTpgCrtzbmVFjDWIidUz4FGfGcy+DJC7q/oEXbg9P/Rw
whwsov9tfQ9IR3OHg0Qv+xBhZWg0KUjBH+28+UBX528BR3/ZnopJT0K7bRPq
2+dIkNwVlNNz0A39sKh52fqbBPeVLtf5c2mj3OsjnIKblmC3tfhghMYZrS4+
1R2gsIJsf+V+q5fK6I6036foy1bg3JvWn/6BBr1MQeGDNNbA9mB7+y7QoIl1
raZ9qdYwID7tkn1/R+PbHD1r6D4bENtI82PtP0Eg5U1qt1+1BZHnZW4hq2KE
1IRRioDsP/h36vOOqUda8WvJk8nHPfeDMq3qT13NEryvjjl7ZJgdRFxcjXUO
X8dXkelH+sd8oPn85PzUrfO44itL7r01AfBvXffs6D+NCxxTLtfcEYFI86tU
Odz6OAUuJM0iLgmkkiy9GE053Bn3kL15RhYuWOeer83nwLf1vivyVsnDvkcd
IZuFFLhWT3mn4psTMEGxy82V/wVzYqdnqepUgtM//bwZxAax4s/tKvsaVWAz
/hhb8cBLDJvu5yf1nIQswxbtuZj7WNXdYLTbpw4u0e2NdvypGO/JOxbBnwmA
Ipt/ZjuEYB26/iL4Lg43SzZ8x+LtsI6G8v3XdgGK5bN9nNi1MQ3G2iPdbDrA
RvGDvIX2CEK8Q5SnBXWBIXXYSIv2C+Gg9aJfhLQ+yAjLRA7KjyL/nqJXkvIG
4OjqTT1c8wE5Xz682SVpCB61v+CJbS9Kvkx5UErJGNzzJey+Jw2hq8JFlOO4
CRzw8RI5QGxCtTu3b1abmgJrywmBzR/VyGvxrkqWqRmw6XY7n26oRnYUkSKl
UkSYfbVh/Fk6Bq3h+bdzHxGBcz3TYVneD+nRBfSlP7EAoWrC50eXkpHnSSm+
FAUSHPwgof/J8hyK4xuwvVtLAjkK6jvRMiQ0d783VVjBEoYs3lkPh+mgrjUS
GKqQYZ6iqluCzRnV+n0zG28gQ972tyP95w1QZVz5kz1tK6ASn/Ye+MGPalWO
qrQZWMPq0Yv9z85To+WYCv3ILmv4uN/Q70OLKMHfI8e8RN8GOB9osti+VSG0
6W2nHB6xhamb1JmE14KESraby0sRS7jjgT7zns5u3KnrsFfpT31wClaqYzf6
iGaS7VYTvxtAntMw35HGATSVz0L9/ochBF+gX7nf1Yv6zi1UCqwbwZpcSZvH
9hDyeEmxWcBoCtszXWTm6Rq0XJo0cz2FCJBeKVfY6Y2On33oLXjNAl4QTZuM
nyUjCLYy6skiAZtdqLEz2COicsVfDwEr4J2j92yulUJJXwK3houtYGyt2if0
+r+G2frN+JYya5hZbGJ/I8JFkBIhyl8Ts4EMnVDPh6rKBKZfPh7Llbaw0Fq+
KzkjTKCLDlyPDdKHEfzUJdFfkyimT+ng00ADOH/kdBpdwTiaudSxd9DbEH5n
t/tslncjrVu8PkrnjEGo47v9q28DyDvgyDeJBBOI49mL/iJWj3oZNAsE08zg
hpsa9ZZfFZLg3ksgmRDhmfmk0l/Pc2j7fS3pQQcRAjs3P2syBSOfEhWF3R4L
GHf84FWdl4haOow8fEdIcGdSjKd5ygA5dU+sWVlYwmOVUXjSaojcX1jn743+
x8PwM5fiQx2RsygLecqaDLUGYam+ww6Iub/sTPZZK+iT+FsuWMqCnO0vSMV+
tgJ+Ldbxje29BoVrfTebvK2hcnLfekPoSkPzhYpiMc+n+NWR2uIj58bwNP7T
fYIaFJDCJ0lRQqjHx66siNy+ywx0nh/xZMECvJFWMHPnKidA8m1FzvbLuGIS
O4FClh82PG7TsdEH4qUrly8XuwtBhr35I8pj1ni+srH5PzsxmF3iqxn1ANz5
55XPtwqlIOzgmZv7EkXwpnu630w95cDwwbOZ5TuHcOZXgUWPnBSA5WV5yULV
BvaTlqX9fZQieN0i054NmMHcqwaUXoUpA1+O0e8Q3w6MCXS09b1VYXnEKIUx
sBrTopB5VZygBq+siuhmh29j9XItIfbRGrCcSMGebROP5XBn3qjNxwAr/rEv
/LU3xlI6K3PPDuAyJd9jPTIJW98iPmOz0oRzLVwMh9PVsWL3rdVAVy3Q7PuY
zh8ihAUz/u3566wNPy49P21tfghbY7w1ORysAwIxV+/fc1gnmFS5ahpc1IVl
5qe1I1wzhLfXE/UXI/RgU5a5xMDoCaJLmqMKSdEHXr1KBlrlDyhdZtDVNtkA
uvR43tCEjaJafc6sw3GGQHfan4M3oxsVhOeG11w3AnKUw6IU1RAq/xy0dynN
GDKmfLxujvUj77Ivz+/lmUA/89e0oy4vEMWRa+32j0xh8YV4Zfd8NXr+6kqY
VrEZJDI9F6LKrUSK4U1Th2vMQegGr49/VBa6+Lsl5KgtEbxz5YRO/QhCcp+H
ZvJGiLBKhc9efx+C9E7Hlh62twBHfQWWYotE5BNm28QzbgGJD/MS9xwuozWe
PgXuMyRglaPTuiHshe5z7r+gN0MCwVH2qzi7Hjpq9qfYwdESJL4Ka6g7mqDt
Q9FSHnOWwPfYXYn5jhNa8pOlYXAlA6Vvtq1Rjx1S3JILHVwiQ2uOaeYhdwJy
H+86/DXIClodeJ0knx5Alg4/dZXDrIHikMHVXLuyBv4GbmfadWtIcWvaNz4o
T3DyvHdgIMwWYsnkUVlFSYKUJ2534ukqHlna2Xz9egfOgLPqyLxmgPWAp7lN
h8rxLKPzVEs32UBXUOVxklEGHue31WQXwgvtRbc/eblG4GI03TK/jQUg94TJ
jF6nC26vk9EefVoEuLMjwzk3jHAXq+PLnDESEP1v6nOp8gm89HOSeP1LGUgz
6q4sfsiN69AHPhtykYfU96krls7UOB8txV2i/wnQvfy378rMEubY0PMlMEIJ
at68uO1NP4LdtDgmaRuoAogprkXvYCPW2rpxlu7CSVj4oxZT3P0As/vJRGcT
qw7MElba63nXMTs6umXZawRY0bVIWBg5j22Ne0+nluHw9M70wDzpNPaL4OaV
OArA7vsskEtDD+OlvkyXOqUJNoxoLK1XDvsrpPN+5pMWtDrOhIk6c2FM6n86
tha1YQVXwlseU2H+clOa1n914Img53oZ1RJhPe21euu2Lszsf6Nv1d9NqP2Q
ZbtOpQ+qXrn7jR2HkVAqQ2n4PgPQcM2RdN2cRumieMgMjSHw+QXJPLzTjzar
Vesl9hnBz7aNZ7de9CHWOwMn2g4Yg9PjYFLDyhByC7MbEGY3gfuaDKsRTG8Q
zQnu322CphARrYDjbjWoy+r9rQgRMwgyR6ufD9egt74Jjr8kzcHgdaNszGgR
Ug30fPyBkwh5UBw57x2Prh9fnT9xiwgdOUsc/mw+yMXo1PS/IxbA/ktkQFAz
Fo3+0mLF+UjAKbEpptP1334JfmNeX0iCHPrlaid/G5TtJWxUctQSWke0OOLi
tJDimCcV3T1LoBb4vKPwyRqpM9E/P19ChuoW5RCczQytZUh/oJC1AqlI9pbZ
JlHkPufQW1JnDToHR8ITPvETfk66FZwr+4TvpOfsXRAfxLtkT+57FE0Dqhbt
MkocNbiydzWb1m0W8BDn1X0ymI0XvjwoUR7LDWb1+esHVmPxCt3VZ+GE49B0
Sq9sKMkTv8SVJnLYVRiU2s+aXdOxwCn06vhSfMWhwFpwyYLnJJ7wLTlR7Ik0
9N23ivpBOo7HjUyYp8/LwZmW3dGg1wz4EecXx6IXFOBmuY+wi9QvTPRMwGjR
P0WwXDme9Lp0HIumvFC5b00Zrt6V7b4r9gbzpMofEfmpCo9JVBv7bMqxvgt5
M9FU6mDIyxG3oZiFpfDMPX1LS4CpTcEC7bpozGBrV6VQCIekfrm1s9FuGFWI
YA/rbYAEC+mWTGcTjHHv24f+HE3gkdduPnJZCevXLfhcel8L/HY9+f3E+LEm
gQ2jiAfakHHW1nxRigGr5Ht5kL1OB/I+sQsN2/4koLoDAz8a9SDy196heYkO
xHblvNV6pz4YD2JiO/wfEbNkss/dNgPQuxLlU985hKyjxUTz3xrC7SSuK+3x
PehUyrP3DX1GwNXx6MJU/xBKll74MjxgDMzueQqMzR3oxUqsVe2MCTS3kkyW
2Z4hjC3Jm/urGbyK/KQ751OGRso7brv/Moe03Ntn+++loaSUcSnxMCJsrVC/
PCbqi7wpDJkuRlhAkswQI9vxZHT8Wa+X9D8L8NO8dsmoIBYptxtofoslQaFj
pXia5xnESrjjiF8mQ8e3TNWjZ8goS09W3OuAFUg4G/Ok6pxA5Mw+CeyGFRDr
7OjtTSiRI3/Jvh+3/vMdkcpztN6MBJcdl5ij7DYgylnJIzCkSHBy6XGTurKJ
kyu0bNimmnEds+8xNRQHQGGgdUG2rxhflVBbzczgADf3Q3PGMmm42/EFm33u
R+HPkLMK9c1z+P2WByW+doLQ6ULnnNjlgCsG5Hh+sxeFY0d0PNcP6uKtbmn1
Jdcl4aT50ZPTntJ431Rij3GbLLS8PuMjGcuG06nR8rDPyIOXxb/0AJtdTKhM
n1No/QRYRwrM3Dozj2Uekdgk/FUC+xRjzgGdfuzJvbYHD5dV4Gf0j+de6c8x
/iSLjfndkxBKe4NF8E0htuAVJSm3qw7uMb73X79OxiR+iorlcmBwK1/q7nRc
EPbOwadF/T89uVJ8e6pw1wabZXhA95f9P57v46St89HEZOiPyerxa8FIf6Cu
ma4EJkSz2MF5XBsMgh0zd3dZMVvL33qRMjrwWy97fHlimyClvXcuW0MX/uRM
plEtzxEe0/B+2tbTh1Ip0wO/ZMeR4tXcrW0DAzjqcD/bNWIKMQutRwhpG0Lg
8a026vM9KPergelzEyPwSdrS+/NqAPU+XkrvNzEGmhXvwnbWIQTNfS/CTpmA
6MsIToo5hGRCI3b3nTWD5Hj1Ru/T1cj47SViV4A5zOZncG/+uosYUpb3JasS
QeyT3lnvtXCU/YKuAHtGhJympKPdlgHoXzTvh7KTFiBSn7nMlxKParUmhk+9
tIB6EUN0nToZXZV70UDESUB3/2xKwH+8r+fDnFKaSFBv7Xu8i2SODmco1IwS
LKHxTALpprQegjyBqokWS5i40plYKGaPWkrqPfi1yVAo3vn9lMpp5DoZzFzR
SYaTAfb+paG6iPxR8kYO0Qoosg1+RurzICHbQu4EK2uY+1zxnkeNAgWUstxw
fm8NibxeeuoBkoToN0/FNwqHceKs0c1Q/xFcSqzHZ8GBCmrX2RWHbetwhoiB
vO+5h4HpsNCkuOodvPJjzjQ5iQsKaPX3FZyIx4MW4v5q1PFDF30Q/nWfL85O
s+/9xJwQeHdeClFItcSJjbF1LZ/FwPsp7c6VGg38l9++Zx9FpOHdISJ3W5Mg
zq/OzPfjoRwwsZ5SYNrcj2dezTqTeV8BHAf2FTs1rWFpNWBOaFCEsJSd9xmM
H7CgU0Ihs8+U4cIhrYGGe28xteM/9z9/ogqN8bad4eKVWHLDvOmNN2pw7sSy
ecXoLUxAZPhGerMGMP9xV5dci8UCGTvdyj5hsBZirLPj44GlHIiPnggH6Oo5
Mx/NR8QodvIMf53XhHDfuRwRZ1Usl2Nzv1q8NhyYi9ulCtuPkbJEQ89l6cBZ
ujDzI09/Ew4mS5crlunDQqyJ0HuvaTRQXu4RUWIAQ041nhXzIyjRpfDDaJEh
5JqneVCPdqNtSjEx1adG0LpkZicTMYRyGcp/l1UYw2G9J+plhj0ow6SWTvK1
CbBsCfBKxdahahHnyKIOU0ilowlxoq5B/zwYgmTazMDdXDbntUAFosnUUW8Y
MgfD0w5J4/I3EDVuNlDkQYSR5YHGY1wBaOypm12HtwU8XM3OOGCbhG4F7Vej
CiLBCw4DsV7Ks2im/Eb59A8SuC1d/jKRqo30HznUya1ZQnVZfpV592nkxi9n
nrdNBpoq5s7jGyrIqc7kaOYlK+ihqMyiJtMhW6ynhCXJGmgme2zSm9Y1Wqaw
HzdobeC1VZMFT9AJgs785Y+yd9fwRQdd0dEvrfhDvdhPqx8Z4c3uK9VngqV4
i5ptOJ7FDiEX4ridmNLxc/qBTNW+fDB7480xLdMw/Kyr6pLssACcWXusUBLk
jLefYues/SAC+Xy9O8oj+jjbB2vLGTpJkF7yXlrrk8PVnen/UYEsHDzreKX9
9BG8uFDta8g9eSC4rstLC1Pi90YF7dSqT0BmzjGKTa+v2KbXv4DRl0pwaTlq
Vap6ELv3N9jzYoUKvNJ8qtl3tx779aBY5vyrk9DQNGs1yVuEFT5pzdF+rQ7r
CwN6opupWMs/DsazQwRg+iVueZP6HDbUWVh16TsOw0qZt9O87DHvIpoXr9cB
tlMsjm2Xa2MK+mv06kw6cA0MV8CNAut56Te5zq0Lrywf/Ftt/UKQDRq7+1dI
HxipR9+SIt4jqTQLC3VxA+jeVjWRZvgvXw+4w6iEDKGHdeXg7ee9KFDQadZI
2hj0i03HpB4MIdLdrFYzFROQp7VaODLQhCTeLWgSdU1BU7Lm1NO9ahR1ttSo
RdcMKlKY9j4PVyP+6LKJajEizGo1x3U2XkSsfpFaLx8QgUOoYZG+wBcd//Xn
eNEjC+CP5iX/yUtG1Td2vujLkCA9N+fqnmQo4u6gsLlfSYJvr8WWWo5ZopGd
xzG4jCXciDp7IeWfNgo7OPnB+gQZBjk/keIVnFH/7wbXjBdkSLqfFhKaaIiE