-
Notifications
You must be signed in to change notification settings - Fork 0
/
illustrate.nb
3167 lines (3129 loc) · 156 KB
/
illustrate.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 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 159112, 3159]
NotebookOptionsPosition[ 156468, 3116]
NotebookOutlinePosition[ 156843, 3132]
CellTagsIndexPosition[ 156800, 3129]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
The following show that the two positions of the same energy will have equal \
mutual transfer prob.\
\>", "Text",
CellChangeTimes->{{3.817630795038725*^9,
3.8176308752600327`*^9}},ExpressionUUID->"a0992156-0e6e-4832-b380-\
1d6325cffa13"],
Cell[BoxData[
RowBox[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}], ";"}]], "Input",
CellChangeTimes->{{3.808249437390602*^9, 3.808249468789117*^9}, {
3.808249632175363*^9, 3.808249635636956*^9}, {3.808249765110931*^9,
3.808249766532976*^9}, 3.808254098303404*^9, {3.8082563792955236`*^9,
3.808256381597042*^9}, 3.808256704271912*^9, 3.8082567676558867`*^9, {
3.808257135602641*^9, 3.808257140940946*^9}, {3.80825720411176*^9,
3.808257204717432*^9}, 3.808257874664433*^9, 3.808293009163722*^9, {
3.808293569355691*^9, 3.80829356980971*^9}, 3.8083195815073442`*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[184]:=",ExpressionUUID->"71fe55c0-6ac5-47c4-9f1d-f5b5c4efbd6d"],
Cell[BoxData[
RowBox[{"SeedRandom", "[", "0", "]"}]], "Input",
CellChangeTimes->{{3.8108453161099243`*^9, 3.810845327307097*^9},
3.811854242216763*^9, {3.817630310701293*^9, 3.817630312491869*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[185]:=",ExpressionUUID->"dcfeec0c-7fef-4bcc-91f1-b2d08ac21153"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "\"\<Sampler1`\>\""}], ";"}]}], "Input",
CellChangeTimes->{{3.8168124021590147`*^9, 3.816812402159195*^9}, {
3.816820205778687*^9, 3.81682020614793*^9}, 3.8172501919117117`*^9,
3.817285432667527*^9, 3.817435500804758*^9, 3.817435604876099*^9, {
3.81744024760004*^9, 3.817440253709037*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[186]:=",ExpressionUUID->"7a36d817-a42e-4aac-a83c-01bf9f19da01"],
Cell[BoxData[
RowBox[{
RowBox[{"HessianH", "[",
RowBox[{"f_", ",",
RowBox[{"x_List", "?", "VectorQ"}]}], "]"}], ":=",
RowBox[{"D", "[",
RowBox[{"f", ",",
RowBox[{"{",
RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.808247361688581*^9, 3.80824736168983*^9}, {
3.808255923249894*^9, 3.8082559242774067`*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[188]:=",ExpressionUUID->"5e089dda-ccea-477e-90c8-a57f6452c50b"],
Cell[BoxData[
RowBox[{
RowBox[{"GradientG", "[",
RowBox[{"f_", ",",
RowBox[{"x_List", "?", "VectorQ"}]}], "]"}], ":=",
RowBox[{"D", "[",
RowBox[{"f", ",",
RowBox[{"{",
RowBox[{"x", ",", "1"}], "}"}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.808248108301323*^9, 3.808248127514201*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[189]:=",ExpressionUUID->"b8801622-cfe9-4518-ade3-cf0a58e217ce"],
Cell[BoxData[
RowBox[{
RowBox[{"\[Rho]", "=", "0.9"}], ";"}]], "Input",
CellChangeTimes->{{3.811090304764001*^9, 3.811090309476204*^9},
3.8110905932548122`*^9, {3.8110936697104073`*^9, 3.811093669891964*^9}, {
3.811165658596994*^9, 3.8111656589701023`*^9}, 3.811166009060788*^9,
3.811462116928081*^9, {3.8114622534166803`*^9, 3.811462260894311*^9}, {
3.8114626082384863`*^9, 3.811462608414331*^9}, 3.811462667488598*^9,
3.811465341096808*^9, {3.811699221000449*^9, 3.811699221542183*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[190]:=",ExpressionUUID->"160a082b-82ea-43e4-9174-d2f623473fdd"],
Cell[BoxData[
RowBox[{
RowBox[{"SIGMA", "=",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
FractionBox["1", "\[Rho]"], ",", "\[Rho]"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\[Rho]", ",",
FractionBox["1", "\[Rho]"]}], "}"}]}], "}"}]}], ";"}]], "Input",
CellChangeTimes->{{3.808319417169717*^9, 3.808319426953169*^9}, {
3.8083195116541348`*^9, 3.808319512297307*^9}, 3.8083195999310637`*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[191]:=",ExpressionUUID->"5885a730-e8d4-4c93-a53f-0dabb5c31c37"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"U", "[",
RowBox[{"x_", ",", "y_"}], "]"}], "=",
RowBox[{
FractionBox["1", "2"],
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}], ".",
RowBox[{"LinearSolve", "[",
RowBox[{"SIGMA", ",",
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}]}], "]"}]}], "]"}]}]}], ";"}]], "Input",\
CellChangeTimes->{{3.808319896929839*^9, 3.808319909777285*^9},
3.8083201224381247`*^9, {3.808327222657301*^9, 3.8083272234173803`*^9},
3.817630036909416*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[192]:=",ExpressionUUID->"4ae4ac9b-3bd5-472f-b335-c86f09d90f8f"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"dU", "[",
RowBox[{"x_", ",", "y_"}], "]"}], "=",
RowBox[{"GradientG", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{"x", ",", "y"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.808319929482926*^9, 3.808319942972332*^9},
3.808320132836433*^9, 3.817630047740937*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[193]:=",ExpressionUUID->"e58f78fc-c868-4f7a-b5a8-0faefba57448"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"ddU", "[",
RowBox[{"x_", ",", "y_"}], "]"}], "=",
RowBox[{"HessianH", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{"x", ",", "y"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.808319946906897*^9, 3.808319956562401*^9},
3.808320141572105*^9, 3.8176300525308104`*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[194]:=",ExpressionUUID->"8a90332f-b6d9-4177-9b2a-6d963c44a290"],
Cell[BoxData[
RowBox[{
RowBox[{"K", "[",
RowBox[{"p1_", ",", "p2_", ",", "q1_", ",", "q2_"}], "]"}], ":=",
RowBox[{
FractionBox["1", "2"],
RowBox[{
RowBox[{"{",
RowBox[{"p1", ",", "p2"}], "}"}], ".", " ",
RowBox[{"LinearSolve", "[",
RowBox[{
RowBox[{"ddU", "[",
RowBox[{"q1", ",", "q2"}], "]"}], ",",
RowBox[{"{",
RowBox[{"p1", ",", "p2"}], "}"}]}], "]"}]}]}]}]], "Input",
CellChangeTimes->{{3.811699597302446*^9, 3.811699609942079*^9},
3.811758183453342*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[195]:=",ExpressionUUID->"2c8f8ed1-91f4-4f87-95ba-c1cf83261a27"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"dK", "[",
RowBox[{"x_", ",", "y_"}], "]"}], "=",
RowBox[{"GradientG", "[",
RowBox[{
RowBox[{"K", "[",
RowBox[{"x", ",", "y"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y"}], "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.8116996218701143`*^9, 3.8116996246631517`*^9},
3.817630298534161*^9},
CellLabel->
"(\:8c03\:8bd5) \
In[196]:=",ExpressionUUID->"3e20f0e6-6359-4f00-99b2-e054a8feb711"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Dim", "=", "2"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"vanilla", "=", "True"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"angles", "=", "1000"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"q0", "=",
RowBox[{"RandomVariate", "[",
RowBox[{
RowBox[{"NormalDistribution", "[",
RowBox[{"0", ",", "1"}], "]"}], ",", "Dim"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"H", "=",
RowBox[{
RowBox[{"Apply", "[",
RowBox[{"U", ",", "q0"}], "]"}], "+", "1"}]}], ";",
RowBox[{
RowBox[{"traj", "[",
RowBox[{
"U_", ",", "dU_", ",", "ddU_", ",", "q0_", ",", "k0_", ",", "qq_"}],
"]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"ND", "=",
RowBox[{"NormalDistribution", "[",
RowBox[{"0", ",", "1"}], "]"}]}], ",", "q", ",",
RowBox[{"dt", "=", "0.01"}], ",", "p0", ",", "p", ",",
RowBox[{"p00", "=",
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]}], ",",
RowBox[{"dist", "=", "Infinity"}], ",", "\[IndentingNewLine]",
RowBox[{"QS", "=",
RowBox[{"List", "[", "]"}]}]}], "}"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "\[LessEqual]", "angles"}], ",",
RowBox[{"i", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"q", "=", "q0"}], ";", "\[IndentingNewLine]",
RowBox[{"p0", "=",
RowBox[{"{",
RowBox[{
RowBox[{"Cos", "[",
FractionBox[
RowBox[{"2", "\[Pi]", " ", "i"}], "angles"], "]"}], ",",
RowBox[{"Sin", "[",
FractionBox[
RowBox[{"2", "\[Pi]", " ", "i"}], "angles"], "]"}]}], "}"}]}],
";", "\[IndentingNewLine]",
RowBox[{"p0", "=",
RowBox[{"If", "[",
RowBox[{"vanilla", ",",
RowBox[{
SqrtBox["2"],
FractionBox["p0",
RowBox[{"Norm", "[", "p0", "]"}]],
SqrtBox["k0"]}], ",",
RowBox[{
FractionBox["p0",
SqrtBox[
RowBox[{"K", "[",
RowBox[{
RowBox[{"p0", "[",
RowBox[{"[", "1", "]"}], "]"}], ",",
RowBox[{"p0", "[",
RowBox[{"[", "2", "]"}], "]"}], ",",
RowBox[{"q", "[",
RowBox[{"[", "1", "]"}], "]"}], ",",
RowBox[{"q", "[",
RowBox[{"[", "2", "]"}], "]"}]}], "]"}]]],
SqrtBox["k0"]}]}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"p", "=", "p0"}], ";", "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", "=", "1"}], ",",
RowBox[{"j", "\[LessSlantEqual]", "60"}], ",",
RowBox[{"j", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"p", "=",
RowBox[{"p", "-",
RowBox[{"dt", " ",
RowBox[{"Apply", "[",
RowBox[{"dU", ",", "q"}], "]"}]}]}]}], ";",
"\[IndentingNewLine]",
RowBox[{"q", "=",
RowBox[{"q", "+",
RowBox[{"dt", " ",
RowBox[{"If", "[",
RowBox[{"vanilla", ",", " ", "p", ",", " ",
RowBox[{"LinearSolve", "[",
RowBox[{
RowBox[{"Apply", "[",
RowBox[{"ddU", ",", "q"}], "]"}], ",", "p"}], "]"}]}],
"]"}]}]}]}]}]}], "]"}], ";", "\[IndentingNewLine]",
RowBox[{
StyleBox["If",
FontSlant->"Italic"],
StyleBox["[",
FontSlant->"Italic"],
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Length", "[",
StyleBox["qq",
FontSlant->"Italic"],
StyleBox["]",
FontSlant->"Italic"]}],
StyleBox["\[NotEqual]",
FontSlant->"Italic"], "0"}], " ", "&&",
RowBox[{
RowBox[{"Norm", "[",
RowBox[{"qq", "-", "q"}], "]"}], "<", "dist"}]}], ",",
RowBox[{
RowBox[{"dist", "=",
RowBox[{"Norm", "[",
RowBox[{"qq", "-", "q"}], "]"}]}], ";",
RowBox[{"p00", "=", "p0"}]}]}],
StyleBox["]",
FontSlant->"Italic"]}],
StyleBox[";",
FontSlant->"Italic"], "\[IndentingNewLine]",
RowBox[{"QS", "=",
RowBox[{"Append", "[",
RowBox[{"QS", ",", "q"}], "]"}]}]}]}], "]"}], ";",
"\[IndentingNewLine]",
RowBox[{"{",
RowBox[{"QS", ",", "p0", ",", "p00"}], "}"}]}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"traj1", "[",
RowBox[{"U_", ",", "dU_", ",", "ddU_", ",", "p0_", ",", "q0_"}], "]"}], ":=",
RowBox[{"Module", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"ND", "=",
RowBox[{"NormalDistribution", "[",
RowBox[{"0", ",", "1"}], "]"}]}], ",", "p", ",", "q", ",",
RowBox[{"dt", "=", "0.01"}], ",", "\[IndentingNewLine]",
RowBox[{"QS", "=",
RowBox[{"List", "[", "]"}]}]}], "}"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"q", "=", "q0"}], ";",
RowBox[{"p", "=", "p0"}], ";", "\[IndentingNewLine]",
RowBox[{"For", "[",
RowBox[{
RowBox[{"j", "=", "1"}], ",",
RowBox[{"j", "\[LessSlantEqual]", "60"}], ",",
RowBox[{"j", "++"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"p", "=",
RowBox[{"p", "-",
RowBox[{"dt", " ",
RowBox[{"Apply", "[",
RowBox[{"dU", ",", "q"}], "]"}]}]}]}], ";",
"\[IndentingNewLine]",
RowBox[{"q", "=",
RowBox[{"q", "+",
RowBox[{"dt", " ",
RowBox[{"If", "[",
RowBox[{"vanilla", ",", " ", "p", ",", " ",
RowBox[{"LinearSolve", "[",
RowBox[{
RowBox[{"Apply", "[",
RowBox[{"ddU", ",", "q"}], "]"}], ",", "p"}], "]"}]}],
"]"}]}]}]}], ";",
RowBox[{"QS", "=",
RowBox[{"Append", "[",
RowBox[{"QS", ",", "q"}], "]"}]}]}]}], "]"}], ";",
"\[IndentingNewLine]", "QS"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"t", "=",
RowBox[{"traj", "[",
RowBox[{"U", ",", "dU", ",", "ddU", ",", "q0", ",", "1", ",",
RowBox[{"{", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS1", "=",
RowBox[{"t", "[",
RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"p0", "=",
RowBox[{"t", "[",
RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"t", "=",
RowBox[{"traj", "[",
RowBox[{"U", ",", "dU", ",", "ddU", ",",
RowBox[{"QS1", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}], ",",
RowBox[{"H", "-",
RowBox[{"Apply", "[",
RowBox[{"U", ",",
RowBox[{"QS1", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}], ",", "q0"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"q1", "=",
RowBox[{"QS1", "[",
RowBox[{"[",
RowBox[{"-", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS2", "=",
RowBox[{"t", "[",
RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"p00", "=",
RowBox[{"t", "[",
RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS3", "=",
RowBox[{"traj1", "[",
RowBox[{"U", ",", "dU", ",", "ddU", ",", "p0", ",", "q0"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS4", "=",
RowBox[{"traj1", "[",
RowBox[{"U", ",", "dU", ",", "ddU", ",", "p00", ",", "q1"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS1", "=",
RowBox[{"Append", "[",
RowBox[{"QS1", ",",
RowBox[{"QS1", "[",
RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS2", "=",
RowBox[{"Append", "[",
RowBox[{"QS2", ",",
RowBox[{"QS2", "[",
RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"Print", "[",
"\"\<Equal prob. tranfer of two points of same total energy\>\"",
"]"}], "\[IndentingNewLine]",
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"{",
RowBox[{"QS3", ",", "QS4", ",",
RowBox[{"{", "q0", "}"}], ",",
RowBox[{"{", "q1", "}"}], ",", "QS1", ",", "QS2"}], "}"}], ",",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<\!\(\*SubscriptBox[\(x\), \(1\)]\)\>\"", ",",
"\"\<\!\(\*SubscriptBox[\(x\), \(2\)]\)\>\""}], "}"}]}], ",",
RowBox[{"PlotMarkers", "\[Rule]",
RowBox[{"{",
RowBox[{
"None", ",", "None", ",", "\"\<OpenMarkers\>\"", ",",
"\"\<OpenMarkers\>\"", ",", " ", "None", ",", "None", " ", ",",
"None"}], "}"}]}], ",",
RowBox[{"PlotLegends", "\[Rule]",
RowBox[{"{",
RowBox[{
"\"\<\!\(\*SubscriptBox[\(q\), \(0\)]\)\[Rule]\!\(\*SubscriptBox[\(q\), \
\(1\)]\)\>\"", ",",
"\"\<\!\(\*SubscriptBox[\(q\), \(1\)]\)\[Rule]\!\(\*SubscriptBox[\(q\), \
\(0\)]\)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(q\), \(0\)]\)\>\"", ",",
"\"\<\!\(\*SubscriptBox[\(q\), \(1\)]\)\>\"", ",",
"\"\<\!\(\*SubscriptBox[\(q\), \(0\)]\)\[Rule]\!\(\*SubscriptBox[\(q\), \
\(1\)]\) terminals\>\"", ",",
"\"\<\!\(\*SubscriptBox[\(q\), \(1\)]\)\[Rule]\!\(\*SubscriptBox[\(q\), \
\(0\)]\) terminals\>\""}], "}"}]}], ",",
RowBox[{"Joined", "\[Rule]",
RowBox[{"{",
RowBox[{
"True", ",", "True", ",", "False", ",", "False", ",", "True", ",",
"True"}], "}"}]}], ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"{",
RowBox[{
"Dashed", ",", "Dotted", ",", "Orange", ",", "Blue", ",", "Blue", ",",
"Orange"}], "}"}]}]}], "]"}]}], "Input",
CellChangeTimes->{{3.811366354088255*^9, 3.811366354679247*^9}, {
3.811366395463512*^9, 3.811366397232502*^9}, 3.811462273027076*^9,
3.8114626387790613`*^9, 3.811465409373406*^9, {3.8114656604274197`*^9,
3.811465661321134*^9}, {3.8116986162435017`*^9, 3.811698737631402*^9}, {
3.811698811167563*^9, 3.81169888045022*^9}, {3.811698991040303*^9,
3.8116990928880672`*^9}, {3.8116991918428307`*^9, 3.8116992043113937`*^9},
3.8116992344116087`*^9, {3.8116994735485497`*^9, 3.811699515087852*^9},
3.811699550796742*^9, {3.81169980843515*^9, 3.811699869144209*^9}, {
3.8116999241201677`*^9, 3.811699980087982*^9}, {3.811700050276331*^9,
3.811700141098805*^9}, {3.811700182556608*^9, 3.8117001829596663`*^9}, {
3.811700225533052*^9, 3.811700248865831*^9}, {3.8117003847902927`*^9,
3.811700407735458*^9}, {3.811700629632187*^9, 3.811700669116385*^9}, {
3.811700706768877*^9, 3.811700769079571*^9}, 3.8117008029579144`*^9, {
3.811700914933988*^9, 3.811700941767864*^9}, {3.811700976956438*^9,
3.81170099387178*^9}, {3.8117010644826107`*^9, 3.811701103255178*^9}, {
3.811701146933124*^9, 3.811701271439931*^9}, {3.811701337780672*^9,
3.8117013649371347`*^9}, {3.811701496531475*^9, 3.811701590960135*^9}, {
3.811701627803259*^9, 3.811701668544454*^9}, 3.811701705778117*^9,
3.8117017424690866`*^9, {3.811701773025908*^9, 3.811701800210608*^9}, {
3.811701986563479*^9, 3.8117020522365713`*^9}, {3.8117020834261217`*^9,
3.811702141108369*^9}, {3.811702813499928*^9, 3.811702827761104*^9}, {
3.811702891929553*^9, 3.811702899785233*^9}, {3.811702968200354*^9,
3.811702988610875*^9}, {3.8117031262041187`*^9, 3.811703140339157*^9}, {
3.8117031986380253`*^9, 3.811703222481554*^9}, {3.8117072695087347`*^9,
3.8117072974269447`*^9}, {3.811707330550557*^9, 3.8117075975225363`*^9}, {
3.811707673747344*^9, 3.811707823821927*^9}, {3.811707880126397*^9,
3.811707921478196*^9}, {3.811708024226787*^9, 3.8117080338986273`*^9}, {
3.8117080657162857`*^9, 3.811708132100494*^9}, {3.811755690000998*^9,
3.811755698584117*^9}, {3.811756300217101*^9, 3.811756324200494*^9}, {
3.811756384672604*^9, 3.8117564394089527`*^9}, {3.811756471520866*^9,
3.811756691649251*^9}, {3.81175672313801*^9, 3.811756723960552*^9}, {
3.8117567560260687`*^9, 3.811756783988184*^9}, {3.811756822059173*^9,
3.811756843882867*^9}, {3.811756916266765*^9, 3.811756940377589*^9}, {
3.811757010403778*^9, 3.811757033265319*^9}, {3.811757076028442*^9,
3.811757085275816*^9}, {3.81175713961887*^9, 3.811757153920624*^9}, {
3.811757195179111*^9, 3.8117572082333727`*^9}, {3.811757277577427*^9,
3.811757314512705*^9}, {3.811757522800967*^9, 3.811757538145836*^9}, {
3.81175757377864*^9, 3.811757612034491*^9}, {3.8117577337939453`*^9,
3.811757734641171*^9}, {3.811757780067772*^9, 3.8117579336941853`*^9}, {
3.811758310894926*^9, 3.811758318108982*^9}, {3.811758351764329*^9,
3.811758386913609*^9}, {3.811758424895781*^9, 3.8117584251099463`*^9}, {
3.811758457026251*^9, 3.811758633525455*^9}, {3.811758687860889*^9,
3.8117587230272293`*^9}, {3.811758753038617*^9, 3.811758772389324*^9}, {
3.811758853315568*^9, 3.811758936804216*^9}, 3.811758967625676*^9, {
3.811759026654052*^9, 3.811759027762918*^9}, {3.81175907832528*^9,
3.8117590784432287`*^9}, {3.811759111324388*^9, 3.8117592068813066`*^9}, {
3.811759238130941*^9, 3.811759269247634*^9}, {3.811759392011499*^9,
3.811759393082686*^9}, {3.811759598420388*^9, 3.811759601615836*^9}, {
3.8117708787443943`*^9, 3.8117708793308153`*^9}, {3.8119609091277323`*^9,
3.811960966711893*^9}, {3.811961010487213*^9, 3.811961160566895*^9}, {
3.811961236561183*^9, 3.811961237997033*^9}, {3.815280786376183*^9,
3.815280788167563*^9}, {3.81528091432285*^9, 3.815280926079616*^9}, {
3.8176302164980087`*^9, 3.817630217036269*^9}, {3.817630258056705*^9,
3.817630266706142*^9}, {3.817630371434236*^9, 3.817630391730562*^9}, {
3.8176305479534607`*^9, 3.8176306103306637`*^9}, {3.8176307438880453`*^9,
3.817630756281435*^9}},
CellLabel->
"(\:8c03\:8bd5) \
In[197]:=",ExpressionUUID->"ddce6090-e086-4b54-8a86-dd09ab449256"],
Cell[BoxData["\<\"Equal prob. tranfer of two points of same total \
energy\"\>"], "Print",
CellChangeTimes->{{3.8176307587307158`*^9, 3.8176307718300343`*^9}},
CellLabel->
"(\:8c03\:8bd5) \
\:6b63\:5728\:8ba1\:7b97In[197]:=",ExpressionUUID->"06b90ded-debf-4953-9c47-\
d038b44fd180"],
Cell[BoxData[
TemplateBox[{GraphicsBox[{{}, {{{}, {}, {
Hue[0.67, 0.6, 0.6],
Directive[
PointSize[0.004583333333333334],
RGBColor[0.368417, 0.506779, 0.709798],
AbsoluteThickness[1.6],
Dashing[{Small, Small}]],
LineBox[CompressedData["
1:eJwVxXs4lXcAB3Az5PKUTUSeMc42kcllIk/xhlEpax3m7uQxbKeSYaiI8Gg2
RiuXPMpio2Rz6ZjLcfkewuN6jpi7OOc976+5JfYUPZNztv3xeT7GoVHscEUF
BYUj//n/jIjBvyc5UpjcD53Zs05wr27DUrBEozz9QLX7GkFBlEnnrXga23f3
qym/IFg7ERESqEjji/hN3snnBNnmmq2GORKsu+Vd1V4mkJxaTSF7JUjaZevo
v0CQyLpwurFcjOqOd0d3EwL+sfruYmsxgr5yNXOSENyQty7aJs5DzTTbdmSW
QH/Iuj/u1hzcWRav+iYIysYr9qn99hSO6/UJ748QqAdPke2uWbzpHxobHCDY
2yQQhc7NwMVA3X6wi8DPbjXf/fU0qnNrnui2ERgHrA02ak0jhn1tubKeIKS8
ebnLYgoT99Kmk6oIDD7S8Eg6MQmDsz3KWaUEEc9pw5XwCVxpMX0mLCAQ/LGa
zEofhyVfZ4drFkFDt2KgRekY5HnHdtDJBHUh+cn6gj9xfTy9qCKaoHjeKWdz
fhQH9X0ysr4k8KhwUhlRGMWugN77Wd4E+1jKDo2sEVzh6b0s/ZSgS9bpWe/2
BNLsC9/2fUJwp3XbaZw7DG2VVI9tYwLpQ69HXoYi1Nz1qz2kSXDEJj8m4fwQ
zm6Q1/FbDFamLAutWgYwVMG7wXvGwOeAd17hzn78Ls+ULQ4zeHjVN6ErtBdl
PA1Gh8/A5KcPGSG/B7qeHJlNGQNhZLOWUKcbGZwWE8fvGXzOfTQ4GfsYmgE7
vayiGNjM0eeUxzrQvp1WoOrNoKF3a4TrIICmixGr054B1zR1yc6lDbKmk27+
+gxetOXlcfT4eOeU38X+LSlq13mlsa8awD2uYa41K4Vde7RuyRQPgQ8MzSxa
pIhaYcuVumvh9wPs9Iqk8HyruPtOVBVKulVdh+OkuH1QlOMUVo6UHy9z2Gek
OKwwVXVxuAS52l+3FJtLsSBflxdt3kSK8TfVlUpSrNg2lOn/kwTZd0+rL8/S
OK4lOm3EjaaCLmn5qvFoFG55ZVX8mksdLVoM88mkkT3L5hiZFVPO87v5nCAa
/M6lkvcsf6Gq9sQEGVnSaKprMyxoekCd0T4XX6RAQxbFMeixqqYep9n79g5L
cNs8xLR1oI4SOh+KrfpZgmCD2GDH+Hqq9TPbMOdICZwiHSaEuo2Uy8eli5kO
EljHXVr0r22majrEC9eUJEjQ5Wr3Ua1UeGVipKlQjI7imyqqPe3U4b98KpML
xFBnHxWllgqoN3LmZVqwGKFv79/wnumgBOfZItsPxPgXxKkhlQ==
"]]}, {
Hue[0.9060679774997897, 0.6, 0.6],
Directive[
PointSize[0.004583333333333334],
RGBColor[0.880722, 0.611041, 0.142051],
AbsoluteThickness[1.6],
Dashing[{0, Small}]],
LineBox[CompressedData["
1:eJwB0QMu/CFib1JlAgAAADwAAAACAAAA72gJai+iwD99COdrsl3gv/LI2EW+
yL4/YfTRkd2W4L/UPqwF2kO8P5deXry2zuC/IGkzCcW1uT+4v5KNPgXhv9rb
KHaTHrc/WBIFznU64b826KqgWn60PyRmqmxdbuG/hBBmCTHVsT/P06F+9qDh
v46PNbdcRq4/e9j6PkLS4b++b/3V1tCoP/sfdw5CAuK/UpvzYgRKoz+6xUdz
9zDiv1I7/TI3ZJs/rxTGGGRe4r9SBg3vqRKQP1vPJ8+JiuK/7inHwFl/cj88
CS+LarXiv0pcyZElz3u/ypvVZQjf4r9qaoVGt6eSv4RB9JtlB+O/FU/OlnB7
nr8lYuWNhC7jv8i5CjY3N6W/l5sjv2dU47/sD0bQEkCrv88T5NUReeO/KC7u
twGssL8hoKyahZzjvyKVOa1fv7O/Qs/l98W+478ye4Vm/dm2v4rjaPnV3+O/
ytjbAbT7ub+kywnMuP/jv9EhzKBbJL2/PigdvXEe5L/HBkO25SnAv9Rt+jkE
POS/DHQJzezEwb8nMnrPc1jkv0l2VrctY8O/ZLVwKcRz5L/jERqikgTFv3i3
JBL5jeS/lx/gRgWpxr95qsJxFqfkv21yJ+5uUMi/jFPMTSC/5L+q7sZxuPrJ
vw3shMga1uS/h0hgP8qny78z1VkgCuzkv+cf4FqMV82/3PBHr/IA5b+KLAth
5gnPv4GyPerYFOW/XxYMxV9f0L/K+3lgwSflv3QjK1b/OtG/l9jnurA55b9u
xm8fxRfSv64ud7urSuW/MacrLqT10r+rdXI8t1rlv3Pz2WKP1NO/F4zRL9hp
5b912IlyebTUv/W+iZ4TeOW/QGRO6FSV1b9pGdunboXlv1KisyYUd9a/YBKb
gO6R5b8ZxzhpqVnXv5CufHKYneW/5DzPxQY92L9SLVbbcajlv3dkXi4eIdm/
R1dkLICy5b+92kty4QXav+eFi+nIu+W/pBQIQELr2r9qe5aoUcTlv5AhnyYy
0du/yiJzECDM5b9fZE2Xorfcv8RPbdg50+W/dxIY54Se3b8Rl2fHpNnlv+VH
aVDKhd6/SFcSs2bf5b8df6/0Y23fvwkLIX+F5OW/RB0Ab6Eq4L9J/X0cB+nl
v9bWXgCsnuC/0Xh8iPHs5b+lmZ8eyhLhvyiMCcxK8OW/zS6+MPSG4b9Oe9v6
GPPlv6Y2XZUi++G/2vifMmL15b/s/Z2jTW/ivypAKZos9+W/plL5q23j4r+F
Kppgfvjlv7k+Gfl6V+O/IFmRvF355b8ZjrPQbcvjvxqOU+vQ+eW/efsfog==
"]]}, {}, {}, {
Hue[0.6142719099991583, 0.6, 0.6],
Directive[
PointSize[0.004583333333333334],
AbsoluteThickness[1.6],
RGBColor[0, 0, 1]],
LineBox[CompressedData["
1:eJxdm3c81u/3x4tSqTRIIiENys7KukQiq4yUQmZJGeEThWSkba9I9iZ7j3O7
7ey99973RGX83t9/f3+9H+dxv8d1rut1znk9/7i5TG21LWj27NmzQLtnz/+u
+1sEq0+sVaHFa2bf6/+Nwc6d+Suje/Ho2fyHEPGcMZBKyd4qYcQj+Z0dnaIn
YyDwRJc+8QIeaX71CVE9OwbGdAw80aJ4pFqTwUrsGIVK0w98KYp45LmnNT79
wyiosg8SQAuPrDpPEt/IjsLWPv+7s0Z4dOWEDf4xeQSOOsbzsj3HI46FiL6H
qSOwlexuYPQKj7qWWAkvHo/AGca9zVnv8Og9p2hT6KkRUEn5z/DYZzwSObV2
eaBpGN5zRG26BuKRHYGpVtRzGCzEUj9sfsejOi9+kzTJYSielN3wiMGjLCO3
Uqm1Idjn91L8TBIemUzS+s0mDsFTmRbxyjQ8UtJ7/eOXwRDEZqgv2/3CI9pt
9ZggxiEIemekIZyDRy5X4u+H/h6ExLvT2ru5eKT9UMe9yGMQAkFkczAPj9yU
Z7s2JAdh8f0J6RosNt0XfOIBYQB+OJexlWH39zIqUPqTB+DNh9EfkI3l7+Uh
4fJ4AASLIqs7MvEoaN+3UHR6AK4zfggjpuIRz2Hl5ktt/VAfkXyeIxGPNmYl
goQ/9IPjJaeXD6PxKG3yS6Ix6oc7KpOB8eF4NGhyoiJ/ow8YL6qE//XHnt8j
ESSQ3Qc3lN+EPP6IRxk5hyZaLfsg53l5TOdbPKKZFnkcytUHMio1nTqOeFQW
s177frAXotVKbk1Z4pFTtnBLVFAvaNwXZfUwwKPLcJ9vTL0XRu2RNf8dPIr1
eBmhStcL2873TOfl8WhpiLV1HHrgTaXJuVxhPMo/6OEW97oHgnKh4TMXHknh
jY38r/UAz2upL/bH8Wjy+b5L6SvdUP72SOCznSokJ6vwiZzcDX/3jx22W6pC
IjZvZKxMu+HXr0Zen74qpGgeQj7G3g2CG1ocGfgqxDGx13SmrwsGC63OT2ZU
IWORPKmlwC7wCHR4zhNahdJMOYTOa3ZBjZWkmPvbKpTzXJby+VAXfKgmlM9a
VKHtZ16Sl2o7YeciUnisXoXaOqYzKe86AVyEDiwIV6Fzie37NmQ64aK59F1v
5iqkki64wv+nAwo8kbnwXxxaTTVljijogOcN3l8Jwzh0PFpeUMa+A94VlDLh
K3FoeotplkmwA9ZDNlSTonGoYzZq/sJyO3wuKnKPcsehHlv66Rep7aD0fXh/
ihEOBQ5ueq88aYdbpzJE6mRwKHXQ0Sn5QjtserUZbp7BIYGspRbx2jYozU3E
uy0D4jQt3Qq60QaazsrXROMBTR5aaHta2Qp2fFxNTA8ASaz0ReRJt8IpeT91
zsOAlteys5xLWqC+If3t/fJKxMBj/rhYogWuuKvzlz2vRMT/Ev44FDaD8K89
vOpnKtEMx3RzhlgzJLMcUz5cW4ES7LQumRc0wUa7g8kfmwr0Qb9LK1KsCZL3
jeqdPl2Buk/TRegU/oY3FZn0TyrKEXVYQsVX4jfkBhHNZ03KEb+LRYJ6SSM4
6Zlrf99XjrSYJKh+0o1Ql2xV455YhoRbttz1Kxtg6PPxnIibZei9V/rHtBsN
IDZPx7g0UYoOfhzQ8q6tB639sQsv3UrRh0dlZ+dv18OttQ5JodOlyMSOk7ev
tQ6sB3kPXc4qQdWrtgXaunWAf9X19L5SCQp0Tdx5NFgLT8UDnuAHitH4IU15
inEtBH8qZXv2ohhZxXjnn52vgb3HrsepbxehZ+JrQb22NfCm6P0Zm69FSNqn
7BT3ZjXsGU1MbD1ThDQWIm1pPaohXZfH1jaxEMmI1/Q401cDxfx8vJZgIRJI
mH3jHYwHzVo+h1dFBQiXLeR/mQMPax7Gp8ZkC1AdvYqSWVoVxA/q1HzF5yPc
O+UGSfEqEDt1Gt4p5aMWUplSKh4HOorCt8vr8lDM7dj5gjs4kN/5kSZ/Kw/9
Gq/Mh1KAwlO1gsdqcpFRsu/7u32VoCbqwcsrn4vkGjj8OakVIOB0csWvNAdt
XE9iF2GqAKMtju6b13KQwTlmxw+i5WBslidzOy0bXZ5g7z2vVwZzimKffnJk
owAjTi+616WwLtIl//VSFuo5KVgkGlUCBjf7Bz3pf6GXWvvSs/DF8PepPnAv
ZSDhD9tBrxeKwPBIX8/LhnT0xeNUQvCJIqgasWKwi0tDkvjRQ7TShVDLhzPh
ep2KaAbcCHUWBaCIjnR8VE9Byg3x7pMB+ZD/ke5J6tlkNG73ceg+5EFK27iI
z0Iiunpw5gLvai50Xkfa5/MS0FjoqvODc7mQryQ09+51PDr77i5x7k4O6GVc
YE+UiUN0jJcTBjyz4SH/3OXArRi0cZui/NgmC1a+8XwItY9G809td/sGMiEt
k/tZ8swP1M1j+J5LJQOOqx/1M9SOQB80Ti6gkjRog/g4KA5D4/z+IhL8qbA0
t8+omyUEVam+MT2ckAw7VT++xNkHIknyYXscexI45F0+YvrVDz3cOqrzMCIB
rNlc7I5e/ozEvzWvDZ+Jh3XtI1L3CrzRytaRq+o/YmFaz2vb1+Ytqq2lHU/n
ioFO9575aX0H5HMqpGZXJwoIjTb4AH1jNGac8S/MKgLslt863ahVARzz4Mrl
D2Fw9sJnmbjlF8CJpJiKUoNh1LrGm2n5NcSRzPvVOwPgcxOd0WK5J/xVLOSZ
2OsLHWL9ayV+H2Gv2NwFnecf4Q9XxewzG194f9ufdWfcE+h27T79xxcILEpf
Ksdr34AsB+Fk2OUQKLZ8FD9eYwuT2/YDATzhIMmq/ezsY124sWRz2kgwEsZt
Zdb+PdVD+RFqm2TpnzCdb6tUJ2KHnLOlmy8mx0Ac49J9h9w3iNJW4SlwLQ6a
4hq8//vqiWzSD3Icx8fDjjlV8lDeB/R5kfV7nU4iyFzcdZxr/opuNEdPaMwn
Qevji157V/xRccqfiVj3FHA1jLG+bB2Eyj3VXSvPpEGzf9HesdkQdCn745fo
wnSoce7p43oYjmYD/5uS18sEsQSW7wt1EUjdvPzuz81fcJfZfFeQPwq5fXrd
PKCZDbletj2rX6ORdISObOV5TEdZz/8dI8cgI698vb8bORBgzJxidDYOdSss
Eb605sK6r6etsUI8slTpGrdLygPfxX36Jy0SUPnJwY0k93zorvNyfOmdiDa7
9u3je1gAzgVvB9/EJKFijvW2bbFCiKdnzrhakoys48jnWRmLgP32SUaP1hQU
N7aEe0ssAlWXpRvO46lI/J+LFW9HMdC9bHdgWEtDnz5Rh0/nlsD2Yuay/N90
lGFp3qkcXAq5mWJ/jtBkouXY4L1FTmXQ+XZm1ZbuFzLVZTlnZlAOzIcY+J4e
yELNSRr9txUqwKx3+bAzTTZqxJuPWPJWwnwk7wlv92w01mvRV3ECQFRR3LT3
bzYqOy9hfesfwHcfGpy1Yw6q8tbi2uOAg7uXNtdUlnPQFLOgK2UeB6Ondsqe
mOYiUiMds/vjKrgZ7UKq7c1Fy08ulf/srYJ850dPLG7nocOnH0ppa+LBtJ1h
TbE0Dx3xWHsYXoeHNlV/S2PefBS+WDxri6rh08aTvJLQfCRueDyxr7gajkqm
5mvSFCCTtsFHrSI1UB+5e+/ciwJ0RutLm3ZmDVhlJ7+/3F2A+GUuVZvx1EJe
3SvxZ1KFKMry5uZufC3YpQ4/nvxZiHjNmQUvcdYBu/rsdsDeItSZ+fRSz486
OPdM/oSTaRH6Q4mNPMxaDx2npaL9q4pQ3aN21YaweljWrc8ZOVeM/CriVg8x
N0DpW6/7hm+KkcmmnXxbcAMcqcuLZ+guRgJScpvMTI3AGXj7F+FqCVLRKm8e
C2qE5YHpIBrPEuS/544DD9NvMLft/U+xtwRF6x2LWQv+DXfkz7/J5SlFKQfb
9osyN8Hlzx5Nmq+xOdfOpkcJawLNvRbfTzeUorYWhRuCrM0AK+unj5wqQ18W
dV5N/WgG5+vs764al6E6hYvxxzhb4ICjy7H/UsvQmJ+eY258C5xpdDw4TyhD
OyG6Ib8vt0LVEn+ul3g5aggxyNbLaIXfZN8Hym/KkSrO3eaecBtoCympi5SX
I8VfSU61hW0gsFIzqbiFzXFr1RLtc+1gyKf2ykWqApkce0XzxqQdDmiPWfS/
qkCyrOELfQntUFyne+phTgWKYLPAPZ9vB7dP6iu7CxXI5hRdvyhfB+CLxR81
clYiLgVWDxG7Drh5VLej4F4lUh8sPmye3wGPVYrqaj5Wop0h+/qGzQ4YCb1V
RC2pRFudurTGsp1wf7acW2WhEpXGzR/l9eyEIYH0gDJmQC+0JMS46zshIev1
lzsKgIK6f9dpHOmCz+eexdC+AOQ+mnUsRasLJFqCzPuCAPllOhkIh2E+snLV
7ncJILHh63/mh7sgKE/VbHgEULBIqHDr+W7gC7piS7cHhw6/Lrw7btkNx6Kb
a/iP4xBSvxt7LqsbfujjDppz4NAGMcn6I7UbvumkMaby45Cz4LfdszI9QCsZ
jduSwqFazqcxQ549IMV3YcVIGYcYTvB9r27sgT4ewddt2jikq2sp0n28Fx6q
X1bWNMShLS6Uf+RBL0xenpEeeoJDLnYvnryI7oWi/Xlyr2xxSKjzlj9pthfm
eTSuczjh0B/tKuefAn3AcWDydK8bDindt35i+6oPnricbIzwwiEnF9do88o+
6CATZa0/4tCUI/6ZG10/jIbRWWl8xSF7ciF9iWY/8FreVpT2wyEeb6d65rB+
ENe/hBcPwKFCx8bpkLF+8HriOoECMV9bxJYsxjMA12KifuphcdLvj9pUuwHQ
EL7z5w12/++I/ivdJQMgYNJIzsDep2nt6dBBMwjK34oCl7HvBdf9sl1WG4TN
s5Ij1z/hUK72ddMLIYOQYDs8HOyN5dOa+815dBBuV6dEbb/FodMUU9H5y0NA
x8wk4OiMQzhj1XD7l0MguOdg9KYdDr19f+Efc9kQzH8cp/1iiUOd/s9+du8b
hhLyc/urj3GIy+BudZrmMAR9XKQZ1MUhG/32tPDwYcjXZWgPvY1D/3nEh0RN
DkPmzN91E1kcqldQqSnmGwHC9s0waSEc6rejdVp4NQIvFEmd3OdxKEJ0ckWg
agSOzzm3sjLiUMxWmeXHw6PQfAuyOWlxKNRG9DL13ijEnZC49XwaEO5BhaVj
zCjomCTF5eEASSWEv9i/NApPHol1sUYC6lAefJksNgYMDbzXYx0ANUlZFem/
GwOcQaKUsiqg35Z0rmxNYzCbl6BwlAPQgPb5XSr9ODhJ7uDIxEp0UEtVzkJo
HK7tM7nyD1+JVtZiWIfujcNmjs4md2AlGuuSunbfZRzEuK+l2jyuRCruZy8P
xIzD3evCuKErleidhnSIcd04vPLQqX1OqUC0Nm/PriyNw1jPSyWO8grkMAaP
3p6YgBsE0c51jwok97yFyCgxASFEJhqyUgWqiXnmlmkwAcHD9+2YDlYg1g71
KFXPCeDed7ziUUM5wtFfIiwmT8Bha9ZvjT7lSP3UDybflgl4wM3+2UixHD22
1ky9Rp6AJq/3pmw7ZchDvld9iGUSDjv7NO4pKkNyz+qjvOUmITSZ2emYTRk6
+TxHStB8EhQjVy+rcJehGBtO/OCnSbDq8YlK7y1FONmfnT5Zk/B3+NRPqY+l
CFAe/bWeSdgj60uhSpQi8//mDo79nQRBpx92ozMlSO2/6uefOaeAdSlhlhJQ
gpafdnSI3pqCzawqQSmZEjRA+3Ng5PkU9JrheLOmi1HzpxKu9wFTELrMUaT7
pRj144pVrhRNgbYgA1wVwuYDSWyleXgKSMUW/GJdRejUbEfmC5ppGDRIJ710
LEJLI6K8h3gwYdQd3jfFWISoAcd2YjWm4ef2g0efcrB5xXe2SdxhGpimDmxY
aBSidUUOwfrwabgdyT/4eq4AUQbmG7Qrp+E70fVInXsBQjr3JAampiE0UzFC
jbkAdTpqnnl4aAY4Rp4GH07PR3ij5qM9AjNQc8Li+AG5fMT/oTfvtu4M+Afd
5kFteUiuwDut6PUMbPnk/c03ykOt/3b9z0Vjzx9oTXu2nIt28rzo3tbMwHdR
TWMj51zk7WSU1rMwA6uSdsrBNLlokZXMdvHYLETLUPyPfMlBnsklf5+LzsKE
Dpt9w4kcdPye/Gyq/iyQDRMuNIZmo2defG9H384Cx0DD7DHWbFSeJ216KGEW
jB44rEtLZKH3HjdFrjTOglWthPaC+C/0IVLCV351Fo5X0jtfuZaJmMo5mdUZ
56B3D/kADV8G8o+9bKUuOQfZLk1p9lzpSIEu4MoNwzlAz2ROujOmIfpw2OD1
nAPWh/mbvDSp6PeaqD1d8hx4cZtwvVpJRvYfA6V6m+ZAsILnwZOeJFRSEEEM
J8zBlb1xL7dLEtH3CGFJzVPzUGZwVE0uMgHd5DDKol6fh3Bzjd9XMQ7ZmFXd
9DWahya9sZEGnTikoaScf8ZrHi55d7myXI1FQtcavEOS58FG/YuG9E40clN8
dpimeR6MxdPqC6OjkHLYyoIRYR5m2KQ+pUpFou8d8vh0pgWQZ9KvOdMWjrrp
ncXnJRcgbXmr7KhhKBqo7pxnNFyAdOpk9ZfpIHSvv8xFwGMB7ia5y0dYBCDO
l78DJBIXoO38sIY77htqpOJbBBoXIOlR5Z99qR9R9XZYJ+PKAkhunMjp9PZC
KjeeK84eXwTTtnueS8mu6P0n3+pE0UXQtFAS2uCzR0Yx3utaDxbBWSTzoq+Z
ERpENO8WXRbBf4fmZ0ysCjTuWaO1jV4Eg8Ezxy3zn8O5208ujeEXocGm32FP
iTPEXy/+IjW7COLr6ZFy/R4gF2nV/+7QErgw2P25ZfwBYpqu5+TwLUHSizDa
L5SvIJ96pqz5zhJ0KltffKsaAFzu4pEd9kugYS/JWccZDElr1HUIWYIiSTY2
5u1QYB7bcAgrXgKFacusJ2Pf4ZHh6eL7Q0vwUmTlSXb9D1D3e/5sz84S6Efs
li0VRIN4mIpAAOcycJo3fT14LxZ2XeUzDykuw9TZX9ZHOuPg8akL355ZLMP8
3vGkba0EwM/T/pf9YRnuGtPd7+9OBKkzwkdGU5fBuE02JPxhMnCryP0lNi3D
Qe0LL6WmUiD/zZ65pZVluPJ+ilRikwYaO35vmo6tAOO0FBvjdjpUqN029Bde
ATXOz1M3fDOBY7FFQkpnBba+DskqnM+C8bv78xocV8BIw4n+o0E2/Gtgd5YO
XYGya1xJOrw5UKrxTSiwaAU+iRwf98A45fuLkz7N/StQ874r7nh9LmheKDi/
+GcFNu5a9WyH5cGj6vrsRdZV4Dp6+5GyVT5k1s4NtkivwvqOEdOSbAH45Oyo
BRmsQlAh28jsyUL4SVbIknJbhaYwR1/JhUJwKZXMro5aBY0Xi9sjuCJg68j9
y1+JPb/tu939vRgSnLukXo+uwmgpvQGbYwnY3TjKnrizCsHGcdO/7pSCzBS3
Zs65NfDMlXn8ja8MnjPxu0TJrUGEuKV/Pn05pBkHaD0zWoOz28bXzy+Wg+fO
EXemt2tQZc3M0v27AhQSsnN+Rq1BCJfwXH1GJdAaTIYdqFiDk1UsRv/8AGL1
jMY1h9dgdN2VRVwKBwKVElKO/9bgvzNPDh7JxEFY9Ys06aME4N1Nv/eWswp8
Y2gJ6WcJsEfrZ8+X4CrYQ16tpecjQLPAv0fXDuEhStdq3wNpAshLkkod3uLB
b3LUK0iVAKbfL1erUfAwHz91tkKfALFsRYqlVtVQviab0m1JADppCfrKiWqw
6pzYGXAiwGW58sl7+jXQZk/60+xDgECdobB3HTVQ/gpvnxlCgKu5mSQZ1Vqg
fdqt4JJAgKMVXaWfqmuBh3FIUTyPAHF8zoUWsnWgSym5M1pFgLfXddLaiurg
uP97eYd2AvC3ZOrXidSD2qAXgTJKgMF8cx+VX/VAG8YnZbZCgMb+JYLulQbQ
LePdrPhHgL+vsCpOagDO1vv799MT4Y3JjTcHLjTCaZU6GUkWIljJH2RMj20E
LXLR6/uXiPC6L3Ski+M3PFaNDzMVJcKY5vugdz9/Q22viL2+AhGa3JlXstib
gN5PaUb6LhH8b0rHGkc1AY/oWvMhIyLgmkIsAtibwfOX8/7q50RYXXVdlf/Z
DHetjJ8/fU0Exm+7q9YcLWCT93Bo3YcIX+fmWE7HtkDecW2Rl8FEoKWFs9Lc
rVD786lGbywRdA+3p/YntoKcFfnQpSwiFJZetVnlaQPDRm8l43IiHK1JZ3bJ
aIPzBZMj3o1EuMLpkizH0g49XEdygnqJ4OeHG3/9EOOW+5JRX6eIcKPl95fV
H+2wfqvvnS2BCDW9qVuRY+1gNBMgJbNNBJr6QR+P8x1Q56aaSjpEAtfg224/
LTrg6iH2jEBmEih2aIuQUjqA0evx+bPcJOB7UbvhstwB++XilnwFSaDrMH1O
WqgTYiSy5xalSfC633JTwLETSp77rQuqkEBV3Ktbr7gThC57bBrqkkAubIQm
d6sTOB7StDoak+CvbFQnutEFKgtN9x1fkKDSxjxw930XnNjtfGXgTAItKV1/
wu8uqH90+wi/Nwn4VRlYTx7vBovEqD8zfiQ48KvQ1exeN/wsdOH8EEmCgDTq
vumIbri70mx6PJkEjp/+7ASOd8MT5+I4j1wSvLNyLrS91APHE+PxgxUkkGfO
ffX6RQ/cvyqSyNZIgs8dLz9n5fbAyTe3+G92k+CGyHGZE3964CYNjZzuGLY/
VMeeKNQLVR61pWqL2P29LDG6Pr1wU0P4wxUqCUSfHpgRbemF9u8XPxB2SXB0
YXxZlqkP3GrxP6PoybCHZ3jX5lEfBEc/yxU+RYbyew0W9XF9sCe7LCmDgwyC
swxPbi32AfvbSGOGK2QQty5QXhLuh4Vevlo9UTIIuz3WKnjdD/+d7Cl5L0eG
1z53G+Or+qH5gfSFSBUynCm2phYcGgC6WZXFUG0yME4gtiWtARjT4CU5G5Dh
uA7VXSFiANwrVc/feEKG4Fozy7LJAfipYGSyZksGMSn6M/euDgIN7YNQr9dk
KJ042cvgOAjmH2xTtj3JYElzb2qufBC01zncDb+SgXNUz3ts/xDM6vLQxIeQ
oTHizhZFcwi0a7+caf2JrfcKr+/l8CFgrg3IHE0mg/81T9dXE0MgblWX1J1N
Bp3trEMTV4aBJebvSlYJGabu45yeOA6DhtnNF3Z4Mlz95HKQrnIY9m7v0DM3
kSFD8sl25YERcJVtyovpIoO0QXWYvxbGKYV+d48Pk+HnxVCWN5EjIEZybjKf
JkMIIR73ZmYEzmb/Yfi5TIZfwq/aAwRHQZ+laX85hQxpFgYfK1+PApWyHFG5
RQaVvTY8e2pGoXYusSRxHwXmTL4S7jGMwYkYL42XRygQU2nOCQ/GwO7YIwUu
JgrAGfsNqfgx2IxS8S5go8CyyG5+w8oYVBhm7/JzU2BUhqn52tlxbL1ysZ+v
UMB0z5rtlxvjECc1e69JmAK+xfYzi0/GISpwmp4gSQGO3ZUmra/jYEmuyNxE
FBDi/4GHnHHQSEd8M7coEH0vyVK0bxz29rs7FmhQ4Au6YJi1NQ6HQkYcn+lS
gNuEHQmdnwBV3r1n9z6igNXgp8xC5QmYvTb8wNWEAtZ2OsoK1hNg1J/KOvSU
ApIf9dI7Ayfg0OLqXXYbCjw7b+b8tHgCgpjfLCk4UuC/XlmNPaMTsH3MfFDt
DQU0W2MaftBOQlXBMQbJdxTYrpL5JMM7CRdlrzgd8KGAD7mEf0xzEjYepNOV
fqFA8JUSXW/HSbA07sm8E0CB99ndP65GYFzyn5hRQygFXqrHJnTDJEymxB+6
8IMCLCiL3n1mEhK5FONMYynwR+6L9ZXDU7At68nqnUSBBEd43iuExUWsFh/T
KfDzVXe8p94UeNef+M82mwJ4Je18AdcpGFMslr1eQAGN21XSg7FTkFqTljtT
QoELaun93vVTkF3VUe9YSYExhx+8AitT8CLX5OUsngJ7lS+09Jychqd8Q/FS
9RQ4Fd1r4iI5DXv2HtZ92YTtf6q6zzmjaXB4SO/0uY0Ct6YuDlR6TcPfPjka
ny7sd8ZDfx6lToNa6sMlsz4KnL+d6E1tnQbRVAOeC0MUWOzOZf1CmYavTv5Z
daMUMLpF0T/LOgPKgb5OapMUmHoiMZGCMI6xlnHKnaEAq5jeHSGLGTD72pz4
b54CbuMSirmfZ0A7pPLvpWUKONR0PBLIngE7qzFr0TUK6JvxSCb0zEC57O0N
bhK23yZ3nRj/zYB1opkflUKBv+62xS6cs/Ana0ogdYMCm6EVEcNKs5C73tQq
95cCTv7vK8Sez8Js1gvLwi0KyJtewvn4z4ICzox6bJcCdAeZFNoKZkGt6Yet
2l4qMLr1TzIMzcJhRa0OS1oqMBQyKSjtmQNmp4OnLPdTgSnCY/flxTmIigsT
u32ACqR/hNJg1TlY+cMlcOQQFfYEaTJl2s7BXCHL32x6KsSSY9JLg+eg4YJ0
iMQRKiR/oOeqKJmDm+mh/6KPUkHD5zdT/ugcpPX78C4yUIG8JHs5mnYeVJ8X
sDEfp8KPl+ObbjzzcO0haxv3CSrseGlcuasxD2vD/+QZT1LhntQNg1P28xBp
v2ozjcVj6t8FW0Ln4eJ8t14YIxUOHtZTciqbhw7D8fWrTNj6tSLlGMfnoU/N
XyMOi4vccspj9i2Ank+u/l8s3n5CZ8LBuwAP5M1YRU9h+d8RaPHVWABRhqN+
d7DYuX86bvXlAtDMWRapYzH936PecqHY/RffBl3FYnVS4cm3pQtw/UMc9zL2
vhSp08SM0QW4yv/G0heLV18HVP2mWYSM8A3Tk1ic/Cr/Ut+lRYhIfcTkjK0X
DspUd6guAmOk2FvA8pmgPylVarMId3ZN4+ax/DXU7gj4By7CTOkxLzK2Pxdd
ZRR1Chfh6m7B+fFjVDDh1jlLO7gIJYrnvH9h+ynyWc4weht7Hj1OM8L2uz99
OPUy1xIUj42EEw5j+XKjpMibS8B/Q0PHHDuvhLMZ63+fLkF0wq+RkoNUMJeO
RkpflsDMQkuGSEeFm76/L775tQSMq50O9Jgehp1Oa37vWIKQjHOf9mN6YY5M
NoyjLIHNfhuPyT1UUK3r+RN8ehnYb6mYxu5QIHREvMFOahlowjyFFTA9stdG
2IgZLgPLMqLW/KEAk+rRoEn3ZeChEvJ4MT2nH+4bcorDuOeorqMdpndas7dr
lJplCOb5T/Y7kQLTqg2O+nPLEFWdyZyw+r/6KzmZeGgF8loTDgQsUUBcOdSu
7+oK7I3UZDHG6k07aEaMoLEClbXLd09g9Ziaaku/ZrsC9ecZShMmKBBwST60
K2AFfiT5GrBi9SyXnOQRlbcCWnfCbtgPUiAnQ+I/9Z4VoBz5aZ/RS4Eqt5Yz
4+srsC/V4k9DJwVaD35iecCyCmKZ22P1rRS4U9PJmn99FV7cxgul/sbmie+j
3+SHq7AoabD+oo4Ck77sHcyuq4BPwkkzYv3q2XHDgbMYt9idED7yo4ICmWap
3/Zj3FLPL+1wAOt3lSUj4e0Yt9A3Frvey8fm2Z7zSW67q3AjrU3UJ4sCsrgl
4yOca9BrT0z+nkYBEdKOo6v8Gki8HZr0S6RA2aiNb4vxGlypCiVZxlAAUX7q
775bg6Zkp+kLkRQ4EqPodSJ2De6R2euqQyjQddGsmLZqDQJCHROV/Cmg9Eov
qWd8DdgDc/zTPmP9i8bCTXEP5vPb7oWSvClQqDmWc/g4ATL+8dWyu1MgiK1W
uZqDADhjDT6+1xRwNH/BbCVIAELl6SF2BwpkfNOj7soRoO6CxATpBQXOtAiV
e2sS4EddjGb6EwqkBXYrbBoSQKIn6JqyMQWOEZRNHloTYFHne0ydPgUU91Su
pbkS4NeKVxKvDgU+a4bVLXwhQHqjrZ6dOgWGZUwrmSMJMGSEr4pQosBRq6Es
4TQCiCpf+ZsiRwH106sOUiUEKD0RfSJSAtPHt/wpoQYCcDGZsL4UogDhnXM/
Ux8BPkn/uXSVlwK6TnlXZ2cIEDosqdrIRYEw3mNFCRQCRHheD1NjpcCBeQvd
u7REsKg04Mo9ieljCz8yfwLjCLLa1hY9pl8GM34bTiKEXN6U5qPFztM9n31c
gAguv79T0D8yzHXf/yQnSwQen2wpSTIZThznVPukRoToA7Tcp5bIkMDIKIfT
J8Lnxp/43kkyhO21kp16SoSSkT4Bt0EyuEQEXCD9R4SDhSpuBzsxv/njStuy
F8ZNyx3Vzo1kOMJJf6YrgAie9BEsrTgyiGxeH02IJkL/naOBB4vJsHrkz1+T
TCKYO4/c4skiA21JmRp9GRG8PfZq8ieRIfUiX2ZUAxGSbAoqT0dh/jfpyhYb
xinXjr7Jngsiw9aJ1lPvMU75VcQkHfWZDF85zfv6MU65zuL3StID86dVL8+d
3iGC76K4S4kTlp9TXgM6TAJceZQJpw0ZwjPZ87VZSHAqO1zB2pwMWUwjlXcu
ksBeUl4i5iEZWNb16sVFSJA2tPag5C4Zlr7yZhxAJMhZU6oqvkWG34/41PFq
JHB5wP35pwwZhjztfJ48IMG/m+ktViJk6GK6IUUxx963pJDEzkOG9Ofcys9f
kkBJKUCgiB3zo659X5vcSDDd+NZDgpEMdVcfLDF9JsFtXrHSnwfJ4PlzSEE5
lASp+zkoy9skuO7b5WAah30vY071PJkEud9XjJ78IsFC65tBuXkSiIQyzmiX
kiDq+j5QGCHB8tjzuct1JKj7ssYk0IldU/WVZzpIkLwyt7mnngTh6m6zn7H7
LOr6PpSVkaA6JD+JZYEEgi/nRw2zScB+is3Kl0KCHqIH23wCCeNSesZljEvS
U1QfPfxOgsAAVm/hw2To33+5PP8bth7B8zFGzNg5yVw32PAgwacRcVUHLswX
k8b1z7/C3mOfaWHLh+VJLe0SsyLB07usvdoSZFBbjJkXNiJBSKqWO7sCdi5/
U4qYtUkg2fJUqU2dDIMfLhpPK5GgwzrniOV9MlT2J+yNvE4C4YXqkjkTMtTW
E8qk+Elg5PefsOYLMizaMWRVc5JgaX7kQeQrMqyNyv4TY8K+t3XvVNs7MkwU
/cEFHMDyD3JWX8B0szu/ztn3lwiB7uenF4IxjtK+IEO3SoTHw/b4doxDcHZp
1zgmiLCj09QSlYLp8tASz/luImQaty/fzcU4puud7LF6Ijj6PqNbLCNDNjsK
mi3B+Lq8au+zWowrXIXVUzKIEPpIua6llQxHbSI/62J1cDYnQ5qlH9MZzabd
ElYnwSfnb6tMkAE8/nBZexPB+qLogOEiGVhPfcINvML4/uzxzgdYne5/Vmcr
/IwI9WTEfv1/3KF7QMfxERFsfjXH72Dccc6171ucBhGe+O+7k3qUAiqbG4rl
iAhUfZeD15mxubNu9L1KGKsrD/rSzHPY3DOUyM3lJoKDdYfGgcsUKLEdz/U7
RYSP6WK/lAQpwPNiBR4cIALTlalaS6xvbQtd3Dn6hwAmNs3udhhn9A55BGYt
EsCj6Uq9gTLmw62E4mSHCTAS9vYD/x0K9Dzw1y5uIYCqqV/GpB7Wd29ea+QA
AnhdenHe1Qj7/vvhq47Z2PPvnvf8scDmUte7mIJYAjiPXc19ZI31uREljelA
AjAUByXGYlxhvuDxYI83ARILOn40uVDAzKFs/tB/WB9dMng/4In55rqUC3uf
EOA3DVW96RMFqOGC12b1CBBsHTcRg82VqDRz6RJlrI8GFAo/DKPAPYMv5m8k
CdCz0XR1I4oCnOqnenl5sfe9rqt0SsD65mpsUf0ZAnhbdQ0MYXNOiHdHSJee
AHkNz15w51DAvUu5+yd1DUZMVZ9oFmG+wsOB2aF3DTL2vSsywOaq2tqF1xeL
1iAk+vFdzWqsj9/T1ygMW4OnDh+5uBux9VaYT190XoMD4tNsQ9gctxcg+jk+
WIPXUslCr7opsDWyHRsnuQY3wu5qUQYoIPZg9EU2yxqUqaXZ3hujgHHGdcWf
m6twwqLR7fs0BSZCDzm/6F8FrvMMNrgFCjxM+27EVrwK5cwXrv3GfIu5s49q
etgqCGjS4ArJFGAsKfc767QKGy+WDvhsUmD+pbOnrd4q5HZY70hsU6Cf7aZT
gtgqPC4wC2vGfHtQwwa+jGkVslTJdYqYTzukXFxaQF6BwWFfjyjM17Eor+cF
dq7A/RVK2QDmCx1drffp5KwAu5bQ803MZ67juGk2/FZgRT3c+w8z5lPFLGnd
bVZgJLt9e5iVCgPmnUZL6isgTW5sjjtHhWePbzyXwXzVB9UqrPFRMV8x6G7/
P99VJ6/Uc5EK3Sb0lG+YL0t9xjF8gxfzzWrvOPxql4Eg5xHvz0cF012hh6/i
lyHRvi+wSpAKH3d0thQ9luHD0/HYbhEqHKY1M9o0WoZ3zckNjWJU4GM7MxAk
swz0OT67MZJUaK9+WcvMugzM37nRQ2kqyNOkuHpsLMEBf/t3VFkqKLmJPeju
XgKCVEG5vTwV8CEdOcdyl+DI4puVDgUqWL21m7rmtwS6qTFHGJWocENHQfnG
iyXwvqvFJK5MBTqhXAnJ20vwL2D/rtRtKljesjh25tISaBzyrudSo8JX8nXe
aZol+KoyYbmoToW27lvj38cWgX3GfTRYE/PlwwN+UuWLEOt1/ALXXYyLzJ2C
6sIXQWoqQMxPiwqGCgO35P5bhOS5u8fGtKnwbWhlf5zWIpyOn888rkuFO6J9
fCT+Reg1UmC4cI8K9rSLpwXpF0EyPoT/jB52vtUZ3PqzC7AtZ3mIgMULuN4M
O/wCfNDwjU65T4UWgTvHX/1cgB8yV1YVHlDhlAhf+rM3C3CqI4NcicWtpJxB
Nb0F4Lra8+usPqaHNK2tsyILcKMx+5wBFgvmpliPHF2A1c+RCm5YHOfv/vXb
wjyw5yLWd1gsfqCxkL92HpSnQxLNsfg8R6xqZcw8KFyZmb6CxWNi6knIdR7q
1vD9ndj3pA6J8+Ten4eyvYyeBlgsUE0QYL42D8e/Hpisx9Z7YL/omRcM8+DJ
zPWXCYvdSENK+QtzMPxOoFkRy49HZp5xrWYOuHfUH+pi+1HL6zJ5NmYOxI/N
J93C9os9of+4nMscBLqW5LPqUMFBzmJHW28OBO40+rRj+92BEukNhOcgu8Ga
0xI7j9Z06RD9I3PQ+VDt/QR2XlfVFrbV5mahL0q3REaDCmz5XxKu4Wfh063f
uNfYeV/XPrtwPGoW/rFKx4RjepCwybo67TQLG7NeBt8xvciKOTRnas+CmG3b
XxdMT1dwv67a8M+CiHnoO3lFKjxnrUq9eHAWOnWiybOYHqNVOJO7J2fgR8SO
vq0cppc7D7xdKmYwX8tcOoDpmVuSPo41fAYMboeyXrpOhV9nloxz7WcgdsLF
U1ecCp5s1VyKGjPwmOHkH/NrVAh3faHUcnkGpK4T398XwjguvJNPi2YGNIV2
hfj4MX2vECTbhqch94QtzRRWj3Qcnb+Vi6ZBWLyKxvXS/+rPV7EsYBqc2aql
NrB6PvlNY9+VF9PAoeGXp8uBcdxeeZvgW9MgqfHtVTAbFRqa4if/ck7DgsvL
b4WnqWAXicMb/psCwx0ZmjKMY+tN5mwreqZgJUF2LBbjUqn7L9VYs6fgTvtR
wecYd/LuVhY6fp6CB3rBW8wHsPPnk2BsMZ+CCIt9Ook0VFCjvYXnRlMQajki
y4xx45u9Bhdfn5mCzrnFhmcYJ85vL0e1kCfh5fObe2IwLqwTsP3G1ToJC480
CMVrFPjN9/DefymTsDiRn5i/SIGLGxfvN3pOQmDtf/xBGPfJbdymsBtOQot6
b7DuOAWi98n4OkhMAvsZy4VNjPPyLG19fp+YBN19B6+/68HmzaSOHPfyBMxu
jQQstmFckq57+m3dBDBMGNBIYlzH+6bPfihmAtL0Fb5b1lDAwUc6V9plAuby
Qi1cKykQ69svHn1vAt6WZrxyKKaA7+C2CZ3QBKx/4xrWyKXAY3li8Ev6CdBw
HU48nEGBQJ5VoYnpcWjrY5jOwLgtJPBC2D0Yh2st15KFozFuOVPM0/Z9HNYl
j++JCMfWu1F5S9NxHFZdcndnArD8j97Q7tIcB7E79kWnvlBAf0Y8w4h3HO6F
ESR5MU5TKvpYQKAdB1/x06GcbhSwm2D9yVgxBl9GWSb+/UcB5peMebr2Y/B0
wo6vzIYC+xcNnydeHoP0vVOfjZ5SQCpe9CPdyChM3XA5NvuYAnSfk387B47C
g7xvLVoPKLAiEEHeUh6FLcXC/ti7FJDprKoO2h6BzC/9N3tVKBD+Mv6fbN4I
jIZzXFqTp8BOT8CjLcsRUC7f9l+WpABNXHp827kRoNfli2jFOGzceCCouHsY
uELT9UJ4KLAa0rRU9HkYtvXxQ/KcFPjavGnVKj8Mzya25TpOY/M5JbJ9a30I
lmKEQ28dw/zAn+whxcwhkDjgtBFDRwHhBqpqgtkQ0Nzkfz2+TYZI65nhc6xD
wM3wR3g/lQyj8UV38tsHYe+zGsljy2T4EGFkYvFhEGSeVSbsYNzFcvh9m4jc
ILRvPvzWNUAGJ82v99mpA4AitOl928mQY/+58GLGAKhfjbsmUE8GJm2VVHWz
AaAkfebKryDD2zS0HMg6AJ9m/Zc58skwK/Zd+k9HP1ibumW+TCODafUVnXef
+qGZ+MgtNYYMZlr/rfLc6AfCJ5f/6kPJ8IQnvoG82QcFW9ZZTV8x7ll//2Mi
uw9MCg8oFXqSwXsl+ALZsg963jy++8GZjPn+RVZerj44eL51SQ7jLHe9m9e9
B3qxfrMqPGJGhsN125L7A3vBQ/GHrKk+GfTCT3RlqPbCs2dU3lZNMtzXZ8a5
0PZCmeyDk1w3yeDRJ1NsX94DgkSzk4+uk2HklqZD4H89GB9uKrkKkKHFIi1l
SKAHhJJEG725yeAXp3Tsznw37LC/LnZkIUPvjWTDldhuuNlmLqR+FONCn7nb
RY+6YUulV+MwDRk6BjM/ZZ7qBu2zS/y56ySo/yPQ2d7WBTIBuyvySySYYL3a
z/W5Cz4+/JVeNEaCh+eW78Te7AKd2+3ezN0kcC1VXLmz2wn5uy5hhg0kiEv7
ZCZQ2gl/4/bTfCsnAcvYdwvZ/zpB5NvOZALGVTq/OhPdhDrhyI6ofiLGVb4H
GHvXljrgbFOel284CUKTipODkzsg5hjew/grCVpVfja9MOuAL0OT9mzvSJA0
dYbkzNEBhInuV+BAAqfhv79Lhtrh36XATPWnGNfc+LhfKLwd2DLvSFQ/JEH3
Pve7w7rtkMflJsqtSYJ7Escfwol2uP2ppvbFDRLI+FZG50W1wdDioSMxoiS4
5LKn8gFLGxD/WvFUXMa4cHfpqW1QK+Bq5WRqWElA4DShbjK0wmH9U88KjpIg
P4G5Yv1zC6hbF7X77yGBWWQok9WBFrjCcClEj0yEtSWpu7rezZCtENFPN0uE
yCne6sK9zRCV710R24/xfIM2xl5NcO/+tyeXm4hQaG0VsrD9G6QYgtbDKojg
3zYXXO76G15qp/uTs4hgkrHnyNGtRmDq53ogFUcEsYtJysMujbAu4mz3IpgI
r409gq5uNUDOG4vtTz5EcBEkXN9ybYDbqz3ng5wx7ire631npx7kwxb3fbQi
Ah2XQ8mld/VAs7xZb2lAhAsKGZxeNPVwxP6zv5gmEWLul+9YvK8D9NjAaxXj
KMcSvbD2g3WwJvW4NgjjqP9C2vhqvtaC5NtCl4sYR+EzWfYonKiFnWvFTfFM
ROBmrbyrFloD8uHLgwx0ROA7OmUwwVYDJ4uOtj3ZIEBRlJzJ/rhqWKM/0Jo+
j3FLIMv3Yp5qSGWlbo8MEEBO9JbcVhYeKpgvvt/+TYCac5c/dUjgQf8V02f6
coyD2OKyr+OqID3I7tr+TAKob3Etit2ugn5l5pSVKALseRP4rL4TB84+J47V
+BJgkMXQZM0AB583f/h9dCfA/LnmfUZtACuhXFrSdgSQMju5yTBbCa9uirsM
GxMgEvg3ydsV4Pk0WuyFFgFOUGINj52ugL1emimLNwjAltXmZiZcDjJtJv/u
ixDguudqGkm9DBY6uO/nnydA/ySdKDwrhaPr9NO7JwlAFVg2b/pQAvzcJb+l
aAmgVinucya5GNxS1p/+nl0D6obsSHp9EcgbBAbs1K1BocPJDO+FQjhx9ZP5
yeQ1GONyFk46Uojp1/4a3Yc1uOQ6nnJSuADe4iS1Bp+sgQaYqzbr5QOJheWE
3601+M6PU+1xy4P7ND/rLlxag60HqkOCSbnwSjqlNnL/Gsw5EXlH23LA/I+Z
KWV6FRRxoiaTf7Phy1ENokDNKrxsPZcidzkbhA7z16vFr8L9C+9uLqb+Au4X
rNfVPFdh+Jv22LGDmfAy3+kHv8kq1M816ZMs0yEkhHSLjFbh/OLKTEhzKrQf
9/L9cW4VCGcjfu2IpEBBuFUj7/YKZB770Cf0Iwk2ha9pRg6twGScmJvQwUQI
CvrlQCxZgeQns/l7neLh/C2qP1/4CmQd8PuaMh8Ln/EW2+qvViCtu42e1zAG
Nvzij+vorkBgfKchQ2oUnDQJPY9EVsCh1mziW1QEhAk9+3Li+AqcP3Taoz8s
DKQMP8X/XlkGXZTHsB4WDD01idVWTctg0rKivhoVAEfHMm+SUpZBRy5zNEDn
G5zra3c38VmGENuPomxOH0A00bW12GwZpjvdaz8keQDxkuHnTfllEKM0XG+p
dgI3iUdHuc4tg2+pyw/WREvYqIvMEfm3BPKvLplm8AsCC/FtB38/xkmXFvUW
fc2RHWN0/cmCJWhqq0bbtv8hlqm21YmAJWjbTdymk3+HnD3Y4yJtluB4pcXj
2rL3iOjfJCOvtgRc7mohV69+QbEquqztl5fAebK4/pSyP8r98DRUbd8SRLXS
Zj0iBSKx+YUjeeOLIKZbGm4cFoIedJAH6CoW4ZytXCa7eDg696ldX/n7Iiit
C1/82hqBIpS25x0xblrj+Xg70TgK/ZTlmfbDuKnhcJiD6XI0+lssmBKOcVPM
3WCfUclYJMzJFeh3aBEubSqp3HeOQ5mHJXcdZxbgE1n7Znh+PMrRHlG/XbUA
Dgc8VH6uJCDbM0n9R6IWQP0L8coT7iTEPqa9B+e8AB7hvL+X7yUjnS+tXGa6
CxDNzXpU8H0KonPhj1wXXICKvTxTvDmp6MK3jN9vDi/Avu5rNwcH0lCHQDYD
YXYeuCZjhWR309HZp9W1evh5kB40ybnDlYlsJ0M1s6LmIfvteh+D/C8Usr17
4K/zPPhLuVU6PcpCKu0bdyV058GqfbC6+no2SpCI+PhUcB7guXcWe002mhi7
wPCFfh4cHVRkGtRzEIvmunzczBzEb/kdxHfloCcLX99l4uZAGriLDunnojVa
H77MyDmsvwmuRQ3nIuOpHwmxr+bgAdctAzejPGTRmKryRWsONhfNfiWN5qGg
fzHOVnxz0H+a/J3FMB/dagoIRgfmIO5Ha1/HQD6a+5G4RT85C6GBShc7dAvQ
w7VzJ1rKZ8HU5dzNU60FqMb0gYhP2CwkOEX8/aFUiAygGiduj/3u4Mz4tKwQ
9fybPzumPgv1SpuG9oJFKEXRItn98izk1o5kVsQWoTPbCzgWmllQ63QYUj5Z
jA60Sf9OG54B37t+1cc8ilFokvAF8aIZKGELVGZaLUaZ5XmCZQEzYLhsraGr
X4KYI0s1pF7MwHOhZGjFl6AmLs2xvFszsDlw/5vnlVLU/MVbiYdrBvq4RNJs
/EsRraPHRNi/aTCQ8WIIpJSirgFHJtreaRjdJ160pFeGqi9GIKvsaQiK7El8
U1SGlvglalo+T0Pfht+4LHM5kpCKpBWwmIZTNpL2og7lKO/gucdf0DTIXFW1
MmwtRwkd5yRnz0zDt7jF4aLLFShVjZYgR5mC6tWlSQX3ClSXJT4f0joF07xt
37e7KxBnBY/nUsoULBXUcc/xVCLpCnUOea8pSHApTt55U4mCPrBwBRtOgS+b
vJZiUyUKe0+7MCcxBaeYcu7mnwH0qvB9v/TJKTg4yDJ09wmg6VxmLf/lSbg1
K3+FJQcQcw5f/kzdJEBhtuuhv4D0P0qbycROgvJEH30gJw6JFuRBsMsk7MuJ
/3lNEYcuG+69unYP46x75zkXzHEoXPsLnbrQJNhzsUvmvschMktuXjr9JOQd
rQrwTcSh7/OtQQwzE3A5PHrQrQaHGPNEDv4HGAeFZeHdJ3FI6/JF+7HvE6BF
r7kVtItD/QVTIpqOE3BiXE6mlK0K6YfEpFdpTsBGWeZ1kngVYqv3kJXinYA/
2QK/pLWqUDlbzO1i2gnQ3Ky8F2pVhRok6IRlRsfhbWzt6l6vKjQSnqVbVzwO
X6WZpN5GVKHCiA+n7geNg9Lrub/0OVVIztWye9V6HEpa5taT66oQsmAlfVMZ
B1PaK+u6Q1UI91y7TZR7HP7//8//D+cF8ms=
"]]}, {
Hue[0.8503398874989481, 0.6, 0.6],
Directive[
PointSize[0.004583333333333334],
AbsoluteThickness[1.6],
RGBColor[1, 0.5, 0]],
LineBox[CompressedData["
1:eJxdW2k81O/XFhIVJZJEUimFimghnUJSpLJEtrSI7CEVopK1bJVCtrJkJ2vW
Y9+3YczYDbNYJlF+hVDP/N8+r+Zzfe7tuq9zvudc94uRuOWoa8nOxsY2xcHG
9r/f3gRXmWIyBf5poywE1mH6FdrHzEEKHKzZIhi/oQ6v2THtHEcpwObUTBcN
qcUnCx0+/8YpcGJnzETJ5los3Gs4bMqgAIn60+zhmxrUse908p+igNbo2Ktr
QjXY+JnrmPc3CpTNSX41iqrG36KXJdVnKdDvIrHHW7QaRWO8jpN+UMBmuW9z
a3wVng7+bqn4HwUe887Vquypwjvsnuk3f1Mgt1prNfcSol+j/z/jRQrsO5a0
9qRBJa67//uGxB8KfBGv3SNkVoGL1kcaC5YpwMmR0n/RshzZlyXktq5SoLDO
KrXfvgw7zZWjVf9S4O8hpc31bqVoU+3079Q/Cmx7fCNs09MS5L6kZM7JNgaP
CHs+lwV+xXeH3xZEsvClaRKh9U0xrsb/ZfvDws4B9a9U44pwceyJqtSaMXi/
9aje7rRCjNKMfrSHhS0U9OLtCwow+rhFIpM1P9hQbXpfVT4K7LGofcrCdBPO
zCtteSgbkt0/zjrfXNLVZKbvC57lpk3ws/DX0+cd/tFzMV1KbI6XxffQImfe
lokc9A788x9phQKXVafEdnBkY1BVwpId674gaMEptisTDQy/rulaooBpDvd5
sdPpGF9tz8/G0uuMylzdLrNU9NW9L7v6iwIeir7GMk9S0OvIgnHjPAU6m4PG