-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathSurfaceIntersection.nb
More file actions
10741 lines (10697 loc) · 613 KB
/
SurfaceIntersection.nb
File metadata and controls
10741 lines (10697 loc) · 613 KB
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[ 627397, 10731]
NotebookOptionsPosition[ 624599, 10676]
NotebookOutlinePosition[ 624935, 10691]
CellTagsIndexPosition[ 624892, 10688]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
RowBox[{
RowBox[{"M", "=", "0"}], ";",
RowBox[{"\[Phi]", "=",
RowBox[{"\[Pi]", "/", "2"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"b1", "=",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-",
SqrtBox["3"]}], "/", "2"}], ",",
RowBox[{"3", "/", "2"}]}], "}"}]}], ";",
RowBox[{"b2", "=",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-",
SqrtBox["3"]}], "/", "2"}], ",",
RowBox[{
RowBox[{"-", "3"}], "/", "2"}]}], "}"}]}], ";",
RowBox[{"b3", "=",
RowBox[{"{",
RowBox[{
SqrtBox["3"], ",", "0"}], "}"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.878776161513973*^9, 3.878776244915511*^9}},
CellLabel->"In[73]:=",ExpressionUUID->"68804dc3-a75a-4d19-9ccd-71ca82c0c716"],
Cell[BoxData[
RowBox[{
RowBox[{"eqs", "=",
RowBox[{"Thread", "[",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Equal]",
RowBox[{"{",
RowBox[{
RowBox[{"1", "+",
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "+",
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}]}], ",",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "-",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}]}], ",",
RowBox[{"M", "-",
RowBox[{"2", " ",
RowBox[{"Sin", "[", "\[Phi]", "]"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "+",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}], "+",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b3"}], "]"}]}],
")"}]}]}]}], "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.878776255785833*^9, 3.878776255795537*^9}},
CellLabel->"In[75]:=",ExpressionUUID->"54c04de4-0bf1-4798-950c-386eb3107a94"],
Cell[CellGroupData[{
Cell[BoxData["eqs"], "Input",
CellChangeTimes->{{3.878776259319687*^9, 3.87877625974728*^9}},
CellLabel->"In[76]:=",ExpressionUUID->"4b4c50f1-9148-49dd-89c9-b2ba011b4c10"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"x", "\[Equal]",
RowBox[{"1", "+",
RowBox[{"Cos", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "-",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}], "+",
RowBox[{"Cos", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "+",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}]}]}], ",",
RowBox[{"y", "\[Equal]",
RowBox[{
RowBox[{"-",
RowBox[{"Sin", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "-",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}]}], "+",
RowBox[{"Sin", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "+",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}]}]}], ",",
RowBox[{"z", "\[Equal]",
RowBox[{
RowBox[{"-", "2"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
SqrtBox["3"], " ", "kx"}], "]"}], "-",
RowBox[{"Sin", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "-",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}], "-",
RowBox[{"Sin", "[",
RowBox[{
FractionBox[
RowBox[{
SqrtBox["3"], " ", "kx"}], "2"], "+",
FractionBox[
RowBox[{"3", " ", "ky"}], "2"]}], "]"}]}], ")"}]}]}]}],
"}"}]], "Output",
CellChangeTimes->{3.8787762601353197`*^9},
CellLabel->"Out[76]=",ExpressionUUID->"711b8541-4b69-42c6-be47-51dd3b276985"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"impl", "=",
RowBox[{
RowBox[{
RowBox[{"GroebnerBasis", "[",
RowBox[{
RowBox[{"Join", "[",
RowBox[{
RowBox[{"TrigExpand", "[", "eqs", "]"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[", "3", "]"}], " ", "kx"}], ")"}], "/",
"2"}], "]"}], "^", "2"}], "+",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[", "3", "]"}], " ", "kx"}], ")"}], "/",
"2"}], "]"}], "^", "2"}]}], "\[Equal]", "1"}], ",",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"Cos", "[",
RowBox[{"ky", "/", "2"}], "]"}], "^", "2"}], "+",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{"ky", "/", "2"}], "]"}], "^", "2"}]}], "\[Equal]",
"1"}]}], "}"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[", "3", "]"}], " ", "kx"}], ")"}], "/", "2"}],
"]"}], ",",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[", "3", "]"}], " ", "kx"}], ")"}], "/", "2"}],
"]"}], ",",
RowBox[{"Cos", "[",
RowBox[{"ky", "/", "2"}], "]"}], ",",
RowBox[{"Sin", "[",
RowBox[{"ky", "/", "2"}], "]"}]}], "}"}]}], "]"}], "[",
RowBox[{"[", "1", "]"}], "]"}], "//", "FullSimplify"}]}]], "Input",
CellChangeTimes->{{3.878776294724728*^9, 3.878776294731744*^9}},
CellLabel->
"In[105]:=",ExpressionUUID->"dbc3e449-1aba-4459-8727-e1d3af016780"],
Cell[BoxData[
RowBox[{
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"3", "+",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-", "4"}], "+", "x"}], ")"}], " ", "x"}], "+",
SuperscriptBox["y", "2"]}], ")"}], "2"], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "3"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-", "2"}], "+", "x"}], ")"}], " ", "x"}], "+",
SuperscriptBox["y", "2"]}], ")"}]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+", "x"}], ")"}], "2"], "+",
SuperscriptBox["y", "2"]}], ")"}], " ",
SuperscriptBox["z", "2"]}]}]], "Output",
CellChangeTimes->{3.878776311458891*^9, 3.878777215076786*^9},
CellLabel->
"Out[105]=",ExpressionUUID->"4c929ac1-88b6-49f2-b35a-7db38ed6003c"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"sol", "=",
RowBox[{
RowBox[{"Solve", "[",
RowBox[{
RowBox[{"impl", "\[Equal]", "0"}], ",",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
RowBox[{"InfiniteLine", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"10", ",", "10", ",", "10"}], "}"}]}], "}"}], "]"}]}]}],
"]"}], "//", "FullSimplify"}]}]], "Input",
CellChangeTimes->{{3.8787763493262367`*^9, 3.8787763493340187`*^9}, {
3.8787764519164124`*^9, 3.8787764525123243`*^9}},
CellLabel->
"In[106]:=",ExpressionUUID->"f97e2b78-14e9-4aa2-82a6-05233daa5a90"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"-0.814\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "1"}], "]"}],
Short[#, 7]& ],-0.8142309185120151},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 1,
0]]}], ",",
RowBox[{"y", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"-0.814\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "1"}], "]"}],
Short[#, 7]& ],-0.8142309185120151},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 1,
0]]}], ",",
RowBox[{"z", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"-0.814\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "1"}], "]"}],
Short[#, 7]& ],-0.8142309185120151},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 1,
0]]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"x", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"1.31\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "2"}], "]"}],
Short[#, 7]& ],1.3099636521649802`},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 2,
0]]}], ",",
RowBox[{"y", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"1.31\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "2"}], "]"}],
Short[#, 7]& ],1.3099636521649802`},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 2,
0]]}], ",",
RowBox[{"z", "\[Rule]",
InterpretationBox[
TemplateBox[{"Root",TemplateBox[{"\"1.31\"",
DynamicBox[
FEPrivate`FrontEndResource[
"FEExpressions", "NumericalApproximationElider"]]}, "RowDefault"],
TagBox[
RowBox[{"Root", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "27"}], "+",
RowBox[{"54", " ", "#1"}], "-",
RowBox[{"17", " ",
SuperscriptBox["#1", "2"]}], "-",
RowBox[{"58", " ",
SuperscriptBox["#1", "3"]}], "+",
RowBox[{"78", " ",
SuperscriptBox["#1", "4"]}], "-",
RowBox[{"40", " ",
SuperscriptBox["#1", "5"]}], "+",
RowBox[{"8", " ",
SuperscriptBox["#1", "6"]}]}], "&"}], ",", "2"}], "]"}],
Short[#, 7]& ],1.3099636521649802`},
"NumericalApproximation"],
Root[-27 + 54 # - 17 #^2 - 58 #^3 + 78 #^4 - 40 #^5 + 8 #^6& , 2,
0]]}]}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{3.87877635044613*^9, 3.878776453785042*^9,
3.878777216117948*^9},
CellLabel->
"Out[106]=",ExpressionUUID->"ba3d560a-6b8c-40ed-bd51-d119efd01173"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Show", "[",
RowBox[{
RowBox[{"ParametricPlot3D", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"1", "+",
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "+",
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}]}], ",",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "-",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}]}], ",",
RowBox[{"M", "-",
RowBox[{"2", " ",
RowBox[{"Sin", "[", "\[Phi]", "]"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b1"}], "]"}], "+",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b2"}], "]"}], "+",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"{",
RowBox[{"kx", ",", "ky"}], "}"}], ".", "b3"}], "]"}]}],
")"}]}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"kx", ",", "1.8", ",", "5.5"}], "}"}], ",",
RowBox[{"{",
RowBox[{"ky", ",",
RowBox[{"-", "2.2"}], ",", "2.2"}], "}"}], ",",
RowBox[{"Mesh", "\[Rule]", "False"}], ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"Opacity", "[",
RowBox[{"1", "/", "2"}], "]"}]}]}], "]"}], ",",
RowBox[{"Graphics3D", "[",
RowBox[{"{",
RowBox[{
RowBox[{"Tube", "[",
RowBox[{"Line", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"10", ",", "10", ",", "10"}], "}"}]}], "}"}], "]"}], "]"}],
",",
RowBox[{"{",
RowBox[{"Red", ",",
RowBox[{
RowBox[{"Sphere", "[",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",",
RowBox[{"1", "/", "12"}]}], "]"}], "/.", "sol"}]}], "}"}]}], "}"}],
"]"}]}], "]"}]], "Input",
CellChangeTimes->{{3.878776360997806*^9, 3.878776361009035*^9}, {
3.878776456795087*^9, 3.878776457356591*^9}},
CellLabel->
"In[107]:=",ExpressionUUID->"5901cca3-956d-42b7-85bb-c925b9f2ec93"],
Cell[BoxData[
Graphics3DBox[{{GraphicsComplex3DBox[CompressedData["
1:eJx0nHlYTnv3/+/mlOaQkEQTSUSIrC2ZQ+apSUIlyZTMiQghEaGSQgkhlaRa
u3nUPM/zPNzdzfOv59n79lz1+zr/nOv6XJ3ss9da78/r/V77nDkmp3cd42Qw
GBt4GQyusb9Hv5BZ83FNOzD8V754WGUFu8IWZkYOiJBbc4pfTvvTA2YmRXem
m3xH/kFVNYYcP5k72LhS0GsY/NattHtuHoVWn7ZkxYgwyGMqXZ2OZRyEafuG
hLcrk3De9A3SynIsNOmc93S9OTfxTXXj59LQdDxwbN2z2OFydFp5+PeW7bzE
lXuD6us3ZSGjbWTzi+woVN/75pDueT5CMXhy5ihmo92Wr71uWTcwvIjM7Erh
I7hlXSp51HLQ7pSBmPThHfisv/bHdz0+YvbISp8bLmM/b2sV3ajtAivEtbOY
/LzEurQ443ftmRh64Ghr/c4E0PtSbJrMw01YTxK79jIsDfN+NX04mFEFXb/v
J8804SBsm7nFPAcT8E/anbolOl1QJuf9eKv9EDirXKjO6kM0ebZrkdfYz22X
OdG9/2U3VDzdyjXpzWdMOnd8afXtSYSvnddKnZQWkLX4KX3xqDPwcai1874Q
JRbbymz4tb4d7Ky+HFKXv44G5/xatLOFSRm1EZdrD3tA7zdY3+kl4VGZeOlJ
Dj7SiHnX2Et2GOY1FtnljmRAspHuuU0CI+hWtyBWy56DmGqx+NeAZwFcVVxp
6LOHicfNLFlLu7mIJUnlN4/0lAHzkWJcyYYiPPbzyxH7RB5inZprj2pBJYhe
UQu80f8D300OfBJfzUuUL0Oepg3VwCjyNe2JfQ5LZZdsWbWCj5A/IL1pbVA1
GB9eopcRa4XiuoHltlG8xLnET3/spKuBdPKawRR0wVo8wP3dgYdghcrbf7lZ
Cd94l+TwfImASINzEvquXISPUKv60pFSeJdZu7PEuxxuB/cqGghwENc/Xv31
zCwfjr06KxxkxgLre64tekJDYDvD/VXa0nRIP7Hvo/ZuBjH0sqdi0LQbQmsU
Ix63/4a2SwtrhU34iR2iHObQ2AJeQhu3Dwy/wbaRtCEbLlEisiXopeWFdjDO
vGnKu+UJDj/qYD3j4yMTmfYZreVd0PbnvYe2TTZMCzUvxsQRbOsIiROO7Yd5
eZtz3l2vhgbpDTvudbSjpPut+8nyI6CTrmBdotUCxpdWBal6ZuBm1rT4HbUM
4onq89KRqA5QK/D05vJJBOM7Cqs6CU5il9gqz66LXVDiGCznfzkXusxVm75I
cBHNwYdfe8zsAbefnPemJ6XCwrX8cva7uIgZOo9nTHrVA5YjHM57l77A4kt3
B5UGOAne8BnXPzd2g/hyg6W65zNQBj9IrxHhJLyIcpcwyS5o+cCz7dLKXDRv
jFPk9GQQ3sK2iX7qHXAghW+v2scoFD3UZ2l1bhhK4xUXZxY2QU63h1rn+XwY
nhKvA719sOt51oXU2kow9M/u6tzDAvvMo88PcnWB55CciwdPBoz4R6ZkBDMI
4HWaFfGjFTIeOHaoHAxE6ZvDh75o8ROeN74sOxzTDrJSlbwq619g3YfZ8+c/
78KhlZtSZORYoDs8IpS7uAjCLX+cv+mVj/oXGVKSG7pAY4Xf71laTfBt0sbg
NUIlsEGFQ6JfrQfEPz2wvsrRCdpDeY/VLrdDiBnPKvWkXuiq3675vq4XfD8+
hZ58Fki+OlRwsb8PRqRs8j47DYJa15rDIapM+K0Xt7Hhez8oXCRZxdHDUHB+
iHiQXAmyGQ7ltQ39oEneiPc/OQIaxj6D5yw+4dfCVUdcXveD1WEpw63Xh4GL
IZ20Fmtwpfr79T1/+sCsc9782XyD8F7nqYIJRwfaf66sSr7fC5mLXL6JKfaC
8vrGBvNwFmaIa69uruqGlivNIxViLHiuu3/GwnWtqBBws/NsSyckBXxxzZ3a
ADrytx2nQw5WmllKvRjrm6/1r+XvOeeC/sA9h+DACqjKCLZ1HmmD1KGUvtmB
oej/WXnrZoFecJTJPKAtwgSymndmsIob3rlQ9eDnySqw7jm2eX9RK+w0cVh2
UaMYuq68Wei7uBcOSJdz5wY2g+01RWU17WY4nGkrIt/JIC7ujQiLim2EphWt
R19EdkKBoN6luTe5iMQloseqlBsg/Zd65/mcPhiOKm3e+Y6TaIjdJ/S+tw7y
zr/hT/Ycgh95fG6fcoagdB2XP+eaOtiytfnlQNMIMG7ttFu7qA1ajHoLP/LW
ARx0Ejz4ahRmPJe8PNIVgqFTH3tc21MH3kzjdPQZgYzY6x9WB3bg8PnDKVZz
6+Fh3R0HQckhuIEJuvnfRrDu55J1gWYNUFL8iTgt1AdxBgcKR/i5yNs+9YI7
JZrgNEdR1ac4FsSu9w103sZFrh80zRcQagFviZCE27yNQCZ+NXXVYJArFwfO
vLGqDXa8t4+afiIPGrJ8+e0HutC20+ZdPMGEZxuSNhf4/sKMt4pv8iqLUFNE
q8p+JROM4zs+y6Q9xS+rYjRMLTgJm2v85KGGeigIZfJt25oHt/WUhj5PmUQc
Kr515NLTclC2VnGpH/v3s+0Tv7nFXIhQe5WyfgFfJrxf9TJhWIUJSkFLXi7e
IkyAqELV/NAfsOisjbPT/i6Q3Vfiu0FNkEg04fO7cvAbbuOYuenG2V74vjJw
kcRGbmLJ9EaffTNi8Oxg3AejPf0wcDHxmWxeF3hptvlFGMSj3J55zwRr+uH3
k8NceXcRbVuOpJ1KjELzuovh3NL9sHSg5U/7lT40Tnxy4JjJB3z49ZhnDqMX
HmUOTlr+jYcUPfX5lrDMb9BQdi68EdQJgS+TygVGBUm9R8PdB/Zng7e2gLJs
dBtMfTm7RMNCmPyuW38tfFYlDM6fN/DYrgaK7svaXKydTB4uXZl4d18jiLdo
rLA7lAVOz51Eeafwk3qXK03ndjFBT9o/Iow7GKuLZhw1VeQg+SLkry3fP9bP
Sd1lFZp3cH3l7yRO2cnEoZJLHDvDxvSX+4mLOG8SLC/kNXi9SYx4MuvRBi2/
P+DZZD+7SLQQel+e0RabKkF0aS1kHGsPR6NNxZNGFKvg/u6wuwHt4kRJxQf+
788K8fC5fWvkcuvgq1mO4bqXokTTb5H8n17VeGWl+RrznU2wPWVhs8pxfkKl
6MwejKtHW5429c1WLSBetMi1t7wPpBuUWNM4G7HBu1XHSaQVBPh8Qm0SolH+
tlDi2vX1eCB3eE2SSgskCi3j9F47jHKPzhhaf6jC7KbptdmZjZC05ZlBjrwA
qXRjKPFjfB7yO57nni1TB2H6kYFMeTHSbsbJ0yfSvmCYcqaT1LpKeGfeXyAx
S4JUi9SZd7A4C2Y+FldYwJUP8/QtHKc3iJPHkkYTvGprwdpLwS/6ZSxM3VBv
F5kpSvovKdmjdbYDjAXn2tx+4YvnXw2KDBcKkIdUW+VrjzCBQf+1pJhRe8dD
hEiVbN7jv7KKPvcK32H+vL1lSJwYOeWyo8gj6u/P37TWeqj2UZJwk5m7eHpj
OrLPu6rvDj3TlSREo1iWKzdU/T3X6fC1iPwiRkxVJZXrWY1/z7m/3uk+XzmJ
4O6PSc572fb33NJhmvTr3QNwXW9ueNNQ+99zvm38dzJKYrBnr3ubx6T//bzp
/V7xP30jmCbIkhksbPh7vmoqo3YZKUg2eN2VWfuj4u/5L99Zv7I3ipN68pGj
UpMS/57PNxq4df+2JBloWX8t2y/5778vZ8LASlN7SbItzPBYr1/N33MBQoNp
c0+cnFdXXPTkTcff8zffGvos9gmTvHQ/21H9DDp0Px+g+7mF6mfUoPv5Gd3P
76h+xh66nxvofjah+hnZ/VxA97Mx1c8YQPdzPd3PN6l+RnY/K9H9fJ3qZxSj
+3kq3c9Mqp+R3c9ydD8bUf2ME/s5n+pnTKD7eT7dz5OpfsZfdD+70f0cSfUz
+tD9PJ/uZzmqn3Eu3c/Hx/cz/qOfgd3PbH0mKH0Gf1qfL9L6XEzpM9rT+qxP
67MSpc/I1ueltD77UvqMCrQ+H6D1eTGlzyhD63Morc/bKX3GQFqfVcbrMw7S
+mxL67M8pc9YTOvzLlqfLSl9RrY+JyZQ+vyY0mdk6/NhWp9XUPqM32h9NqD1
2YfSZ5Sg9TlkvD5jHq3PRrQ+C1H6jHdpfd5F67M1pc9QRuszmzcYNf/lDbCf
wBvbKN5AFs0bB2nesKF4Aw/RvHGZ5o0Wijcwj+aNVJo3MijewAGaN5rG8way
eaN8PG8gmzd6aN4gKN5ANm+EjecNZPPGpAsUbzyieAMv07zRSPNGKcUbGEvz
xkOaN6wo3kCS5o2tNG+8oXgDwybwxn6KN7CG5g0bmjeCKd4Ab5o3PGh+9qL4
GZom8LMOxc+YSvOz4Xh+RrcJ/CxK8TNq0vz8ezw/41uan6fQ/DxE8TOq0PxM
0vwsT/Ez5tH8rEDz8wqKn3Eqzc/fx/Mzsvl5Dc3Pxyl+Rk+an2/T/JxF8TMq
0fycQ/NzI8XPeIvmZ0Wan2MpfkYpmp+raH4OofgZTWh+Lqf5uYHiZwih+TmM
9oMVlB+EUdoPxtN+sIHygzid9oOttB+UpfwgMif4QYLyg2hN+8GttB98RPlB
lJrgB/UoP4gZtB/spv1gPeUHke0H1Wg/OI3yg2g5TPnBEtoP8lB+EKVpPyhL
+0F3yg9iDu0HLWg/+Jbyg2g8wQ/mUX4Qu2g/OEL7QV3KD6IZ7Qdv037Qm/KD
yPaDq2g/WED5QZhF+8El4/MNYOcbs8bnG/iPfAP/kW/gP/IN/Ee+gf/IN/Af
+Qb+I9/Af+Qb+I98A/+Rb+A/8g38R76B7Hxj+/h8A9j5xoS8Dv+R18E/8jpg
53Wm4/M6+EdeB//I6+AfeR38I6+Df+R18I+8Dv6R18E/8jpg53Xl4/M6+Ede
B+y87sP4vA7ZeR0ref8Zj3WdYHezruaVqjN0dUTs6Lo8mXQ4bLttXyUfsTY+
Ee8GxAAjMMRwx6A2xhVM+romnY94b6NjPMNt7Hyr5CMzQQvY05rjlc3HAvIz
i5STuw0qOw2evU8TIo4uf5cmlsICM0XDP6G7fsMb3jtna+r4yDu/mQ9KK/ug
IU6zwj80F3gyvhemdIziapReK2o/AnatxFa9iEo4rrlt4aqzHRikMmnKhm8c
REWvUJzV93pwy1/tjmWFOLUpiEPRk4uomvXcZYl5CzCMn59x3nAWj/4Ycmpg
jvlLxj1Lj/ftkNF/QtS/LhKs8nSflz/nITbVD0lEeDAhw4DR2OX7CerhzWzT
uzyE/beLGWn3x7hAOy4uLy8QVTckVurEcBOj+NGIQ6MdjN/fUej8EYFYteF2
+B4uYtuPAM2Pwc1QYZ3GHXDSH+QlB0/fCemA8J2fVjftjIcNK4PyZ6QO4Nn4
OSmXznaDxvor585FVkKydFBpYGc93m5bXPG9tg9KuSrj7QZbgFyfqaq93x0W
5b1KmVU1CInXOrKTDnWCjInJwmOZFRD/ZUpZGYyA1Rsz6ePLeiHjsA2/I38d
WFtw9uiXj4K++53s2UX9oHPnrtac7ZWAkwp9TXxHwX90oWfn9X7wrnp3KyWw
Cj/NP9ibVTsMCvNtYsK9e4D4LFDOXV+Lan/megRaDMKpm8KvN7uxYKjlygII
LkXPefca5+zvg8tzTrjOXNAMxG23Ne83hoCd8YvhHdpMiF+vJNM6Ixnm2Z3/
UOn5E0LmnthdlsUEyxVmui+iaiBH7MoasT+tsOW48sxTBzog4oNSwN1OJkxt
m6ayyqEP3OOThXlSO2BF691nrT49sOhYWcZjqyE4zu8b6anIgtRYYiOfwiBU
4Jf2eGY/vNu8aPaQAQvip9o6DGUNA9+XWQWzTnfiZ+3v9QxpFkgoaM2xvTwA
84ZsGg79GUDxi6kqPb864HLP/itmZDfYRP7x8Fk+hJ25Prpl2h3gfmEZf25C
OygkP9cKJ3oxLuy0VXsYE05VXJxhf7wKSm7M3P3yQyPmc9u5EMpj/THHjBzR
IoEQe3Zu2OsrSDO5FXVutEKTPv/cBPVECJ1e+3xz4QgknH+gdVugAV7tLUtw
0KoGVoBJynxzHuJC3ZD8JEY18CzhNFY/1Q4u4k7X98byE/dd9y1at6EUymrk
hbyHuiBc2Vg/24OfuPPikxRxNB94XPzvWDzog1ML+5Ze8OYh68z7OzR9yuDg
J4GEMOUucH72sNO8gZ/M6LBJ3ptRDQZFD+3eX2gF8Vflu11P8ZN7uL2OxXY0
gNnsI7OvClfC6g7fT4+cuMltlbqyZmfawN1l6of3MyMhjPXZzv/yECZuvX/j
G08z7GlsjenhioLies3huPOCxBTPSyHJCypgy+5UkQyhYlh+9co5J0tRQvVG
QNTAs0QIsHgsukyqDmwCvyl8uSlGfLg3UBJEBOG+Owm2N+xbQXobX8fjJ6KE
1xE7j5mZ3qjayNHVcbAFeGw8XvtViJKpm9UPigmnQ8jy5FuqITUQ4h4hqHlB
jAzhNY0O7q0EW/+MpsjBfBAreDo5QFSUlGkJvzLLuxGYvUtVT130g9TaBu9A
aTHCefW8jEX2+bDU+MPmbedi4CrnpRAFbwkiI2WpnP/0t/huTWiaQFw6eC2+
WfzziCTBGEpSXM7tBRX+L62GE/6Av3VS7ZvjkmTq3sihr9IlUFI2IKC8Jgqk
zEtF3+pJkI6DnUpiAi1gJ6ZhmGHtCTdlmGuz7/9/z4Ps53lEP89q6nnwyoTn
CaKeB9nPUzBIPU8D9Tz48f9+HvzH8yD7eSTpeulR9cJldL1U6Hp9p+qF7Hq5
ja8XsuvlN75eyK5XCl2vMKpeyK5X2Ph6oShdL9TkE9CpbwEiJMT5xYJgjC7/
88dqhwA5YY7wGz1HSfQceVJzhC30HF2i50iQmiN8RM+R8/g5wt/0HIV9OLJv
ln4eWF17I3N2fj/K6vWHVodyExPmCNlzlE/P0RFqjnAyPUcH6DmyoOYINSbM
kRc1R+hPzxFbJ/9QOombblA6+YvWSStKJ/EXrZPbx+skCtA6+Xa8TqI8rZNm
tE4mUzqJpbROdnu7zTh3nAXEtqv8PlojyJuYN7s7sAMCaJ2UpHQS5WidnErr
pC2lk3iS1sne8TqJUrROJtI6aUnpJCKtkwXjdRIJCUonJ9x3+I/7Dtn3Xdz4
+w7/cd/hP+47ZN93n8ffd/iP+w7/cd8h+74LOnvdeWlOF+gNNW5MvFCO0h7h
3RcSWqA/VW3JH9EOyPlcHD877jc6Sw98vXFoGNh8UkbxCTJMKT4xpflEleIT
dKb5xJrmE22KT5AwpPikgeYTO4pPUI3mEzWaT7opPkEGzSckzSe6FJ8gSfOJ
ffYj/sbtHERf5qe9e2PqME19sYukwJgueNaJvnw5DCPBJQ1ibyswuUD6jEhQ
J0i05y1KKu+F7Ra78f7abIy/sqKcd+yfL5I4uyQvvAP82qZa7hv9T27kcyVZ
Y9JEnsR/8CT+gyfxHzyJbJ78PvtW75tVnWB8s9Lk+x0v9F58X0Vn4WSy9+HU
faazOsG6Q+WC0Zxv8EurOm8HnyDJSo85kbJ6zEdwvxls1TLEA91CK07OEyV5
r8clqB0ZAC+xDLdnYR9BJ+umheBubvJbYEj9u/p+UAu5FXT/fdLY+7qjGNvB
Seaf+GE19TKD8Mtw2aHMFQzHl+iW267px0r9C85ZUxmEW3fZ74GOHAgpNc+9
HteDm7WZ04+N+Vzj73IltqLhEHo6ZO/t+kZ8AozgjD+chPoJC7E7/cXgKNdz
lkyuxQdCF5ruPuUh3E8tFXPMJsHrnYGpkWc2enMIlZjO5iOMLJyylCKjgbGv
OfBl51eMaPiw4mEXN6G7kRkbF1cB5AqbFPGBZEyDjsXGA7xEleqCkkXzstEu
3bYoIvwj6pZ4Hrz+io8YtpDXCdsXA2S5l8dJCWeUvKbsut2Kl1jI9W1vbkAV
MBwrWWuTz+CvZQI1NmF8hIalmPan6Tlo92SGjkadGawhj9qMsvgIhv21FSdj
x/jf9tnpmzYqyHA/xDRW5COSM0yDT3tWg558OVNV7gkYTV36+IEfH7G3yyEr
ljMHGeTVGPlNRpgq3OAi0sRLyMSclBa5Gg0ZyuU+fQ9DgBwBm4dxvIRud9TO
gIqssf4/uF9r2g/om1adY76Kl5AOdfbvsKwCYtabbQuYL8DZ/5FEhBYPoXNM
6rm6OQkHNGMkFuXlQ4iUAPdXb24i6Lf66sTDFcB/bFHn4d2Z4MpKXuFZwUkY
Q3LBmZ4w4LybPe3yzVa4Sc/FKDUXEE3PBbPssb5DzSg0/DRvuRT4A/znr79j
9HXw77xwh/x3XiCanpcK7aqgV+r9wLRy5OVO/wCSQvKwNISHEKPnSJeaI2DP
0VVUVC7mZAHj44onvx8/Q2/JP6LBdkLElR+jCftjxu7xpVImMyfbYUPSuee8
kaJECT13JDV3sIWeu/jKd97bK7pBzfqY/KWIVGCs+KDt085FasSkGqrvZIF+
YbK/h100PPjEZ9mznotcXWKx9UvcIDg57nIsSioB/6nRvJsN+9CkqsSi7Q6D
YEqevrqlqAb4G3ZPrcqvw92uIbsW/OIkjh++1rlyZRPIVtb9arGKw6YyxnpN
fW7iceL8rTyabZBBxqy3/xkCytM3NsMhHiJxZeM7991MqHh2offt6zC4v1v/
2ZZwTqJwR+m3BVbdEOantTF0zB8ftYm9M/Seh7ii4xCz/BsTCozUCuduu4tp
7yJVbPYyiIbS1DWf5AfBzJfz8BRWPLzafT//6ioeYsaZwVujC5lgbbli0g7N
CLzL1HHs3cogrrdYOkTwDoK+KIfHpZ50jBbOV7B6zkloWyk9eiHXDSrfS70W
eeegpbCcq8FxDkK78dpOjxgWVIhadvefz0RZT1dv1zOjsGRZR45aTht88z21
4uvxHyA68n1GctrAmL+OlppuVQT73vJw/KgYgIM741RXbRqEDvUfKg4OdRC+
yKzfS6cebHWjzMzzumBo+734Eb8EKHCpbpbfwkNMuEfgH/cIsO+RurtRtY6v
OiG/75j8msv5cORa+7026ETtWd1uKit64YG7zqLDt+uh68s5TWvtXHwndGP+
yQsD4KDAAyn8HXDgce2V9xtywMrR2ZRbcRiYz08W9pzshss+lVbd8TUQ2uD0
jqE8Cokhv97N+N4H/rNbbeVja8Cls+nBhdN9oCd5ewdPbT+I9axty1NjjflR
G6dNagxCT4ef9fjzAHzjWL/h3fICaBvPGyBA80bPHKPCMEMG0VnHnFuybhA0
Dkc8Vmr3RkG8tXqVAwsauMSULkaNQmiS0pYB7Tx4esNdVt+eBRqfo7Zu9hgF
7cnmVsMK5VgSZL/LLWIIXN+mOGmTXfDhycK2rohqVE07kObc0Q/k2pgrAkuY
YPdwy5OPM/7grcP8q1Sie2C2gHq3BnZBgWe1Zc7JDjwg98w21qsHjuQYPlF8
UQuWK++9kOoohfNePk8727ugP/zLgJBeG/QcWn+9rqEGie8F7ra9LFDibz1Y
uzsbLIPvf0hO6YZZLCOrtOKx+Yp8dJ3pVAnWoy58WxqDwLf+5ZTrb5kQ5qyn
VTMYCZ+msFo9TrdDl6d5lN0DJpy+qLnqsU0xuN4uy7vFVQkKjLLAS9M6IOFE
T/JR+2ZQTNyjeWFPFzjsUxEof9YBRTy+x6VaO0F5Mver72PvL/1Voa3CcAcI
ZQhpMoXH/L525FmZgUFY1Mz1vnwLCzJ2vw5dUjUEsna+fZfO9kDgeO6FOTT3
6mhtSRy6xIJH+UXvUtXHeOo3b8lwaz10Rj0W7fLIBOfbP120+IfgMp/ghRKp
QXiW3Dpd43kaTLIy/Bh5eoxvysqOzwmvgCBu41eLedIhdY7RymXbh+HGpq1n
Cs7Xo5Kv/1HHtVkQcr9l8pH3g7BcoTZwYeMwNuUcTWuz6wBrU1vPqg4WPD+W
xnC4OYBKihL2zxgd0PrlVu+FtY3Ayqv1e+jBwmdPl+0oOMuEOxcePi59nQe6
B+7LLA4uwPfZO6ZLhDdBisRTS6kLhaCZbnY3/RMnsblqYAsotgC5oDqxYXEc
qO2R553WzUMsSNWqdz9XC7q3c1JF7BrBqCd10oGrfMSvrwpznWOrQZ91vv9a
YQUcktGs7PMUJGxHLoSN/KyAGSW/0h7YsuDETU/fTcqTiJnFb+rNBAqAKy5s
012HZlCB2GnP9ggTrg4DbjnOhcAIkTY+rdkL02bZnKocu79vdHJcuqQdBw9C
d0ULHmGBbJRn6Mo1QsR203UDhYHZkLfskfzC+AE4cU1h2YVrHMTQkNucHxlZ
GH/A+bLJ9w4Qj+GoCRQXJHK8+aMuQAnqfD8nsym4E7J+zujtPMZB5Fy+d3zt
6Qp8Fimtxx/YBY5+T1Qv7q4F6WOKw+IrKnDyoi1bzr/qgtN6j6c/aWzFDQ+U
S2UHinCl+UhZ/OlO4B5QiFsFXKTTppPChnEZGK9QLfhsbQeUlq54te7IZHJ2
llBErl4dvLw+6KBjXg+mN48fvV/KS44YvJfZ498MQbuk4pW6c2BhYOKjw54c
ZDTt75QofwextL9jXNnyZundWpDaLjy5rDYTkocV73teESYOVg9ufz2/CdRm
yL3QTP4JK80aUvM9RIgjrDlrEwTywO74ycSnrZXAaWeRLqguRmjWehoSdUWQ
F6bQ796YBWV686tjj4sTfivF80bGeM7zcVN4fl8jjOiqzni7TYxQ0vZctm3Z
J7BQrtaRrymD+Ut3DCzcKEFs7dJpL6tIxDJOX4+Pe5hgMu1SVO9rYUKK4wLP
nb0FeOj50Jaty3Ohiit8n2+3OPFD97X74p4avOuywjvJuhAyPua5zPERJrZe
NmPajjShuObveW77SuCIss5zwcdcxJ+Y93eEOlvRqMpMije6FOY47rebUl8P
ikZbTC59bsWlSR1NMl6lYNW1bd7DWCY65C6usVBtQj2RBVt8WcVwfWpaUrov
D7m9XO/dLd9qXBeQfpcnowA6vjMOqE0RJcfe/dvs0hxUUG+8NdUkBzxfl6Yu
05cg+Uf7DYWO1MNnx7Qlwxf/gJP0+aoIByFyu9YGiQHDFtBbYLh198svYPxA
Xq3zijC5x9du7sLXVVAgIOlGhP+CzHN6gjsfiRPbr5VeiDnX+Pd7gE8sl3aD
B2KEO0eVgprRmM4yPSRCS2zwqAfnXO7fY+9t5qY04lMsVMzc+2ZfQhI8cFnO
PaotSegfXBB97HDe3+8xggVcYwkhSeLDZ8dMqWspKCksEKHyOWtsjrt3yHVJ
EMzkHdy54t/+fs9QGeHm5R8lScz8UXt3b1UFhs5Ks9W/kg/hRduf+p4QI6bS
9TKl6oXsen2n6/WUqhem0/XaTNdLlqoXGtP1SqHrdZqqF7LrpUDXS4eqF7Lr
dYuulyFVL2TXS2lF8ieu4DLcxwgvmYd5oKcfWLxvRIzUpeu4naojMuk62q2R
ycnPiUGbg4GZsWQmqD8/tcNCWZLcStdXg6ovsuu7QjvDQMPmD/gdGXII8k+A
htzjZkekJUk7P7nOSa5Wf+ul9oGjNDdRkjRmZagdnlL897zm9a8TkhES5FJm
9WTJsfuFscCneP06N4welIlQPS9ETg6vtF7V2vz35735vzYXD4qSSXQeqEfl
gcjOAyf0A7D7QZLuhxaqH/A+3Q9edN1lqLoju+5qc0uOprZn4Nlb8/fPeFuD
VTF6Dd8FxIm+8XqIbD1MpfVwI6WHyNbDdFoP3Sg9xH/oIbL1cC2th6soPURO
Wg9V6DoaUHXEHXQdJ+gksnXSka7jVaqOuISuY6LElI2Vt1NQNcczPyWzGn1n
/OxbYyY+sY7IruNUZbe21/HhcH9A9cKC7hKcl2FpY7legoy7N7NCq78GrOtF
ZCMWBaNuzYGNI9PFyTzbsxc1Z5bDHcOSA9NjM7CgN+PKBjHxifUFdn0n6ACy
dWD0MqXPMyl9xkRan7fQ92ksdZ/iQvo+ZTCSfwjY84DXU2+91kk9qJCtcP/T
An6CP5riBxeKH/ACzQ/RND9MpvgBJf9vfkA2P+jQ/PCL4gdcTPNDrQrPzqvT
YpFTVdH/rk87knJ3uJa4ipCO43NyZOfk17V8B4LzQmGdQe1J8+8NmL1222pD
KzGywX0ZzDuXBIeVOjYc9OvAd0cGHi4LESJfW9yx+HW5EGq+y719fqECl8c9
qr86JEreGZn8/ev9IliVc0L/QWojLtuuetp8tjApROutH6W36EjrbRtR861f
sRa+lJYY5pWX4tm62TNN8gTIHc7Hrxx73QqysRZrUwzDUI/H8/DiFG7Sl+aZ
VIpnUIPmmYU0t2ymuAUNaW4xPlpYNXy/DpZI9F0V6R3GQwP+WSZXe0Cc5up6
iqsxiObq5+O5GtlcHR1ycq7HvRxYMfMaHtEZwNdGvvemfOYk/cfvWZCT3rOo
7igPqY0tAs+0BQsJux7cHbNCVj2cj3S/ta91SWg9uM2+bv159QAyO17pa7tx
kLmHfF3f21bCszvMuIDHY/xdkjPwcM4kch7NFe4UV6AhzRUchhRXBFNcgQo0
V7D52YLiZ7Sn+XnC3g3ZezdFmqvjKa5GGZqrJ+zjkL2Puz+et1GG5m1XnfMm
B3S64fFL5iOpBUwMIkMH8vkboGiqT2Apfy/IlhvlrnTrRjWN0No4fha83XSg
/eSxfmBNVxwRLRjCvF2fU74zmmCCn0W2n7UT1FfddH7MR+xWmtrUPII3/FTu
AaMFHcf7WWT72YL2A2+igQWZp7ROLdw3hPf8Ig98SOxFDnAN3qTTD6km2+7u
PzuEoQzLpDbTFnwk7WexrHWMt7uPxk/Z0Idxod3FHTuGsG08h+N9msNVaA5v
oTgca2gOdx3P4TiZ5vAJ/hTZ/tR0/P4Xfej97wTfimzf6kjvheuovTBOovfC
E/wssv3sKnpffI3aF6MFvS+e4HOR7XOD6T1yMbVHRi96j9y4SJOp48NB/LBS
61HT6sTYiNJ5N7WzJ/pWZPtWxUTZp33ruYiwKlcFzps9qKes4are/RviL4SK
nJzFIOR4fkjMPD2A9p6qJ4+LluBbv+mvZFdyEUGGFh8NT/Tg1W1N/t6HYjBL
MOeSffkIbOB5ZWW2sg+/PQoe5bapnehz8R8+F9k+d4KfRbafneBbke1b17Qf
ZPCe6YDEFIM/ajw/saq76TGvCzcRS+dISlSOhOwcaUK+iux8dUKOiuwctWB8
XorsvHTT+FwU2bnohPwT/5F/Ijv/NKXzIhsqL8ICAyovek3nQtJULoTOdC40
crKiaLssN1Fo2Kp5zbcViZvuv2Ruf5uYQyI7hxzp8PJafpSTMNBY4cZ/ohFD
X2SuevMubWLeiOy88euhwilWMxhEafcng+PZ1VhUsLe/6UfTxLwR2XkjO0fK
onIk1KdzpAl5I7Lzxst0jsTY8d8cCdk50oS8Edl5I7P/luVkZgeQix1UP4u9
w9XNU+X/HJw8MT8Hdn4+IScHdk4+ITdGdm48IWcGds48IR9Gdj48IU8Gdp7M
XM8b9MCGj1h374WHyvwY9OqUvmRe/Hpibgzs3FjjnEGC1wJewjFYo3OScxTa
pry3in+cNDE3BnZufHV8jgrsHPXa+NwV2Lkr6/qGaR/O8BL8JvMNs+qiwG/k
iltfbQxOeE5gP+d3JQflNUc4Cf55gZ+OVRTB/MH0h/qrG8Fwv9BPi/mjoD88
4mvonQ1uRlo3bBX74RAyThUl9EFJnLKlQXg8ZK2QEL6sxz2xXsCu1zxaZyIp
nQGC1hkvWk9CKD0BS1pPHkZfc13t2AMyS4Y/bJAsAoeLHw+1BvYD0PPeQM07
lNHz7knfU13UPQUl9D2Vt0vpd0tpP9whjpgMbh4BWxEdvUi/DIi0Fe0OyOqH
w60WOy4pjoB4g+5+jQ35OHYPbnKyH7v3HQ6rf6obgO3ZUUkLmxmE4Xg+ATaf
HL2hsSRYrA4CX2RcO3RxFFLXRKdb9hWDpewkT55ZdcBt1Cd6c/co/JreuZf7
RTUGyWRKHvhcBzmNnRVHHg4D+H9ZvzSuDxluSf/lzzcUfwKbP+3O5n99JB+B
WW1fm+479cGCkelMqY2jcOD3S06tn3E4e91vwjihH+o6q2boaVbDvFkq/f6T
49BN4p0Hh38/bM7KmqHV04QNcpGnW5f/Qkvt24q9y/qgKDV415zNHCSpPGJ0
audLsPP6cNLvejd4+vqUS+hOIrfSXKdGcR3spblOnfY1pyhfA9W0r2kXMZyR
fKcCY11Gsr6HN8DWxTN5q3SECPtTt4oihepwXUfljOGHzWCtb9zyO4yTCA6+
FtQb3YAWd7sPGba3QCep9qGysQ5Mn1S/07VvwKgdMxaeSGuB575N6+6KM7FL
2nfOrB216H3S73e6SjMkdtScmfObm3RK4J23wa8MV8psTjCSbQDZ5KR5RiHC
ZAXtU5QpnwJsn5Kh+63cs6Dorx/X//4i6Mx6CcJy/u/pcRfr/p4v3/cQxaVF
CL3PbzR6P7X8PX/08YaYfD4XYWV+pVkk5H//HYRbTmajmFYDSMa1d1w69b/z
WRtGIzg4OjBgrdfFkdzmv+fbnyZVLerjIXeJPnniMrv277l/VnrgoT2i5FL9
kdQvG/P/9/ttgkVFvkqQ+uPzJWTnSy30+0+m3j+y3/91+v1vpd4/st9/IP3+
ran3jyz6/RvT7z+Rev/4jH7/HfT796PeP7Lf/0P6/WtR7x9l6PdvTPdtDtW3
qET3LUH3rRzVt1jzf/ctsvvWme5bK6pvkd23oeP7Ftl9e4qexx/UPGIsPY9W
4+cRf9LzSI6fR1xFz+PyH101blwdYLh61GVZfSK2KPL/1DdvQm3joxE3VnbC
me2aMfNTq5Hhfe5QxNIAKKZ15jalM8jWGZLWmUOUzuBUWmdW0PssU2qfhY/p
fVbavmudewx7oaA0IJf/ZCkG2zRLCcj34oaorEaJSUOw/3nJyyu8DZjjyZ8p
L1uN4hxP6h9PYhCyRSJr9wgykTFy+sEMdxuYsP9F9v53wp4X2XveCftcZO9z
J+xnkb2fnbBvRfa+dcJ+E9n7zeDxe0xk7zEn3Dv4j3sH2ffO4fH3DrLvnVGH
B5tbezgIa4/jmgt/hMLMRcMn5dxZONQbupR1l5swvnrhd4IIguhOoQNLyIqJ
9y+w799UEx2lm3O5CccnJ/PTpkbCmm2WepNeV8OViEse2y5xELZHSl8mXvkJ
pSE3XpZGdYEI/f7lqPcPDE7r/77/BprDgygOhxyaw9n8Vk3xG3jR/DZMc5oh
xWnA5rQAmsdqKB6DNJrHJvA5sPk8ezyHA5vD7433TfCH9k254/cvqEnvX36O
36Mhe482YV+G7H3ZhP0j1tD7Rzl6/zif2j+imh+1f9Sk97Z3qL0tfp6wt22h
9rYoSu9t99J7WyNqb4tqE/a2D6i9LTrTe9v59N42mtrbIsOV2ttO6BP8R58g
u08mcBSyOWpCn+A/+gTZfdJt0nDZeioH0XdWxH1DZCoUTGpXSxvowlgunakz
oriIvobf9eHPs2BP9CqrW9+r8eWqdkfVeXyEVN32330d+UDUpxz3+fQEA3Zv
5zf8ykcsrBLknppWAHZSze8feh8HOaaz/qy9vMQ3Dhm3m/55wBD7em/ywXDQ
inC0CL/JTWjKNa4zSsiGvtj5c0qul4Cx3WjpkPpYv/IEqGY8NIWMdxv5di0S
IaLDVCNrRjtBJ3bttQUqSfDki1SOZhoX+efIsvhaES5iY+wqi/wTHWCXPyl9
YU48zB5Wm240h5s4vpjn2NypnXBgyY8oy1mJMIen5PGrA9zEy4ZNUoKnOiH0
HNcZXcEY3Ju0WHhP2iBsLYloF33UD9J6YXUfuFvgHZhZzl04DD1c6ZyWDkNg
cdImftf+Ymzr9vmTz9cMSbdtCzEmExzETOz5znETywUOjCoq10DU9dZLR5Rr
oWZuzXDLl0mE6aKEydMzi8Hx0UGv+QuYYLLmeJFegSDhH+pn0HCvHoZUFFbu
9IqHp2Z+Vwo2ixD2uTna6QH5EJQdOYXZWQRCvyq4rjaLEbxDsScZc3yAX+TU
8uYXNVCSu2dJzzFxYlQo5fydmlpg2Ky63iXvAwbrcvmal4gTStpxpvh8zHfs
KbJPbTkAKtIWhiODImT0dvtZnraxsPTy9qqm/Gwk9jarJWRLkLZqz8becyVY
zrN/uu1TLIZmJvBdWSlOflsb7vTEfez3tD75YHzwLD747vh4mqYw+bSrX5HF
GQVSfR36JpPasP+TSfem5aLk3nB1uYiZhRCal7H48ORaZJ1v1RKcJUoWuJ4j
PctrQWvBOtenLbkoFVT/WcRoMqnyyaN02vcuMPRY1tFxoxezN245v+9CL6z8
czx3Q9EIiObtdtVKY6Hd8NXVFrW5MPLSV3vzLwYRIv1z1r5pPfitABudfcrh
jLDf54A/o5BXkuZwirMLHQdOuMw6UowBqr1pUp+HwCzmRKzd2H1kvLCab4pc
AA7DSQedDTzEcpVjw3fNGpDxVbaiVvk6Rt4w9PFM6IFNLBPB4LZE5N8acsuO
m4fwf9NVGNTPRyx9vWNDQEQQMBgOAQd7mGvi7Xav61z6n/9OdsfKkEVrIXkb
3wzpWBFC4s6XUvOzHcAwWqDjJGyFyg2lP+9fFybXm8ksDTjZAcbrXpOy/U8g
2+hyaUOjEDnh9yP79x8wetVP3uIjsiKGJWfdKwDG0raNxKmzMOHPRfafuyvE
+r1keS/ISjQbnmSlwqO48n07P3CRr+51HX+ynZMwreZLeL+cCeQSSXIg7gfc
OrdO88l1HsLsxBoNoYwGYPA/upLS6gFcD939JvOO+dAzN/kWJDUCYdShu5cI
AInPVXKqn7mIHIlLNpEGLJC60XbbVDoFzuRa8j8x5CU03yU8aT3ZBNbH19ou
TPWD+FTNhu47vMSeSAFSPaQJZJM+FJVOdYLZYSJqzFPchLp52mf3R52gxAp5
Q5xDSDvsJrCtlIeocmke0breCIzn1x/cTPiEQ3S9VlL1AsYXql4LeYLTZ49w
EeteWId9KWNBCW/LqG9uEioE8ucedeUimrJSa57p1gL/83PLFm3NhE0+7/ab
ZHESGvt1PwrxjfnC6Q+eProYg6/uLXSYv3YYnvu0SEEiExKNJhe7X4yGFXQf
ilN9CM50H77ak/4pS5ZBuOxoV1/r2gWTW9ZcTdpeCkN0f/6i+hPC6f5U2nnI
ZXR0AGZcris9m94L/LrPDW8ebAWFmaOujks5iOYjBvtO7+gFUVScK2tUClKP
3Gs6nTmIoZn2Hy6O3ZtSeq1tM3XyQU+UN0HmwhAwji35OL13AKY6SZTs1WyA
864xPYoXOIiyuTp8kp29oF/nPftScCE++pHZJSvMQUSveBuj2tYDsWuVFs30
KMNDS1c6xTsyiF39A5uvu3fDUrx3osq0HK3oOSqi5gj86Dn6TM/RCWqOwI6e
o2TT1ctdxnzwc0PRxaHDXXA+72buUHYTRqW8X76jsg9Wb99e/fVqHVQ8XiYc
+L0aBtY+mTS4rxeMoywCjse3g/WJxiSdPcV4S3dE0kOuG94kvX9yL7Qf1qTx
o315DzC3Cvm323dD160dHsOvBuG1YnuQbUYXmP9cGz25vxVk5FqybJcNgMC5
xZfrxvzut2pjse47XdCQ051f+7kHflhMNrHb3ou8n/mM4hhdcMjBzi9lches
ELxuoJbZjQtMWSdmjtX16sC8HJ/pLaAsvXI3K7cFdzselR690AFOpeSOQZ1S
MDMBN2HPaPREpe5Zq/Pg2FMDi5CPLGA6eOiP3hUk3rsm2iquagH9lOlXdrOS
QbXzWfTR/VwkxB5YumxNDVTkZK3cfL4Q5rX+LJ69RojQyjgmzTknD9YsZ07/
/bEO5PZrPJRcL0qQbQN5OqdSwCw/UOzP8maQ4Bw6/OuEKLFnklaZ3IkkyN41
bbqffCX4vNbec8hcnPjS5mJyMcAZza2iS5a96gCz3geVjuuFCUt7O/fuc+F4
Maw8SCOuHppbBD2PF4oRLj66B09tagXnTylmzoJeMMV2z+dfp0WIeqn73UGz
W4Es7rvonuIAI9ZuZ5+7iRI5lUK/3K6Vg9mSdMNNownA3Le5xMNTnFjtK3LN
50wo6PnXwFuDXHhzTrnjbYME4ewyz+22VAiaSLyJ7fQrBD89bj7fZAni6/h7
BNj3yIQ/F9l/Lvs5vajnREn6OQvSPeytE17hJoOtN5W1ClCfucPzbI0EaR9h
XO+ytwHsslySXwt6YH6By4qd38RIxfH3ILLvwRW8n377NH7F0LDd9hW/6jBv
8Niur1LipHR70qSg7KdQyWSayw4x0ZarNvP1iDDJ7cuS2KbxE2wYz093n6tG
e66l+9RtxEn3LM5rzjcy4PKNWrWW4Ea0D8tr2GIhSq6JWyB6w7EYVhtIqwwZ
F6DLZp6aD+/FSJtGXf75ro1gq6vzYopyFK4esH93XkuYbKe5Io7iCrxNc8WB
k9GBuWIFsGV/6bLM/A58sHP4y0YUJCXlow9JZpWC+IbI2xUP2lB1ecazY78F
yeo9pYKvt9TCXX31l1zzqrHs7my9xgF+0ofuTwuqP1GF7k/eu7JqywM74PHT
xXa/hyswh2XAHXcxC3jtvT+o83SC/yxnA/u1bShjuEzpkRATVuQlLC5Z1AXl
SyxZH/q6sFOvdlvF9R5Y6Rzjbc/XBr4zOuxuXe5H+3foaa7IIF23OkXOjeuF
vEaVgfeZTGxp8Na3960AR7UHT4Mf98PvE9rPdmqM3b8VPbZYXYa9Ljy7cyYP
wDSz8gXTorrRtNTQu4DVDK9pfX5G6TP60frM1mFnSoeRQevwPlonR0z/q5Mo
SuvkQnupe/Jrh6DAzs7y5Z0BbLP/eI0V04j6tB5up/QQZ9J6OEcid1P1nhGI
+0bGHZRnobWTvOAvtSycwJ/I5k/fVCmbo01cxLvWl0vKG2rx20BwT860ZHSg
71lD6p5FO/qenXAPIvselBvPq8jmVfb9WEbdj0jQ9+PW6z/PCfJwEj/yG0ta
7rbjmGB0BNS+wiscCwYe8neCbnumo+20GJRZFi1zU46HmMAhyOaQCdwCbG6Z
wOfI5nO7hIJ6w//k28brOh7hDpilZdewaaEoOffqxVPuZ8Z+T5Flgc3bq8Cw
d/PPthQmA2t5dtmXchPufXqvbh7IAUbC0eZI/3zk1OSTmpLIQ5yJ1BHW3N8I
5JMnSSsCPsJ10Y2QncRHtKjP+Ki7MAcZa9+tfUUYQObePIsZXnzEzXajps7g
Md9hztvLumEDmjH4KduWlzj/Y5jUcB/jEOOA+tYqL2hYJiEj+JKXkMpNX5DZ
GYheMd/WTemNAj9HIWNBex7CJEwjLlg/F4zlba5uO5EOgm3KJff3cxEaF789
TGnPhNSv/d6JpXXA5LWLLtjMQRw4lsO1RrISPO1PXfrxsAlsDR20FMb81BbL
o0/k3VrAbn1RzNqT8bBl0DP6pSuD0Pd7lRenkALz7NfsstPuhU73RRyqUcPw
8ln30NF3tWC27NVUbY9GYPOnPsWfwEPzp6VlsryZZx/MixU02jO7GLZEWUut
/jAAhcuCHu/5Dweq9olIhBrjmafzy9abihDTYYubT8bYc9zmGE6Rf4Z96vdf
ul4VJczO7KrN/c/PH1pwuNtQD/aMXLxhGChCTOgTYPfJtN1K89NNuIj90x+f
1W3qAOetP++u6E+C7rdt1Wm63IRDMfefK4c6oeX1gltFNrGg9u36srhmLiLv
8odfqVEssLXnsW03SwbJIbPHAWN+fLW0f7FlWi9MjpqR8SeiCBwtYh7bqnET
lRlwo3ZhJyT2bXDRkUhA1ff51XfWcxAWcj9cL17uBcbtbUca3pTgq6Mvubve
dIHwstzBVZb5MJCf4hIn0gdKbxccO/qyC9YQPFLTG6vBOHiWS/WZLNi4/bDY
Uc8O0F5Sk5rIlwoVX8/knh5kwdKYgcgMhSGoXzjnwqttA+Bc7DR7l1kTBLTI
FbFODoJpfP2LWb/6YKDthvFkrlYIP2ugw7m3G9pDGEJdQQPQkxcbN5unB2K4
bp8jegpAaXOhA471u9mjw6Mn0gWJ87TOn6d0Hgha5++XNtyEaB88b/Z8pHx5
HQjs3CorvViccIyPPPwmIwKUet41BdVVwb7fmbusbcUJ8hxTr0bdDFq2r7mr
9DAfDs72TedrliAIXh8tZ+nWv98RTXZ6P7U/SJTY/H5ywVf3sblTiXP4aOWI
i80cbn+cJUyqC7DOdbiNzam4dL4jaQazePycdXeJkJ9oH8pY+F8fis60D3Wa
k7J1+X/+/z9zLe694jkBH78cwh9jfWVx1FLrzO9vEDlQOiNZIg8rbM4cndsq
QS4d//uR/fvrYifxiJulQo05LL9C9KDZkXdR7QX8pGxTlHNGbiJEWWxLPNjc
jCsW3RuCI6IkY0WLgJL0D+Aw3VQ+m4+JWdz3Nvbli5Cr9m4zMZJjgY33luX+
l7Mw8vaMJNtzTGwuScz7MaUbLFq0p23gqkO3hKSUvsUxKPv7ttlzw25I9J3d
+XNxMSoMSaz2v9mNmXVVeg3rhmFBYq50cMgQruD4Kq4ulQcTcgZk5wy6M3Le
G97th728ndbWQ6W4DWx12mb04Q7aZ0lSPgudaZ9VUMoR4PF4BJ4uesm9rLEO
ZxLrt9u/rEXuQlseTe4RyDugxpEZV4hmunIDV5p6UUFZR+7as1FYtCZClSus
E0l3daf9a4tgd2fm5VALDqLg8b3nm7xbkcGxO9ZE6RceTX1ybOZPDuKL1fbn
CqWVGL9PlmiorMNHr/1VX1gziPOLR9rAqBt15lga5Uwpn5jzIDvniRuf8yA7
55mgz8jW50fne19np/EQH67qCszNy0U77/YXWz4lYt73jYl1T3iJy+mX5voI
5CPzyJK5I9OD0G18XoTsvEhPdapPVx4v4X6k32+Jzg+oMBh6zMnzGyfkQsjO
hSboMLJ1eIKuIltX1fZLfJKw+s+9w7twp9ZDmMUImJ4uI0xmS91c3Qpj/bl6
hpnw0FlI/3T5Y0qcCDlv/D2F7Htqp0Pa4gdRLcCw/ZkVXmCNw/matyvj/z/9
RLZ+7i/T3zHjMjfho/skBRQrUA3efZy9OxtDreI8N1TyEs3+fiXKjvlonM16
J37nPU54TmQ/5/t7eRo/xvTN6dlOD36OKFgoejJ1VwAfmTP++ZH9/P5Xz7hO
5x6EgPfNe63W5kAztPoVG42iekgpGWTBILwOD3tWd5YB/8r0dnP1DtzxfvDn
QBUnIfmpqf4udw1YB55LSJ5ePLG+wK7v2yuctjVJ3MRaHxWzvFtjOmEU8Ge7
wmdkJdpN3fyHj2jamb2e83YQMlY77AkMuozrjN2XeTH5iJM3B1Vrvo+dawZc
4VZbD9ZaelmFJB+xyGW1T86cHLRbnbj8uKAZzq011DGu4yFMjxQHFY5UApGp
/2pZ/1coeHnpivQcbmKQK6xQwKwcZNNvnj/SlA/MQ9yHDfg4CXedtdVbQ9NB
ncdY6kk+E2b2yOZV/OIgMozJrQKLCsH/7NFtSbFtEC+86bEz9zBYd946F9qe
D2liecVRRgMwob7Aru/qTXM+RirzELZt3LFTljYAw6nio0ztBWyZe3R03Qpe
4o/Aej87+SYgjuvNC/D4Auz5XUzNL2TQ8+uzgBz+NMAg3FXP32gM7wEz6X6D
iknloHjjo4dLUiXcP1+fsR/GfO7HmednjXAR5upCmwKyfkLcnqMNPy3awX9n
nV05nyjhlaezU2CJF4rkLto0JacCsjrk3VgqEsTxs1+qaozSkbN8UmGddxOU
3KtOuB4lSmzb/UsvsaUQ5T88MtNxqwDLhRzCJ6vFCA2hud8uBZSj3/X1Siv9
KmGKpeejjjhRsiXW8K7b2STUuTxNfNXOEggfObptubkE+Vf/5/1X/9Gf1n9l
yXtvMu2yIOOxfB+ZnAHXoDVtn50EYUF0vs4/kwiRZl8TlszMhln707e5xkoQ
Be/VeZhf0qGv0dphZ3oImKlVtznoSRIl9Llo03/PkX2O5T9250IG/nA8VrfP
sxSOc+6emScjQZwVjV46Vy0GC87f8nQ6FQUPiocKxOwlifPjz5F9Hnr1s81F
tXLUfV41O+J+PBSBwUO7sfu0h9vRYPm7KpxcdklNaSgB5vtICZZvEyO76XMl
6hyV6XOltIDiucszkfQyH3RdFQNvObbo6UyVnHiO7HPGIvHFyw+Fg/Wi3zZW
XWHAnfxzdqu7JDmkOu4c2edb9r97Iq5aCc4Gm667SHxGybRSmdWSEmRLnuMC
rv/4VvovX+lJsubSoiTJKooJjX8PTC9rlzSeCDSQq2GGfpEkDPU9OW+xctBp
10bImRWLjhLTrqnGSxDz6bp/o+qO7Lqnl1lfKk3Kwcsle9ec3N+C8qJGkhpX
RMhnoUe10sRDsCA+P/SiDYlMWctzCV6SZBfdJ1uoPsHfdJ/seSOzbuuvKCw9
/VFE06UK5x/m87XMEScfDzsomlnngOPcG/2+b4NQ8Us+Z6+EJDmBT4DNJ2qN
p7f6htwDj6+GH2Jn12KA/4NdYQfEyc5pEjw787NgIDzFmbmwC0PTPo7WJwqQ
PdWKBaURreDfrjD0fskwVil8WfU4tw/v9TYWzC5vBec7ppwfbw/i6vKoWpcd
w/iO1tXblK4iW1c/0/rpR+knttP6uZTWT1dKP1GA1s83g1znVmYwCG2G5CGB
rDy8zDGrb8PpTtxO66oopatoR+vqoPH+65uPcxHJxxV9dt4uRZ1NXTqOqpXo
Q+vqCkpXscKQ0tUWp0uS80x5iF+d3AVa+ypRz0C+IP3C2POWOg1WOHMTAxYP
8p7HdeIKV13p4w5e0LTo+cXwe2P3dVl7zMExfninKjxy3+kbTti/IHv/ws5F
K6hcFGXpXHSC3iJbbztov7Of8jvI9jsTdBjZOpw4Xm+RrbcTfCKyfeIEfwds
f+cUG+CiHMhHHNhsfiXFOAjsVt/l1uB7AFLjfQ2wfc0EHwRsH9SSbhZSV9oD
xn3odIbnLRR4HfXXseInV8XOZgqMdME3bZB+8+wrFPaNNHP2TSJvRu5u+znm
y0SjtdzjJ3+H3msBmZr7GOSBnz1hpTWDoGRoeHlLejSMMOMXMXM4yLPXr6jO
mspH5PzK/DTFsRoq3NZ5ZX9yxVR6b55P7c3hC70316D3+HrUHh/y6D0+Qe/x
z1J7fKig9/js72OfUN/HQjj9fWzu+O9jYTn9fezQGuo71TTqO1WIpb9Tvf6t
fbeiZh+EBMcGmqztB43hUwYqLiz8efbe5dCBNnjVQoqPlCUBr7SrYfFpJkgK
Hq9TGLsfy75H6Cs8qYEQnVcMjdnD4PrlYUdTYBMorX7N05fChMWsnlKTt5zE
rzmcAzbNDXB/2bnuqad7QGoR/4VTt7iIS+O/v4W79Pe3cpdmHjV2aQXN8++D
W01H4cJFtYU5fxowi/7+Nov6/hac6e9vO8fPO1TT8z7hu2hgfxe9RHNRiMeS
RtDLdSh897AbWHH/r64zj6ppf//4aZ7nEJJknuMackvPRrkKyZQmSSjzGEWG
EDIWEqVIQtIkSRLPLpVMadBAg+ZOaR6USn6t+/ls33XOb91/99prH3bPeab3
6/05ziktN8XYEMpjzyA8NtyiPPaNQ8napxJfgIH2yi/S7j9BJOWd7dZFoqwG
5bEDCI8NlpTH5gvONcDNNe8pj32F8NhgR3lsIR4eOB7eVZB7hzmUexefTDj8
X5P/5fAhmXL4JTRv7yN5G7i87U45fCAcPnAcfhDdc8aSPSfk0z1nEOXwywiH
D3qUwx9H/REexB8BnD8iTdAfAZw/Qqh+AVe/GOoL20t8Ycj5woZSn1EK8Rkh
5zMypX3LRNK34Fbat2ykPrWjxKeGnE+N8xXaEV8hWlJfIUt9asHEp4a/qE9N
i/oKJYmvECdQX6ER9amVEp8a2lOfmqegrxA5X+HCzcteBw4twFcz05y/Zrdi
dobPhvjbkkzKlIC1nZ/SMbri0ckjEe3IJE7VWRMnxbjWdKXX67dAqdfdfY+/
RcPcgB/fIsIl2SDnqBz1+P58GBxV43vgJL5a6hwc+F2CPUp9kVrEF4kO1Bd5
g37vxpPvHU6h37sQ6osUI75IVKO+yET6fTxLvo/IfR9XCfoikfNFFlHO8Drh
DHEp5Qzlp4phodR3uGIw4sTDI03oHpEUcH20KGu7oCpo3/AGeGd54YG3fjka
TladkjOnB+uOHTrx5VYzhJpNENflxUNE1qT7V3U7MDNTQurA2iYICFR4PVCD
xdP1wxV8ZetxOvXVZhBfLQ6jvtqr1B9tRvzRqED90Yupr9aL+GrRgvpqnx0c
tFlX7gfcGJwXsse4D51dlKI+NNWCO817z0jewxk0752m/mVN4l/GTOpfPkB9
ym3Ep4zl1Kc8nPqR1xE/MgZRP7Lqy27bua+aIEb+9R7jyufwOeq01g6/0j8+
5QDiU8Yk6lM+T336+cSnj+HUp/+a+ujnEh89qlMfvTX1v9cR/zt+oP73S5RH
1SE8Kp6lPCpjssBCe2kjuC6yNbbQiYOHg/ON09eIMpqbJh8q0G4G96C1rVPe
e6BV4cszy4LEGKHzJZA7X0KoriFX134fX7q//O8GMM/bG3RI4xG8jlSpWVEh
z7w+579zjG7/8++3JS3UOYB1bIaNRom8cP1Frv4u1To9O3ZAB7gvOTRRavM5
2L5JUaXmnCx7QrD+Ild/Nxr5rUz/pxtKzdp+By9MAOnrj8Sny4iz7scijZfs
lWAKOjZZak/NAunVEVHLLbJxqmnj7pXDJBiHf0YrnE3OBGnfCgv2YB44XgvI
/5UqymyWcQmfffY9nLNMbuwKqIf9VoEys+70gnbad5nNZlG4dcFdEcswEaZk
18+xt17+BssjwbUxw1KADbT1bhrbA556A1ynPOgA+3OpNockA7BOToePWTLM
yqvyx9SH/QR2varE5MXPYOPYrSnbt0kwQv9f5P6/ngNW3tt3X4rxSDl6vhcL
gHn+fv5vbWcUej5wz8+8bGW67G0b6H4MiAtpewAaGvtFOlbJMYX/iOxst2wC
Nkrug5LnWZyYoGMZESnHHv5yhBd7VJLROlOdI/6gDnTjonxvJ97EwKGPuqqN
JZmOTXfvXzWsg8zrUwatdQ1HVZv0Z5f8mqA0Z6zH3EIfdOrrmfpluAjL7dM+
kX0a6NB92j+UA7QhHCB8oBzgKarvJBJ9B/jfiL6juXNpFK+rC9619Nwq3dIG
1r62q3Se1MDi1/rRERd64NbFg3N//O6C5ti4ZQqv63HGEU/TE839/54BuSb3
Nt3A4s9WDr3L8lGK6laXiW4FRUK6VTjRrUCb6lZbtp0p8Q9q689jo9dqT2yF
BtnihRl9baDYamGdltEIk7RX5Lt+boSYTM1DT/r7C07nKiU6F3RQnYvjVdYR
XgU+U14lKKszqdSyEQJWsuDR/570b10dqf6Zx6h4DEwM+tkBBml5SocbfkHM
gpLlQZMbMNZ1Z7bJxg648fZJm5tID8jzbf3XPepAT9NxHUdnNUIX31hHdmgn
9E3Yfy47j8dK0jzsQ/IwnBLKw2kkD8N8mofVBPMSxtK8tGVkhY69bDNoB+3Y
0dBzFq4MEM86oM+HDJqfw0l+hls0P/OdvE9sG90M0Q/kbh0P9MUtrlP2HY7p
hN2Dp4/RN/sOF5JU2h5FlILaMOdHAVE85vBQ9deSITVQ+/ng7NaYBrDdmn/t
mqYEE/7w7q3vHVWQ4fj6kUduC9hk7HDKmSTJyA8zfJosVgExg/ZazlbpBIU+
j1EXgySY+ccV12aolcHyslH6Bie64UrJ3xu7douySwKvDFrvWgG9jm8iFeb/
gNA0j2um3RLsB/XjJxiPasgzW5oaUdAMGlk7XojXSrBfx8qJuBvwQd6yXsbg
2neo1jGZYOwtzhY7NfKr2r/D1idZ8ioviuB+yCaHj/1zH5/6CqOJrxAOUF8h
V9dCSF3DRFrXdHwPPFpU0QTMIpeRR4+fBQmjNXPa2BZ8p7BYfLtpfz60HnA7
7IMXmg0oXODYIsXwqL78jejLcIjqy95UR3YlOjKcpDryV8pNzSLcFAym3FQQ
1cEZooMDp4OvoH5zPvGbI+c35ziEKsIhoA7lECwpV5BHuAIMplyBJeUHDhF+
AOsoP8DFA0viARxoPHBzRDCZI5BH5wg9ysWlEC4OyykXN5DOF2VkvsDHdL6w
35iaFtRUCpWL7vq01Neh/2gn6YUz5P7EFZ/EFXJxFUR5jI2Ex0COxzBgJ1kZ
XXwLHid69F5X/ECt1d89O35KMUJzAXJzQQ6Nn1wSP6hG42d3ov/NZTeboD0y
XXTmms/4ZGTKufZuPnB5SXPgv3kJxHNJXppvd0bk2qJmaO+bbRgTH4cr8nee
DIEK9JFNmDhMuQnm8S+Kt/zNx7e3J7g8S+gCVcH8g+E0//iUnXqx43IjrLt+
ebP9yVbsC0xzZYbxmBDBPIN6NM9spzzMMMLDoDTlYRZ99yoxiW2AY8kFTm6R
v9DM6ciBEPFOeCqYf5DLP140//wg+Qe7af7JSjDb77KtETyiClZO3dmO8W4v
41sieGws5W2qCW+DjyhvE2g08WXCvUZwnJgdyj/cghnsomrs/34pn7/QIzG2
CazF7I9drKvCuB2bt11q+IGBVg9WBz9pAu86tQvBXzNQ8tGGr/sNKnE85aaG
EG4KRSk3ZU25TVPCbaI65TaFOCi0pxyUGa0v/qS+YDmtL4NmvG8+87gb/u4c
fsuvfz7nbbE+fUe3ATMo13SVcE24kXJNilGn1+lc7YJMMfmZ57xa0ZzlidcH
VmEv5ZocCNeEmynXxOmSc4guiReoLrmI6pJmRJfEcqpLCvF1yPF1Qlwrclzr
UVrfVUl9R92npL4LcXHIcXG8NUr6bm5iTJrNqz1aKq3Y1RM01/HCOxTSl5HT
l4X2JMjtSYS4ROS4RCEeGDkeWKjPQa7PeapmGz58ixQzI3VXxW+7/uesskku
DViAWYJ9DnJ9zoqB11bNK2wD3vltU6SkPODMYcXMB3w5RmVW1NVc8xbItB5u
kmF9F26mbvP/x09OuG9Erm/89F1yYueDLuDdTun9ddAX1+9cZXNntySruetN
kcGwXuAlly+Jcb6NaW2D7e7YibLRDT4no451AX9EqMH+sgQIOB1Z6FgmwT6t
+WYkYdgF3ooPfQ+uSQUDBTVXfqc4ez0ufo3FJh6z6qie8vn9qagTteP+2IQu
PKM3/Ezzlt9gr+k4UEosBC22DL0a09eL0XXxsbN0fvc/Z2B7YeMbqB7BDv31
swdFqG5YQnRD4HRDnX3dMyNniTHqvr/r5j78gJ5b9rntlatFmUTrUwpRogxT
ujm08OUj9PW7NLfseBPOzEhcZ3FElJHvGqX3zDcTrJV0d4x904AbqZ4YTfRE
eEX1RKG+Grm+2theqcptiATjWcM/OVLsMfLi3APurPqKHPcSSrgX8KbcS66g
PgicPtioP7h2fJ0UozI0wdH9YSzy5lcPVmi3wPbQU9KHH0sxKZXjx+/KKgB2
7GwDi6o1KLQPRG4fKHnQzEC6QIrZA9KDf1zuf05Dr838pD0g1P8j1/+/PRjx
WjlZgsmU/VjytOgxyntaLihWyYajlYulxV6KMwWP1HzTz0ajxtxEFzWzbzBg
6Va9NyfEGdu/PKLURGqg+ejd+smyz4XnCPyPOQL+Y45Abo4QmheQmxdGae39