-
Notifications
You must be signed in to change notification settings - Fork 0
/
SIRS with vd.nb
12161 lines (12076 loc) · 606 KB
/
SIRS with vd.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[ 620465, 12152]
NotebookOptionsPosition[ 616940, 12047]
NotebookOutlinePosition[ 617398, 12064]
CellTagsIndexPosition[ 617355, 12061]
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[{"\[Nu]", "*",
RowBox[{"r", "[", "t", "]"}]}], "-",
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[{"\[Mu]", "*",
RowBox[{"i", "[", "t", "]"}]}], "-",
RowBox[{"\[Gamma]", "*",
RowBox[{"i", "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"eq3", "=",
RowBox[{
RowBox[{"\[Gamma]", "*",
RowBox[{"i", "[", "t", "]"}]}], "-",
RowBox[{"\[Mu]", "*",
RowBox[{"r", "[", "t", "]"}]}], "-",
RowBox[{"\[Nu]", "*",
RowBox[{"r", "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Beta]", "=", "0.04"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Gamma]", "=", "0.1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Mu]", "=", "0.043"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Nu]", "=", "0.01"}], ";"}], "\[IndentingNewLine]",
RowBox[{"R0", "=",
RowBox[{"N", "[",
FractionBox["\[Beta]",
RowBox[{"(",
RowBox[{"\[Gamma]", "+", "\[Mu]"}], ")"}]],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"tf", "=", "100"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"szero", "=", "0.6"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"izero", "=", "0.2"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"rzero", "=", "0.2"}], ";"}],
"\[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.5638355904558997`*^9, 3.563835704468421*^9}, {
3.563835950077469*^9, 3.56383606512805*^9}, {3.56383613822723*^9,
3.563836367357336*^9}, {3.563836441114555*^9, 3.563836444201731*^9}, {
3.563836535058928*^9, 3.563836589524043*^9}, {3.5638366644783306`*^9,
3.563836665858409*^9}, 3.563836695909128*^9, {3.563836848841875*^9,
3.5638368525680885`*^9}, {3.5638372255084195`*^9, 3.563837265661716*^9}, {
3.5638383151507435`*^9, 3.5638383180759106`*^9}, 3.563838362903475*^9, {
3.5638384377487555`*^9, 3.563838575206618*^9}, {3.5638386112096767`*^9,
3.563838657701336*^9}, 3.5638386889021206`*^9, {3.5638388049467583`*^9,
3.5638388069938755`*^9}, {3.5638388495773106`*^9,
3.5638389667440124`*^9}, {3.563839052506918*^9, 3.563839071571008*^9}, {
3.563839147339342*^9, 3.563839232596218*^9}, {3.563839314060878*^9,
3.56383942477421*^9}, {3.5638394696337757`*^9, 3.563839480051372*^9}, {
3.5638395362795877`*^9, 3.5638395453271055`*^9}, {3.5638395986821575`*^9,
3.5638396665290375`*^9}, {3.563839704558213*^9, 3.5638397404972687`*^9}, {
3.5638398239280405`*^9, 3.5638399312031765`*^9}, {3.563840014577945*^9,
3.5638400676949835`*^9}, {3.5638401105414343`*^9,
3.5638401221330967`*^9}, {3.5638401613373394`*^9, 3.563840168990777*^9}, {
3.5638402066729326`*^9, 3.5638402331264453`*^9}, {3.563840302771429*^9,
3.563840436440074*^9}, {3.563840514500539*^9, 3.563840531504512*^9}, {
3.5638405784992*^9, 3.5638405821144066`*^9}, {3.563840618338478*^9,
3.5638406387506456`*^9}, {3.5638407249235744`*^9,
3.5638408190029554`*^9}, {3.563840856229085*^9, 3.5638408646885686`*^9}, {
3.5638410103198986`*^9, 3.5638410835630875`*^9}, {3.563904889849905*^9,
3.5639048947221837`*^9}, 3.563905146271571*^9, {3.563905177983385*^9,
3.56390518226363*^9}, {3.56390525411874*^9, 3.563905299470334*^9}, {
3.563905472683241*^9, 3.563905489465201*^9}, {3.563905647569244*^9,
3.5639056597809424`*^9}, 3.563905703646451*^9, {3.563905742419669*^9,
3.5639057793877835`*^9}, {3.563905869972965*^9, 3.563905882011653*^9}, {
3.563905928065287*^9, 3.5639059436801805`*^9}, {3.563906210805459*^9,
3.563906215908751*^9}, {3.563906296232345*^9, 3.5639063025807085`*^9}, {
3.563906365490307*^9, 3.5639063658563275`*^9}, {3.563906424897704*^9,
3.563906433242182*^9}, {3.563906471104347*^9, 3.5639064759836264`*^9}, {
3.5639065288356495`*^9, 3.563906529453685*^9}, {3.5639068329810457`*^9,
3.5639068636147976`*^9}, {3.563906895527623*^9, 3.563906897037709*^9}, {
3.563907110560922*^9, 3.5639071121890154`*^9}, {3.563907274969326*^9,
3.56390743274035*^9}, {3.5639074760618277`*^9, 3.5639074793070135`*^9}, {
3.5639075838239913`*^9, 3.5639075917484446`*^9}, {3.5639076912791376`*^9,
3.563907827329919*^9}, {3.56390786003179*^9, 3.56390792333441*^9}, {
3.563908035557829*^9, 3.5639080787833014`*^9}, 3.5639084082481456`*^9, {
3.563908463246291*^9, 3.563908492127943*^9}, 3.563908626086605*^9, {
3.563908825166992*^9, 3.5639088270190983`*^9}, {3.563909238651642*^9,
3.563909240037721*^9}, {3.563909322232423*^9, 3.5639093344691224`*^9}, {
3.563909382587875*^9, 3.563909606090658*^9}, {3.5639097046552963`*^9,
3.563909707405453*^9}, {3.5639097449105988`*^9, 3.563909747819765*^9}, {
3.5639097912782507`*^9, 3.5639098366698465`*^9}, {3.5639098902499113`*^9,
3.563909890652934*^9}, {3.563909925842947*^9, 3.563909927602048*^9}, {
3.5639099589818425`*^9, 3.563909959344863*^9}, {3.5639099995731645`*^9,
3.563910046916872*^9}, {3.563910148071658*^9, 3.5639101691978664`*^9},
3.5639102340005727`*^9, {3.5639102721207533`*^9, 3.563910351650302*^9}, {
3.563910421532299*^9, 3.563910510495387*^9}, {3.5639105664275866`*^9,
3.5639105760811386`*^9}, {3.5639106580418262`*^9,
3.5639106824562225`*^9}, {3.5639107247046394`*^9, 3.563910728398851*^9}, {
3.563910840744276*^9, 3.563910862496521*^9}, {3.563910926792198*^9,
3.5639109604401226`*^9}, {3.563911036075449*^9, 3.563911037511531*^9}, {
3.563911431340057*^9, 3.563911437320399*^9}, {3.5639114789257784`*^9,
3.563911493954638*^9}, {3.563911525254428*^9, 3.5639115374281244`*^9}, {
3.5639115750762777`*^9, 3.5639115880030174`*^9}, 3.5639116582250338`*^9, {
3.5639116984503345`*^9, 3.56391176008786*^9}, {3.5639119457344785`*^9,
3.563911983828657*^9}, {3.5639125210403843`*^9, 3.5639125490409856`*^9}, {
3.5639125979987855`*^9, 3.5639126611083956`*^9}, {3.5639128670021715`*^9,
3.563912872138466*^9}, {3.5639130702837987`*^9, 3.5639131278820934`*^9}, {
3.5639132526262283`*^9, 3.563913434192613*^9}, {3.5639134979902625`*^9,
3.5639135222086477`*^9}, {3.563913559982808*^9, 3.563913587937407*^9}, {
3.563913643129564*^9, 3.5639136454726977`*^9}, {3.563913893811902*^9,
3.5639139485590334`*^9}, {3.5639139840870657`*^9, 3.563914064888687*^9}, {
3.5639141429901543`*^9, 3.5639142000144157`*^9}, {3.563914269810408*^9,
3.5639142840542226`*^9}, {3.5639143487949257`*^9, 3.563914404738125*^9},
3.563914455294017*^9, {3.56391451899066*^9, 3.5639146165662413`*^9}, {
3.5639146469249773`*^9, 3.56391465099021*^9}, {3.563914684638135*^9,
3.5639146860572157`*^9}, {3.56391471638395*^9, 3.5639148044359865`*^9}, {
3.5639148449423037`*^9, 3.5639148622662945`*^9}, {3.5639150099297404`*^9,
3.563915010347764*^9}, {3.5639151494727216`*^9, 3.5639152487053976`*^9}, {
3.563915326603853*^9, 3.563915346531993*^9}, {3.5639157813918657`*^9,
3.563915788211255*^9}, {3.5639160749616566`*^9, 3.563916084421198*^9}, {
3.563916863145738*^9, 3.563916867956013*^9}, {3.563916974604113*^9,
3.563916998251466*^9}, {3.5639170471642637`*^9, 3.5639170731837516`*^9}, {
3.563917169124239*^9, 3.5639171692732477`*^9}, {3.563917225221448*^9,
3.563917281815685*^9}, {3.5639173866606817`*^9, 3.56391743454642*^9},
3.5639174785899396`*^9, {3.563917530621916*^9, 3.563917555013311*^9}, {
3.563917726810137*^9, 3.5639177280732093`*^9}, {3.563917953474101*^9,
3.5639180072371764`*^9}, {3.563918038582969*^9, 3.5639180432792377`*^9}, {
3.563918079533312*^9, 3.5639182234085407`*^9}, {3.563918264612898*^9,
3.563918265321938*^9}, {3.563918306125272*^9, 3.5639183204030886`*^9}, {
3.5639183512768545`*^9, 3.563918397598504*^9}, {3.5639186015911713`*^9,
3.5639186261025734`*^9}, {3.5639188701865344`*^9, 3.563918974285488*^9}, {
3.5639192569106536`*^9, 3.5639192590437756`*^9}, {3.56392175996782*^9,
3.563921763652031*^9}, {3.5649408605169983`*^9, 3.564940939182498*^9}, {
3.565699720038359*^9, 3.565699860512394*^9}, {3.5656998929582496`*^9,
3.565700399603228*^9}, {3.5657016419202843`*^9, 3.5657016452584753`*^9}, {
3.567773137024293*^9, 3.5677731407058997`*^9}}],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"General", "::", "obspkg"}], "MessageName"],
RowBox[{
":", " "}], "\<\"\[NoBreak]\\!\\(\\\"PlotLegends`\\\"\\)\[NoBreak] is now \
obsolete. The legacy version being loaded may conflict with current \
Mathematica functionality. See the Compatibility Guide for updating \
information.\"\>"}]], "Message", "MSG",
CellChangeTimes->{3.57109286031627*^9}],
Cell[BoxData["0.2797202797202797`"], "Output",
CellChangeTimes->{
3.565699864028595*^9, {3.565699900984709*^9, 3.5656999321214895`*^9}, {
3.565699969247613*^9, 3.565700399952248*^9}, 3.5710928606222873`*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwVkXk01WkYx5FhWlRKWZItUyg3InSl5yH7eunipxLZ7pUoN0tKRRyTLpUW
LfaxNUxIoUiWKWuNaDGRSvbtvo0yypL5zR/vec/nPM/3833PeZW9jjj5iggJ
CbnS5/9b5eJ8drpC3a5+gxkTHZ0+sK4Wtx/kLECaJ2evsGw/PJfQWyzqvhJL
jOvbMpYOgMlFs9p8Bxk08K4zYA0MgO/WGQltCSWUTjzVWFUyCHbqvv+stVBF
huSU+eLoISh9dSq02UAN+dJS11B/GAxXaRSLntmCNafrmQZkGDizl9NyKS2E
jRc2myePwJGGIZ1iu234YyCZb7NxFM7vFftl/rAujo64HZAqHwXvxpN+2lw9
tLGUb2Vqj0GduKxjgJsBPmcVnFDNHANtyajkkRAmGrZIcTWWj4PnF6Mt9od3
Ynxlx6OmgHHISn8qP524C72dI2wVm8bh3gVmd4cTooy46AJTbgJCKo1YabbG
GHQtPrbZYwJqWnIn4/eZoHs+o0knZwL8LXZU+1O7UZcn1SnyfgJSb2uXBRwy
RafJcXb2WgEkfamOPhpuhgpFS6/etRTAYEFoGJdnjhmXHSz7eQJ4Yn75xtmz
Frh/+6OdWqkCEHbo+nQxyhIDB5REJx8LQMi4XDgzwgrDxupPNfcI4Jz0tF3H
r9Z4z8Os4N9ZARzd610/HWODPienF9hrCPRG6/rANVvUYCTbSzEISF252RWf
ZYeBre1ve3YTyFupEdGYYo97kjxkNrsQqOMpmm0qdMBvYR+TxLkE0pUy5Y3Y
LPQLahkSDSdQHTkjxH3BworYNzLWsQSUbcTnNrg6Ijtl0SgkEbjfps5/+8oR
58My/BRTCSQ0/NBScnfC5vC53vI8AleOjRzyfOeEz2quzl4vJhC5RmGFitse
fNNTI8arIPC7+k8tIR/3YFLSunnZGgJiX4Vl6w6w8VmzOol9SiBf6+7p9CE2
1gZUSv/ZSkDdq8B+PtAZLVxMr+9qJ9Bt9bDUecwZ1xv3ynS+JtCRfUITjrlg
t8pAG6uLwLd7jjtVJl1QsxMLVr4nsPCAVxYX5Ioxb29KH+wlYKNlbhY55Ypm
VtsdU/sJNLpv9Q84RmHt6uIgxyECRSVHpxcWKNTXES4UGyXQktQdPMwhwDwn
2ZeW1wdMcVZJ2UMK73gqrysfJtDMjHk4zvkAIsv9NogE9kFUaU3CpWQKi75z
9onR86/7Pr7sekfhavXyl54jBJoSOiynOG0Q774gp8TqA9Vz0zPzsRTeUl2c
+Afdf11bxsSylfYXTK5ZRu87xt32DRRQOKy/ye0+zQ1awRwxbg2UnhQaLjDu
g395q0t1Iynk7NikdobOc+V1v5+up7A3JBx76X6V4CfyOQMUFltqVp+n81Ep
CUV60xSil+HtHprrjhflqHGLoV42wuMIsw8OcLMkXoTT7/ERd+HSvnbJerbc
YwptjdorWmkf38/pOPRSmPtaoyOUzg9fklvSMUnhtwumS5ppls774ZQyS+HW
DAWdSZorflZXDORmQG2WlHe7Xh+s8uZWpYVSGHdi0cB+2h+3LF5kuIrC0NLZ
sVrajwdlN0p8oNBnysM5gM6n8nsrMj7T//Np3f3HND/LfCrk/53CQtPln4do
9irUz30+T2FEW1DwHM05f5MuPkMFYzqn7lTJDsH+lNIHfMY2tDLO1qTWj4Kf
9URMu5khKl+95Sy6fRwezF3J5zNMsLjgtVLnjQkwYmsnMirNMWVqYrDplADq
C//itZvZYFPcJROpFQQsFh2mQtodMHI+qtbmEAG/O6IZfIYTZvoPrZLJJqDm
9lLSPYSN02XRu6MbCIyK/hbLqHRBJc9Urv87Av8BHsenZQ==
"], {{{}, {},
{RGBColor[0, 0, 1],
LineBox[{1, 71, 65, 60, 56, 53, 51, 2, 3, 4, 5, 78, 6, 7, 8, 9, 79, 10,
11, 12, 13, 14, 80, 15, 16, 17, 18, 81, 19, 20, 21, 22, 23, 82, 24, 25,
26, 27, 83, 28, 29, 30, 31, 32, 84, 33, 34, 35, 36, 85, 37, 38, 39,
40, 41, 86, 42, 43, 44, 45, 87, 46, 47, 48, 49, 72, 66, 61, 57, 54, 52,
73, 67, 62, 58, 55, 74, 68, 63, 59, 75, 69, 64, 76, 70, 77, 50}]}}, {
{RGBColor[0, 0,
NCache[
Rational[2, 3], 0.6666666666666666]],
PointBox[{78, 79, 80, 81, 82, 83, 84, 85, 86, 87}]}, {}, {}}}],
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox[
RowBox[{"\"S\"", "(", "t", ")"}], TraditionalForm]},
AxesOrigin->{0, 0.6},
Method->{},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.565699864028595*^9, {3.565699900984709*^9, 3.5656999321214895`*^9}, {
3.565699969247613*^9, 3.565700399952248*^9}, 3.5710928617113495`*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwVl2c81Q8bxg8hSsooI5VUQrI6Sf787tve+zhOKHtLtogQLQqhkFBWRFZS
smWHzJBRNAgZFaHh6Xl1v7jeXt/r873325wzsqcnkUiRdCTS/69gzJ/MtL31
hJSktlRychtoVW/W++y4AVyK27Yc1WqEzm0yLAynd6CsawT9TGotKMWo1j3U
58HFV8t0NbwVYC+xvk1qmwAOfp/d/5qtDHRF7Jd2qR9EwS4WN/bLRVDaH+zb
JiuMRbdV9jKdyIP/OESLGELEUM/L9/63gUxw/BWXmk2TRPpEnbpmh3Q41zx1
rEhXGjOXHsaf0kyGSDOmQ3/cyHi21liDlSkBbFsuOEg5yeDtsr/HpvfFQv1m
XkPXU7I44cJ8hZk/EqTYQ+988ZHD28++DKu2RYDVdwUxPTd5DGuq/FvYEwoP
0pr4f94ksF7Wzzn3ZSA8iZYb6TVCzL15aFbZ0Rd8XigYpOooIruzXcoFW0+o
bc/+dt1cCcsK+YoYW13BWf1ktTNNGRtKtqhdinSAe7lST11dVPBD+2TX01/W
cOt7dZiHvypqYVXH7z4L+PzI18/JSw2nBD2OnPEzhUa1uKRLl9RxXGLD0avZ
COj0307GhGpgGlP05a9hekBSLKe7H6CJw+XiGyyfNOAa90/d3qtaePF1/CvR
NWXwMLNt+BmujQqePI3qKggTYWQ7uK2DGQljW0jb5IArPvnt9Qe6+N7utDX/
dzLk7BANaEnRQz3Vix0BvFJQ77VP9XC+Ph5WSHBv8DwCaQL3+RUoBkhReqEu
YCYE1UHrJKduA9w16+4zri4I+7U3/z5gaoik94ux2xz2QNlrkajhfkP8/d57
2PUPN9xo/ispcNoIQ4RM7pOBE+K9v7hYjRphpCpVZpm6HYJ27t0ueMoYozvP
JBemsUCeCGO7z3tjzNMg9TsmMgDTDzre+jMU3Js1QDN4RYKHkiUX06YouDLa
KP1yfp0QsXmk9+esCfaUxiXwFi0TI5oVpSazJugsFF3J/myJ6M0MPAreVFx0
bnCQ850jVp8Yygt+o6JZEFdRktQUsfHc6+kVd1Mc1WZn6VufJLQl1VSDlk0x
41+R+SLHiZbTEs6u3jTMfWiUoTY/RBQWe/zc2KAhnZNUdff7PqL91ojntOMC
5OL9E72tzSDw59m9PzLMqLv7Sp7VVD1c66r+aELlwrGctO5hkSpg+sk/z3Cc
Hyvy7O2r3MvhVEHRxgFRQXw1QR9s1FwCpBDPY4sghM+jbAwuSBSA33djjilz
UXzUaBA1Yp4DEQzCy4eTxDFra+XDjqwHMGtV7SRqLIWS1hO7A8XuQdAOfzHq
6WO406NGzirkDqzsqFdNdDmOgXyXWbV/3oI4x3bzLqsTWP433iNk000IuNHh
Lep4EreeeWJMr30VNOgtJ1I8/kPdUOd09tJLQGoLyxIMVsCsx+X7LEgXgeFM
qnNXMuDhxaOaV+f9IeFg1c7nrxEfvVBt0tHyhrsNng84+xVRxjYl5wbXOeC6
vjC4NqyEus/obzcxOgP9+eMVLhPKmKYVpdux3w52dfnvnl9UwZ3+aa5UXiuQ
f80bNPdDFf12ll4KpTODgdXwq77ralhxrNq84jcFJrXKFa7TaWBnFWvJgTkD
IGlofOP+o4HO+W23q6K1QaXgzeEKkhYWcqbGntRRg9+yR5k/MWpjsSejBjIp
gfC8dnwKqw6O1u24EcivAIFbHNPtduuikvbfLzSLEzDWM5GXu1cP3X8Gb5f4
Kg1l9Qct/Q/oY+6nhO05GeKwMZdoOc1lgJs9RN87NgpDoV5ntFGcAUav0yU8
ljkICx9C/zZxGqI9b3eZk50AuNxMevn3tiHSUu8eT4jmAwFtrvHfPEaYX2O+
8DZoJ+zY80O37p4Rrux7Yl1xmh2uNW5SqOczxoZcYZEJE1aIS7/l9uqeMXIs
s5K//2SCoejwpbv7KDhbwn7uThI9UOOGX0lkUnAoRbU7kvMv4dT2+HOnqAnK
ANctf9FVgv5NSWhegQl2U1Q4fZO/E29WeD20JKi4raOgfKxugfh05HNZcSkV
T75TFf1U/IUIkpk8T5E0xeUMt2mBXx+JbB+LYLpSUzzYa3NHlGmCYGmRG/Qm
05C8FCMwnDhCyG02KH5aQcMiNXb1CaU3RJtceMWc4ztgLU+v0RRqhcmvY9bn
lzehxmSb9iDbSxCdX1spOMWBQbVflhPfVMNBZSlTqhIfcuyVDVd78Qzefe2/
8mj3ftQZs2h2migFIYprXLjyIVR9bCj/tfMxcH2cVqZpiWB+2uL6rpqHYLVV
xOhO1FH0DbidMFWWASLiwpumByTx75HUY7tqUsHPVaL40oA0bt8nHZ01nAh1
T22Pt8+SMdVyncXsZxwEch+YCPkog31eQ2SPO9Gwp91rUuu9LKr9J0RlyLsG
VdcXKhu+yWH0TL2/R3k4SI5bOj3/IY/fwve80mwMAaOLUi8nOQEdywskWvQC
IOvUhfOqtxAjd4yFPzT0gas1g+dpsYpo2NT67oCHB8zkaX2qv6uEkUPvr107
7QLllif+m7mnjIu97BPGuvZwqpIyXfBYBQt6Tn1u3GQNDhOzF8XLVfGjSeq8
ibU5jPf8/uteroYzWiFs9yWp0KdxcJaxQR2VStWPcPYZQqj2rZm+Gg00tSS+
Zf/VgQy/5gSxak10lwsTdIxRByVhjvyRZi20sBnjV1pTAp2swUeardq4Hnd2
m0ozAeLNnXRGQzqoEX/nyn/5snC+VEGJ970uvjYde3DtyzE46nzYLf2dHtZ4
/c3glpAEq5zD9EEz+hjVID1FcIhCzrPyAXVvA2SWcrsYrHgIqMXVG7u+GaBa
xbvt0fr7QW3uteUtP0P023K9YPvSbnDZlk6eXTHET84bNmplu6DLMwAWg4xQ
5J2GLWpwANu9lb7YX0ZI7+wcXdOxDa6+HaxeCTTGsn0Gu1mVmAEtGsc2kyh4
TiJY1O/WJhB/edZrNZSC/waCmtG6QTAHF0nLsJhgUZDTXqbwNWIXf6Gj8U0T
5DPRkohO/UHQ7juzG2yjIsfthsuaTYvEqKPk0s84Kjr/V8q7xjJLyC0nm/3H
boqrtyRuKG//TPTeceRViTXFt+0XyOVrE8SjC8i2ykVDC9+hZdvNY0Roae2N
2Ds0vLVZv2Q8epDgeRmo/PXmT5A2OjUwFtIEp1njBi+dZ8WiEpfOoN46mH7E
fPnJ4b145VPJh6XjT2E4z7/fQv4w+ksr2OnlPoIh+tr3vmePILk1+vEcWzb8
MH/f93aUhrafj7CfZBkgWm/0aiw7vobnTu12AmOtoPLZo/FpPx1eee0Wv5j3
En7/4Q6lmbNj5GJVRV95DRxzrty6ocGLdyRuSnk6PoeMaMf2swkCqCyj91+t
7xN4QUwr6jYdRImwy+6GOYVga/FkNLZeGHlYDKfT/XOhKzovVJfxKCb3Rtkp
KGdCa8LpbrVESbwUeFjEPyYNviVt/ZZ0SxpTMw72v4hOgkPMxx/bFpCxTJyn
qScoHtzPvc6k5chgKv3GScIpBhamm3hH7sli2u3G+W7L6+DXeULkzRM55JQ1
3UciIqB1Wk8rvFge3ydINYnSQiEw+tfCXDeBVy/dehDDFAhiP96KsHohWsy+
Lrz/1gcuVkV+DHRXRL8Wx2CuaQ+gyK3OHrqghJsqdu232uEKscayv/MCldGU
nC6f88Ue1qVe6rHeUMFNc27zpIvWEDRDCX19RxVjTro6AI8F2BAVL1wT1DDS
jZ2J3EOFDvURRZUMddQVSFUQ1TYC7Hz6uSRdA3dxM37prNaFHnNC2yFFEzd3
BXsPn9CAT+kRv6kPtZB+NmAy2l8ZzoSklgvkaCNSPAqUswCYD9I9aS3XQZc5
w8XiUyehOcM1uqdOF3kw2Pm+GxlmqnK9nlTr4badTfN5GZIQeyMzs6ddH3+t
riRV1YqCoNuRgAZbA7y2M3mKpeMQRO6yCIUJA2QrkR9gbN8PZSmzj986GGKL
Hu3u0Tx+iHH/aDn5yRDZ4lprb1G54Si/vCzzWSP0lrq+KtjAAd2bPdJXZ41w
9ThPn1UMG9CvBjqwuhijk3yG+PpXZvjR8NGYWDTG7Qe4Kuo5GODuyLOm0XMU
pNwvf3GCSoIhajHNZZWCJOtvHedk1gkBEDPPDTbBO/p7Q5MPLhN8vuyTRr9N
MEWKY1OE/BIxeNkkZnc4Fcdbiq1Vq2eJ5H7fV4V0pnhx1dkgNe0zMX34rKN9
iCmuHU9xv3x+klg3eSkix0jDR79JoprNY8TBaz/X/0TQsEm5ZxOX/BDxV5Vc
s3j/OzgVc+8f3d8MNQyH5OpMt+C911OWxWb18LM77qSF+B48KVhtyd79FNwS
foQL1Qih/V7pjxzO+XD4wvWAwo+iKBLsZNBdmg2JUjxKGq9oaG3z4m1f5gBh
eCXX/uw8Dd8J8uC8Vz/RLOnpyORUC9J6SVXZcm3gY/PCl9JIQr1fgsmOay/h
I+NDFfHRHRhkKS0dKlILckEkB+MOHkwcH2xYrH0ODO8k3V4bCOBM26XLj+uf
AENZDq/0lYMoQtdbOclVBJFJPCJHg4WxzCZZNuFNLkR0f+EWaRDDgYGwQZYb
mcCjsSqz4i+JLZXq186sp0Fx0RpH1jlpnBA6ZkmdTQJayR7n4htk9JO5/ulH
XTw4pRxiVL0qg3nArq7TFwO08+eOH7ggi7uMXmfJV18HysLIO8sEOeT+YZcm
Gh4BXO5qlU7R8rghJXLxQkIoRN+xYvIqIrAu+V7gObtA8KiPvHTOBvGxpLzF
BQFfCPjPju/GaUW0U3uXYXrIEyLp9RdMXZWw5OIr6URXVwjqfD563VkZRbjp
WMNOOoD22XU6+QsqyHmzo9alyhqYzrHQCVxVRcblzIwwXwvI8knkcQlXwyGh
Pl3yAVPQaJE5axunjvOVJmMO14xAgK31o0CsBpJjd9TL79MDJ71AofBITexK
DPkweFMDykcUshMStfD3QHPK0jNlmHILlWe9rY3e+s75XVsQfvTe3cqQo4N3
3baTK0tPwte7l8qOlujiYNvOyx5FZGh1+OJbXaCHI9bBrbljkpAbYu3z+4U+
CuiyFYnxHAGr5k/tG+YG2H52MCKDQwiyNPPr+oYMUGjpk20EsyC4eK9RNs4Y
Yk9VscKmzXtAcN/V5dNjhuiyE5Zc8rjBwjOVzGVvhMmDk49VSJxQWi12Weyj
EcbwjldGfmOD2isdconWxtgW/+5UhSoL6K38iD41bYwLd79x5lIYoHvKYkHP
kYLT7Xk+FrEkqHtZ7MqwSMGpvk+h67fWiSyD2f5eXxOs989wcHZeJgIt7bX6
v5sgwyn79TMBS0Q/l9OZoxeoKJwgv5LLM0doNH4In16jos5IS0j21GfCZ/Ec
24qfKe4XvtUgVzFJCL6MOGH+xxTZmvRE4w6PEytenKXkIBrWHv96MSphiPDa
fOuC16MlCN5iWvDUpxkWQt1nz6izoCa7oqff/Xr4+2EwvriEHzUD6vX79pZD
OcsT860hQuhLymiJeZUPes/mm7Tui+JAzy6lQqYccOInr11soOHluSn29m8D
hKBnI3/WJxrGXzvRH1XdT4Sm3CiU+UlDk43wDXORfqL+fGGWsFMR7Fxxdlpy
boMbApZdB+RJOCY5fvj+/kYYuj4vlHxvBz51o9IH29ZCDeOB+D83eXBJ2OvU
47nnQL66U4EkIYDtWgwJ3NNPIHf+8uUsh4MYpc2/0aRTBGky2vq/zIRxtTaa
K5UpD6yWrn+4c18MPT5GOb8tzoS6B6ozuk6SOJkhJmMjmQ47Kj0y885Io0xG
bL4oXzIsMXK0vAkio/kkX6/1h3hwKOk+Xukvg+WzshlFDLHAiirK6i6yyPtS
96Tf3HVQIolXZkXIYbrE4xnenAh4IdngbR4sjwxbz5auloZCEl/8rbI0Ar2L
tELexwcCR+6k+AMzxMHlwpJ6ZV9YWTd4ymWiiHz+jDq6Kp6Q5bC+6GGrhPzm
nSoPk13Ba8tax5qVMnrNVKeMWjvAjy13Rvq8VPDhmtDC9xFr0CmxVdQIUcWT
8t1vE7MsoCnmivqXQDUcfwFtmTqmwHT1I713pDqOhfy5a/TQCGLFe2xp1zSw
40dstpuBHpSp8ySwX9LE/PPscloVGpAecDfgSYwWeqlukvMeUoaCD54bYTe0
MbkqLtzzMILLo6lnD1J10C4jZ0v5u5NA4o5uMc/TRY6m4deiPWR4VnndXylL
D2/SDQ35/ZIEckDdKPsTfRScS0nvVTkCId8bvPfSDPBDz85EFRkhkPzQ+z61
3wDtcjKbdgoLgpplaC67uSGmhp2PtpbZA67+tLrdw4aoRmURV+vjhh+7u6T5
rI3w282lZxn7OCGDZ+t5tfdG6Hrc7l7iwe2wV285y+K0MTKGSi6werHA7+3B
Yo4fjdFvTiplzpcBZs9KMDDbUvCr1A7TqsckIK9L+vbMUnBsB12fTN064TDc
zj7taYJ2fSrKh24sE48oMX8TFkxw3ZpNek/SEmFssaQq40/F5C/k56/U5giB
Kj4rxhUqsqiPxoaxThEvKEvWNV6maJE/ZDgzPEmUB9jXuK6Z4jhPCEchbZw4
4/RgW7f/Px89oETF+iFCzAnMjhUtQtNCIaSkNQMzcKqI1zJjeGn28+aX9XDp
7HqdmTc/Vrhezs9QKwdzlbiWYEsh9GsqtDD9kw82JvvneC6K4uRgXhafVA70
sDdQ+GpoeDLCSb2E7w0R5WB0HiZoyB2sJIsf+onpWL4tvd/+5Sol65v0+gnu
nL9GKb/+8btn6uUyQz/xjFlk31mndCjUGO2qCGsDqQ2b5V2LG2DX9Jt+9Hgj
9GpzJHUF7UCz7g9ns0NrIVSrLK/XnQdfXdm+/8imCtAZuZrNs1cAG03/mfuv
J7Ae5Vs5Rj2IYfsMMsf/cZ972c5mUlMYM/cdudq0Nw+66MzaOmPEcKrC/erp
l5mQoDf3MsBKElmZLUyuaadDHRa4ilGlUWnr9vpyqWQgJFlXvviQsfi0J6Xo
RzzYyQvprZyTwdigok1qHLHwOVNca9haFkWvdCiQGSLhBfOgcEvQv/9Q/7Za
7rMI2MZB6u73kcc3bZJFOg2hcM9CbNDkDoHRwsU3K/IDwSaFjfEFFTHjooT2
BRNfeGrWQj9soIgNYZ2mY8aeUMrT7RlsqYRkjR0PO4pdYV6Yt6nEQhnJC/GH
hPwcoHXYtjXdXQWHA740dM5ag3tvrLXFBVXc3Zbng1UWQN3zdzzXTw0PZ7nc
9rc2hU8OlZ2BV9TxN6G5PlJuBK01Jgd8IjRQ6hPvvKCzHpipTLUPBWtiUn6C
g8BrDUhobf3QEKWFZU2M/9xbGdjSZ97LX9NGP0UjqTfHEYR2sUV7JOtgmwvu
MVw5CVsGOHu8snXRxsbsEf8kGcQ/n1pluK+Hm7nLiz23SMH68kIsf5E+Pu1k
UVE0PwJskVKtR6kGWBegLXNEQwhWij4vXuo1wHr5NLHvJwSBPKDcy3fKEPuk
D/Qc1dsDsQfqs9ffGOLHKyOP9ae5AT0/S+2xNEKvXscj9hKcwKcW6Gk8boT9
re80qhS2Q2cJixmYG+MusaJV42ss0Muv+OvspDG6MX4wz7zCAKvFSdEfrSiY
/6uiK6SKBO47l9aefaGgbcihn16D60Rrqlzxt3Mm+EhLxDQpdZlQ0s+i+X41
wWzT1JzYh0sEB6m8Qs6Xip+SYyLXTs8RrLIyXZw/qMh7ZmsV474pYpSpNCrD
wxRn5uvDTGcmCdLAAJ3PT1O8QHLmM/McJzhsnSpTfWmoYsmsWjE4RFxLutMW
L7UIY8yvIgqeNANXmv3izgvMmG414xk5UA8qpVGDEhR+LPBiPZZjVg6FE2F/
yvWFcOPumdrTuwrALUtuY8NNFIPqE4M81HLgCut1+ulKGubr6B6qE39DoDWv
0LZ3NHzy2mMbeaWfuBc18Sx9kYazc16vPKz6iY77TSTnNRp673QvFeDuJ2zy
T2R3/qFh3mKy4Mi3PiJraOFtlLggjpIp3PsaSsAipfR5lLg0XjkkelU2Oxme
/45/GCWuhCcP+AVdfe8GDfldXj2q2niacSCLV1ARHB4zpEeJG+HOOu73Fi1c
MMOQESH+DwAhJz8GJ/svhIPW1/Ae1f8wPjPbKnlLOChQpG6Kv1DDbXN0o7n7
qKC+yY3m06OPryy2plmniILwqT720z4UzKy25UJZEvwPPHpGIQ==
"], {{{}, {},
{RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwl1HV0lmUcBuBRskmNkaM3YIMBAtItJaCEYtFKw4jRtpSEmIQBSgliF0qX
IKAo2IQSBhiEYiuY1+/wx/U993NvO+c77+7zpg3I6T46V0JCwhwfcV7uoyAf
y48ziHrcrGtBXva6P0xuOrgX4hP5CQZTn1t0LcnHPvdHyENH98LslxczhAbk
pZO+CAfkJQylIbfqWnER77k/Sj6ucE/moLyUYTRiavyu80pnUQ7JyxhOY27T
XUZ+3nd/jPx0dk/hU3k52TThdl1rEvnAfSGJdHEvxmfyk4ygKUl01RfnsLyC
kTRjmu5iZzdnCY7IKxlFcwpwlb4kR+WnGE0LCnK1vhTH5FXk0JJCdNeX5nP5
acbQisJco0/lC/kZxjIrnoOzCNfKZfhSfpZxtCaZ6/Rl+Up+jvG0oSjX68tx
XH6eCbQlhRv05Tkhv8BE2jE9npuzh7MCX8svMon2FKenviLfyC/F/mKblKCX
vhLfyi/H1mKHlKS3Po3v5FdiO7E5StFHn85J+dXYAbPdOzlL01euzCl5dfzf
Y2Ok0k9fhdPya9wR+6IMN+qrckZ+nTtjS5TlJn0G38truCt2w93x3Jz9nZn8
IK9lcmyH8gzQV+OsvI4psRsqMFBfnR/l9UyN3VCRQfosfpI3xOZiN1RisL4G
P8sbuUeeHtshjSHuNflF3hTfM7ZDOkP1tfhV3syM2A2VGaa/hN/kLcyM3VCF
4fra/C5vjd3FbuLv47k5s511+EPeFv+T2A4ZjNDX5U/5jfi+sRsyGam/lHPy
dubEbqjGKH09zss7uDd2Q3XihVefv+Q3Ey68++6L7ZBFjnsD/pZ3cn9shxqM
0TfkH3kXD8RuqMlYfSP+lXfzYOyGWozTN+Y/+S0eit3Ec4rn5hzvbBIvYN52
zI3tUJsJuqbxfmaP+7zYDXWYqGtGbt5xnx+7oS6TdM3Jw7vuC2I3sV33NiTx
ofui2KzcNt5NfBQ7dbajAJNpz5R43n72P+01p8E=
"]]}}, {
{RGBColor[
NCache[
Rational[2, 3], 0.6666666666666666], 0, 0],
PointBox[{385, 386, 387, 388, 389, 390, 391, 392, 393, 394}]}, {}, {}}}],
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox[
RowBox[{"\"I\"", "(", "t", ")"}], TraditionalForm]},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.565699864028595*^9, {3.565699900984709*^9, 3.5656999321214895`*^9}, {
3.565699969247613*^9, 3.565700399952248*^9}, 3.5710928617423515`*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwVVnk41N8XHrKV3diTbCmRJUsRzklEKkVo8JUszYwlIe0KkaWJiiJZsxZJ
EYUQETOfaSFpQUp7KVT27ff5/XGf+5zn3nvOve/73nOOqs8Bp328FAqlixz/
n9XOzxXkKDdb4hDRkpHBBvsGQYcvjAU4au5y7okjAU9ETRbzeUpg0aO7TSxe
Llidt3lYskMedZIWlr7vJ2Cf3rSogagK3u+kOt4wJmC71r5RWVsNbLHVSkxT
5EBl98lD7PWrMCBCWcP8ZQdskFpdwRepg0HXXwflb28HxkxKdhFNHzW6G46J
1bXBgcdfDSu2r8Uef73cYotWOOsusGIuyAjD3ifNDRW3gG/7CboB0wQ3C4k/
kVvXDM2CCo6BbuvRImhIwaO+CQwko9K+h5uhwRElvtixB7D3r4WOQ5A5ClWl
mMbG18G1nDaliSRLdB3NliGk7kNVsllvlxOid/AadtWragivs9iZvW0jTovO
qUb8rIImTtGfRA8rvDIseVGVUgn+tqYN/rRN2CWvY3zZ4zZkXTeoDgywRotb
Vs7ih8vh4t+G6JAjNmhoPvB6c3IpfCk9dJgZthmFXcoShWSvQ+vmlCunT9ti
9+3ZvI7iIuDZ8XbwfJQdaoT2Cr5OKwDKxhqevGNbsCepe+ZF0zVIkJvY3hVv
j+aBE3u9P+RCiLtvy0TMVlz2fLjF4kUWfIg28oPL2zCnVWyN1csMkE7NeJt4
bTvyLtvofFQ/HYolVh9rz3TA4r5ekXXjl6A5bLnNyrIdWDNJhdUdKZCjkqdk
4bwTBxrVHt/6cwEaIqYpzOc7sUc2aad7YTKobhWcVd/tiAVbJnxND7Hg7jMt
1ptuR3RzOxCl3ZEA5x7P66t4OmGK7zU5tnYcpB78HrC3zwmTazJWlE7HQISM
sria2y4s0z5wbWDxabihxc8Jf78LgzNtHb9rRIHAPx6F5j3OWBvszxldehJK
9O+cyvnqjGs80r3OLDoOWj6lDnP7XZCzBa35tx6B3i21lS4/XTDmY8HJB7yH
oKvg+Bo46IqWBhtCPNlhMFnlaK72xxU735lTpYZCYOF+WHVc8G60C1yW+osV
DFv1N9tEjO3Gir+pDf7OgdDuqecfeJCGhrKZaxTzmXDrdsjEwgIN3+sFp/kN
7wPOxd7Qb4xhCM1RMK/J4oDK3L2sORMhhPX1IUoEAQlPGz65uErjGcVycQYf
FwQmlH7zGSvh0fYRrl06AW43KxbUV6vhdnfr2BuVHKBEhhqOgCaKSxHC6SfZ
YCa483Z1LQ2dODP7H1kygG0WUzvEGIBRyY6q8hUcGPzV7310bBFuLuBoBZP+
Vv+eGr/pJoX5puIWU0Jc0NhksNvVShG/3+urtm0kYOBXd1zpUlXUNgm1VJzg
gKZzYErMphX4fXJi6vYTNkRVNp27kEZDkdLyVf2/GCD/6PimX0kTIBARKmhN
IcBTJOXV6aMiOLzkwyG33wR8vhXdNe8ki2G8ms4Pxwj4Vip0pmqlMvJm+K5l
+RPwz+P9i7d9NBxzdSgvuUSHjnNddmOMZ6Clul7Pa5AN1l9CWqu7edBaavd/
Y4cJmJ2Ti6J5SGLDuvP3DQS5YOhfL7xgp4CXN4wmCj8jID+Zwdl/SQWjrkop
rJIjoM7y28btbRp4mQgOu/CPDRoJE9NzsTRMjngz5GzLhHkbo8aRvL+wS/ts
UPdzDjTyrTB7uHsJXliiSal8R0BNY65HqbMMlpSOetjNEzDxPMX0P91luHx4
vmVdNAHpBvJWdgQNFcYY6vLf6ZCe7sAYTvoKjneNc2t9ORA9Zr9gJiuAestD
qgOqSD7f9e7noVExr6lvkELiL7OsJL98w1J85pYTXlpGgAE2d3w/PgOPI1d1
r9Um30tdCLfJEMOP1KEWy2kCpHN9t+xylEMjNabrjx8EOMZd37f/Nw0/Cl/3
3r+RDo/1QxkCzCYQfLSv4FYjG8J96g45t1JwORF5ddqPgE/8Jda6fRL488S7
cHN+LphFUOi7uPJIbU7Xqe4hgG9AP+jZThWMHZ/Oz9Qk7bvFCmvjNPC5S2hM
hSAHxsOolUYRNDRTE0xNCmdCmODFE2Glo+C53HbZnnscGI4K/rnHdjHmp9v1
Ml4SkPtYVyXqpTS+JT5fKebhwvzHV6m37yjh1rM8n8PPEcBUMpo61ULDy6tP
qz+mMuCUotDgn7yPcNPIvpFmQ+b/0Eu33/DzY2KdY9dkEQFWXwLWvnothT0P
stRekPgp6Q/N8jQpYtjsjj9rqwko5FGOC4ydAqGcDXftZAnQuWCn+SlBFC25
F1qa/5F4Fv8R1eqSxbTgIwHZIwSohbYqFX6m4TbjJyKfAuhADVhJlfvMi518
lJDWeAKo4wn1p99KYrjcS68DZLzkhN4EVaYwbq7JLKR8JeBpjUR69wsZ3HNi
ir2D5GfnMR6HLE1BXNNuPMNuIiCkZ9sN7VdUXNelvCNOgAtEevWj4hviaN5m
UvGS1JN2enfsyedyWLql07D0EwFRmedumUzQ0GI+W3aTPB2aj94qXMWsgM/5
dxQcb7LhnIrXU3VzCn4puNxx1oOA14m/NTOyJFDQtu7xzCIuNPKrp84lyePP
T4E33N8SYBQvY0HRU0G/HUw/PV0Crv8+c6aQroGBfi7qqyU5sId5TfT5ERpe
tw6hbWExQYcJ7oYVI5C4MrpJv5QDQkC11m0SwhcMkTKf5wRcsj/M+/OyNHZk
qLYqkvFO759+6H5QCY1YsVpDFwnolGxxVmykoe+Cco+ZFgNGe/1yw0s/gJWu
4Jebxhzg6JkuvnGSDxf+3JzpySHAJKBS2ipDCtdmF43OkPjm1YqtvhmliNUJ
vmNELQF7fLh+OomTcMSuf3GAMAHWDr9OVVFEkafW1q1slICR1WYjqSmy+Khd
U27jXwJYdKej8IGGO03FjB9E0kFnFTfwy3+82Mc7OhYdSYDQ8efZv65KouHF
DKuvZP4o3hz1YWRgCS64KtAnBwloMaMdg0syeCohRtNrltTnuuLiwBwB9Iqo
cHhA3sfP/RIrLZ2KgQENdBp5/pOLu9JhVXFU3yh+vZjcr7JUifHtghyO9tGq
e0l9fLuguKTrDw19Kg7WPtWhw0i0l8+jpZKIKwNsHUg9lFp9tZCRkUF7tvLs
EgoXHHzUdwbJU8l4aB9J4iHdVfBEWlIMJZpbadkTBCiHbXwjpyiFu24lneWQ
6xFO/O0PZcn6YewryiL/85ughOqf8hKYoRQYySX9yRXPO2XO0PB33yk6ZxEd
7glpLd/PzAVn/c5CSgEbDBZ8xmRHFmCJu9HnalcCurZKXXkaIYE6AzdqtpD8
RtnfvdEVLI9OL7v9VfsI2NYbXySvrIIHo0wyhw0ImGYdqu931cAuwzGVQRkO
SPky67MP0fDF2ODyijQmJFxJY6cajEDe31TdkwUckM7ZNyJzQgi3Ue7UtT8h
IFnk6kGVw9LoUt001EbGs65kvdJzVsJvHxu6sy8RECeSyPutnoYen8wKZw0Z
UGZOPWVc8R7+xqq03VvDgZ5kVvUFUz5suXP3lHYmAfo6PWeOHJNCDamD4pTF
XPCZz8Vob0VcccooRvQBAfUp6t3nhCdhMHB8foCf/K/loJ35QAQVsu0+LCbz
wY8hT292mCzWl7jZnyHzBXoraIoOkPX946I4vbN0WB+vLP9UnhcXsounrE8Q
wGv+T0z9hCQunetam0LqobKoK+l15hK8t+fKov73BDyYPfEjPFwGjzMpSbQ5
AsLTJp12uQvgnigD4awaUt9K92eJI1R8xW4auEOeH7l7s8OzXwwPlJT3Hpsh
QLEfXi8LlcMCNTUm4zsBWawP93JHaBgZd69s6To6/NgU82JyRgJ1P44ZECT/
hXzCFfZj0pjt1Xmnn+TfTsO2+dqkFO61/3f+FqkXFYPxY7NcUSwfzOo3GSdA
3mif/fZpSUzl7Wflk+uHpbXWj49TUb8onytK6rNxva9c0y9xrH/aSH29QAA3
r43iP0X2G8t8vPnE6PCunVPt+0AC58PfUvLJ/sX8enKBSK0UWu8e+TNE+hO7
f1Pnfp0kngzRmdAn7YArIYre96gYr0CN/0f6H2/aMJq0QRKzWmNb20jbOW/h
ncl6Km59xtlr8//+5PYy5qCpFOql6wd+JW2fsnVFT+ZomGmaq355ah/8l1l5
n6W7Fr1r4nrtXVqBbv8rptNmA4aJjZ7w/VsP92dTS1i6VmjBr+BV71gJFs4G
Sbp1m3HHeYGq/NASaCl7GtZpsxUf+np7ClzKBttFQbTwzh14zpA3Vzc6Bejl
fLksXScUEc18YikRD6vcXkh6hjvjKn/vgu6cCPjBlx+rW+eKm5wd2xIbQ6Hw
9fBblq4axtPawoLucOB/4FEQsQ==
"], {{{}, {},
{RGBColor[0, 1, 0],
LineBox[{1, 156, 119, 90, 70, 58, 168, 131, 102, 82, 51, 163, 126, 97,
77, 65, 172, 135, 106, 86, 2, 157, 120, 91, 71, 176, 139, 110, 59, 169,
132, 103, 83, 180, 143, 114, 52, 164, 127, 98, 78, 178, 141, 112, 66,
173, 136, 107, 188, 151, 87, 182, 145, 116, 191, 154, 3, 158, 121, 193,
92, 185, 148, 197, 72, 177, 140, 195, 111, 189, 152, 199, 60, 170,
133, 194, 104, 187, 150, 198, 84, 181, 144, 196, 115, 190, 153, 53,
165, 128, 99, 186, 149, 79, 179, 142, 113, 67, 174, 137, 108, 88, 183,
146, 117, 4, 159, 122, 93, 73, 61, 171, 134, 105, 85, 54, 166, 129,
100, 80, 68, 5, 160, 123, 94, 74, 62, 55, 210, 6, 161, 124, 95, 75, 63,
56, 7, 8, 9, 201, 10, 11, 12, 13, 14, 202, 15, 16, 17, 18, 203, 19,
20, 21, 22, 23, 204, 24, 25, 26, 27, 205, 28, 29, 30, 31, 32, 206, 33,
34, 35, 36, 207, 37, 38, 39, 40, 41, 208, 42, 43, 44, 45, 209, 46, 47,
48, 49, 162, 125, 96, 76, 64, 57, 167, 130, 101, 81, 69, 175, 138, 109,
89, 184, 147, 118, 192, 155, 200, 50}]}}, {
{RGBColor[0,
NCache[
Rational[2, 3], 0.6666666666666666], 0],
PointBox[{201, 202, 203, 204, 205, 206, 207, 208, 209, 210}]}, {}, {}}}],
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox[
RowBox[{"\"R\"", "(", "t", ")"}], TraditionalForm]},
AxesOrigin->{0, 0},
Method->{},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.565699864028595*^9, {3.565699900984709*^9, 3.5656999321214895`*^9}, {
3.565699969247613*^9, 3.565700399952248*^9}, 3.5710928617633524`*^9}],
Cell[BoxData[
GraphicsBox[{InsetBox[
GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJwVkXk01WkYx5FhWlRKWZItUyg3InSl5yH7eunipxLZ7pUoN0tKRRyTLpUW
LfaxNUxIoUiWKWuNaDGRSvbtvo0yypL5zR/vec/nPM/3833PeZW9jjj5iggJ
CbnS5/9b5eJ8drpC3a5+gxkTHZ0+sK4Wtx/kLECaJ2evsGw/PJfQWyzqvhJL
jOvbMpYOgMlFs9p8Bxk08K4zYA0MgO/WGQltCSWUTjzVWFUyCHbqvv+stVBF
huSU+eLoISh9dSq02UAN+dJS11B/GAxXaRSLntmCNafrmQZkGDizl9NyKS2E
jRc2myePwJGGIZ1iu234YyCZb7NxFM7vFftl/rAujo64HZAqHwXvxpN+2lw9
tLGUb2Vqj0GduKxjgJsBPmcVnFDNHANtyajkkRAmGrZIcTWWj4PnF6Mt9od3
Ynxlx6OmgHHISn8qP524C72dI2wVm8bh3gVmd4cTooy46AJTbgJCKo1YabbG
GHQtPrbZYwJqWnIn4/eZoHs+o0knZwL8LXZU+1O7UZcn1SnyfgJSb2uXBRwy
RafJcXb2WgEkfamOPhpuhgpFS6/etRTAYEFoGJdnjhmXHSz7eQJ4Yn75xtmz
Frh/+6OdWqkCEHbo+nQxyhIDB5REJx8LQMi4XDgzwgrDxupPNfcI4Jz0tF3H
r9Z4z8Os4N9ZARzd610/HWODPienF9hrCPRG6/rANVvUYCTbSzEISF252RWf
ZYeBre1ve3YTyFupEdGYYo97kjxkNrsQqOMpmm0qdMBvYR+TxLkE0pUy5Y3Y
LPQLahkSDSdQHTkjxH3BworYNzLWsQSUbcTnNrg6Ijtl0SgkEbjfps5/+8oR
58My/BRTCSQ0/NBScnfC5vC53vI8AleOjRzyfOeEz2quzl4vJhC5RmGFitse
fNNTI8arIPC7+k8tIR/3YFLSunnZGgJiX4Vl6w6w8VmzOol9SiBf6+7p9CE2
1gZUSv/ZSkDdq8B+PtAZLVxMr+9qJ9Bt9bDUecwZ1xv3ynS+JtCRfUITjrlg
t8pAG6uLwLd7jjtVJl1QsxMLVr4nsPCAVxYX5Ioxb29KH+wlYKNlbhY55Ypm
VtsdU/sJNLpv9Q84RmHt6uIgxyECRSVHpxcWKNTXES4UGyXQktQdPMwhwDwn
2ZeW1wdMcVZJ2UMK73gqrysfJtDMjHk4zvkAIsv9NogE9kFUaU3CpWQKi75z
9onR86/7Pr7sekfhavXyl54jBJoSOiynOG0Q774gp8TqA9Vz0zPzsRTeUl2c
+Afdf11bxsSylfYXTK5ZRu87xt32DRRQOKy/ye0+zQ1awRwxbg2UnhQaLjDu
g395q0t1Iynk7NikdobOc+V1v5+up7A3JBx76X6V4CfyOQMUFltqVp+n81Ep
CUV60xSil+HtHprrjhflqHGLoV42wuMIsw8OcLMkXoTT7/ERd+HSvnbJerbc
YwptjdorWmkf38/pOPRSmPtaoyOUzg9fklvSMUnhtwumS5ppls774ZQyS+HW
DAWdSZorflZXDORmQG2WlHe7Xh+s8uZWpYVSGHdi0cB+2h+3LF5kuIrC0NLZ
sVrajwdlN0p8oNBnysM5gM6n8nsrMj7T//Np3f3HND/LfCrk/53CQtPln4do
9irUz30+T2FEW1DwHM05f5MuPkMFYzqn7lTJDsH+lNIHfMY2tDLO1qTWj4Kf
9URMu5khKl+95Sy6fRwezF3J5zNMsLjgtVLnjQkwYmsnMirNMWVqYrDplADq
C//itZvZYFPcJROpFQQsFh2mQtodMHI+qtbmEAG/O6IZfIYTZvoPrZLJJqDm
9lLSPYSN02XRu6MbCIyK/hbLqHRBJc9Urv87Av8BHsenZQ==
"], {{{}, {},
{RGBColor[0, 0, 1],
LineBox[{1, 71, 65, 60, 56, 53, 51, 2, 3, 4, 5, 78, 6, 7, 8, 9, 79,
10, 11, 12, 13, 14, 80, 15, 16, 17, 18, 81, 19, 20, 21, 22, 23, 82,
24, 25, 26, 27, 83, 28, 29, 30, 31, 32, 84, 33, 34, 35, 36, 85,
37, 38, 39, 40, 41, 86, 42, 43, 44, 45, 87, 46, 47, 48, 49, 72, 66,
61, 57, 54, 52, 73, 67, 62, 58, 55, 74, 68, 63, 59, 75, 69, 64,
76, 70, 77, 50}]}}, {
{RGBColor[0, 0,
NCache[
Rational[2, 3], 0.6666666666666666]],
PointBox[{78, 79, 80, 81, 82, 83, 84, 85, 86, 87}]}, {}, {}}}],
GraphicsComplexBox[CompressedData["
1:eJwVl2c81Q8bxg8hSsooI5VUQrI6Sf787tve+zhOKHtLtogQLQqhkFBWRFZS
smWHzJBRNAgZFaHh6Xl1v7jeXt/r873325wzsqcnkUiRdCTS/69gzJ/MtL31
hJSktlRychtoVW/W++y4AVyK27Yc1WqEzm0yLAynd6CsawT9TGotKMWo1j3U
58HFV8t0NbwVYC+xvk1qmwAOfp/d/5qtDHRF7Jd2qR9EwS4WN/bLRVDaH+zb
JiuMRbdV9jKdyIP/OESLGELEUM/L9/63gUxw/BWXmk2TRPpEnbpmh3Q41zx1
rEhXGjOXHsaf0kyGSDOmQ3/cyHi21liDlSkBbFsuOEg5yeDtsr/HpvfFQv1m
XkPXU7I44cJ8hZk/EqTYQ+988ZHD28++DKu2RYDVdwUxPTd5DGuq/FvYEwoP
0pr4f94ksF7Wzzn3ZSA8iZYb6TVCzL15aFbZ0Rd8XigYpOooIruzXcoFW0+o
bc/+dt1cCcsK+YoYW13BWf1ktTNNGRtKtqhdinSAe7lST11dVPBD+2TX01/W
cOt7dZiHvypqYVXH7z4L+PzI18/JSw2nBD2OnPEzhUa1uKRLl9RxXGLD0avZ
COj0307GhGpgGlP05a9hekBSLKe7H6CJw+XiGyyfNOAa90/d3qtaePF1/CvR
NWXwMLNt+BmujQqePI3qKggTYWQ7uK2DGQljW0jb5IArPvnt9Qe6+N7utDX/
dzLk7BANaEnRQz3Vix0BvFJQ77VP9XC+Ph5WSHBv8DwCaQL3+RUoBkhReqEu
YCYE1UHrJKduA9w16+4zri4I+7U3/z5gaoik94ux2xz2QNlrkajhfkP8/d57
2PUPN9xo/ispcNoIQ4RM7pOBE+K9v7hYjRphpCpVZpm6HYJ27t0ueMoYozvP
JBemsUCeCGO7z3tjzNMg9TsmMgDTDzre+jMU3Js1QDN4RYKHkiUX06YouDLa
KP1yfp0QsXmk9+esCfaUxiXwFi0TI5oVpSazJugsFF3J/myJ6M0MPAreVFx0
bnCQ850jVp8Yygt+o6JZEFdRktQUsfHc6+kVd1Mc1WZn6VufJLQl1VSDlk0x
41+R+SLHiZbTEs6u3jTMfWiUoTY/RBQWe/zc2KAhnZNUdff7PqL91ojntOMC
5OL9E72tzSDw59m9PzLMqLv7Sp7VVD1c66r+aELlwrGctO5hkSpg+sk/z3Cc
Hyvy7O2r3MvhVEHRxgFRQXw1QR9s1FwCpBDPY4sghM+jbAwuSBSA33djjilz
UXzUaBA1Yp4DEQzCy4eTxDFra+XDjqwHMGtV7SRqLIWS1hO7A8XuQdAOfzHq
6WO406NGzirkDqzsqFdNdDmOgXyXWbV/3oI4x3bzLqsTWP433iNk000IuNHh
Lep4EreeeWJMr30VNOgtJ1I8/kPdUOd09tJLQGoLyxIMVsCsx+X7LEgXgeFM
qnNXMuDhxaOaV+f9IeFg1c7nrxEfvVBt0tHyhrsNng84+xVRxjYl5wbXOeC6
vjC4NqyEus/obzcxOgP9+eMVLhPKmKYVpdux3w52dfnvnl9UwZ3+aa5UXiuQ
f80bNPdDFf12ll4KpTODgdXwq77ralhxrNq84jcFJrXKFa7TaWBnFWvJgTkD
IGlofOP+o4HO+W23q6K1QaXgzeEKkhYWcqbGntRRg9+yR5k/MWpjsSejBjIp
gfC8dnwKqw6O1u24EcivAIFbHNPtduuikvbfLzSLEzDWM5GXu1cP3X8Gb5f4
Kg1l9Qct/Q/oY+6nhO05GeKwMZdoOc1lgJs9RN87NgpDoV5ntFGcAUav0yU8
ljkICx9C/zZxGqI9b3eZk50AuNxMevn3tiHSUu8eT4jmAwFtrvHfPEaYX2O+
8DZoJ+zY80O37p4Rrux7Yl1xmh2uNW5SqOczxoZcYZEJE1aIS7/l9uqeMXIs
s5K//2SCoejwpbv7KDhbwn7uThI9UOOGX0lkUnAoRbU7kvMv4dT2+HOnqAnK
ANctf9FVgv5NSWhegQl2U1Q4fZO/E29WeD20JKi4raOgfKxugfh05HNZcSkV
T75TFf1U/IUIkpk8T5E0xeUMt2mBXx+JbB+LYLpSUzzYa3NHlGmCYGmRG/Qm
05C8FCMwnDhCyG02KH5aQcMiNXb1CaU3RJtceMWc4ztgLU+v0RRqhcmvY9bn
lzehxmSb9iDbSxCdX1spOMWBQbVflhPfVMNBZSlTqhIfcuyVDVd78Qzefe2/
8mj3ftQZs2h2migFIYprXLjyIVR9bCj/tfMxcH2cVqZpiWB+2uL6rpqHYLVV
xOhO1FH0DbidMFWWASLiwpumByTx75HUY7tqUsHPVaL40oA0bt8nHZ01nAh1
T22Pt8+SMdVyncXsZxwEch+YCPkog31eQ2SPO9Gwp91rUuu9LKr9J0RlyLsG
VdcXKhu+yWH0TL2/R3k4SI5bOj3/IY/fwve80mwMAaOLUi8nOQEdywskWvQC
IOvUhfOqtxAjd4yFPzT0gas1g+dpsYpo2NT67oCHB8zkaX2qv6uEkUPvr107
7QLllif+m7mnjIu97BPGuvZwqpIyXfBYBQt6Tn1u3GQNDhOzF8XLVfGjSeq8