-
Notifications
You must be signed in to change notification settings - Fork 0
/
SIS without vd.nb
5605 lines (5560 loc) · 272 KB
/
SIS without 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[ 278675, 5596]
NotebookOptionsPosition[ 276843, 5541]
NotebookOutlinePosition[ 277322, 5559]
CellTagsIndexPosition[ 277279, 5556]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], "\[IndentingNewLine]",
RowBox[{"Needs", "[", "\"\<PlotLegends`\>\"", "]"}], "\n",
RowBox[{
RowBox[{"eq1", "=",
RowBox[{
RowBox[{
RowBox[{"-", "\[Beta]"}], " ", "*", " ",
RowBox[{"s", "[", "t", "]"}], "*",
FractionBox[
RowBox[{"i", "[", "t", "]"}], "n"]}], "+",
RowBox[{"\[Gamma]", "*",
RowBox[{"i", "[", "t", "]"}]}]}]}], ";"}], "\n",
RowBox[{
RowBox[{"eq2", "=",
RowBox[{
RowBox[{"\[Beta]", "*",
RowBox[{"s", "[", "t", "]"}], "*",
FractionBox[
RowBox[{"i", "[", "t", "]"}], "n"]}], "-",
RowBox[{"\[Gamma]", "*",
RowBox[{"i", "[", "t", "]"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"n", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Beta]", "=", "0.7"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Gamma]", "=", "0.9"}], ";"}], "\[IndentingNewLine]",
RowBox[{"R0", "=",
RowBox[{"N", "[",
FractionBox["\[Beta]", "\[Gamma]"], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"tf", "=", "100"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"szero", "=", "0.90"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"izero", "=", "0.1"}], ";"}], "\[IndentingNewLine]",
"\n"}], "\[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[{"s", "[", "0", "]"}], "\[Equal]", "szero"}], ",",
RowBox[{
RowBox[{"i", "[", "0", "]"}], "\[Equal]", "izero"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"s", ",", "i"}], "}"}], ",",
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[{"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]", "\n"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"ShowLegend", "[",
RowBox[{
RowBox[{"Show", "[",
RowBox[{"plot1a", " ", ",", "plot2a"}], "]"}], ",",
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[{"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[{"{",
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)\>\""}], "}"}]}], ",", " ",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"Thick", ",", "Blue"}], "}"}], ",",
RowBox[{"{",
RowBox[{"Thick", ",", "Red"}], "}"}]}], "}"}]}], ",",
RowBox[{"LegendShadow", "\[Rule]", "None"}], ",",
RowBox[{"LegendSpacing", "\[Rule]", "0"}], ",",
RowBox[{"LegendPosition", "\[Rule]",
RowBox[{"{",
RowBox[{"0.2", ",", "0"}], "}"}]}], ",",
RowBox[{"LegendSize", "\[Rule]", " ",
RowBox[{"{",
RowBox[{"0.2", ",", "0.2"}], "}"}]}], ",",
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]"}], "\[IndentingNewLine]",
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}, {3.565205413893314*^9, 3.565205456988779*^9}, {
3.5652054965510416`*^9, 3.56520559460565*^9}, {3.565205626407469*^9,
3.5652056380641356`*^9}, {3.5652056687088885`*^9, 3.565205776364046*^9},
3.565205847304104*^9, {3.565211009007336*^9, 3.565211175587864*^9}, {
3.5652680442214127`*^9, 3.565268078220357*^9}, {3.565268111401255*^9,
3.5652683648857536`*^9}, {3.565268399068709*^9, 3.5652684114434166`*^9}, {
3.565268456068969*^9, 3.565268492771068*^9}, {3.565268570906537*^9,
3.565268636284277*^9}, {3.5652694812356052`*^9, 3.565269547792412*^9}, {
3.565269593761041*^9, 3.5652696563066187`*^9}, {3.5652696908595953`*^9,
3.5652697156060104`*^9}, {3.5652697915233526`*^9,
3.5652697923624005`*^9}, {3.565269868986783*^9, 3.5652698836706233`*^9}, {
3.565269983819351*^9, 3.565270014109084*^9}, {3.565270046430932*^9,
3.5652701314807973`*^9}, {3.565270162578576*^9, 3.5652701900931497`*^9}, {
3.565389617101556*^9, 3.565389739684567*^9}, {3.565389820249175*^9,
3.5653900798300223`*^9}, 3.565390515397935*^9, {3.5684976211641216`*^9,
3.5684976258693905`*^9}, 3.568497659288302*^9, {3.568497711848308*^9,
3.5684977400989237`*^9}, {3.5684977859845486`*^9, 3.5684977880826683`*^9}}],
Cell[BoxData["0.7777777777777778`"], "Output",
CellChangeTimes->{{3.5652685850513463`*^9, 3.5652686377083583`*^9},
3.5652696587827606`*^9, {3.565269697437971*^9, 3.565269722190387*^9},
3.565269886796802*^9, 3.565270017329268*^9, {3.565270077489709*^9,
3.5652700917365236`*^9}, 3.565270134431966*^9, {3.5652701801975837`*^9,
3.565270192356279*^9}, {3.5653897122129955`*^9, 3.565389740563617*^9}, {
3.565389821115225*^9, 3.565389879871585*^9}, 3.5653899112033777`*^9, {
3.5653899431612053`*^9, 3.565390081153098*^9}, 3.5653905161639795`*^9,
3.5684976274334803`*^9, 3.5684976598273325`*^9, {3.568497694503316*^9,
3.568497741366997*^9}, 3.568497789145729*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwdmHc4lW8Yx5FEVJJUUoSEn052st7bPBx7i6QyQyUjI6LsjJBKKSoaRvYO
URKKJKJSWck8xw4Zv+fxl+tzua7Due/v9/M+z3vgzAVjBwY6Oroqejo6/JP/
xkp66v465fbZ06MfPowTlOpN+kNOa4T3TtUwhW1UomWLLAujDTucdItoyfxE
JVRvaNQ+M9gN9O36b9iVaITDkaUtElv4QOlt7x3XIBqhJ+IwxUUWhF9vEmv0
XtCIwo5A7yY5YXA8Y8U/2EgjFDhE8xiDxKBCf/u34z00wulf4oMnluLw+PXj
vPxhGnGh4Y9Unp4k5DqpxlMmacR1K6aDK27SwN9UlfhqjkbYvbvsKOEsC3Qe
Wyp4FmlE3aY9Rq7H5aD2xhNPs380QmJ78O0RL3l4X/2hfXKZRpyaURLTd1OE
oJDNxbQVGvEo9S3P31hleDna+N+zVRpRFCf/vd0YIJcxlX/fGo3wqlQyfKCr
An92Zfg4In7V/GQ6yloVNltJX0pBfJZ8rPqspRpM9XaSqxDffy5R4uqiDjoE
88/PiBNmqq+6+2gA2d10wwDioSzvS84emmBafMNlHHG9ZmLytWtksDbZIzWD
mN7gW/+NYC1g8m2y/4uYTqWU/qGfNlDqOf9bQhy5669eewQFXiUqHF5G7G5l
9/pviA7YHR9zXUHcd1XanrilCxKBWUariDlv3v0W9UgPev799xvzU3ZRv3cp
+tAeEQlriOs8eDUOZRvAlseJ1zCn8j3kUTI1hMHTodmYqwOW6JzbDOHTkcgG
zAd0Ni0LWBhB7OLWb5iLP4pEf+0wggFezyHMMQ2r4nw2xuBQ7UvDfNNzxOVU
jzGULdLNYQ7YuX8b/3ETeGKSv4A5U2Rjs1evCWQ5y//DzDRLv6fupClUbX+1
jPmZeMGV1D+m0NjXvoJZ5EyW/so5Mzi5GrGK+bt2RaHZmBlwuK2sc3u6/2HC
0xz4pczWMC8UGSnyT5vD7vnodV4r9ygJP28B5TL566wjrqkRMGcB7wea1/md
zZGzrp6WsI+rZ51z893/rq1ZQp7KyDo3J3y/OOyEPjdN+Z9J+ATBt1J2f0WW
GWp7jl8UukIlIlurB83MOaFSL8qmcQv6Pn95qIwyPPDdhk926DSNOJ6TtyYg
yg/TvO0GK8lov0EXpSYJIVgodFuJqKIRl2ZMOP5Yi4Jw5wuJws80IpRReO5Q
MgkOR6+dJw3QiLFT1c6iJhIgpqO94+c4mie7j5i5jRTwkppLDWZoxDx7ncYd
FxlY5hJeSP5LIxKdmq1bTx0F7g7xWtMlGuEX88FT1OkYHB0qa5FDfdBisO1L
cVeAt7yTDQGoD3RNVzP4A5XA5NZfcTPUB8aTD8623iVgk1uT/jDiJMGqneUf
AQ5kspzUQPO49/riox0dKmBv0vo6BOctita1+FUVsq68CXiBmMFXpsKlTw2E
14pNmxBztfrspU6qQ3TMe4UexIof9wSMz2pA9sFMsWHEnQshEd5LmhAzUyw6
ibifUqoURa8FJUazwXO4D1pa07tWtMDpidHoAmL1nC+HKugosI9bhvwP8bLc
YebfG3Vgr6dWD+6HMFXnZgqbLsSyjY7ifvhvdkqz36sHbw33p+M+/PjUl/l8
vz5k6OrvWs9znaCtj4ABCDmv2q/nZfyO7TCnIWy/wXl3PQ/6LXHGiYYg7m1R
iZk2ELz6docRCFuZfMLsEpv8ZvWWEahPdvZi5tPh/Lm82xguFLmMYWbfN6tX
e98Y5pzDpzFH1m9QquM2gS8djH8xJ6YluL2/bwIWBM8S5u64kKl7vKagqsG5
3gfzxK/vj6SbAj2r+XofnJteDLWImgGfpsx6/hm+FARn5pjBmF3jOn+Z3+NO
OWIOB34cXs/v7/+GivMLzWFhxmOdA2T7fU3FLaB3MG2dn3idCKQvtAAxllfr
zPJOvstT2hJk+TrWWX6TYX5JhSVY+Aysc5N8SMW40y/iSskwZVh8guif+HHa
d24DDDPbBlZTqIQodXE+5zgHtPanGdyYohKCahIW5qrcoP5kxiDEkEb8mugI
z9p7ACbDzMWmY2iEkKlrYojaQSjLpdIblaA8DQ6rWVJEQC6kbWdyK/I3q4jx
7ejDMK4mKTfzC/mAJLxhuFMcfoaJrPWNov64Hsm/1ikJ+x5KDvRN0YjaEjuZ
5jFpyP3PjbN8Hu1/l0Bf0KAsXMvtM69Fz4t9zR79lF458NnOwdeCnhdVUbSX
r6flQfHrNpBHfRD/aetcPqsI+0UFDx1E+Te+IvGmfwcBfSP17k2IM45f9tVI
AKjf1q0gieYRUdPlaxmvAotRf/u9EY9mUn7X3VOFk3TPxjIQl9oeVRi9rwYF
SiWH3iA+/tJ0OOeFOpwoWlDvRuzYN3aFVKoBIXF+5b8R//y0vHq+VBOSVuTt
qYg/awmObXxNhsrWbOdZxME6CaOfa7TgqumyOO7D40sNSWLV2kC6nVyLnxeq
whzZ3xsosMy9NQr3QTejK0u7UQc6uF7m4j6QGlrojbt14VdQZADug2+hkuqe
Xj1wOBCwhvnw2UNuab/0wXP3SxO871NPDzEEjBpAlqVdHOanZaWdZE9D2HZR
p2g9n/nVa1zThhDfLfcBs+b4R9uES0Zw9+7tH+v92JImPTZvBOUvto1g3v3G
X20i9i/xmMfFsObjBGHDlth1zZcN3k050t5mUYnfuVfbV425gO6SXFGbEI0Y
zmIOKzq0H76cuC/scxHNw+W67wM5ARg+lnqrIp1GfM306TiheAj47vCyMb9B
/WF41et97j+QooQeVutGz5sVoeSKjCMg4yJqvHOIRry4adzq1iYBnc+kd2ZR
aURISc1zkd9SkGscaJk8SyPE3Pab24/IwFGDL0H2CzQiq+ziF54fRyGxTfH4
d+TTTn2hjZupx6Bui2z2XeTTDv9WJm+qAoiZJnFUoPy0FaX3xmxUhjQ2Nptg
lJdvqYzMP7YByCVJcjOi7x9KWy4U3aoCqaI9FeaIcxpYm1t3qUJD916leMTK
zknNEnvUYOuH00QJYoN4qaNjB9Xh4LMXfa2IO3rv6R+R1oBlJ+6WX4i3V5P9
sqQ0oX83q+Mo4g+8rOwBymRo0fN+M4V4Wlcrnk1FC35ZD16Zx/lRX1u0VtCG
+ubwxEV8PukbqtRQp0DwLoNa7NNSbS0HqroOJOs/P4TzwihXe3nCUBe8B9kO
4Hy806dT8jmpB8TQwFvMqlNh2+Vs9GE1YVxs3Rfc/MYGTigvorFemFuZeGoS
pQ3h3KW6x5hLfw+QbQoMQWigsnbdh7vUV7xkjKBjfqQT8w3S6UNuJUYgL9g4
sO6zyKiIrQrG8MskdgJzY0y71pzTR2JXdvx0P90EoT7kXl/SQQ8UAfHyXSQq
sbyyK9jSejscNKp9wd1PJaTOvmRd09oDlDB9OXYy6kucU/O5JD64v/+Yxelw
GlGpPKyi91YQyBYZ1i4F6Hx5oqgnvk4YUrgs/Nvfo/8/LjNYb+NhoHOLWZT8
if5+kk2b5h1xYHVmMW0bQfNNZp1OTpAEI96EvL3IPweZZV7Y5UjDxrzpMWXk
n/MXPqZbPpWFql8Xw92Rf2jDb/d8vy8HrBFivxKRfy61HBX5UiQP/Il+F7lQ
fhqH9Skh+YqwVn5bgBXlxz/uH228TRl4zx+sK0csNvtNhM0DQL1Oe+oQmseV
quuD/udV4Fa47+nziE3lF8YOXlaFc9Od+g8Rx5vILWf6q0GYrtPHV4iXJN7o
s8WoQ8ysAt8X/DwYNQ3+eFsD/EjOjoOIzyhXVLomaUJstMD5CZwn8ncV9cdk
MHCo+YTPq9BSMlSQpgV6H21b8Hn1k7WyjmOKNtyyb/XB/vmdFrps/owCnpVq
J7B/TgY9KOV7qgOMj6XicJ6YBemLGkt1oSoy2wnnp+Gxa9ynWj2IEH4+jXm0
6rlHUbU+RGnZUfC+42PS0z81G0DPnYeRmPnd/vN7bWcIwh7BeZivc50IJvoM
4YqLVxPmVQ3pmsmHM8T72I/L9jkTRA3jQflai83gvs35P6VbVKK0Js06y3Qn
+KmVZ2dz04i/bYnHTpD2wWmOBrffZ2nEbdZ9VacK+EFbtb3qWiqNcEuaDRGq
EYKGP7rPfryiEYcuR/nlDorC4DR/1+Yv6Hxz9VmD+/4j0L5hsTt7EPl2KSLl
fYIEHOmKbWqaQOeXgXujxCMpWLsR7baAzmsJNbSWlUwZEMxcUaNHfkndtv/u
XOpRyGr6mn8X+eXZOT/Fe1nHYH9veagL8ktcDZdXZJ4C2FM+ZKWgfDjXJL62
qFOCDG0OK1eUhx69z4r1PQT8/vghYR4x2/UT9lNUAFPSXlV9NI/0+hqZqTkV
SGw+JXodn78XVD7dWlSFPlropQJ8XgqVXsjaoA7KwTKkD4i9LVhK4jg0QEIm
gP0nYqGGfH8dTk24aWZkOILYMObde18eMhDJPsv4vDaqszZEOaAFDd7bN2C/
QIPA/DEebfhR4wvYLzvfnbwYxE8BZx6bMOyX43WVpWcP6kBe5x8GnAeuHeER
9SRdMIorZsL7f5gjRZ1Q1APV0blizPF9RTaTSvqQGHqKH+/377l0SzsNA0hf
W3DDPLT4Na5a0BBa81lSMUvesxzremQI6bseVGOuCuwPfCBoBG50uZ8x37mj
70SL/UOUSVkvq9tMEFfnKGvyXEzApf7A5qQDlTj+8/s5essdsI2uiHcPHcqn
+EUnJudXRK4Kjev0n3HC60ylt2k9HQQNmb2a2U8lBjc+Uyf1sENw5I8phm9U
Qj6AztHkw25QpM+9xa+K/PlL3O2jIR9wJ5D9Zq4hLn66RzJcEASi3VfO5KH8
Ju8WORwoDOVNN2NPNaPnQ9vILpHXYhBhkpqi/QM9L7UWZOd9xKF08amuKvJN
ft4iR8YFSTgy7R9ejO7HlgX7zubHSAME222mQ75xTjm4USNCFhh0avnJyDeW
vhdkBC7LgTo1OsEL+caU9v2XbZI83OmJ/bEB5YnzvOZL5zhFoCSG3VtBHHf7
FJNHnjJ0pCzvyEd5cq+7fu3CGQCy2o97/Gh+fgr23DE2KiC95Z+fC+4fgwHN
wlUVEv6Ijj7APmkp74k6qwaOdXnba/B969wSveJldQj/N87Qie97F1jo+SI0
YGx0cxm+H2d43dntEqIJDw5zUvH9WOud7Dm7RDJsETwbgn3Dt7VxkC9eC07v
LL+LfeOs7y8Ucl0bysgzRtg3pd+VniTdoUDBZhYN7Js/bsGKbLd0YMVfOxDn
a7b9HivjU114Fdp2Audp4t614sMFepD6YG0Mc6PjiHd1jj4Ylodr4Hw8Dzrt
tVxpAINy6WGYPTYlXPbImiK4tXkeld2dIGjB58dOklkgaadC9u8oKpHWQOIL
7uSE4K/D9405kX8Gum7mF/BA/PjKyIwDujdaD9m8PscPwyzfhJtT0P/LUmTN
GiQEwXuHbeVqaIR+GfUt5aEoeC5lZPd30IiW139sdL+TgE5To5iG7n/3X/Az
XPaQgMqtdvZmyCc1i5GsOkFScFmCLPkE+aTeqTuKcl0Griqfrm9B979kBa+l
vqCjsFhByfBBPmFf4bfsjzwG1fuurhghn7j2VopyxCuAiJrf2Ri072K161JK
aUpwbK+16ym070kzORPtlwS8dCAiaYgNLdiqTX8B7LQe+aGN5pGqluzK8FsF
UkQv00fg82ptcse7EVWw8x0PysN9Xy3uvj6tBjaWpg3N2O+J4FxDrwFSuxVE
fiAuoDj8S2HUhAUrci2+/yW+fiLiw0qGM+8XkrBPHN9bcbhs1wLK3bbX+P7H
t0l29jmbNqySx3Zin+RU3NO9v50CSeGFjtgnoiuets926sB0QCcN7z+N7uyW
/H268PDh6gLef6x1i6aImB6wzR3Kwfsu/docokvSB/OIYZ7188WuQetz0gYQ
9Yb57Pq+Z+T+tvMYApn9TgrmK9zM/dMPB4iu5a7aMPIEUXgxKf/rxo3w7dy1
fB1LKqE65CLZ1c0Bl8+vXZtcoBJ1vrkZws55RCx//VH5b+NEDJ9tq4AiHYQp
aaSI7qIS3VFUobv32SFya1nRz04qUbNR4OZK7G64z86vkUrQCOmInUp0R/gg
VLv3i+RVlEdqWFiGoyDcCybzWeWi+cvqGPyzEoauyqqc9CZ0Pp+KGrj9UAys
tRiut/eg+88jjVE9Z3FopjpMLw6j/b90T888KQmiXUZDQcgXUxs53n0JkIbZ
ICaG3jk074I2mZc+spArx7lNEvmCDdTVyC5ycEb/H48D8oUqHellRqg8VMxL
Oy2h/FSKv/a0DlSEhvsHiuZRfpK5byYUpyoDrdvmWTbKC8fzftIjK4BThYwD
vGh+80uGJZxmKnCJ8WiAM+6749Kku50qCEbsL7yP57158cPiKTVQldIlVeO+
br79/bOHOgwPuF/swPeXAjsVrSAN8D5ZZYN98fZGOHnEXxNCZ+tysC+YIgYZ
PK+T4a/FVWvsi3jSJzvLSC1gF3UJwr4oJu9O2n5NG4LHMpSxL9L87vkV3aDA
8Lj0MZyXnIGLa1djdKD3YKYXzotL1p+yRw90wZiebIHzQrcr7p11ph4chZRh
zGUvo3xUM/TBup9VDedDzJmwksqbJJij+wT2xU8QzMQOddIrZthn3UnZEUol
kiiXGMZucUJEX8OrHdtpxLVzS7VWnjzQ3xSaPmmHnr+M70kzOvzwzcXCwPEe
jbBWT3wXaCsEEYaFt09Xo/OY2YHx3VdEoZsSG6CD/PB8IEKksoIE/5qWym4g
P6hvci9pPyMBF98EcnMjP+zbSHff8LwUzPK+qvZHfrCp+vDH3V8GrnQ6RxYh
P9wy4v3P0v0o3CgL13NEfqiv/+vA5HsMasvn/DXQfq2m2JgsghVAIfxKRija
rxUT0/iROCWYLKvrskL7Xfrq8jM6iwC92FtTY4inlezPhncBVHz1fURG8+DZ
EMYU3aMCZS+H7MMQLwqqf+ntU4WPV2lBuYjZFOYal0bU4K5U7Vbsh/PiPSrm
i+pArtYpwu+H5mNeKdQva8Ajo/FL2A/FP5Is5xnQfnvFT2E/CEYzP/dj0QI1
R67H2A/xhwjPXkZtyLUCOuyHhUK5SlEWCvgHb7XAftiR0ibVsEUHbgXuH8L7
tvexaju4Uxf26tFNr99vpLhnGgT0QPNF1jO832azL7f9hfShu4qXG+936rt9
mldWH7GN9XSIieIE0XzkGEtmICN09brbJhlRCVmXQk7VuxzAbTWzjzxHJcqY
RXjPOacRVsdm/GI+jxMSa2fmuCbXCLeKhz1Hd1CJdh2O5NYAdjBYvWLV8JlK
BFOKM9vP74ZLvn05lsoo798jnuzezwezfa/YDIPRvKO9X/4wF4SK7PpTHS/Q
/sPsz/RrCwNdh7WSOup/K71VU8sNMbih+J/jbdT/JP3xN36nxGHkqOb5TtT/
WshxFTOXhOeBksYOqP/K4mzzI17SEPbJrOsj6r+9opD+/AVZ4Go/63YI9X8o
nUT5eloOavRm9W1Q/yuZu4TfBciDuH3jwzmUjy0cdG0dXorQzpIoO4Pycf+E
WJfZbWXwO+ntmonycCZl68ZKc4DRksGQ/Wh+JVbvGL4aqsD8dKO/E+LC3W0X
A21VYaL7CSfuP1V4z9uCE2oQu9HrJX6f3vjVrjHtvDosJy0dwv0/3x5/+sRl
DeiozxbH/Tfft/rz+SVNMKF/G437/9vxZYt/OBkYL9mQcf8ba8wEvEK1oFZe
0AP330r9T3N3oDZY8/2Wxv1PamwceB1NAZZJRSmch61po72KkTrAPa/njvMg
xLU1zv2uLpgJXDbFeYhMvt10U2KSKDDI81GOniA4Ux0md15mRuenUFm+q1Qi
ju2eJ98lTniUNPT4wjbUx8LoriOmPHCr2s9k/Aw6D1sq3WFV4IdupkfehXdp
RG7f1ZVSAyHQjHshOFeF7hcZ8mtrbqLgaMu9d+0zOg/PWmybek6Cg5tOKNmi
fu8u5rqxaiUBQnwVpOVxdB+7rmD/2UEK+DeU6DqgfmcNNXCFesiA8c0BylPU
76QeL7Wbzkchd7Y81gb1O1W43y35wjGQWC0RIdD+Pi2WePr5KsBx0dbyYLS/
F3/UtDLDlEBgzzsmS7S/tb2MgocfE8CpVXpgFLGaRDyR9hmglfI4WBPN4/uZ
fyYG3cjnIjFJodi/XXssJn6oQhx77Sp+/xtxfGdd4W81iJlhKcTvf1XXMkn8
c+qguaPGB/d7T+64dfeCBhyb+mOE+63uyjZZsqoJT29HaeB+b3D95XF1oxY8
OZmYgPvdWPZGcplOG76FGs3h912peZumphgpsBoibID7zTy6+qCTRQcybY/2
4X02SH4u6GXXhRySLxXv81UqNN3i1QMxe8cneJ/ZijuuyOT1Ei4nPfY/kZ0g
vsRFl8QfY4SC3qAxMX0qIS72JczHjwMKRI8wn5uhEsLHP2+38TIFuRrq+vvY
jG7at2gSP1g9rwsTQ/t0fMGYFk0yBrHb9eP496OMj0NJqACTf7zX35c6UiZC
PmkoQLG8sLAjmreSqUQsqVITDGo1AvH7vRMpheXRJEmIrbi2aYBGI8qXbz6L
JqnCrS0P6ZPQ719nt3p80tABM1KqMP4+5A1ull6fDMDnUWcU/vz/AYqwhH8=
"], {{{}, {},
{RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwl1WWUlWUYhtGRRlpppLsEpVsaaUbpHrpBwFZU7G7sbkpFUMLu7u7u7nbf
yx973vf6FotzDud5hoZFywqX7lZQULDLj5wH+9GZH903MJtb9Uo68J3uwHnu
Y9mTZ/TJFOMQ3YWf3Dcyh9v0Kjryve7I+e7jqMqz+hSKc6juys/um5jL7Xo1
nfhBd+IC9/FU4zl9KiU4THfjF/fNzKMzF3o2geo8r0+jJIfr7vyaz8l8unCH
Z+uYSA1e8Ox0SnGE7sFv+WwsoCsXeTaJmryoz6A0R+qe/J7PwkK6cbFnk6nF
S/pMynCU7sUf7ltYRHcu8WwKtXlZn0VZjta9+TPvm8X04FLPplKHV/TZbHff
3bnG2Ye/3LeyhJ5c5tk09uJVfQ7lOEbvx9/u21hKLy73bDp1eU2fS3mO1X35
x/1OltGbKzybQT1ezxxRgeN0P/51v4vl9OFKz2ZSnzcyN1Rkre6fgWW7Y0Ve
J+/ReZWziAa8mVmhEsfrAZlzduiD6MvVehYNeSuzQmVO0AMpxk69kn5co2fT
iLf1Oqpwoh5EcXbpVfTnWj2HxryTOWEPTtKDKcHdejUDuE7PpQnvZk7YkV3L
njmHUJJ7sq8M5Ho9j6a8l1mhavZK708p7s1+Mogb9Hya8X5mhWrZIz2U0tyX
fWQwN+oFNOeDzAnVs0d6GGW4P/vHEG7SC2nBh5kTamSH9HDK8kD2L9933qPz
ZuciWvJRZoWa2SM9InPLg9k/hnKLXkwrPs6sUCt7pEdSjoeyfwxjvV5Caz7J
nFA7e6RHUZ6Hs38MZ4POL8c2fJo5oU72SI+mAo9k/xjBRr2MtnyWOWFndil7
5BxDRR7N/jGSTXo5e/N5ZoW62SNdSCUey/4xis16Be34IrNCPQ7QlXk8u8fo
zLduz5eZD+pzoK7CE9k5xrCP/iozQYPMpn6y4P//E9Y6C9nX/evMAQ15KvuU
1+WbfN804unsTV6Hb/Od0ji7wdh8ZzTJ7DOOrf7e9c6mjKcZE2jORFowiZZs
8WcnO1sxhdZMpQ3TaMt0tuXf0zmDdsykPUX5rMzK5+E/9H3IFA==
"]]}}, {
{RGBColor[0, 0,
NCache[
Rational[2, 3], 0.6666666666666666]],
PointBox[{435, 436, 437, 438, 439, 440, 441, 442, 443, 444}]}, {}, {}}}],
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{
FormBox["t", TraditionalForm],
FormBox[
RowBox[{"\"S\"", "(", "t", ")"}], TraditionalForm]},
AxesOrigin->{0, 0.9},
Method->{},
PlotRange->{All, All},
PlotRangeClipping->True,
PlotRangePadding->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.5652685850513463`*^9, 3.5652686377083583`*^9},
3.5652696587827606`*^9, {3.565269697437971*^9, 3.565269722190387*^9},
3.565269886796802*^9, 3.565270017329268*^9, {3.565270077489709*^9,
3.5652700917365236`*^9}, 3.565270134431966*^9, {3.5652701801975837`*^9,
3.565270192356279*^9}, {3.5653897122129955`*^9, 3.565389740563617*^9}, {
3.565389821115225*^9, 3.565389879871585*^9}, 3.5653899112033777`*^9, {
3.5653899431612053`*^9, 3.565390081153098*^9}, 3.5653905161639795`*^9,
3.5684976274334803`*^9, 3.5684976598273325`*^9, {3.568497694503316*^9,
3.568497741366997*^9}, 3.568497789169731*^9}],
Cell[BoxData[
GraphicsBox[GraphicsComplexBox[CompressedData["
1:eJwVl/c/Fm4bhpFVJDuJEpnZWzz3fdmyx4PHTJRNiKIoRYOyhZAts5IZIiuk
yCoNyciWUXbi/b6/XOfn+gPOz3GcJxwvml6gICMjC/jv/D95Y/7lZhxrRr5W
y6Upj2qxbgON4ZTLHm4rHCnRpn+Ouw/K76e0YwTB0xGCTpM5WC1Gs6nAiAOo
m/iTpu1S8AWJ7YNSB3nAhDS9L4o5FhsIX1hh1z4JqzRGNrGl4bh8MCTgraIQ
nHrf078YdQ0rM4s8p7whCkk5pA3mAV/s8jf+cT5JEuRtCHbGRi74Yvu0zHMD
afDKkLz2w/0sjrSm5v/nKQsZlV53voYSsVPHNWcpV3k4v6LQdTRKDzfTHDHx
sFKEudkmN6liNSzFFJo0638a3gUxpboFKWKHPwRRQ08V4Pk8UXL2liTOznjD
tRGFoH/VzFc9ih9XRJ/+1m8KMOXthkttuLB/HcH4sb4qeGhShNBJseDXXfm/
I2zUoOjZtVk/3v3YTVupwY2kDl7cLXcPLJLh9EKpKg93DViNtxz9y7iO4v40
3PS5ognelJUxJTsLaKo44LKrnxYMHi3Q83k4gdq04lNu3dKG+CvYqvvRF0Ru
9HU8JlQH7mbNHvKS70NkqtXkWUFnQM2vdbHUrAPdO7xh0H9XF17vG3LjS2pE
PtZOLRthelDwrEI2P6sSjd2UPY8f6gPd2B1R0mQxYk149DUi2wBOpo+owtNs
9IRRJKgjzRCEnixRLi0mo2a/45qCJUaQti1xzdImBmXwZHERiMaQwx95e/rY
bdQQvE3m2msMFHNkbDP1weiEHs0On6UJdO+bkg/b9EOVH4Tvfxk0AZa2Io7R
B27oQfuuJI+dKaxK8SWp3XVACZdm3R2GTWHrokgE/a45CmY7dojXygwS2JUt
XiADVCRM1eU/agaHHJfAjl0dUa+SH2m2J8L74I1lX5rTqEDyxfWMaSIYspES
VlylkLBjseE/L3PgkORvOxIqiL6dqS03nzeH9RW760iYB/XnXhXDlywg4Zfz
TqoVG9qsMFHh/W0BEfMWaYEUB9HeS7+qO96WkCfb3BD7ex/Sk9TSDF6zBAMG
21Mc0duEDjsJN49LJCg8n+a/tLpCeFbms7G3R4IBOtML+o1ThK64b74zLkvY
kIxOwKG7AvP8q0n/J08LjLKOcynchfheT8NPcwtWmNL5ZLiqlYGpN7gWKeW4
oERjVVmpNQFblT7f4xPhBbl9gXSlTRGY7IavzDIWgFTVgTyykVB8+Y8Z87SN
CJz3YBOeFriMwymF1gRTxGGs158srs8Tzzs0uIqYScGfvsrqhRNOOJjxiqiF
nQyUvMr1iHewwuuMzZrJ7nJwn6z/+aEeYxzv0mXT46AAo3efq+yVauGgB+8v
ibgoQSpnD8QrIaxDcXYszUcZ1IrIb3QekcVkb2/m8YYQYNij+M5yqQimtH/s
1vMIQ7YAdmO5dgInnnzF9vIDwAGab12IeBintvhmswyqQrsIv4f+oUOYNWJp
aOuLGriOmvF8yKLEFIFyte5j6vCY0i003/wvYu+5cnRxWQOkU2qUOMdXkMqH
I8ELq5ow5mb01WJgGn3cDLsbsK0Fo4LfOZKKf6Bx3WpCBLkOMLLctz28O4jI
dHR+H/6nA/1+Ju61Fu+RRuknwVoyXVj5eyn51EQL2lEUo52k0gMhkT/Or7hr
kdCiXkIavT7cinZIe15Yhq4ecMk8f9QAvFsIfzh78tH3vrGiwmOGwN732l2e
/TGqbD559gqfEbxScN6WdkhAewvJZ2dYjaFTSOL0cF0EembYHW0abwzpb5+X
87rfREsTobtvWEwg5fK9WVOLK8g9KqV196EJfCWzrzkkchHx6LGO7HCYgs0F
qij9h+cRI/eqQVO6KZx4vk0n/sca3WvbR2jmNAOvKpMzFusmKD4zzvNduhmE
6nTFUMRoo8/RYSupx4ngUH6s/aMdRhbxX95J5BKhJpOzs8BIDrm+fTrVLWIO
flthBY9LRBHFpxehRaXmsMrt9KspiA99Wj/ioythAdn9H4VtNDjR5KmpyrJy
C/jc/c86NJsJBcuPBxIlLcFRa/dLOzUtyve3DSEvt4SvjINvKghkaH/H6aFL
siSgyyweGl9aI5ymMS6rqiVBUqPIyfx7C4S3p8NqF1x+YDZiXWFDYDUe//X9
XODaPojwPSQiRCzBIotb66VWzMA3QE5NN52FT6pLWVqoccI9y/vkSUsP8Y9f
g3eKj56AxjD7Of/WB1iA6BEfps4PVKbE8NcltzDrzxl1kq4wEJwaD69VBGIH
OmHTpPti0EE7MuVE7YOFxYX2zXyUhPHjVz+RHbuAL3tIlN36KA0FTa9TlVRs
cVOVk1zXvCyspRUKDsWY4quH+cZu/JQHiVnNoLxdHczd5TeuO6oI2Qx1NTfe
AX4VsVTf8vs0EJflbu6JyWPJkbOuL1dVQKDcriHqqRg2vS7VOs6C4UDsuRd/
XPlwntW1QM04gMpn1LvhL4/gu41DgaRYVehcf82Qrc6E54p0J5tT1UCxdgqz
v6HG1WcVlOfS1aFG9eQ5fs9dZFVPnCl9qgGH+I2ehH3/g5zH5q+LV2uCbu4V
3zvac2ikb2fXu1oLeg+7ULxJG0MDOifnqVq0IUkv1LFQfAiF6sXNDTT+xwuJ
macMKh9QzuX2RNGGM1DCIl3w+F8bUhNiLvnWrgtlfkuCMTX1SD9vqPhMpx6I
qomeEYwsR+Lt3eSmn/VhiC6LPFChEAWWE9SOjBrA4Iz2LIE2C4m5CXpm/jAE
CG5eTr/3EDk8EaQInjMCwhOmQWGXB+hJTfVH7UvGIMC5XccnFoYsyhr22H8b
Q2GUbVXCRBDSWvhwNu6yCfhdG5askPVF7gczZefXTaDRW0LddNkZ9fgG4eVg
U0jYjlDPGLJDDOnrA7F/TYHSrCBh/08zdPfrUMP6VTOoPazAuCGti8C27TsN
GRFkJnTHx3RVkXirl99mKBHmWqns3eMUEG3Ic2n5/ebATy+acz5JArFzPXMx
izKHVzJ5Tdrn+REpy43J+KAFdK49GtkT4EbDLpIrG/EWEMWURDBwZ0Gn1x5Z
KzNZQhZNFLmJwwHUn+RyRCPWEt4pDqvKklGg4mvAsMlKAkbGi+oPRTYJoeWv
H8QmkSBUdO239L5lAkfrVfVfURuY/Ojz4IO3XmA7+vihW4H08HAidn+2WT6e
fHazf9eUHVhhJWWWPw3PFNPerhA8BuIX/Bc4++LwgHtk4GNFPsi/7ZuSynQX
fym6MmirIgiyF2WXt65dx58pXo8GeJ2C8smvBLeRS7jon0BKbZ4EUL19+PHp
nBt+mmDa49krBZNCpHdz5udwWFVjofCkDPx5bE1l9MECi3oeszg/KwfDzK1P
yLwNcXGN7yeu7wqwmdtAnxavgT8aClAdWFQC+w9xRk+klfHg1R7qgEVlYPGw
t/NQkMa9FbmjD6gQjOXLYYltQfw1g5L2+yGAbNnAp7XBx3H40k65CIMqXLYu
udomxIZL2+m6eg6rwUzt+/21y3QYuSZ2SR1Rh1tvth8IllFgo1gZhXl+Dagv
chcVjtlEg6OphhKymlB/cu+d35FlxNSgHVQsowU1uRzTq48m0fvjdIzBSBtu
0NqcZwsaRr/1dWLpVXXgY6iCjgLFAArV2NuyUT4DB3a/DP6seoumxqbqNDV0
QaN7Ort1twlVn9G5sKihBxy5B+83tFQjSsWma7+M9cFC/MxLqaWnqMOQjHDF
3gB+EIrSmrNzkdrKbSZFO0Mgt/F87eWYik5z8poauRiBx8y3ImvvONRDzdUY
L2sMAatmTjlGd1H15IS23Qtj6A/gcY/OufGfr2n885czgWqxqwmUdAEoRvyc
oGeVCVyF+83fFDzRp3sRdxmUTeFn5juuUj9HZNrpbRv7yhRYVLuU9UgkRPGL
8cmKkhl89ev5+0bWGGWc1zua1WgGC8p87m4KmuhisxqVOCaC6r2gO8vhKuj6
OU1+lzdEkGaUlxBRk0GdbKWpF/TNIWT/ye57W8JoMabq0d47c9gUee1R9/oE
Ihp09WiYWEDs1Ysj0cWH0eJSWqBbvwUMy5z1WNs5hBqviCR7GVnCPUmjQd0Z
KhRoajLWOGAJ+UEmCccs/hFcmeWylsxIgOoZZshG/hBWbUYHvg6TQHp7Izfd
Z5bQ+aBfZ83lAzYquOoxs1mDNaZ82qoGycH/3b/g8GNP8c6/w6EkGyao+bXN
Zz6SjWXc6un2dI7AoIXE+iZ3Ms6JdunySuSBDtJLqWKXaFyHZlQN3pyEYbE5
drf8MOxkWzEc2ywExS/8B1vYr+Ke6KJQAyoxGNt9qxQ25YM7E+16tZIlYWMl
eb1Ewhn/TqH7nRInDaWk5u9ZnXaYn1buqVOpLDxdyuYUumWGvS9+yCU9kQfJ
vVDFKYIuXpp5c+RbuiIcapgwCKtUxZe7FYQ/VZyGujfeMv76CrhzxlA3rEwF
2kouJrYbSOCr0X+XFnoRXKnkLvp16SQWXf0qTO8HsP3sU8cn8qP4+qvIn1e9
/+uHn+hSrwwzJp7enOe/pgbirE5hugs0ONZMcafoqjq07fsuKkBPhrelWg3p
H2iAg3IN1fjSKgqeI4Z+SNKEh3w3+5IfzSNHVFvnkagFmZte/sNx4+i99jdV
jRxtkC6wkAuz/Yygu2rqRaYOKCxan4MLvajPBuk5p52BY4JJXu087WgyM3zH
okAXTmZIiY5zNyD7G4+reZ7ogXh7/pJjdAWiPUle0VmtD3e2Pl3Svl2E2nM8
ovuaDKCuUPxKyMcsNPeq0K+iwRA00xkXkyqSUOyD3Ny+LiMgXn1A83UoCvF6
ngpqcTKG1Zc3/CeOh6NIdttQPGYM8l46jV2a11Bl2vzTr84mMGRzPHGcyQ/F
eP88Oz5pArM8FTu9cq5IjEtFkdbLFAIL+A1kXc+iXhqfzM15U6CgzLb0mSci
is2rzvTuZvAy4O6zzyV6aLXlP2gsmwE9ZRnVcxE1lPqt5s3wRSK0d3+/O/lG
EX22KCO5bxLhL3vWqzcJkogHi9oUhpgDU5WtLY2nAOIMYBo33TEHT5kbDrkV
x9DQbfOYo2EWYLv1/cKILSt6NBjw7hm5JeCHIjwmhXRoRtDL5cINS1DjGthz
/kOBts1bhU9TkSAuZLRAxXKLcPLexva/cBKIXXv+vGptmbCrKdu4nPUHN8SF
HfmcXo4bKflPN1kegAP7b5wRK32CqxszbYqJbEAVt3Evsjgdb/TGK9mKc8Nw
BcHFtDUeJ9Fxv3J4wQvhXhffGZ+9hz0TV8MEGgVAjqFdNIk2FAteiwh69lME
ZKPueCdu+2P2mwXtPsckYPzUDrUg0QPrb99NexcnBRs9Qn7y8o6YbCJ1DmfL
AMvfHzMCh0g4rnGp+1+RHPietdm/vG2EMw4de7SWoQB1Gl9SVB5p4gKvIJXU
YiUYndpV/FqigqMb2f3vPVcG7ZLSx0W8Mti1Mb7FspkA2mYXvtnrCuNhgwGV
tmEMjDrkMhNJPJg+0vb8yiJAVr5s6isFdpzb1ii3sqYKeMCq2yXgIC7aVO17
uKUGlStUbfY5+7BQuOxm8T4NmO35jidntlCA5f6qaGZN+JBxfrJ5YxkJtJdd
1WPVgop9rsXPnk8h4wcd7wK5tEGLFPiBkWYEzentTeme0AGGNFfW71sDCNr5
1pW4zkBbSgAl924XYuuw973BqwvU4VrR9+eakVVzXbUbvx6cVbUwbZ6pQews
d+62iesD/eBdzprE5yirVGbxl4oByPUF22SX5aHYsQq7ZYIhhJOPey5VpKEN
r1ySk6YRfKZrPjDqEI+mtr5EN5w0hgr7A+8SQu4h6VTS/FC2MTyyomSvSwlF
r0LGQx6fNAF9s6iOPaHLaOd7cEBVvgloSPLRphR4ob+zPeVUIqbAaXrim2uI
E3rhJ2VeWGoKl0MT7h95bIXekCJ9rUTMoLVzrWU/owniuLXYWvXUDMZsXO3O
22qhk0JO185LEIFu1Jfh1iQBzQXFBOpVEsH+No2tg5osSlkqbXqkbA55I4Kh
1/VOoc1ymSjiK3MIjfhatZXEi9h4/TK1wQLkf62vWy1yoLmDzexRLRYw4fs0
S+MdI+K2pKYdw5YwfO8ni8tvahQ9YSEx2mIJjSm3D70f3iW8XKipl9IgwZ00
y96hsVVCshSHms47EuQ1//g2uzRHSE42dFmKmsYGjgLcn9mq8M013b3T7NRQ
3nhKgLO4CFuNfPMiJ7GAwm7m/uGdTPxgLdel7xAREvdef1q7DOj78N1jfIX/
7ZOVz8E9pxcJJncKL3gtkkDN9VJww4EZQrukrwu162uc8Pu2zkbbS+zvWBdA
bCMDnpTsA22Mz/BPqgIN8WFGoA0UHO89m4NPB5M5m73ngNTPWPdpfTKm/CHp
+cGYB+xrjz47oR6DKSufHJG+cxKW5b/ngmA4jkzhEBYLEQJJvqPytR+u4vDe
2cPCLaIQp/xJUtrGF3PobMqvX5GERW1qq4oGZ1z2fIs576I0rKQdPfAnyR6T
XnC7lT2QhajAPcX5PTPsmsZPpXlXHkSnB3b9J3QxKfCiHN81RfjwVau3VkwN
E5e+/TibeBr4OY/GsYwpYFZvrXrXaBX4Eug1cIpFEkcnOVD7PUfg4aIm37R7
Evs0R9666AhwUCfJMq/0KA5SPs/5wE4VOJ+QMprKmHEkhdGSpYcaLHJzagXc
psXB3S+HI9zUYag971DiNTKs57VNrnJNAyInlW5FXF9D1Bf3k/Pc1YRo3kNs
azoLKM8/mcM9TAvale3755kmkE6HvJdTvDYYNH+tEdv6jHgYOn/yxOqAJo3M
x5j3vcjV8KpAWOQZCOw4OCb9sh1VfyPkJybrwi8e1eKaoQY07RmqQv9QD6xO
jAlWslWi1f5UOson+uC4k1GxIV2MfqXeqhR7YQDnF6v5qSEbdTrPBjSUGsKq
j+iakFEyKrxxzn+nzgiiI+6/kUqORg7tk117Nsbw+Hh+zcDjcJR3pqRp4LMx
JK/pburtXkPul7aIe/YmcADRjK+F+SHe43fX7L6bgL9AdQrLD1dk6/tYlvWC
KTyZeNlxhssBlTeI3hb9aQpbiSbk2pfM0es7708nnzOD999lT1D76SPD9dVo
qxkz8PYgkh3IVEO907ZLhi5EyPW5lbj/ghJqai3zoFwmgh8pKbVzRxLlGc8P
9geYg9ufRvk7qwLo6tkLuoN/zMHhsanH29DjaJDV1V7smgW0jzQVvPjCinTa
JsJmtizAelr893EzeuS/fJFh/bIlvFd9txvsuw/xtoYr2PyzhBjyhkmnlS3C
uh9LuWwwCR6oHH9Rf2mF4EcTd82veAWr1nBvqElX4KVQ73l77f3Q9cRQP/9c
Ac5sF+cJ/cgK7L0/HUsDHuPdiaGEshdcsO6l1nRLNQFv2kzZtXjxQsntkLcB
ByNw9f4KG7obApB4fItqxjcUG9YsvtHNEoGJN7fEiqMDcHfLtJ3+N3Go/1db
p8foidOf8lJc85OCMvWV9LRiR9y4dY9O74YMvNj32DetioTbXD5H6EbKgfT6
LcKV88Y4Rdl/e+yGAlTScXz0PqGFGf/xksbvKUHkAtMGiw8Be4zWiTDHKoNt
xJRcZboMrlSPlCFkEiCPWt9Qd0sYL5srmp2px6Axn7s/mOYENrakbyD+ALgT
nD6o1saOM9RTPCgmVWEkzlIuGzPgkaaUwY5ZNei5OTlGyUOJ2XYrP0f+VoeX
CT0DH9O3UXs8uDaSawJvx8QlR48V9EL3wt80Si1I9aHIf4mnUXxLvvAVOm1I
uC3txvBpBDm/s2Z2Z9KBHTGHbvurg4iHRn61kP4MxPxb/WFa+A6V1qbqpzPp
ArE+ZVLuQgsS+XfpbAGbHhxlV+Yg3nyJMsncDpZx64OSUOf2AcUyFGXTrSUs
agDNS0PPQgn5qPpLV5i+uCE8eRXKuuSYjj4d/mnjJWsEAslbUdMj8cjvj+JG
P5cxLHfULyhzRiCfVWYCW6oxXF+55CdAdhN9KpNqlOIygU8GHw+VNlxGFVWR
XF8fm0DAgPljL2dvFJ3qq3PzhClUXndfKqU5j0Y74k0q80zBs+G4zKymNXrC
oG0ZxWsG2qKExZ6bJugjI6RAvhmkv6xi2Pujha7bXKq+IUiE6VqTvIoMhHju
yOY0lBBh8rvQ3LkOWRRq+3RARsYczCd0Jzu3TqGQtt/UNyrN4c4g12Y/NR+6
TN9O81bBAiJ8j3azJB9BM+kDkxl1FkAymleXlWNC77X8mzgVLaHgXQuZ5jUa
xFlq9cmnzhIWDz06yFW0Rxi1CnprqkKCmc/He5Nc1giuXLJb11tIEPRls+FI
7TzhOift+O+sCWx4BLb6WqtwuW9i2RcqKrjM+2VD41oxVptylx76zAz2HPt7
N3yy8JRQZc4CDRFeCUfWXvwB6FlH7kerTBI4IvOLntxLBF7fNq68SRIk6gob
3qibIYSmPXgmv0ECE7NP4pbnpgnNgc/yhFyf49z5EVVx8Vr8gOdsD58KGXQI
lupZRD3DnyMWBR6lM0Lk4Zg21pwc3EjFl/AvigMSFhpsFOhTsOxdNgKZBA/M
ZvzS1q2IwYWLt2/nOZ8ExS/HF6+4heMMeT2jv9ZC0JqQky0sfg07rERMJGWJ
wpe+b7t1Sb64KVtzzsBVErxlmFIPUrtgxnqf3CJ7abiUb2IrQ3YWr1Axd3wK
loUwJrI/fzERO7/olau/Ig/9itzPWxX0MD1oqGu7KwK7UzvzAw81rEYmXp8X
fhpsDumS2kUVcZ1kyyWbEBXwzTryx1lXEqdwJsRVZiBIfurDG6nJj5kLx8Wz
rQF+KFe4Ze7jwuvbxlWs5qrQ5uRzPGiZGec5by/7OKkBhb2cfeAALfY7sPV+
y0EdLD+I/fGsJMOrB5K+DfhpwLkj+6O1O9aQ/gsnVZ0bmuC5eDtbPWcBvYm5
oz17VQs8hy7TIbMJRH33J8WlSG34q0Nj/E/nC4oV73Mi3dMBQ+WPJafI+1Cl
Nkci060z0H58xsSTqgNlBqUGVcTowhhHdApJvhGVTvju3XygB/NVb3qrrSuR
e/F0TfZjfdD7fHku5n4xIjsc3WFTZADx1JI7epezUU19xBW1PEPgnh2sCMtI
RrJBTcNMFUbg57gkX7QZjW78abl0jGQM6x0rVJvfw5HkRP/o40FjuPxTh/Ky
bTDSOhtayGRjAuYXBHzmG/yQxxVS09EvJnA3Lxg7y7mh1aM90pznTEH+DNeC
h4kDyuGgC9QaNQWP2txK1Zfm6JjhWp6tnRkwjHruezCkj3YOhYi6/DSDbtaJ
ZOJXNTTvJUFJ60SE6NkWxY9VSkh2WzKgb54I3h2PzQKlpZDzly6mGV9z2Ilj
cvmlJIiKiTG7iUvmwHVbNoxp9jgys13RlL9iAQ81BVMEj7AhnlecDlTrFiB6
aYROuZge1RFXzjX6WULETS2Bghf7UHXQhUaPLUsYkU908cTbBHvX7IO9V0jw
rjXt61LFCkHUFVvLPF/G3QV/RtlvVGBazKIh/poW5A2LTd81FuBE3csU8w9Z
4cwHPpp73x/jW17bTdaXuOBpq3PM/M0ELET5TvyPHi8YnuEX8XSOwDYa8R0h
ZwVg7buHWWpBKHY0P7HAcV0EXBOsa1W/BuDCibvCdbXiwLZPnkrHxxNr0PhU
9TtKgUmUyq3JSUfMTUWWbuwtA77MTTfDD1hhu1fvp32uysHXv0TxpRRj/NDk
+CmSjwIcLD8x+dZZC7e1bVygDlQCv0ElU6FBArZeoae2DFWGXxvddMNfZbA1
NfWCRDQBlDJDf4aZiuDtL+4j94sxiH8QKZiCE/g34bzbnSEA/9WOn710hzHX
vtvU94dVwdae0aT9CQPeOqnxaXRMDb5kRh6ndKDE9Mprnduz6nDSVnIqc3Mb
eUsOq1psaQD1HasT1k9W0PqD18ptO5qwf8DqGUvENKr8nkhap9CGe1q0KAP/
QCfv0xYG7dcBY52LMgP1gyhWEF8apTwDh9ulH+7uvvvP/xXrRPbrwszPvX9v
clsQS1qvTPtBPXD/eJKF/sNLdP6KdS8/mz5wRawm/btYhihlOP+08xlA2wPF
U+a381GX+aekqwKG8OoT33naV+nojWe43e9TRrBDpvuc73gCUvRxLfnOYQzR
RjkmdO4RKObE8qRMkjF0eHd598neRI66ViN/D5tAjBvvWMj+Kyg6m8ln4ZEJ
9Hiw4twObzT0W40Fc5sCfQhpmcHwPAr0azWqyzIFrgXG19tp1ijZjV+34JgZ
nHzfbV5Xb4JkP7tSUGebwXFO7SoWVW10KuCvWRQfEYzj63Lv/ENImZ6m5nIB
EfJfCBpx0cmh1Xix72QS5iBf2rataiqKwKBfPbHMHFyDztpOYz70hHy5PkLG
Am5M5Xfx/zqCnH/adhdU//dvbzUxeDChPFSSriBrCUX9l+cKmmiQA+M3YaVq
S6iPNxZvpiVD/ozWxk6K//FobkHao2SN0MfUQuRsJIFteTQZK88CYeXb+Uz/
4jEcU7Dx04m7GndJKO0vCqGEvN8rtHVTxVjevZxV7REzuNh38USXZOEENRXY
3vdfnytPvXnLqooON0hJraWSgLqcq447bolw39k0EI+R4OW9qHstrLOEmVjO
A/2/SSDYdmISVU0TDj/ZNU37S4Iptjgr6kPThBpa4eNerplYQc5ZMt2pFkvt
Oa6xL+9hU7dVuVvvn+F+PeaUnmBG4DrKSEptzcGhupVF/d4ccGGl12lONgXr
f7ubz3GMB2pu0fTuzMfg7fsB9d8tToID08ee6chwXHj7vOP4GSGw22MzHbO7
hnvIrd92x4jCXM25ufeVvjjRcKE1yEESclg7a+NEXHATlHqIWkhDbTFlmo3c
WYwk6ddn/WXBddGWaOJAxOdVBAzXL8qDTVINydBRD0/liut+OacID353irdF
quE62iGhjuDT0GIx8mrLVBEfZCbrHfRXgWL04qKesyROtxUdMk9C4M3KUazj
xo8d0xio/tML6JxWdQuT4MJV1h0UX4xV4Y5T1tgTRhZcztHrG3JWDYy4WyM+
bdPiRaEjb17YqkOmRNoubx8Z7vzi1JnprQECk9erKBbWkHd/7Dnba5rAuSGg
TPywgCy4d0cKL2uBx0LB7mn/CTTpXN999Y42eFQssGj6fkGdjeZ8/uE6UHU9
ILLmaB+y1pju+hxyBjTjr1/3EO1AiZ2dEy33dcGRGHqDw7kRMWTOjarc0wM9
/X+GsSGVSICdIdrnkT60DOh0GbwsRgc+svT55RvAlY4JibK4bCQ+ZbVJmWUI
OPAQU3J7MtpeW4rlem4EcqPnz3wRjkEMkVKdYhbGYNpafNpiNxytP59avtVv
DB6v7SxpHgQj2Y/q/ZxWJtCyta+P95sfiuVrzt/+ZAK2LyW2/zm4IfCdkuI+
awrr/O0PwNMBcWpd9TUbMQVCb8oj+m/mqPvFfmtsYwbvZJsb5w4YoH4u1b9e
42bw946SiOqmGtosS4n+6UCEhef2fhEjSsibbWWrZpYIIfF5LWkGUqjz8emy
3xfNgUF13/KcrSBSM8ojBfwyh08tY259DDyImay69nSABXDl6NxXVWJD9Iry
PSyrFvA2bixc7hM9GqYuv5/jYwmVC/soKfv3IbKPH8n9NyzBiuH02FPnbQKz
k2v94wAS0EvdV7IYXCHcS0l6myC1jMNv3Vajf1KBWTMuLLNdo4UM+fLknV8F
OJo+9RLPZVbY9BDLpWPOwBrl94ckiFxgZtAZG5CTgIFESKZT/m8vnWlMdkyI
wM/Gbv6rNhIAFe/0weqWUOyZd3pvz1MELuj3+vmRXcbsq5aHVgrFwf3zZm5L
lifmqGSP2bWWArLPP7T69jnh95HK5wcuyAB71dZDZzkrXDzVzh7uJwfe5l8U
3SuMceKwv3qCqwKMy4txvY3QwhlC454pF5Xgui6DcTk1wn1bVZeCApXh6P6n
0UrbMvjptLpO0W0CxLh7RSfcEMF7RylPiuVgyLC8nvbT7gRWl4rFmQMA0R75
1fMSh/E3x79mRp9VYbc1a234KwOmHjpi+eu7GkxVfv1rFEqJ71qxNZdPqgPj
QQXqvwJ/kdpekTjvmgacMLVe52lbQUeeLdh83tSEvsZmBo5n00jDg365alcL
9nKK+xf8f6B9Hj/8blLpwKOoYrVvw4Oos6ZVeofsDPy6ELZQK/IeZTynWVmh
1AW6wpUlu9YWRDu3+/jjfj1o8Amn0Fx/idqlB16MMupDyi/3SPSgDL3OgLcP
jxvA+4DeV4+K8tEgi21WHK8hoJLvuTsz6UglmThSKWgEzxJwQyxOQOnioVQT
7MYQUn/SwSw+Ah17aX9aPdEYnoeN8V4wuomEDnstjrKZQMJFraoOiSso+PDH
R2vJJsAewnyhfNMb7WSUbHIfNQU5z6y9ON/zqK3jm2FDhinkeHJNXeywRozJ
ieqVXGaQaMmEfn82QRPKITUzGWZQfeXA8JKLNqK49Zo38QQRLJSiyq6LYaRj
5KNqlk8Ej3zeWW0xOST+7Ne5aVFz0MlRvj5wXRRlM6mNZzwzh/pTZ27U2/Ih
tEPJ7i5lAcVi+xnSOTiRwL0TqqWVFsCN7Ypl7zCh2Lr2A1rSlrBp3OubPEGD
brWkmXFXWoKLE/MDXj4y1OjRleAhT4JXZQcfx3evEe7QR1DM1JOA4rdytTBx
gVCiwnJd7vkorqH3oh/Srcafou9XxSpRAvkauY4UdwmWFP10+0oQM2wSD4ZI
vM3C7kZkh8gpiHA2+T5To4Qqel+zOdKSQoIDr8f/RL9bIvwn8gIHf5CA7O+J
umC1WUL6/bGazGUSDHwoIr77MU14n/WGzG3rP14tMs/8VZgmOJYo5Hf/I8G1
LuOwy1NThLzPS1/vi/OC2OzDmJHKCOys+yusT1MZWKVFGht7ZDGBKBUlXqcF
6hkluO7WGLJNK395X1wa5H20Qja5HPDLnYSC++JqEL/v0U7y3gHcUtLj16ep
B2Zn+aMY6quR9j5Pkn+fEWxV1WYlR0YhIasBJjt/InhyHD/xyl8ROT+lzLwv
bgrqzCTSQ3knNEeZEy7+H0DYjvTfF1FmQv8DFT/2Og==
"], {{{}, {},
{RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwl1ne4iAUYxuFDsooySoXslZEVIdl7k6yikFUZkSgjyh4pe5U9MsosKrTM
ylZGssqoEEp296s/bt/3/K7rPc5xnYNsrbs07JwoISHh1cQJCfFM5DnRSyE2
2m1po12xx5CXYlpRbbP3Vly2x5M4+GWS9iib7BfiY2hX7XfJR3GtmLbF+3P8
a0/gjqBP1gqz2W4XH0O7Zr/HIzymFde2en+eK/ZEkgR9ilaELXZ7SvCYvk1r
zVV7EncGfapWlK12B0rykl6C7/Q2XNMmkzTo0+LPgW12Rx6PO/17rS3X7Skk
C/r0+Lr5zu5EqbjRf9Be4IY9leRBfz++Tr63X6R03OjbtXbctKeRIugfaCX4
IT53ysSNvkNrzy17Ol20lJ4pPWfoJdluv8wTcafv1Drc/mbweXjcFfSZ9uPs
sDtTNm70XVpHEvGBdnfQZ9ml2Bm/L0/Gjb5b60RiZmipgj7bLs0uuyvl4kbf
o73IHczUUgd9jl2G3Xa3+Bq08nHnuVd/iSTM0u4J+lz7CfbYr1CB8vo+7WXu
ZLZ2b9Dn2WXZa3enYtzoP2rxg5OUOVqaoM+3n2Sf3YNKcaP/pHUhGXO1tEFf
YJfjx/gZpHLc6Pu1riRnXvx5aOk803ku1Mrzk92TKnGnH9C6kYL5Wvqgf2hX
YL/9GlXjRj+ovUJKFmj3BX2RXZEDdi+qxY1+SOvOXSzU7g/6YrsSB+3eVI8b
/WetB3fzoZYh6EvsyhyyX4/vJa1G3Hke1l8lFYu0B4K+1K7Cz/Yb1KSG/ovW
k9Qs1h4M+kd2VQ7bfagVN/oR7TXuYYn2UNA/tqvxi92X2nGjH9V6cS9LtYxB
X2ZX54jdjzpxox/TepOGj+L7UsvkmclzuVaDo3Z/6sadflx7nbR8rGUO+gq7
JsfsN6kXN/oJ7Q3SsUx7OOgr7VoctwdQP270X7U+pGe5liXoq+zanLAH0iBu
9N+0vtzHCi1r0Ffbdfg1fnbstzwbxp33k3o/7melli3on9h1+c1+m0Zxo5/S
+pOBVVr2oH9q1+OkPYin4kY/rb3JA6zWcgR9jV2fU/ZgGseNfkYbwIN8ouUM
+lq7AaftITwdN/rv2kAeiu/Z+Dz0XMH7Z1pDzthDaRJ3+h/aW2RkjZY76J/b
jfjdHkbTuNH/1N4mE2u1PEH/wn6KP+zhNIsb/aw2iMx8puUN+jq7MX/aI2in
XbfHkp/m8TG0c94H8zCfa/mCvt5+mrPxd5o90rMFzb2f14eQhS+0R4K+wW7C
OXsUz8SN/pc2lKys0/IH/Uu7Keft0TwbN/oFbRjZWK8VCPpXdjP+st+hZdzo
F7XhZGeDVjDoX8fXyQV7DK3iRr+U6P//y4zwzMGXeqGgfWO34KL9Ls/Fnf63
NpKcfKU9GvRv7We4ZL/H83Gj/6ONIhdfa4WDvtF+lr/tsbSOG/2yNprcfKMV
CfomuyX/2ONor92wx1GANvExtH+9v0MevtU6aDe9j6cgHe1bnhPodPs/Yf69
pKjX/wARvR4d
"]]}}, {
{RGBColor[
NCache[
Rational[2, 3], 0.6666666666666666], 0, 0],
PointBox[{620, 621, 622, 623, 624, 625, 626, 627, 628, 629}]}, {}, {}}}],
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.5652685850513463`*^9, 3.5652686377083583`*^9},
3.5652696587827606`*^9, {3.565269697437971*^9, 3.565269722190387*^9},
3.565269886796802*^9, 3.565270017329268*^9, {3.565270077489709*^9,
3.5652700917365236`*^9}, 3.565270134431966*^9, {3.5652701801975837`*^9,
3.565270192356279*^9}, {3.5653897122129955`*^9, 3.565389740563617*^9}, {
3.565389821115225*^9, 3.565389879871585*^9}, 3.5653899112033777`*^9, {
3.5653899431612053`*^9, 3.565390081153098*^9}, 3.5653905161639795`*^9,
3.5684976274334803`*^9, 3.5684976598273325`*^9, {3.568497694503316*^9,
3.568497741366997*^9}, 3.568497789194732*^9}],
Cell[BoxData[
GraphicsBox[{InsetBox[
GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJwdmHc4lW8Yx5FEVJJUUoSEn052st7bPBx7i6QyQyUjI6LsjJBKKSoaRvYO
URKKJKJSWck8xw4Zv+fxl+tzua7Due/v9/M+z3vgzAVjBwY6Oroqejo6/JP/
xkp66v465fbZ06MfPowTlOpN+kNOa4T3TtUwhW1UomWLLAujDTucdItoyfxE
JVRvaNQ+M9gN9O36b9iVaITDkaUtElv4QOlt7x3XIBqhJ+IwxUUWhF9vEmv0
XtCIwo5A7yY5YXA8Y8U/2EgjFDhE8xiDxKBCf/u34z00wulf4oMnluLw+PXj
vPxhGnGh4Y9Unp4k5DqpxlMmacR1K6aDK27SwN9UlfhqjkbYvbvsKOEsC3Qe
Wyp4FmlE3aY9Rq7H5aD2xhNPs380QmJ78O0RL3l4X/2hfXKZRpyaURLTd1OE
oJDNxbQVGvEo9S3P31hleDna+N+zVRpRFCf/vd0YIJcxlX/fGo3wqlQyfKCr
An92Zfg4In7V/GQ6yloVNltJX0pBfJZ8rPqspRpM9XaSqxDffy5R4uqiDjoE
88/PiBNmqq+6+2gA2d10wwDioSzvS84emmBafMNlHHG9ZmLytWtksDbZIzWD
mN7gW/+NYC1g8m2y/4uYTqWU/qGfNlDqOf9bQhy5669eewQFXiUqHF5G7G5l
9/pviA7YHR9zXUHcd1XanrilCxKBWUariDlv3v0W9UgPev799xvzU3ZRv3cp
+tAeEQlriOs8eDUOZRvAlseJ1zCn8j3kUTI1hMHTodmYqwOW6JzbDOHTkcgG
zAd0Ni0LWBhB7OLWb5iLP4pEf+0wggFezyHMMQ2r4nw2xuBQ7UvDfNNzxOVU
jzGULdLNYQ7YuX8b/3ETeGKSv4A5U2Rjs1evCWQ5y//DzDRLv6fupClUbX+1
jPmZeMGV1D+m0NjXvoJZ5EyW/so5Mzi5GrGK+bt2RaHZmBlwuK2sc3u6/2HC
0xz4pczWMC8UGSnyT5vD7vnodV4r9ygJP28B5TL566wjrqkRMGcB7wea1/md
zZGzrp6WsI+rZ51z893/rq1ZQp7KyDo3J3y/OOyEPjdN+Z9J+ATBt1J2f0WW
GWp7jl8UukIlIlurB83MOaFSL8qmcQv6Pn95qIwyPPDdhk926DSNOJ6TtyYg
yg/TvO0GK8lov0EXpSYJIVgodFuJqKIRl2ZMOP5Yi4Jw5wuJws80IpRReO5Q
MgkOR6+dJw3QiLFT1c6iJhIgpqO94+c4mie7j5i5jRTwkppLDWZoxDx7ncYd
FxlY5hJeSP5LIxKdmq1bTx0F7g7xWtMlGuEX88FT1OkYHB0qa5FDfdBisO1L
cVeAt7yTDQGoD3RNVzP4A5XA5NZfcTPUB8aTD8623iVgk1uT/jDiJMGqneUf
AQ5kspzUQPO49/riox0dKmBv0vo6BOctita1+FUVsq68CXiBmMFXpsKlTw2E
14pNmxBztfrspU6qQ3TMe4UexIof9wSMz2pA9sFMsWHEnQshEd5LmhAzUyw6
ibifUqoURa8FJUazwXO4D1pa07tWtMDpidHoAmL1nC+HKugosI9bhvwP8bLc
YebfG3Vgr6dWD+6HMFXnZgqbLsSyjY7ifvhvdkqz36sHbw33p+M+/PjUl/l8
vz5k6OrvWs9znaCtj4ABCDmv2q/nZfyO7TCnIWy/wXl3PQ/6LXHGiYYg7m1R
iZk2ELz6docRCFuZfMLsEpv8ZvWWEahPdvZi5tPh/Lm82xguFLmMYWbfN6tX
e98Y5pzDpzFH1m9QquM2gS8djH8xJ6YluL2/bwIWBM8S5u64kKl7vKagqsG5
3gfzxK/vj6SbAj2r+XofnJteDLWImgGfpsx6/hm+FARn5pjBmF3jOn+Z3+NO
OWIOB34cXs/v7/+GivMLzWFhxmOdA2T7fU3FLaB3MG2dn3idCKQvtAAxllfr
zPJOvstT2hJk+TrWWX6TYX5JhSVY+Aysc5N8SMW40y/iSskwZVh8guif+HHa
d24DDDPbBlZTqIQodXE+5zgHtPanGdyYohKCahIW5qrcoP5kxiDEkEb8mugI
z9p7ACbDzMWmY2iEkKlrYojaQSjLpdIblaA8DQ6rWVJEQC6kbWdyK/I3q4jx
7ejDMK4mKTfzC/mAJLxhuFMcfoaJrPWNov64Hsm/1ikJ+x5KDvRN0YjaEjuZ
5jFpyP3PjbN8Hu1/l0Bf0KAsXMvtM69Fz4t9zR79lF458NnOwdeCnhdVUbSX
r6flQfHrNpBHfRD/aetcPqsI+0UFDx1E+Te+IvGmfwcBfSP17k2IM45f9tVI
AKjf1q0gieYRUdPlaxmvAotRf/u9EY9mUn7X3VOFk3TPxjIQl9oeVRi9rwYF
SiWH3iA+/tJ0OOeFOpwoWlDvRuzYN3aFVKoBIXF+5b8R//y0vHq+VBOSVuTt
qYg/awmObXxNhsrWbOdZxME6CaOfa7TgqumyOO7D40sNSWLV2kC6nVyLnxeq
whzZ3xsosMy9NQr3QTejK0u7UQc6uF7m4j6QGlrojbt14VdQZADug2+hkuqe
Xj1wOBCwhvnw2UNuab/0wXP3SxO871NPDzEEjBpAlqVdHOanZaWdZE9D2HZR
p2g9n/nVa1zThhDfLfcBs+b4R9uES0Zw9+7tH+v92JImPTZvBOUvto1g3v3G
X20i9i/xmMfFsObjBGHDlth1zZcN3k050t5mUYnfuVfbV425gO6SXFGbEI0Y
zmIOKzq0H76cuC/scxHNw+W67wM5ARg+lnqrIp1GfM306TiheAj47vCyMb9B
/WF41et97j+QooQeVutGz5sVoeSKjCMg4yJqvHOIRry4adzq1iYBnc+kd2ZR
aURISc1zkd9SkGscaJk8SyPE3Pab24/IwFGDL0H2CzQiq+ziF54fRyGxTfH4
d+TTTn2hjZupx6Bui2z2XeTTDv9WJm+qAoiZJnFUoPy0FaX3xmxUhjQ2Nptg
lJdvqYzMP7YByCVJcjOi7x9KWy4U3aoCqaI9FeaIcxpYm1t3qUJD916leMTK
zknNEnvUYOuH00QJYoN4qaNjB9Xh4LMXfa2IO3rv6R+R1oBlJ+6WX4i3V5P9
sqQ0oX83q+Mo4g+8rOwBymRo0fN+M4V4Wlcrnk1FC35ZD16Zx/lRX1u0VtCG
+ubwxEV8PukbqtRQp0DwLoNa7NNSbS0HqroOJOs/P4TzwihXe3nCUBe8B9kO
4Hy806dT8jmpB8TQwFvMqlNh2+Vs9GE1YVxs3Rfc/MYGTigvorFemFuZeGoS
pQ3h3KW6x5hLfw+QbQoMQWigsnbdh7vUV7xkjKBjfqQT8w3S6UNuJUYgL9g4
sO6zyKiIrQrG8MskdgJzY0y71pzTR2JXdvx0P90EoT7kXl/SQQ8UAfHyXSQq
sbyyK9jSejscNKp9wd1PJaTOvmRd09oDlDB9OXYy6kucU/O5JD64v/+Yxelw
GlGpPKyi91YQyBYZ1i4F6Hx5oqgnvk4YUrgs/Nvfo/8/LjNYb+NhoHOLWZT8
if5+kk2b5h1xYHVmMW0bQfNNZp1OTpAEI96EvL3IPweZZV7Y5UjDxrzpMWXk
n/MXPqZbPpWFql8Xw92Rf2jDb/d8vy8HrBFivxKRfy61HBX5UiQP/Il+F7lQ
fhqH9Skh+YqwVn5bgBXlxz/uH228TRl4zx+sK0csNvtNhM0DQL1Oe+oQmseV
quuD/udV4Fa47+nziE3lF8YOXlaFc9Od+g8Rx5vILWf6q0GYrtPHV4iXJN7o
s8WoQ8ysAt8X/DwYNQ3+eFsD/EjOjoOIzyhXVLomaUJstMD5CZwn8ncV9cdk
MHCo+YTPq9BSMlSQpgV6H21b8Hn1k7WyjmOKNtyyb/XB/vmdFrps/owCnpVq
J7B/TgY9KOV7qgOMj6XicJ6YBemLGkt1oSoy2wnnp+Gxa9ynWj2IEH4+jXm0
6rlHUbU+RGnZUfC+42PS0z81G0DPnYeRmPnd/vN7bWcIwh7BeZivc50IJvoM
4YqLVxPmVQ3pmsmHM8T72I/L9jkTRA3jQflai83gvs35P6VbVKK0Js06y3Qn
+KmVZ2dz04i/bYnHTpD2wWmOBrffZ2nEbdZ9VacK+EFbtb3qWiqNcEuaDRGq
EYKGP7rPfryiEYcuR/nlDorC4DR/1+Yv6Hxz9VmD+/4j0L5hsTt7EPl2KSLl
fYIEHOmKbWqaQOeXgXujxCMpWLsR7baAzmsJNbSWlUwZEMxcUaNHfkndtv/u
XOpRyGr6mn8X+eXZOT/Fe1nHYH9veagL8ktcDZdXZJ4C2FM+ZKWgfDjXJL62
qFOCDG0OK1eUhx69z4r1PQT8/vghYR4x2/UT9lNUAFPSXlV9NI/0+hqZqTkV
SGw+JXodn78XVD7dWlSFPlropQJ8XgqVXsjaoA7KwTKkD4i9LVhK4jg0QEIm
gP0nYqGGfH8dTk24aWZkOILYMObde18eMhDJPsv4vDaqszZEOaAFDd7bN2C/
QIPA/DEebfhR4wvYLzvfnbwYxE8BZx6bMOyX43WVpWcP6kBe5x8GnAeuHeER
9SRdMIorZsL7f5gjRZ1Q1APV0blizPF9RTaTSvqQGHqKH+/377l0SzsNA0hf
W3DDPLT4Na5a0BBa81lSMUvesxzremQI6bseVGOuCuwPfCBoBG50uZ8x37mj
70SL/UOUSVkvq9tMEFfnKGvyXEzApf7A5qQDlTj+8/s5essdsI2uiHcPHcqn
+EUnJudXRK4Kjev0n3HC60ylt2k9HQQNmb2a2U8lBjc+Uyf1sENw5I8phm9U
Qj6AztHkw25QpM+9xa+K/PlL3O2jIR9wJ5D9Zq4hLn66RzJcEASi3VfO5KH8
Ju8WORwoDOVNN2NPNaPnQ9vILpHXYhBhkpqi/QM9L7UWZOd9xKF08amuKvJN
ft4iR8YFSTgy7R9ejO7HlgX7zubHSAME222mQ75xTjm4USNCFhh0avnJyDeW
vhdkBC7LgTo1OsEL+caU9v2XbZI83OmJ/bEB5YnzvOZL5zhFoCSG3VtBHHf7
FJNHnjJ0pCzvyEd5cq+7fu3CGQCy2o97/Gh+fgr23DE2KiC95Z+fC+4fgwHN
wlUVEv6Ijj7APmkp74k6qwaOdXnba/B969wSveJldQj/N87Qie97F1jo+SI0
YGx0cxm+H2d43dntEqIJDw5zUvH9WOud7Dm7RDJsETwbgn3Dt7VxkC9eC07v
LL+LfeOs7y8Ucl0bysgzRtg3pd+VniTdoUDBZhYN7Js/bsGKbLd0YMVfOxDn
a7b9HivjU114Fdp2Audp4t614sMFepD6YG0Mc6PjiHd1jj4Ylodr4Hw8Dzrt
tVxpAINy6WGYPTYlXPbImiK4tXkeld2dIGjB58dOklkgaadC9u8oKpHWQOIL
7uSE4K/D9405kX8Gum7mF/BA/PjKyIwDujdaD9m8PscPwyzfhJtT0P/LUmTN
GiQEwXuHbeVqaIR+GfUt5aEoeC5lZPd30IiW139sdL+TgE5To5iG7n/3X/Az
XPaQgMqtdvZmyCc1i5GsOkFScFmCLPkE+aTeqTuKcl0Griqfrm9B979kBa+l
vqCjsFhByfBBPmFf4bfsjzwG1fuurhghn7j2VopyxCuAiJrf2Ri072K161JK
aUpwbK+16ym070kzORPtlwS8dCAiaYgNLdiqTX8B7LQe+aGN5pGqluzK8FsF
UkQv00fg82ptcse7EVWw8x0PysN9Xy3uvj6tBjaWpg3N2O+J4FxDrwFSuxVE
fiAuoDj8S2HUhAUrci2+/yW+fiLiw0qGM+8XkrBPHN9bcbhs1wLK3bbX+P7H
t0l29jmbNqySx3Zin+RU3NO9v50CSeGFjtgnoiuets926sB0QCcN7z+N7uyW
/H268PDh6gLef6x1i6aImB6wzR3Kwfsu/docokvSB/OIYZ7188WuQetz0gYQ
9Yb57Pq+Z+T+tvMYApn9TgrmK9zM/dMPB4iu5a7aMPIEUXgxKf/rxo3w7dy1
fB1LKqE65CLZ1c0Bl8+vXZtcoBJ1vrkZws55RCx//VH5b+NEDJ9tq4AiHYQp
aaSI7qIS3VFUobv32SFya1nRz04qUbNR4OZK7G64z86vkUrQCOmInUp0R/gg
VLv3i+RVlEdqWFiGoyDcCybzWeWi+cvqGPyzEoauyqqc9CZ0Pp+KGrj9UAys
tRiut/eg+88jjVE9Z3FopjpMLw6j/b90T888KQmiXUZDQcgXUxs53n0JkIbZ
ICaG3jk074I2mZc+spArx7lNEvmCDdTVyC5ycEb/H48D8oUqHellRqg8VMxL
Oy2h/FSKv/a0DlSEhvsHiuZRfpK5byYUpyoDrdvmWTbKC8fzftIjK4BThYwD
vGh+80uGJZxmKnCJ8WiAM+6749Kku50qCEbsL7yP57158cPiKTVQldIlVeO+
br79/bOHOgwPuF/swPeXAjsVrSAN8D5ZZYN98fZGOHnEXxNCZ+tysC+YIgYZ
PK+T4a/FVWvsi3jSJzvLSC1gF3UJwr4oJu9O2n5NG4LHMpSxL9L87vkV3aDA
8Lj0MZyXnIGLa1djdKD3YKYXzotL1p+yRw90wZiebIHzQrcr7p11ph4chZRh
zGUvo3xUM/TBup9VDedDzJmwksqbJJij+wT2xU8QzMQOddIrZthn3UnZEUol
kiiXGMZucUJEX8OrHdtpxLVzS7VWnjzQ3xSaPmmHnr+M70kzOvzwzcXCwPEe
jbBWT3wXaCsEEYaFt09Xo/OY2YHx3VdEoZsSG6CD/PB8IEKksoIE/5qWym4g
P6hvci9pPyMBF98EcnMjP+zbSHff8LwUzPK+qvZHfrCp+vDH3V8GrnQ6RxYh
P9wy4v3P0v0o3CgL13NEfqiv/+vA5HsMasvn/DXQfq2m2JgsghVAIfxKRija
rxUT0/iROCWYLKvrskL7Xfrq8jM6iwC92FtTY4inlezPhncBVHz1fURG8+DZ
EMYU3aMCZS+H7MMQLwqqf+ntU4WPV2lBuYjZFOYal0bU4K5U7Vbsh/PiPSrm
i+pArtYpwu+H5mNeKdQva8Ajo/FL2A/FP5Is5xnQfnvFT2E/CEYzP/dj0QI1
R67H2A/xhwjPXkZtyLUCOuyHhUK5SlEWCvgHb7XAftiR0ibVsEUHbgXuH8L7
tvexaju4Uxf26tFNr99vpLhnGgT0QPNF1jO832azL7f9hfShu4qXG+936rt9
mldWH7GN9XSIieIE0XzkGEtmICN09brbJhlRCVmXQk7VuxzAbTWzjzxHJcqY
RXjPOacRVsdm/GI+jxMSa2fmuCbXCLeKhz1Hd1CJdh2O5NYAdjBYvWLV8JlK
BFOKM9vP74ZLvn05lsoo798jnuzezwezfa/YDIPRvKO9X/4wF4SK7PpTHS/Q
/sPsz/RrCwNdh7WSOup/K71VU8sNMbih+J/jbdT/JP3xN36nxGHkqOb5TtT/
WshxFTOXhOeBksYOqP/K4mzzI17SEPbJrOsj6r+9opD+/AVZ4Go/63YI9X8o
nUT5eloOavRm9W1Q/yuZu4TfBciDuH3jwzmUjy0cdG0dXorQzpIoO4Pycf+E
WJfZbWXwO+ntmonycCZl68ZKc4DRksGQ/Wh+JVbvGL4aqsD8dKO/E+LC3W0X
A21VYaL7CSfuP1V4z9uCE2oQu9HrJX6f3vjVrjHtvDosJy0dwv0/3x5/+sRl