-
Notifications
You must be signed in to change notification settings - Fork 0
/
sling-slow.nb
1733 lines (1705 loc) · 90.5 KB
/
sling-slow.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 9.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 157, 7]
NotebookDataLength[ 92523, 1723]
NotebookOptionsPosition[ 91467, 1683]
NotebookOutlinePosition[ 91809, 1698]
CellTagsIndexPosition[ 91766, 1695]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
RowBox[{
RowBox[{"f", "[", "h_", "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"2",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "1.1723328*^18"}], "/",
RowBox[{"(",
RowBox[{"l", "+", "h"}], ")"}]}], "+",
RowBox[{"1.1723328*^18", "/", "h"}]}], ")"}]}], "]"}], "/.",
RowBox[{"l", "->", "13599840256"}]}], ")"}], "-",
RowBox[{"Sqrt", "[",
RowBox[{"1.1723328*^18", "/", "h"}], "]"}]}]}]], "Input",
CellChangeTimes->{{3.680585829094075*^9, 3.6805858351312037`*^9}, {
3.6805858782596283`*^9, 3.6805859357091427`*^9}, {3.680586064676009*^9,
3.6805860730090456`*^9}, {3.680586134585082*^9, 3.6805861394061165`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"f", "[", "68773560320", "]"}]], "Input",
CellChangeTimes->{{3.6805860755752125`*^9, 3.6805860910582876`*^9}, {
3.6805861701510763`*^9, 3.6805861767066126`*^9}, {3.6805862085174775`*^9,
3.6805862121335287`*^9}, {3.6805862482940664`*^9, 3.6805862530120745`*^9}}],
Cell[BoxData[
RowBox[{"-", "1756.230095731537`"}]], "Output",
CellChangeTimes->{{3.680586076597076*^9, 3.680586091340562*^9},
3.680586143690735*^9, 3.6805861772280917`*^9, 3.680586212500877*^9,
3.6805862532732687`*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"vel", "[",
RowBox[{
"in_", ",", "th_", ",", "satevel_", ",", "GM_", ",", "close_", ",",
"peri_"}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{"ink", ",", "z", ",", "result", ",", "result2", ",", "varth"}],
"}"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"ink", "=",
RowBox[{"Sqrt", "[",
RowBox[{"in", " ", "+",
RowBox[{"GM", "/", "peri"}]}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"z", "=",
RowBox[{
RowBox[{"satevel", " ",
RowBox[{"Sin", "[", "th", "]"}]}], "-",
RowBox[{"I",
RowBox[{"(",
RowBox[{"ink", "+",
RowBox[{"satevel", " ",
RowBox[{"Cos", "[", "th", "]"}]}]}], ")"}]}]}]}], ";",
RowBox[{"varth", " ", "=",
RowBox[{"2",
RowBox[{"ArcSin", "[",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{
RowBox[{
RowBox[{"ink", "^", "2"}], "/", "GM"}], "*", "close"}]}],
")"}]}], "]"}]}]}], ";",
RowBox[{"result", "=",
RowBox[{
RowBox[{"Abs", "[", "z", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"Arg", "[", "z", "]"}], "+", "varth"}], "]"}], "+",
RowBox[{"I", " ",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"Arg", "[", "z", "]"}], "+", "varth"}], "]"}]}]}],
")"}]}]}], ";", "\[IndentingNewLine]",
RowBox[{"result2", "=",
RowBox[{
RowBox[{"Abs", "[", "z", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{"Cos", "[",
RowBox[{
RowBox[{"Arg", "[", "z", "]"}], "-", "varth"}], "]"}], "+",
RowBox[{"I", " ",
RowBox[{"Sin", "[",
RowBox[{
RowBox[{"Arg", "[", "z", "]"}], "-", "varth"}], "]"}]}]}],
")"}]}]}], ";",
RowBox[{"{",
RowBox[{
RowBox[{"result", "+",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "satevel"}], " ",
RowBox[{"Sin", "[", "th", "]"}]}], "+",
RowBox[{"satevel", " ",
RowBox[{"Cos", "[", "th", "]"}]}]}], ")"}]}], ",",
RowBox[{"result2", "+",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"-", "satevel"}], " ",
RowBox[{"Sin", "[", "th", "]"}]}], "+",
RowBox[{"satevel", " ",
RowBox[{"Cos", "[", "th", "]"}]}]}], ")"}]}]}], "}"}]}]}],
"]"}]}]], "Input",
CellChangeTimes->{{3.680587311169629*^9, 3.6805876296075478`*^9}, {
3.680587713929714*^9, 3.680587777346264*^9}, {3.6805878231169176`*^9,
3.680587826596884*^9}, {3.6805879450709343`*^9, 3.680587985010673*^9}, {
3.680588069272805*^9, 3.6805880877574263`*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"Abs", "[",
RowBox[{"vel", "[",
RowBox[{
"845", ",", "x", ",", "945", ",", "8.1717302*^12", ",", "21000", ",",
"14175000"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "\[Pi]"}], ",", "\[Pi]"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.680587635280297*^9, 3.680587673248394*^9}, {
3.680587879328382*^9, 3.680587889603097*^9}, {3.680588093073917*^9,
3.680588112530334*^9}, {3.6805933097739315`*^9, 3.6805933099471035`*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVz3s0lHkcBnDp/b1G1LhsK4pMrVIsWrVZ2fEtahG5tqlW885VZBDFNiFU
WMmkXLpRzqxoi8OKamr1bSstqzQ0uugiRSrbWrcyzLb2j+c85/nnc87DEcQE
iXV1dHRCJ/N/+3jP8Bjn6ONcmeLvhRdyYUVd0gDLSh8DzTLVvUW5YM95V/y5
hT52Zrx2KtmdCzPGbn9cYqKP76dskIyuygVVRVrVVl19bNssLfVTHYRQ1oh5
RzcLFWHG80v6cyD8dudAtYKFsyoXN36yPAD7Vp8tkdiw0DTYO7hWmgnb72tM
jpjqoWnr/rYAXgoUBvlvymfrIfeHw+csOSlw9Z6itMBAD+Ov7Xr5tjsZ6Du+
Tken6uHCpLrKTEkynLxV4lc8ROOBxDWd7dIk+KN+VVZ5O416qYqDrGQZWB3L
1l45QmN6hpv+gtIE8DR77vmbnMbmFve8NmECRBQ45zQcoNHagG95aEECXMh7
YnF9L42HBYx65fmd4J3t4NIYT6NW7/Sb45d2wI7dbXH3gmmcPmj347yOOGgO
s+h7ZUrj4iG7gku8aPA/5eauYdO4zzF0xOofKdzv2lLINpzc8+Ls5uyVwlOh
wtOVmvQc9lJjFVEwEGFfKh8muFH2zMVGGwk7z637+PMAwcuH3Ktt8yNB0x/r
r+wnmL4nMXaafSRQsXXaV68IPmBsmcItEWCWwN3kqia4o0/Kdr4XDicvMjX+
KoJezTKjrJhwsB5LZ4nvEIzaaediaxQOi5Nu18tvEexd2caSr5cANz3AtKeO
4JTGhb6GQyK4cSMuUlNDcGRmRWBziQi8SMF1dhXBd4OinJ2+IgjMehTjeoag
SBjTNlopBFGuoEVeSDAzbs7s3BwBvGndN7/sMMGU5StuFPkIINq4XKbMJXh3
kddQpIEABoOaVK3ZBH93/95kQSsfEvPf2fZkEHy8hNXwuoAPE+rpqZp0gmbf
ao/d5fEh1czpAXsPQVK6yHzMng/Zx3fsd00kKJ84HnFIxQD7SeET/3iCrtKK
oa/PMpBvedlZHENwoUGKbv4+BmbxOrNl2yb/TfF1uC5koPi09oU8nODLrVNd
OtYwMK977jdlQoLTqqY39H3JwJn5qw4peQSvuNgOG5gzYCcWvW7dTLDMbhsT
wGKg+kwGt2cDwYdrP3EzP/BgWV9FgSaYYKpXp7CrlwfKRX/2swMIrrQWVyif
8YC77S8PG1+C5fm6jZbPeXDjPPuEqxdB6yi/Y9P6eOD1fsmgv+ekJ3Sy3j3G
gxbHEG8xEIz3q/Wu0WcgcHvCaZkbQV3Z8jsdFgx0/Hr0g9yF4FcnDXw5Tgxs
HlauK1tKMKB8F53mzUDXsqdlSieCxmeXZ7yUMCBO/DTRak8wu2Nzl3MWA28v
cUJ6bAnWDVye2FTJQIzG45zmi8m/XEWXewcDwyskukYcgj/pzRy/rcsHbcMv
1a7mBCPcm0SBEj60PHY0HP+MoCW9vqq7lA8nRy+EXzEiGKYf3OT4gg9uDtes
3FgE2y+qWgSRAjD08ZRNTCXoI7qgHK8VQKe4SX31E4VpUeppF3WEICtuz/n2
A4XFHWHFUxRC8Lkc2qcdpHBtv2O990chmKufejS8p9CdQ89eECCCS4Z9Gm4v
hZK+vDgLPTGM7p7YCmoKfSWR/JpcCTQWpd3UUVGovJr3OGhEAgW1tPX1Fgo3
mCaYVISFw9J3Rg9W3qSQ/pepn+m8FeI32az2qKXwlGXz8/vvI2DAZR1nTd6k
F2KccfeEFDCkPYk+OOndlariZ0aDPDb0YWMWhQ+LtlU/lkeDQ4VA/l0qhX6K
uRlNGTEgNUvUesVQ2HBC/MGmIRbejpx6tNaPQiutgmuKcZCccktpMoPCyK55
tzYeT4D/ANC9i58=
"]], LineBox[CompressedData["
1:eJwV13c81d8fB/BrbyIrLamUUbRL6rxpIalEaYgkoxCRlKwSiYiMe83rXqQi
iqKowzezhOyGva+dPe/v/P7h8Xycz/U595z3+3WOdRY3Da9xUiiUreTH/3+3
69xoaObhLvBplDK0XXUHYnnnHlwgVnxNxU833IGzFzhdrnJwF7CnT3U9vnoH
yjiWn709y1WQUJFpnN9+B96c2ilH6+cqyLjyEnO1uoHbsCujrZyr4EllbIX/
33sgrLzw7uYTroKfN11ERls9YXcib91TIa6CcAHLmMl7D0FToWZr4AxHQYSY
ms6jbUGwyl3UI4hFKbBZWC9BOxAG1bHbTKya2XiLp+/tqf4IqLLJ6t2Vt4hP
Na+t9UqgQYrSz72nPs9j56Zmdm1JLKhHVzveKp/Fxq0KF0wCEuCxvivnbMEs
Fp3Kdw4MT4C2pZXhXrmzWPvA/qychAQIuWqdE/RiFotbNtLm3yfAsOrSYsqj
Wfze4x1W70iANKwS+EdrFl/9pLv+xSo6bO72Sz78aQbX7L+IFZzosE79wG/p
19P484tzsmfEEkHfsD9yK3Mai4Tcz8xZnghuLhFnjsZMY+tOitgK2USozB2q
uB04jQUm+VbVyCfCPa14XHdjGnP0WZqv3p4IdYbspDDVaRwyJdUnY5QIfrcL
b4q+mcJNMtk1URGJ8DbKXlUxZQrb0mWKn9MS4e/HFf0H4qewM+OFVVBcImxf
crKwD57CD2NSsp2TE6HVT8G4wmEKW2uW28q9T4S91IcaT9Sm8Mpk6ZeltYnA
+nSEh/ftJC5Kz98eJ8qA16E535xTJ/Gnc9Un7kkwwM5G6Vl7wiQ2Xv+f3Rlp
BgxJiaz8HDyJL5d6Z82uZsCoY526i90kZq+9f1p6CwOmFK9e6tg0ics2PNIM
0WMA53PvrC/xE1h3dbu3+wMGfLX956YaOYEVDPWaGx8xwBcsD0Y/ncDOPEhZ
PYABPMPHym/fn8AlPPh5XQgD+HXFWlUvTOBdu6WjBmMZIMqOE4yRmsCnBz+N
m3xgwMobeVfuBI3joV/nJ9f2MCC6rOyoiO84tmf1bR/sY4CcYoMK030cdzoH
380eIG4fnay8Po7/Mf7K7xtjgOw5xSeKuuN4fkGsSmKBAZKHQt/V84zje/Op
1/nEmSC80pprp9c/vIHZNTS+hwlBbi595a7/sCHnxh+gwQShBp8fZg7/8GXm
yjcBmkwQDI2NCjT9h+8pxrmIaTGBn69GpXP/P/wz9sKyAV0m8IxrGoXNjOHa
7ozrHReYsPRNInnUaQwLbRPY9fMuE05TFd/LXB/DpUUbtnm5MyHpmkbxQYsx
/NM16p2SBxP0KBbdQYZjuDeGccTFmwmRu95uVNoxhh8cckxv9WPCFvrJlCsT
o/hbasAzsedM8La/+uHx4ChWblfuexTOhBqNOyUZXaOYN39DyUwEE9zq43sW
60bxUa8AqToqE4oEhxWj34/itJebeK7EM+HS7aAXNa6j+NqxrkN1qUzI0Kbn
zDqM4vNfb4stf8UEzmXZpfLWozhpjedVg9dMSH31p9fh3Cju12ckf0pnwkSr
8mahvaM4MxtF2b5jgnRE2Gn/fyP45NSXi+vymJCv9smpvmME299/Xrw5nwkW
39pD19eO4N8vV+epfmZCOlu9BmeN4GEPPT1lzITD16sMZ1xG8EiWIGXhPyY4
I2Ej2+lhnKQtyjNXTsYP4UdXe4axy03Jgj/fyPuO3coxrR/Gf98IOX78zoRc
gyY5w6xh3MZrk2nzgwkLpkntGjeH8VEruy+R1Uz4ceXc8l2Xh/Hug76qJj+Z
EH9N8IjaiWF8fdtgmFQNE8DeMXW9yjAe9b+n51vLBN/7mjeFe4dwsmGw07YG
Jhh7j9B5G4Zws+yKjBpiRV9GDaV4CO87/qnvZiMTygL5d08yhrDxqJZhfBMT
aCF51iOhQ7j21paHar+YcP25A63fewjn+ux7l08sElO30Hx5CGsbVQl//82E
lnj/rU0nhvDxXTf26f9hQiZDw7xGcwjrtRy3Kic2fEX/WiI3hAWLxL/m/CX1
8dHuSUbDIH5z1Z9q2cIE689r818WD2KPCKnGYuK9hTVDzOxB/C7PS259KxN+
l+09TQ0bxJVFhu9+Er+uGHgQ5jOIW/ZO8q1vY8L96vjsIMdBbFPUY+lIvKaJ
W9bHYBDvu6u5h0088idH1/3AIPZz7k2HdiYUtF53v606iNdvvrnVkzi0c3X6
zZWDOF/R/WMOsUVvdYut4CBeSs0+OUS8Y+DhMsvZATwZVjK2poMJ3CO7tS/3
DWBvZ1O6PnHDv35nk8YB/NZ30yVX4hdTscmGJQN4d1rdxlhit7mTjfrvB3CR
u/TSZ2K9JU6BY0nk89vo3X+IV3J+0NB6PoCv5io2TxIP8tja7X8wgNNSzbuF
OpnwWWBV/C6nAaytw8teTRwsUlWlZj6A1wX+VFYlNhN/wKF8cgD/t9zTdjex
utSuHRsODuDAHeV5msQcK/os12wZwNsHnBQQcc2qmEjZVQNYfWBv/AFiprxB
mYTQAN4Q3qa2l9hlA8ec8BwL60yrNm0lPrI5W4Wvn4WrGjqj1hFLq1qbcjSx
sGVKxc1lxL1qciHzJSzsH5hpPk/mn7vjR8HkexZmm1640UGsLvNURyKFhdUM
PsaUEqfO6VdvjWTh0Y4K65fE8i3CJsf9WPiI2MFqf2JqYUWrtSsLC4nn1FsQ
L0sOsva1YmF9p043DeLHj/VH6GdZ+DG34ydR4jsnK5Z+7WZh9dCrK9PJfg1v
D/KbUmTh/AxlzTvEVtL6ostlWHizvsTSAeKzzd9X60/34+YIcZtCUh+VBYHJ
Nr392Kr6ubgn8dGk41seNfbjUaaQ2R7i3Te+a37J6cd7AnFfAqm/NwaBRb9f
9OMz8Y6KJ4kVtx/Xn47qx1W6HELzpH6lZ79dVHfrx1t+fvyuQxz890mnvk0/
Tt985SWrmQm8BXo3bE36semyXvUA4km/b/cYe8nzSzGQT/rD7voTTry5H1/k
Nf52irjrhN6TP7L92ItJ4egg/VQn+Y0mOduHkzcaus2R/stmlOf6fezDNxXk
krlIf6r6BQDzZR/mvBmd5kH6OclWtwzT+vBCrYv1BOn3cPXyxpm7ffhMt+9c
I8kD5y9lU9c1+nBk6Z30O3VMYCU+9vRX7sPpyhutG0meWDzS4U2S68OHz/f/
3kFsqF8m1TzXiyfat0x0kDza/rt058m8XjxnvY9XrooJrz775994TcYDN740
q2SCQuKxw49jevEJof0iDJJv4jalZwrde3Ev36NVqyuYMDpZcmuHZi8+edLY
abCM5Ll4SabM5x6c3mVWMlpI3ndZXK/tdQ/mb+DuGS5gwtSrS52p0T1YOF9Y
p5/kr+bhcSkNtx7s3m/cVEPyucxV3v3Sjh48Uf76xN2PTGj7c+9I4stuXJcU
lSmZyYSHm0tabKndmKF2IqzxDdm/2+Ju2/278QGFR98jyflgJ5r66qtlN+Y3
4tMTJOfHDNSL9aztxqEDwfHFyWT+KWq/lSO7sO6bP93B0WT9/91zHvftwla2
AiwZGhNMUIlwvnMXDvUrvRQXRfK66RI6cboLXx5aeB5HzjtloSdJN4W7MHPL
uLNzMBO0Hbscsh904uMjRcauPmQ/NGhcB2924Dab5N5OS7Jf4Sv/nr3cgRMt
KjJ5r5L3Dcdl3zzRgft7E/YrXiHrwWBcS1TpwKvMJmMvmpK8EUor5e5rx65p
Fp8TjJng+fdL4Hfzdjx9+fd48BHSP55dy88btuHYzJhQqfWkH5usBp202rDY
+qfDtfJMKN/eX/REvQ3nCqw/GrKG3B96h27ni7ZhPN0yuLCCCXIJwxkoshXv
/HrPIHsZEzZqSe794d+CPYymZCIWGaDxyEyn7/pfjJo3/lWtZ0D5jEfy5eN/
sZuVzN28GgaY2MVy1qv8xZFfzdOPVTPg9plf+YWDf7B5K0enyXcGZK47sy3G
4Q8uHeLNvlzIgE1fjsoZ3PqNDZ4kec6mkfvQ9JbBrLtNeIuSUJQbub+N2CyE
eD2uw8yiFQcFNjHgYVGMBpopxauCrB/VHUuEhjzNoldcmfhOeKZb2hU6mKke
XTQKfYM03Bn7p4fjoV3HxNghvgwxztqPv12KhYPd96/8aC5DDtkjxz9Nx0K0
T6K96upytEL1yvni0Vgw/MR6xIopR3u6L2d1dsTCfyoeH6xo31Dgpl1HtUtj
gSHKkDELr0CP5UTZT5/FwnbW5OCmI9VINOViBnVjLExkoAzZ/XVo2/iZoEMX
Y6BCZb+Eo2kdulQ1UZhjFAPMF7tvl3rVIReFgsfbDGLAMGHLfteiOqTz59G/
XVox8DZ4ZWntiXqkfqnc8OGmGHCwn24JNm9A99F12X2T0dCrnCHK49eEHhiE
HT4bEQ1fUl45XUptQkuyrYshIdEQqZBSl/WtCW3eYeZbHRANR+Xioi3EfqGD
Vt/XWHlGQ5JAoGIB9RdyunbsWZd1NJj3WR10f/0blQqMee/WjIbfyWscxqr/
IreZi272LBp84E0zTP33F93C51e+7KbBc+t9e8wkm5H5E/PLg2000Fcyovw4
14xeTl3uDGikwZe0gLDUlmY0sUX3y7piGiRmTXwwG2xBXnpX93xPpIGn5IMY
adFWVOfRYeobR4MLt8W8f6i1IimncdoRGg0k9ijr7nduRQN6px3antHA95PZ
H+m5VlSP7h9x86ZBtnk3u+dhGwoJ1LXZfYUGCuKqt/4Lb0Nh+rWcLqY0eFZ4
qysumdhuu9yn8zSwV6CUG5W2ocggmUYzQxoodsmF/SfYjiIfW6hIHKYB1frU
hvjQduSxu0ZSRokGvLJRkXcZ7ahpbd+a1o00cClr5jfOakdzV+IXMhRocFrp
xpBQfTta36Jy1G4VDQQHHn24K9uBDvs7hZ9cRgN3hzwd44QOtCkDTRTNUkEb
TVqKZnQgx6nddbrTVBBcpuZT+qUDpZjujmuYoAL1LePjvpYOZJvKMcU3SoUP
EwHKa1Z3ojV/zol/66GCZ/HXo42qnejAhcePg7uocCRy0eKZZidaHxq5dKGD
CnV7nGI4L3Uio2CRMaEWKvy7ZyLcG92J0m/sUq+sp8Kn42Gb6a86kapDaGtN
LRV8VlUcPv+pE1XwZyS0/qTCsi/I4/uvTrRHeOS8eCUVtnJsGs6Q7UJup7f1
15dQ4br/ZJVbZBcyq4Oaro9U2G6iNrgtpQvdvnfA6EYuFWY32/IPvO9Cze8f
Lc1+oMLjb3/BtL4L/RMNr1LLpkKKSNFbLclu9LZ1593FN1RwaFn8Mbe+G5U6
PjtTnE6FXRl7WFk7ulHYwQ8XI9KoUHTqtYLimW5UFdA8deIVFdqfhz0XDOtG
9IfHzdalUCHVsiLja2I3yqXLZygnU8FxF2/F/bfd6MnTpS2aSVRYarjLPVLd
jS7S1HLvMKhQ8iJLPrWtG+Wd5fmamEiFp25DmldGu9GBtf3ztXQqrJK7crtW
rAfZ/F0/dSqBCl2s6NCgtT0oMtTlQ0w8FV7n1aUfUetBHnUljOE4KmiY6vTk
GPSg29dfLKTFUoFz6wNOp8s9SOZjio0scflS3hplhx4kHhqwGBhDBRO62rnY
pz2oSEI+PjCaChN8LVdexpPno9doihI/uxlk9yGjBz0YybCPplFBtVHjzteC
HmQSdWSbOnH5wX6f6p89qIe671EVlQpWKVFBzR09SNCv0tqNmFP0aBRrvAeZ
P1RrVCaOvz2ROM3di4ROO/7siSLzbWakcUv3Il6jp0ZpxA2HT+eIb+pFMvQb
l92Jb6WxC9fs7UWmnL8GzxCLSr6pUNHtReeEgxZ3Eb92v9S490Ivilq7IXgd
sU6nYMeRG72I755SrDRxl97HQcP7vShNW2SjFLH3O+tps6e9SDxKSXF11P/X
V5rTPr4X7aIVxm4lzvUpEr6XQcY9Ofx0iY36b8n4F/Si5NilFnvi0VPrFMJ/
9qIIj5ykGOKnuVWqiR1kfgPL62qIleU997wZJ9+Hq9t2Ofm+Jf6q2nncfcjy
SLeFKbHFyG9yg+lDw6m5eZnES2cDztUr9qHlaYOOImQ9o7/ssejY04dEbGY8
nYl3K/bYjej0IYskzY4O4pqn4XcWzvehbBZX1EWyPw6T2g8EbvQh3ZgHtGZi
QdOxIOn7fWiKo7HbhuxvSlFC1PqnfSj6hoDbInHr84W0Axl9qJg731Sb1Mv9
+Vc5egV96G5Ma+o/Ytmr5/8797MPjW7brJRG6uvU9g+NTuN96LO3T44GqcdB
mmWHJ3c/WtIKyBcj9fqYY/lQoFQ/kj6W1DVEXFB9kzNlTz9CpbKBxUwqXNq7
RiRLpx8NLCtgfyH1P5NQIVNwvh8FOkn7FpL+2HZTactv9340o+3/ov0FFRgi
nedEC/pRpMBCyVPSbwdvh1qs/NmPKnrUN1WRfvz9F9lv7uhHv0Mc7VZkUGF5
WuwDbW4WWvkbQsreUsFX72y6qw4LdeVtC3uWQwX5d9y5D8+zkNOi0R1Zkg/5
K7L+e3adhVQMBsRffqLCZJ9Y06sgFiqv0zDo+kzq17+cs7WahRj8f7f+V0Tm
L3Yjp7yNhVb9CZR3JHlzOkrYLnuUhX7sqnDfWEYFzRcn6wOWDaDApFcmKd+p
IFFa/2LnqQHkWJ0uTCF59pm3Q/9J1QAqXxzYJ9pJhXfBDzlutw6gj5sGhdS7
SV5Ib/xgNjKANBMu9pn0UiFM0WbtLrFBJO4m8QIPUMHm6PBYq8EgOojX+reR
PJb0m4vcVTlI1pOP7zo/yXfR2OPyLYNIrLzOblaQBuyIAxSh4UGEJtRvhYjQ
gJXsbdsmMoR+B9uhKgkaFBTzagaeGEK7RXJWxa2mgR3P8va2iiFkcE74yY2d
NPjqq6IS9H0Y3Xj/JineigYa7Ai64p9h9P63TF2HLQ3e3aVIF7KGkcdilpuq
PQ0YDvWUKYERFBCzk6vamQYPTLwazHRGUGlTz7FzPjTQVq312VE8gj6dXF5r
Rs7TNtsDdvZmo+iN8J9lCr9ocFm8+z/J82NIRjv4XqppNHBH100EXhlDgX5w
7oZFNLxSKFLkuj6G0o5c0d5O7gPTO5hPxu6NIY99ESsqHKPh2dkrhpWxY8j7
qgFzq280fI1pbvdvG0N5pl/akl9Hw2bFBsq89T+09rThsuWL0TCuUXqw3XUc
7eKTeh7/IgYyborMWjhOosMZLkOq/2KB26vHxyJ8Gl0ZUJzlPxEPDs5cxldi
p9Gdr2XpjafjoclafrN50jR6Gzo6zDwbD2mnLlSaZk+j5ceyl3aax4ORQqXc
+TryPNsgba9zPDCL32edlJxBP1ooEoXUeNAWftStGT6DBFK/UtPbyecpjNz9
sTNoKHpZkEJPPEhPfgnUSJpBWnvBKJwVD6zmmW17s2dQf0vECYfxeAjLsPPc
UTeDcpN0T8xzJ0CHoZGMsuQsSlv89stlUwL40NbryoTPorLov12Ftgmgype8
MBQzixSONilm2SdAk4tiZhFzFtXlLcN0pwTYekpJ5lbWLOK/O9XqeDcBfvOq
dVfUzKKLZ7Mtux4nwE4XDZ8H4nPo6deIZ5IvEqDH4PTHoeA5JHjw29K9tgQI
za+xK4qcQ4elcjyHOxNAU9lIPiZ+DlkFp30z7U2AMJ5z/jpv5pDlq7sf1IYT
AOVfMmb+mENei+N29PkEoCrZjJkIz6PV037DTYJ0OM7tpVQUMI8Gel/0a66n
g/lnjbrDYfNI4YvlYOUGOri6TnkWR8+jW+HcWy8r0oHRZ19X8noe7ThSbemi
RIeZikte5T/m0csY0yteanRIitCorxRfQOseDyWJ7KfDJ4Mpr5NyC6g4uSzg
liYdqvjeKVcrLKB/DzSf1B2gw9xdJe+fOxZQY6VLfQjQ4dRlWZU64wVkfdT8
QucROiwoTnn/oi0gyZvLTy2cpINE21uVC4wF5CBvMbntNB020ewbfr9aQAtW
DyMsDelgKNSt8jePjBtafCkwokPqSG1DS/MCOsQjd9/oPB2Mc9+qdq1bRJJr
6fO95P+D6072jZbKi4hhaGk9aEEHb2WlB93bF1H62eSMoat0SIulN/YcWkQl
Axcre67RgdMn5EH/tUXEt3RaIsuWDm907ZuGXy6iwvwSpXeOdBgNUzBpebuI
ripnvr7tRIftfxsbf3xcRP4f7M7vvkWHDw7ajWnl5P1jo9sznOmAn8s0XO9f
RLWhD3M9Xcnfb64wvjC2iJ6YnXNWuUOHw4oP6nVnF9HhToOT9cRluYN1mwWW
kPnKW9fk79KhurmwtmfzEvIPt+l57k7WZ9OdMw3qS2hKf4uf0n06GDmq1hbv
XUKZR+o084l/cUbVJOksIYvgsv4mDzp0bLL7aWGzhE5vTVk74kWHDU7rThs6
LqFfuRVn7L3pYPWpoVrLbQnd8aVE9BGz9LWq5R8voSat4n0NPnQYd5Kuanmx
hO5dihx4+pAOu/O+G1RmkPEDw1YTxG7cPpWfc5bQ4Wv93ed8yX5GDvyILV1C
tuwdLKlHdDjYRj8RVLWEzlRo2jsSeyud/eHeuIROhV0cKSXmzi+ouNC7hGa+
mQ3a+9HhKI+rvt7IElr3XuhaPvFjA5WKfdNLaOnRhSY+fzqItEd8l+Vjo66P
Xq8iiE8pHz/OL8ZGkzFKfE3EYc6U79PSbGR/ZfUFmcd0kOa98a1BkY3WW8W0
BBGbnJTXK9nKRs7SvKJfiaOp9eXvd7MRZ6C72hTx3/YnuskH2SiTc1RLMYAO
a1SgPPwom+T5+cNniM1dJnV8DdiIvfBmlwcx4/OrMuezbPSYgyWZRNzFa65z
9TIbvRTiaS8l3nRKqszQio18LyxQ+4iv074d03Zgo76iH/t5n5D66/Aq3ebK
Rh9y3crkiYdVdh1b58lG+zWWDuwlVr/NKlnmx0bBneYx+sS3viQcpQSzUW9n
UqspcTafcclIBBvxsqsE7IinTgkebY1jo8Cgftk7xHujcXFlMhsZ35kT8iLu
2pqhafOGjYwK7270JQ79Gp/FkcNGVeO1EY+ID5oEK8dgNqpvqbv4kJg16JG4
s4yNJu7puXoQR/rYy1ZWs9H5Pu4uZ+JD0qYh1r/Y6Jd9A9OKePSVPi9HBxsJ
uz3PNiaORZoe0Sw2sqRKSGkR69SpTOwYZ6NDvqpFm4knbVbe+DHPRiKS+QXC
xD89eZx2cFJAXzteYOj/66froerNS4Evyu/iyok9l0/2VghSwF3glzuDuK7Z
jrlCjAI8GlOMO8Tety7JZclQwCbydak0cYNmXT17JQVcRaW+tpH9VeXTD9WX
p8CCgwzXC+JfMfv5ezZTYNtbuv5mYnWrrK/bt1AAcgSMOkk9+amreHlto8C3
waqYaOLtxXJTshoUkNt7sm/x//UZEvb22kEKxG2dZ70kbj0vaP9OmwJXJFlK
hsRBw7Odx49T4JbhU5MoUu/dMr9+el6iwAqFo6b/7x/NjlNPv5tTwCLgnOID
4rC0Mh3ZaxSwH/m9Xpj4oFYufmtPAaa5a67kAzpQr0eld3lSYLmA0tYZ0s/D
O8Vstz2kgNnObm9r4sNsvw2e/hToKy+m1HrSYey5a4zMMwpMibj4xpM80Pts
HKCXSAH/lbcf8pE8ofv/OByVTObnu3uTyT1SP6ePULpeUiDR89RYEskfZs/O
Ox7vKJB7zWByuxvpbzHJa5lfKZDvsz5F/jbJR4sakO6hwJmaQ2wtBzok294f
NGBRoPqpC995e7L/jopU/2EKNM6MnbCzo8MTz3sjM1Nkf/fsee5/nQ7XohXi
//BxwIrdYkPPrOiwqubWPF2JA4ojloyGTemw/NeqF7+2cMC+83enyy7RQbCt
xFBiOwfwbGBWJlwk59OQ3KuHGhywbuq3rDbJ/1qBryZWxzlA9Iu86HlyPjzR
Wv5BxZ4DElhTQtW6dPDR+XzF0okDektRvr0Oya+T1iJxtzlgKcicxneM5KFp
nqWoJwf82+nUpHaYDlp3ry4fC+aAZ4XrLxseJO/LzHL8kMkBqkdXqbduI8/L
G6loTXCAXmPYEylpOjiP3NPaNsMBq49nnzKRJPX4JfHcugUOOHT15ccoCTrE
XBp5yMHNCSvoJs7CYiTfqU/+Fi7nhE2Bsj4/+OiwT7zwqfYOTnAPP1Z+cSwB
hDm3jh66xQmDfMmclRkJsOKnEe9OV07Yufj5P/fXCaBId1+14R4Zd24uVCD3
Ca2D5TrcDzghUpvicSk+AVzvWSYWhXKCww2xHs+gBGj9F2145C0n6JcLP/Ky
ToB3nXzvj45yQnbhgI+jbAK4NH0vCZzgBKpop1ODRALs+RHSVD3DCe8f57rs
FEmA/A8yCxc4uODF2mPZTRwJUByw6dBNCS5QWkjZ8B+5XzWpHfsZtZMLzH7T
FFbkxcOSu99Qnxt5PiNv/rdxPBQ46rG3eHDBRU+t9YdPxsPDa6Lizj5cUIwp
Dsk68cB/MnLnYgAXSKcIOJ7aHw/iCinu4jFcQP+Pe7+OfDxsKCsW0PjMBUUz
n2NxXxzoSfIoBnJww3Gr0cxa5zhY//MIw5+HG0RozJyvdnGw8NRvja8ANxw7
xLOYdi0OMq9SdWaI+8591m62jIN0Pn4ZD3FuWLfm2rEb5+Jg2bwg/8113LA6
SWQiUiMO6tuXsQy1uEEqaK33a3YsmGeseiPnww03Z7qcqn1jwVVv5+5XFB6w
/Tp3S9c9BkYOvWCEcfOAk2+Hz+dbMWBzQE7MnZ8HHtUcuqByPQYuqlP6ji/j
gejyVrsJkxgA6Qra0Foe2BTeiTR3x4BQx9VFdcQDhzmzG5gj0UC/+7wox4MH
TNoE5/6S+3x56j/D4jke4HzQPmWwmQY2LnNn1wryQUSQk9SvfZEQc2Tn3Bk+
flAcjm1S3PAcXqWfi/oixA8dLkVRaVLPIVfKfafSMn5Y9VHujirfc2joLnRY
kuWH1AM5KXKsMFjmd7IzVZkf+LJ2HXibEQZ+JdcrFk/wQ/uNiJFv+8Lg1rHE
+Bfh/BAo8dC1UjcUfDKK9kvQ+KFSQMeZrhEKz2T6mu7H8cMCPVvOXiUU0nu3
LjdM4Yd5pUH3KeFQ6PPP91/I4YfCW+kC36qeweWyRsfTf/ghVJnfyc3oGdir
z4vkt/KDWLGu3vThZ3Cfuua1Yhc/JFhyyTntegY0m2vd84P8kCXjn2ks/Qzq
+MfPpyzxQ4SqTkBJRgh0OkpPL+MSgIw1qbkuz0LgX9O+cHc+AehvLNiwwikE
xF56V55aJgAXE9v+HN4eAmvEk27kSQrAUarN/gqJEFC9W8qvuEIAgrxeU/TG
g2F/Oyv52WoBWOcjsPtLbTDo6Yoeml8nAGE/Rro3ZwfD/wBNbZjc
"]], LineBox[CompressedData["
1:eJwVkX881Hccx8vu83WXuPxYF9uUZA+aiVTjZrzXHd35defSwhWFo1/ze7Xu
JAmTX+f3zjrK/IqHeeRnuWvt7SEVMaT0qBj1SKv8SCXqaJv98Xo8/3k+nv+8
zIIjRRKtZcuWcZf2P935epwFMwauIyusPv8xC75uiZ+hmzLwXmNo2S5xFlib
TZSsNmHgS8U8SXDOAr33N97ZGTDw+IKmro2WBQMXTtUf0GJg3+hp85GcTPCj
vzUeekxHRaMR90F1BoTfeDhzsZyOlRF9cQ/6z0Cya01pmAUd6crmIcmqVIi+
ozHIN9TGkI8tayaFMigSCQIKmNpoffarREpfBlf6y8sKdbTx1ytey00GpED1
etoqPtLGn90c99v4SEHZWepV8oZC9UiPB9PnONxs3Z5WPUih8VrN5E3BMTAt
Tv+gzqcwT/IPd44fC1zWKPd3OYVs9RZZAiMWDhbaZ17NoPDZE34FrSsGmnOH
TdpPU9ilaK7V4cUAP93G4XoshdP8eT0dt2iIk92O6d+51HdZSDBJjoDuvSbP
nhhSaHCB/aW9NAwE55xcNEwKU33+qlMtSuDOWGARcyWFCTFNyn0nJTASUs5l
0ygUukZvNkwJhZmD1mXyWYI5Ri9mY/KCgXXUOYB9l+DkwEElPXYvKC/taxAM
ENS/lr0NEvfAuvdJdEkvQdHlYlFKthg2xt9olXcSdL0+yWmp9wfnJKHheAvB
H55Pqe21voOOjphDmgaCgbeHFcJPdwGPFLYz6wmWujyKMnb0BZ+0+5HsKoIj
+zdKXeJFEJod3CMvIvjSQ34kM8obnvclm1fmERSozkS+mvKECP1qqSqboJcV
NWwV7wGvRV0DfekEuYJv1TRTdzhWMGE5nkpweUX342tDPFi8q5uoSSIYIW3o
3VC7AxJZtveYJwkap1g4CJRuQPxFNhYyglHR5151NrhC+i9xKexjBH0zItxs
X3KBOVw0LIgl+GfeoKOuPxcKPmuzl0QSTPvNecrzDQfWBD1Mlx4mKL00ETWK
HCg5/+GRPJxg1wZZtU8HB9Y/XutYGULwJzvdsXwNB6rMt+eogpb6bVTgYigX
vpCE/t0nJsh+mv02+V8uXKxKdR7fTbBGq8nT9rorbH12oVCzk+CK7s3T0mY3
UFndmmQKCZr2v7fDzh3gfHiKY+G55LOKm1Le8qCjjnmWzSN4Sx3eVuviDrxp
u9cCLsH13e8aT573gJ5NvnwJEJScSEvsMfKCoUbFvNyBoHdr/72c1UIQz6q8
K7cQvM+Paw57KoSxrSOVKluCHtvNTeVqH3hx2cx33JJgUD6rfk3sTvhwtfYi
25ggb8D1ExZjN/Q82LRywYhgoV9+exndD5RzzeHqVUt/C8SHAlb4g5PNH6ZO
dIIV87t19xiKQVoymPnNPA0rxePeSduCYE62eADu0tDbfU/QwmAIzDh4m7nl
0pBVAYIT3x+BEwmdKgM9GmblCjJtFEfhP2LU9DU=
"]], LineBox[CompressedData["
1:eJwV1nc41e8bB3AkIyMtSiUtuxIloZ47SWSEzL6KUEahkDStIpEicazjnBMa
SLIa9Egqm6xk7733PH7P759zrtd1nevz3Ofzud/3/dlp5WxwmY2FhWU3+fj/
d7vG1brm1ex5FQFmx++K3AIFOkfNU55VeXtusLTUePuCyq6q/YFzrHnrniQ/
SmUNgm13+e8HDbDkPcuSPRO7FAKVMQdNrzSv4MQKmsnVOy+hwi699/DXZfyS
quXYOEOBRMk/inq5i3jjHo4y1a5okI2qvO5SNI/j2pO55hyo8FjbnW0+bx57
DPjn77pJhTbm1jDPT/PY69x6b3VPKjyzts0Oej2Pq5U/rLr9ggojMszlxEfz
mGfbcGfAVyokY+nAxhPzmG9bw3GzNXEg0e2XoPZlDsdopWrRGHGwU/ZYg2DS
LP5PYMAt/SINtA36w/e/msXvDItRlBUNPNxenlOPnsXsb/XK7l+hQfmn4dKb
gbPYNvPPPmUnGtw5QcU1V2dxfIP9ZPB9GtQYrMSHyszi+xGH10RF08Dv5ndn
/vczuGL3WyulWhqkRTjKiCXO4DRWO4tP9TRo+ryl/xh1BssrfeSVb6KBHPOG
lWPwDC6373EX6aRBq98uo1KnGZyXOHW3bIwGihRfpScHZvCai9ubA3noMPDl
1GqOtGl8YeVetd9xOiSFZBe7vpnGLdJ1g40n6HDNTvJ5e9w0LrbGa/edosPw
Jr6tucHT+H2BsccPLTqMXa+Rdbs2jTl/38Z5pnSYEbM27xCfxuM+R74v3aAD
2wuv9G/UKXwx6O6XW3Q6/LCf8JAJn8Kzhhk3dOLp8BBsjkc9ncL6lf/Ed7ym
w+qR00U3701h44GsZ1nJdODSXNsqc34K86TFncrMpgP/Suya6E1TuCtfkn+i
lA5br369dCtoEvcn3vK0naJDVGGhOt/DScywmv6UNkMHYbE66Vd3J/HW51Nr
5+aI28emyx0m8e5TnIK3l+mw2UTsiZjmJF74meOpzcGAjSdDPtaunsQO3xaz
FYUYEEajRlxdnsA1N85jtS0M2MBMusc6PYHNV4t0am9lwPrPv9T3dU3g6iPf
PPV3MEDgwFLDw/wJLBG6lU9GnAG8W21XHfKcwD4b8u+IKTAgyMOtr8h9ArvM
NmlMHWEAT513mYXTBM5xT5b/dpQBa0JiIgIvTOCq8FQL9WMM4OKsku5UnsBO
6b4bJNQYsHpSxTB0bhxLe21gPajPAGbx+oSxG+NYNXnPjJkdA/QpYplCDuP4
c/CrmTB7BsRfVvp53GocV+kUR5Y6MOAMi1V3kME4HqO0FR50ZED44bS9kvLj
mFdHorL6BgP20c4mXpoawx/8YrPc7zDAy9E66/HQGBaR8li6e5cBVUq3fqV2
jeEhsx999+8xwKOW2rNcM4YrH/7MvPWAAQVrRsSiMsewjNt1yZM+DDC/GfS6
yn0MxxZ0fTkawIBUVVr2vNMYXrikHMj9hAFsAhm/RW3HsP3yGEcd8Zt3jb1O
JmM4wqqZ9XIQA6ZapSR4FMfwzY+f1/73jAGnU44rysmO4YYepprAcwZE3THQ
MJMYw2ZdARL5xLDpjt3rzWM4tHNhcWsoAwRfhur7T4ziXa66ra/DGJBz4MuN
2o5R/EQJh6i+ZIBVcXvI7upRrJbzevofccqKbBVOH8XSJs+qmeEMUHOoMJhz
G8Wt0r8nRSMZMMA+66J+eRSH5y0eoxM/jxN5EWY0ip1z9x4TiWJAU41jtawC
sfUBc/5oBrgiXkP72REc6CEn9TuGXO8kfmTdM4KXRXNeHY4l9Z12yb5QO4JV
3PhUacSfdOuFDdJHcFf/s1F7KgMCDAK1tRkjWJjKzVJEfN74+AP1kBH8bmfT
kT1xDFi6EN+u5DyCGzYfEK0gLrtksuHwxREcem91lQiNAdTLa04d0BnB3563
pTgQO9vnukuqjGD6A+eMj8TgeP3NbukRTNlu3TNDvO7G7obtwsSTjmpH6Azo
cKvj2cw9glMOq1W4Ej+8p+LM2zuMr4OMUzuxkdcojaNuGNd1vPddz2CA2ENG
FcvPYZye0FKIiGf9jdgX04fxC1mP4/bEhYFcCtOMYTzWdLQ7mDjy2Vfb0ZBh
nJnVmfuB2OGFU2S/1zDOunekqJxYJWJnSafzMN5Ln+PsJ+aLrllqvjiMt54Z
urNC3EL131+vM4zXGvzbtf4VAz4wlCyrVIaxocZjlp3E3onDIaXSwzjixo+N
MsQG72g/fgkP4wuFxpbyxLvfn5vO4x7GnstruxSIp9I4xL/ODeF0j2+0//tn
5mfTzN4hvGAqHyFHHP752pPUuiGscPNEoRSxbe6OnLc/h3CIee6hHcSK36uG
X2UM4c0m95vXEnP/fLSD+moI77pi9mOJ1NtQqKhPCR3CvWrSvd3ESaWDPqHe
Q3jOrlyzhPheJTUj6PoQzmETHUom1qnR7/GzGMIry8zyJ8Qi9eybvXWHsD2b
9qQN8WhjtubdY0OYz63fTJk4r9Xh7k2ZIdzu+nUVP3FI5/YU561DWLft9XgT
eV5WvZUt9muG8Duj8D1viOUHfQVs5gfxz/3XY52J2UcVVC/2DWK7R6IX5Ilf
z8QkGPwaxBGVnzNSSX94LJz9q505iC0/eevaEZ9hsnGfjh/EJ1a3yG0nHlpt
f03ZZxBPs+1tvUv6MZd7G/XwjUHMkbH0Zi9xMF9FxQHLQYxLVL+XkH6W3XRY
fs/xQaz7krOdl5h1S5+NyL5B7HnZe+A1yUPVtujwzdsGsVSDhToidtvDusC7
MIDFrxrOW5E8fZIvy5vOHMAtpT/sHpD8yQo91VifOIC98MoQO/GbBe3K/eED
OK+BR+0xySvle2mrrfsATvu7xukRhQG3zpYy/ykM4FMaMsVmJO8jckF+M2ID
OMHcbvkXmQdXBLX5NwgNYM43D9ccJDZuLtmuPduP5/tWMhdfMEDhaonKt+x+
LCOfc/VqCAPe6wYWNLzux7qPk/wxmUdiclrasxH9eLyv4+H/55XgfPF/sh79
+OXkrNjbYAZM+xXfYSj2Y8+ppD2ZgQzIYBR98vvchwsWz3UoPmKAjF8AvHrb
h/2Olvv/95DMc3vNQhzZh7/XOojc9SX7Sbbo79ztPvyxQW53mjeZN98KZxyU
+nCOlPWbufsMkGv4fejs117sXiEYMnOTAe9y/XOuJvVijj6dwj43Buyin1Z7
HN2Lx/kHWetdSf7tfp/7frcX2wU1eKWQfTA2/ctFXqUX5yO3skNkX6Su+/VB
KLcHz4k9TMi0Jnm6uO5MW1IPzum5O3jBigEz78w730T14I9/YwvYLpH8qk1u
UvLowWG3X5LGIvl3F71rLt+DofLoKYopA9oa75yiv+3GNerKrgxtcn7igQap
8C58U9MrRUmO3I+JO66TD7uw1n51bw9ZBpiiX7w5rl14wKlse/p+Mt/qzZGO
fhduULL4JCLNACmeJ/HOvF1Yfel7YsVuBqhe73LK8OnE2Ho7d9RGcn+UIlcd
d+7AfyTfnG+ZpENY2NYm44sduDzl5Xz+OB0yRmIznHU6cP/F0i/xo+R9hsG4
TJfuwBaW80kWg3Tw4En+zd7Xjv0cV3mnd9DhQdO3wBLLduxWpf6lppIOjx90
bTAzaMPRHgZLPEl0UHpkodHn0IT/a5R0LjKhQ9Hc/YSLWk1YTfnsuIIRHUyv
xbDVShPH6d5gGNDh5rl/Od+HGnFyHo+rqw4dPuw8dzDaqRGviBzKXlalg/g3
dWFdlwb84UHapsx9dNg4u28o/XY9vhL7Z883VjqM2i0983xcg3dHmepnx9PA
tyBaCc39xkvexZXpIXFQ91Wl4N2qD1j08Vf/Ta2xYCGjvmwY8h6Napjznjka
AyqUxOAZWiryFJTJNZGLAV6fVQLaph/Q+osnh22lY6B/2Jo+bJuGOmuDJp5t
jwHGzz0Fsv7pSLKFK2GFNQY23kzkzv6ZjU4XHtG9UhINS/2GxZd25qGF/FXa
4ZbR0K5hauRELUQsQS7HPoZHwfHue5fKmgvRgkzaxMaQKIjypjvKbC9C9/ew
sd4LjAKDLwOPBqKL0PIZeU0TryjIl76fdSWyGO3xLLc/6RAFDH6GkEVYKUrv
9p2aOx4FcgPTQ+KnKpH+rfIugaFIsFTkcv1rU4lOm04e4e+NhGA/4Xm/h5Uo
Siy7XaAjEgZ2odU9PyrRXq86E7H6SGD893h7vOofZLnKNz60IBI2lAnrikIV
+lx40OIWNRKmUlHqZuUaNEN/LpFgHAml0srrr1+oQXKcc1tYDCLh1WuFm789
a9CXY5qXrXQiwSBun7J7QQ06eSdX5phaJKQFb/1drVOLJCQ4XLTlIsHJcbYl
2LIOvVzU0fi3NhLU+ydO9PjUoWgxWlMKTyRsvzwSfyyhDp3RKNd/wknOM+92
GOqvQzcXcpnnVyggpV09o+n6F80+60DmoxTolUrlX+1Xj3JGVoS6KynwLfHd
DfM39ejgtThz5TIKhO9KrEkvrkf3Kc8XI4sooC4cG2W19h9i00+1uZ5PgXju
QLE8yj9k+wQl+mRQwLLvyvG7SQ2I48nzfNUoCijaWNGryhpQTMKWhp0RFFjb
doFdaqwBfc60z+QOI+f9NSyuP9SIbnDoV0w+pcD236rGCt8a0eWg3/JyPhRo
SBBxGq9sQlOYLdb8GgWyOJIN3kw0IY01Tn519hR4YXv0iMXGZlQVcVbBzJYC
2pKGLGUmzchI02vRw4pcPzkg9E1LMwoyev9ylxkFoviEblmwtCDZ+1f2dRpT
4JZT/H+Cu1qQjnx6SrIhBWRl8Z6Hl1vQ5JX4W+f1KEBPn8qyGGpBL/S/RO7U
oMCDjT7Rgvyt6FzSzVZJdQqcv7nWq+xAK+IaGZZSVqPA+iNSmsqurUhibG+b
B1Dg4ReLRsGFVnTyXpzwpaMU+PB0prWHrw39QrndaUcoIHL8mEG+SBvyz4ou
51GgwCK1WOE2tCH1Nl72RjkKZFh2r/T4tqEfhodUe2UosGudjEt+WBs60WYB
3tIUeP7dpSs2oQ1t+bXKeY8UBRx3sRQZ/m5D78PcL/uJU6CpSl1Ztr4N3eAa
kz0lRoEzvk9TePrb0OcqhUN8eykg1iUcmr+mHcUOP+/7tItC5tsldurWduQd
HhcTs5MCbKfeuN+WaUc5si7BAaIUaE04/J+sbjsatbeV9RShAMVWbw81pB39
1vob+E+YAhybI8JvM9qRwoFubU5it8JmLqP0drT7x4Nzx7dQQF/y6jBPbTuq
4gyVLxSiQN6/NIve7nZ088oE31bi/U/m/uTPtKP3hUkHbglSYM3go6zbmzvQ
/R6WM1qbKOARXSphJNmBWnQqVX9spECP1oZoWaUOtJvX0leN2HDJjI9XqwPx
TobwVGygQH4yzbP3vw5k43Cw0YpY9kLveP61DtTsPDm+sp4CVL79NtT7HehC
hPvZRGLeb251t4M7EAvvyUkj4rtOXzWM4jpQSEhxEz+xKpq24U/tQBkxSf01
60h9Age8f3/rQHqWtzleEf9ps4v1Ku9AwkWBcreJKWmMz0dbOpDxX0lHE2JL
n6baieEOVEMPyDxGLH5OcCJpuQPJmFWukSEe2a3Hf5mvEwVprb+2izhrKkBK
ZHsnSugzqBclfvDzh/pfmU4k+euhrjjxqfBlq+cqnehRfEDZYWJe2yOemtqd
aKObnKE2cc2RG9Fs5p3I7zxqdyCO4UrK/nq1E1V8snANIbb+11XtdrcTmR6U
484jlnonMrYvkJz/Xo0+Qzxxx5S3N6oTldhvVDpM/v8XrVAJ2rtO9G6Cr+ou
sfe2UjWzL53orcJT22JijeHVl9YXd6J1Ac1zouR+C3xD90v+dSK7pKMPPYnr
g29HPuzvRFx1JZw9xHEW6ZnH5juR4vcUb0Py/Pazio+kbu5Cji439TTJ8575
Y7nGXqILhTxQoVcQ5zKixHYpdqGJwU1tF0h/aKvxW4SZdCH/Gv4doaR/Nm7S
uKtj24W0JJ8KK26mQGO3dwTHrS4U/rZxoZvYwX+6wiO8C+kfnLYwI/0pZ3pg
6GBiF6oUoraKbKXAvIQ912BmF/rIWn10iPhxcRNcqO1C6tuGXajbKaAXLXhB
sLsLHS6T1fYl/S90Te92xVQXOkHNnXLeQYFEvoK0Exu7ke/D/NeWJD9OLctl
C7u7kbjn1mQbkq/DqUcG0uW7ka2vp+v13RQo0EvaJXauGyWPex6PI3lsfxH6
Yk1oNzJ4+87WjOT5jU1p6g96NxobEhCPJnm/fpij9F5aN9LctCe4i8wDZt1t
9tHKbuQQ5KsWfIAC24Qv3axe24OMbmv3qhymgCntgEnM0x4kKnD7wSpVCkxx
tlx6S+1BQZ4opPMkmR/OQdeyUntQ/KVm49JTFCg63u9d+acHjZpLdGZrUkCp
mZHMLtiLvmUlmMzo///6gmyOVGKrmoM8NhT45F3Aeye1F9V1HApyvULy1+8i
5J/Xi3yzz8t32VHg6acKGXpHLzLIWqfd4kjqNQ4wqRXrQw3njdsOeZB58mIp
+VhqH5IpPSN5OZgC9xbfZZ/J60PJXewDuiEU2Gxtlm/ypw91xf7yArIv9OSy
/t6Y7EO/dOMFlMl+yat0Zks80o9+8ryoD0qkAIOv04Q/rx91SLpWR+dR4Ip/
EVtr5QCyCXilKL1AAfO1V7OL2gbQvvgaZu4ymUcRvNcyxgaQ08vo18askaDy
+mxtgMAginG28YnkioT1v2tfH9IbRNwdnx7ZCEVCLkeH9pOKQTRWIDDNpRAJ
G/0Wwg+XDyGllse1525Fwo+H0tJBJSNIJ61xYC9XFCitvKSJNY4g1kz5z/94
o+DjbRbB7wMjSN6pXiFkHXkfcaplmeEeRdt06nbwbI0CH1PPOguNUdR9WOy3
5P4oUJWp9pb/OYp+pQ9c4DGKgjb7Y9ccLcbQV87Zk+IJUXBxXXf+RrNx9LXB
U5ZVJxomlX4fb3efROefst0bLYmBVGe+eavr0+jt6uyS/H1UGKqWnZ3wmEbK
24T+PpOngpSi4bSP9zQay12KOX+UCvEs0eOM0Glknb1JqV+NClEhEgMdGdMI
vrL39v9HhYfpJxut5qeR43bdIJcAKpyfvZ1r5TODZK/wazl1UIHds8fbKmwW
OV599fnpkzhwcl1ldClmFum9Wvdd5Vkc1NuKSljGzyJX9bX6fS/iIFnvfPmF
jFl0eWFq5lBsHBjuKhc2q5lFq0bZ/rxKjYNXPzPTz26cQ/dqdntXV8eBKu+j
bpWwOWQvId6dzEuDZBbGJ+WYOST9bYN+IT8NBKe/BSrFz6FPebpZ7QI0GGie
O6iYMYcclpqt+TfRIDT12gP5mjkUbCXIf3Y7DToMDIWkNs6jeorB16v7aOAd
uVtTKGweCR+N/bVKmwYynAlLw9HzKGeqgvuPDg3q3cQ+FLyaR6NOVsIxZ2mw
X09SyCV9HnXMacRIn6NBA8eB7tKqeRT4cXz0wHkaHHJT8vZZt4Bag1ROltrS
oEdX//Nw8AKSSVv732lvGoTkVF0rCF9AYUXeDcU+NFCRMhSNpi4gb7Pfo1oP
Sb2rTfw13i+g2s6zter+NEA55kavyhZQ7hYtix1PaUCRtBs35V1EMSe28Z6n
0ECL3VOyIGAR7dVIWzOYTAPLXKUatdBFdEDaVkzxPQ3c3Wce/IxaRLXvjQ76
pNKA0edY8ytpEa1VTeDj+0iDuVJzz6KyRRQkw/V+OYsG8S+VasvXLSHlsjF2
zTwafNGd8TwrvITia3TA5jsNKjg/SlXuWkKceznP3sunwcJtSa8/8kuIx7Sr
lVFAA72Lm6VrjJbQZydhelEhDZbEZrz+RS6hZthw5F4lDda3pUmfZyyhibOf
uqz/0EA80rGu4d0SGj2zRlujigYGPN3STV+X0Ke6Fyc4a2jwZrS6rqV5CR1a
VxJ74S8Nvr155m3Rs4Sk3Y1HpOppUH1JS6ZthJwn4js0RcyszvduX1lCfbId
gw8aaGD0KU2ma+cyUgrYUu7QTAOHG45/baSWkVEkm/X2Fhp4SUn6dMsto7YK
jqhS4uQY2t+ek+T3nQV9om00YPN+5tN/eRnxj5zuoHXQYLOS1j4Hp2VE55O5
pNhJg32THPUD7stIYk7iVimx6eUH+4b8l9FYznrrgS4avNd0rB95u4weGVA3
LffQYCx0l2lL2jLiCjBZvt5LA7mmv3/LPpN6Ml91thNnOan+TS5aRqn3y19/
6SPPJ3vWOObPMvqcHOu7q58GSqwpdYH/ltHFmTYzf2L8QqjOoX8Z4cTM8dMD
pN7mUqPz48tI9MvJtFfEamI+tZrzy0gZkO0SceGnoRoJbiZqC3z8MX6QBmvY
GIab1zGR7urWU1PE2lomNZxbmGi7aUYJDNGgsvl7dY8EE4VdWX5XSbxe/Na5
OlkmUil5xLZxmAaG12Wqfyoy0c7zWRqGxP/YIqriNZhoNF4ztoR4m7a2QZge
E3ENl75lG6HBxZesVb6mTJTifpSuQNwhfu2PlR0TTb7bofOSeM+NnfoG15nI
KP3nCia+8qWu8oQHE6Vlp0T2Er9ZFaR30IuJ2MLGt/COkrxrn6gUfcxExdHJ
nvuIZcJnzgo8ZyKXyMUiLWKn1qSKlQgmKi9bmr1C/EHi0tnROCZ6OVrB5Uk8
eUOwouU1E5VKhS++IFb4WqJbnspEeqIuFQnEHuze5bnZTHSi2983g/iLjoJu
CmaiB2zLgnnES+GDZTG/mUjcbe5xIfHxNppOUAUT8b5l1JcRe0kal939y0QH
fPauriTOd+HRudrKREcNkvkqiNlz8krP9zIRs8l6vJhYfbW79plRJhowvpP6
g/ixrnTp0VlSr82O05+JiyPatCRXmEhM+GF2EjFf+8uSzZwrqH2hYSWKWE9K
S4tr7Qqizpns8CcOdWUpmRVcQTquspuuE9fkZJ7pFVlBdorRbUbEghxXi+vE
VtAjn0ovRWLTs6Jnfu1fQXiOe1qIOIpSW5SpsIK03vuoTJH739T+RDPh+Arq
23nLuIxYRBqKwtRX0NVje1VfEVu6TWs81F1BR41KmTeJGbnvCl2NV1CxeErg
KWJxvU2FBldWUACH2IZG0h8OkcWnVZ1IvbHqvHTi5A7P3wfdyXnsH6usiWVv
DvwS8FtBYWUBBa2k/1y+xamzBK+g3XMBIxTiDE6jX6MvV5CS91KHLrFiFP5Z
nrCCFmjOEmmkn7v2p6rYvV9BAq1HDYyJQ35Q01mzV9DEk9H0WZKHgaH79EOF
K6hnZ+rWA8Th3o6byytX0GgFh+hPkq+Tghee2f5bQaevd5uaEMcglftRA+T/
cvXZXCf51KiRnpKfXEHPWe0VZkl+p+22Xi1bJPWPSB29QyyueV/Gi4MFjnpD
hRvJv5eLuXC6EAsUXC8QVyHzoU6lpnZlKwtYM/Qd4zv/v3+0Q7RFWSDy2Ugz
F/G/aGWuHgkWiM6wOvarncyLn8Izm5VYoMxxoG19Kw26hf79eWDOAoz1aiKT
ZP6pdOg9LbFkgfSwKbydODS5UGPzZRZgzf/krUbm5/ETn3CaIwuEvVu28qsl
+8QhIqXrAQtItDsENpJ5fCbXKOAMnQVWVxqOFxfRgOZfphaRwAI72+V1ssk8
n9E/xdL1lgVGM1dw3G8avOo5dOv+RxaoFDnXZ/OT5GHtxssffrBAqPoqhTyy
P95bVYFgDzlvRN78C9kvCfb3hnQHWIDLzX6nUya5n9fFKP4jLIAtqcIiGTR4
8uDO6NwMC/hs833tkkaDy1G7qI2crPB6blF0JInMkyqXRZokK5i+3qq+n0aD
Df+2vf63jxXa8tDOBCqZZ22/DNbLsUJQ9TrpzbFkfg4Lv/NVYoXnda+KJiLJ
vuD+YXpFixVOndwofecFOe/EhixpR1bQD2LX/POIvA9o5F6yucEKTwy4x9eQ
fexx1pYv9iYrwMKqAiD7+sqFrzb8D1jBQ2jPAO0BDU7ctt4wHswK5xVs5eRv
kfM+pF/P+sAKjc5VoXFk/18RNZQ+McUKx20lOsNP0cB19M6Jg3OssJt3f+i9
k6QfvtFNdi6xwp2mbf0XT9Ag2nzUl5WdDVp884K3HCPzl/Kk6fsGNojUfPv5
yiEaHF33/amqPBs8EWcPt9xFA162/WMnXdhg3RPZwoTZONjyx5DjkDsb8FVf
MYqeiAMx2t1te+6wwcwq6yNPhuPgxPEiDXYfNuhOaFs41xkH7nds6AUhbMDu
W/M4ojwOWieiDE6lsUFIfC9/66s4+NjJmak+xgbXxqZeu2jGgVt9ya/AKTYw
R43+G0/GwZGyZ/WVc2zgr/2L9l4lDnKyhJbOs66CJaUwr5IDcfAzQPyk8/pV
wCUtyPdmE3kfPHD6T8ShVTBa58oz3EoF5l2/4T6PVbChZOePZCcqnNm4WiyQ
lR1emPyIPnYzFnb/OcXwX80OfAWZOj6OsbD01E/kITc7zKYpu+VdjoUP1hSN
OWK9fS80mm1iIYWTS+j+Onao1H2zW9I4FgQW13A572QHqjHbvmSFWKhtFxgw
OMEONRl6c7dnYsAyddt7YW92mCvMPXvsRgy4nzmk8I5lNdgwd7h9OBcNRW8m
DH4urIaVjz+DU5ooYOe2YLxjDSec7bIvCSt9CdGnDi2c4+SCq2WVL1ljQ+Fd
iknENx4uKFBltO0IDoVPm+4ekhTggvhLUemKnqFQ1/3dibmZC/r7YlnPXwoF
Ab+znW+kuGA0Ji/pv72h4PfLoXRZhwtUu+NTLiaHgMtpOvV1GBcMrHl0zzPr
OXinFiivj+SC07eYvjOJz+G5UF/9vVguaLK9H2YX8RxSevdvMEjkgnsPDrQf
83gOff45/kvZXGCmdZaacvQ5XCz8e12/kQvs/Effnbn1DBxlF/lyWkk9J88W
xZo8g3sUkSSxLi7QPJDi1XvkGUTaXe5eHOKC84Xvmy7OBUMN16RZIpML/rPP
VWrxCIbO64KzAqu4YbNzz7pe02CYqD8adpeTG0r41v7uUQyGtW+9yvUEuEHH
af9I8dxTEFkXf/XrRm7o0Yq4kVr/FGRu/+YS28IN+1m6pp98egrK7QMJz7dz
Q7xgve8FylM4o8l/cnEnN4icU94j7vEU/gdAtE1k
"]]}, {}}, {{}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->
NCache[{{-Pi, Pi}, {0., 2545.623309669457}}, {{-3.141592653589793,
3.141592653589793}, {0., 2545.623309669457}}],
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output",
CellChangeTimes->{{3.6805878908731318`*^9, 3.6805878955310802`*^9}, {
3.680588090829983*^9, 3.680588113441083*^9}, 3.680593310391466*^9}]
}, Open ]],
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.6805878760556183`*^9, 3.680587876942367*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"Arg", "[",
RowBox[{"vel", "[",
RowBox[{
"845", ",", "x", ",", "945", ",", "8.1717302*^12", ",", "8000", ",",
"14175000"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "\[Pi]"}], ",", "\[Pi]"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.6805882185428295`*^9, 3.6805882188861585`*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVk3k01vsTx/F8v5/vo54I3aQdLUKl0CLpFoXI1nJLWSJLaZGKG+5PRRJJ
thJFJSV1UVJSfN72PfsSibq2cEtEJXKf3x9z5syZc+Y985oZebtjFg5iIiIi
vkL7v99qKKH7U14cNKBYIiKGw7p07wH+XHF82Ht/ou06B1X5vpvTZ4pDJHjv
auUoDhI/ir6vkBZHeOhn+xcRHKoTzyY7i4nj7vrV4rnBHHbzh+UaPvCRMoP2
uftwcCpqGUiN52OHcqiJ9n5hXn1rbHYcH+qzO/J0bTgYxGUYlcfwMdXKPdjA
ioOyR+S97nA+ik/YDBvu4fBpkem+uX58OHfwR5TNObj75xYFHeDD532AyL7f
OfhtfhDruJAPNY/bKvNnc0CQp9k8eT6c6jpa/eU4jFUbiTXN5sPaRudwz3QO
J6w/ORhO42Nz2sOsGCkOdh7qS1V5wn6HGtxBOGx4kPXySzuHx7uGTHs/E3h+
vnwk6S0HQUGK8WA/wTNN23n2TRxe8+Sshj4SLMsV862r5JA+5axq+z8Es1v0
tz7P5tCq0OFp20gwKqhp8r7B4aKDlWLTK4JV2+MDNYUcJ8Xm11i/IHC7flL7
UziH5UYqvS3pBB8Xyt6yCeIwVXWuZ0YyQZPOPqdNnhxcFmd3fY0T6rt2DXO7
OcwcNnfRPUMw+Oz5fWzncPD0nVYZb2E/4wF7TptyCGuWDmzwIEi4qJLVu4VD
nE7gmhXHCCLuuPpVaArjovq2FVYEx+tGpcNlOCi6hzqyqwiuWphaRkgK9fK9
BxRXELyqir8dOVmoVx/4l5oqAakwVovicSDhoQVSCgQ3CmK33RwiuOijsqhb
QJCrOxQRK+RW56bx04kj6M7RfxvXRzBc1fuwRpRgZfbAoTtCTnZrH6rZj7Ao
frYp4H4tQf9zYw/rdyw+aVyrTKwkEDQeGJvVxEImrW96UhnBKaOl7VnVLKxT
IhIe5Qm5HLjp+zCfxdD9rtzHaQRyIq/+/u0Bi7nXA8dfhhOYip/+sNGFhZ5s
m15WCIGBs0tilR2Lg5Hql7KDCFbHNxuvtGTxNPTtzBxfgtawQUlnAxaGgcvW
FJ4gON9q+eKoAoujfF+fIiE3pxtmvznOYBHh31hY7EJQKiK7UEGCxbtzZ3aW
2RMcUb6U2vadwUmvGreq7UKuJQ61S8sYRP9YlFltKuQ6q6SATxlQDy/RWiOC
PLc5SyKeMJh0csGVel0CGyeb1JQoBnGH3f9uVidQTD5VGWXLIL+39GvLcoIF
Kqqbr5kz6HWep92qQlDhwruutImBhkNxaZsiQbZ4r8UXeQalVjN7OmQIItfd
uxfUwoNpnPaGUUkCLb2c0qFiHurara9KCveUZH9BUJLOQ6t9vJ4WQxCSIDt2
K5gHu4SCGFMRgiDJ1o7IP3no6uoePDDG4pqdmuVbOx4GDqreDvnKwuvN1ws6
mjycemjy/e4ACzv9XbGGc3gY7Xc1zexn8eTdom9+DA+Ma/p4RweLvftTyxyr
xCDrrmOpVc+iO9J6XcY+Mdx4bvvYVLhXVzmLoCXrxTD/xzm+QwWLeLeSGadm
i0HZu+hZSAGL/T4fw50bRZGS9XFKQg6LicrjMSNPRKExMdkhM4uF7h/nVqgE
i0LnnJlMZzqLcyoH3CzWiyIvz+3Q6GMWf95qtFeWEYUBG5kjmczixWXVtBPd
IjAPeHNM6x6L2LS4XUpBIjhw2a485CqLfy1GlWu9JujHSj/FhDAWlt/VFdYI
JuhRqfuemZdZ2F65MG0o6hf1iOhT6vQX1r8U7nXs/jgNjD55XsuDhaxU3oR8
9E+aes9fp/MPFjevap3+9WuEavYkRo5uF87TuOezt8MIzVxS1i9pxsJ34P3d
f8uGad4jyRgt4R2qaYilmEV/pQ1Por6FrGFxped81gyDQbr3a6ZJgoaQ57Ro
89kvvtB2zdaETDXh/ItdxvpVv9DeDPkdnUos1BsWxdbOGaDj2UmpWnLCP3o/
PynavJ+WNy8X/JzGYmljilzdzj56Y+Sp08upLAbfX9ySbNdLtZfRudp8YT1X
xSDe1R4q2KrnOcZj8XoiYP+GV920xaGk/tUEg7KEw8G3e7uo583aS+u/MSCP
GONf1p1064vdPeODDMrnr1D1udtB5epbdbM/MdgfZpKi+eUfmiHoGdXpYjC0
/UGhRtwHOuI15vx7PQNm8OniSVvaaeG1s/ki1Qwy3o/oJHq10cg0Mj+nnIHS
wG0J31PvqEbf1MaN+QwOT3m9QMP/LWW4qyvFwMCj+fL0JTdaaJ3CrMu5LxkY
r7M5diezmZ6wXLhZN43Bq4RQvzDJN3STe9ItXgoDgy/8LOPNTVQ6bPlYXhKD
PTI7FKXPNdInJWvT9O4wSNmYp31KsoEOrDGR3xLKYJuydc3p47U067W5pXMw
A9EwI0als4aGuO5uKgxgIN+TOPFtYw1dlmgXon+GwfsP1wL4C6vpEVmPcYNj
DEZObt9mu/Y11VYf2y3uwsDR51BGvX4FFZiefVriyKBqZ1PCfKty+sg/yGWr
NQNlk+FcuZhS2jsc98ZoGwMrB/uOIKsimim1UFNgyMBXc7mTX0whvbg06Uq5
HoOe5tq/XN4VUCWHp/rbtBlEOCl5bXHLp451xekmKgwcAiaT4xU5dNXANinJ
xQzuzFBqmLIsh7KC2sOVCgweJ5+4ZBYGele3VdFsJoPTsYH5F0uz6V//K8iU
lmDgdMhztCPiJf0PG4EOKQ==
"]], LineBox[CompressedData["
1:eJwVzms4lGkcgPGZ9x3MZkpSJFuU2oRKVyrH+T+SQ5m6SipUq4NUIpSVtqxm
EikWUZtRpIwa5ZBTDvE85DitGmKSneQQyhYmXCKy7Yf7ur/+Fh/ycz5CMRiM
/T/6/7eUJwTuSizS0au1qOd4Md7tTgUeZrKIKcsuZ+dvRbiOqbH7t3GacE9o
aMwWFuKs7aYLEj/S5PVjrpmlVh4+mmalereHJhm3S3RiUS7W+7ppMqOTJnmp
NZmZvo9xfMqu9qdvaLLa2KBALMrGwQNBdzvqaTJUqPuPWpQYr9kYGv+xmiYO
dT285HkPcP/1iLAvFTRJ+xTvsV6Ujvdb3/RildAkRhwz2N5yD2+MKjY0yKDJ
fTXF93LFLcwxnMz1u0KTsc/suJ/uCHB1CCstOJwmNjrb9eWq5/EfjZzrfMEP
/9bP2tljp3BYvGa66u808VxmP7246225rJZr+cSbJglWm5fVhV2A1CP8Z9+3
0KQo4e+lR/VugntxZL6KA02Smx6DKVMIc2deE822pQnn1ZrjAfIkCM+/G7HE
kiY8/acL8wUp4ENXOdkb0aTnrV1iqeM9WJ+q3BytShOxqbeaeZ8YLF4G1Nuo
0ET/+7dki8oM4E7Jy0cpmhho743tEj4Ee7c88f4JimSuXb5tj10m7FE/ELry
I0X23v7zwyyUA8H8YsOGGorwtzn2mDvlwfnspXr8SuqHlwyFP8iDC29j5q0r
p4jftn1Px1n5cNnci3G7gCJNn40GD5Tmw03FHJlPGkU6raOq7y4qhOKDPnyO
gCLqtv4C9coiKI95HURCKFJmHNERqFEMlWUbfQLPUsT2lBt3nWcxSLS198j9
KWLhpLqilS6BtsZq40ceFJnxoC1JYVkKEza6rU7WFElXVnQI48rAaknTqqtf
meTMTXly5RCBOxqqpjtHmOSdo7She3YF0EqbzHWGmMTsyoUFMpMKkPQV2D7s
YxK+smG1iX8FuGYmuj6XMQkvdLhA+qkCAs0OClQLmMQ5PJi1qb0SHm0barka
wCQbVwWyikVV8PO5WSFR/QzC44w3ubbXwvE0idbtXgbxtPJeNjlSC4UN4bmZ
XQzSeIonj1atgx16030v2hhkt8Ev3BCzOgivHnRWf84gar57F8bF1cGgWqPB
jYcMIjt6IrUA1UNlWnxzig+D5EvtgrQSJCC9tcbV6+00Dv95WccnzQZ4eSyv
b13pFFYb/8bN4TRC+opGs+1l3zAjbV90uNUrMBFK/U/Vj+PpbJdonnMLLDax
btN8OIaTY4X2hxNeQ3+JnZLy41HMu/TCxKTkDeicKD14JmoYB2ydp3N/qRy+
S+aIhgIU+OKBhRsiRe1wGjgux8cGsDD33/qWuA4oWttARgv6cbTH9TEHdhdk
q9fkaJX1YusQjWDLiG44bZFIc/268Gu0YWfjzh6wuOTh+MFbjteXcayDd/TB
4LHJmNDLzbgqrH3yrM9HuFiVZAFfa/He92OH18/4BLJSq6oMOgdnpW7qlsoH
wMPYfsolLgtwooa0Q1kBnY6uu04m10FSuiy2TGUYRrIhe75lMzSf6477bDwK
baJFJxVSOXSv5Pj7SMbg3MlSx10pXdCs7f7ojxsT4Hpn9Z5b0b0wYn6s3M1p
Crwi6ql30n64pv/Frk+LgZ6FGRlFPR8A3au2Wc77mOhX9Z7KuW4KKN+QcjTr
PIWGLWq5nUHDsGjUv8YmlkbZfjPHD/mPwqggQN3Rj4VYob38QwljwFekLu8E
JcRP1N+slTAOJSsarPQ/KCEnVuiKqshvMOHp0hIarIyyNvu2DoinQLcqydOe
UkFmQlz9QjQNQxk5T24EqqCsQ01Is5eBIgwjL89/r4K89FyMbEaYaMssXwl2
YKPcbpUC+yEKjevu6mMms1Fg6/OaqyMU0lyu7OuTwkYbGmJapV8pNMBLE8ru
sNHTQq1JdyaNvH+vGhTfY6PqyOW2fnNoJPF0i+WJ2ah1tUPjX6Y0Esne6Qny
2cjzUvPhyWAaZdQJgookbPQfaTzdGw==
"]], LineBox[CompressedData["
1:eJwBkQFu/iFib1JlAgAAABgAAAACAAAAwHBK/y19A0AQALS3EswIQINlDhFz
gQNAQmUx+XLNCEAIT5Y0/YkDQOiKqagy0AhAESKmexGbA0B5w2HIrtUIQCTI
xQk6vQNAsmnQgZrgCEBKFAUmiwEEQKnJs5lD9ghAI9FFoocFBEATiIpihfcI
QLBjmEj4CQRAyqUPkOv4CECtBwgXfREEQFm55aNJ+whAEPsKCG8hBEAV9vfQ
TAAJQNXhEOpSQQRAIluBBksKCUCunlFmT0UEQKcSkhKKCwlAhluS4ktJBEAt
jq32yAwJQDjVE9tEUQRAGMzYSEYPCUCayBbMNmEEQMR72SE/FAlAc4VXSDNl
BECw2xf7fBUJQExCmMQvaQRAPzwPsLoWCUD9uxm9KHEEQPrtBK81GQlA1nha
OSV1BED74vH5choJQK41m7UheQRAqSJ1IrAbCUBgrxyuGoEEQHXaHg4qHglA
OWxdKheFBEBRNzbSZh8JQBIpnqYTiQRAzarFdaMgCUDusyToTYkEQFFsNYi1
IAlAdGuQjA==
"]], LineBox[CompressedData["
1:eJwVjX1MzHEcx09+v76fr+hOUZ56wJyHGvdHUh66jzXTE5KmyUO4tZGnW0V6
8FCbSvNQXE/koSTXEZ3VLsvsY05lK6SmrNEDrahERdzVTf5477XXP6/3/P3H
QiKtJBKJ+8T+M+Pr6fqULAH13oWlla6cCiL2zuXZAmZEpFrfdOEUO68wV8wX
0Dc3rGOfMyd5zoJMy20BpU3DA1/ncrKEqKaaiwW8/KfO9eEcTi3S4vRRrYDV
um3q6Nmc0s8vSv5RLqCgvKYyO3DqS1oS3UUCNq4Iumhrx8noffD7J6OAl5wU
eU2yif/fpVFtdQI2h3VW50k5bVa7qZrfCnj4ZHngwmmcylXLQ2s6BCzLmXrG
GzidCPDw1ElEhG7Di4MWoEHfe0VXBBF7+mP9PceBDqybI02ECW9Y1241BrRT
IekNlIm4+Z3Jq+AvEDrU5w+4iHjc1Sf3zQiQTZdqXKEU8b3odcmrD+h2/FWj
4ZSI+p8Wo7QV6JV2KOSlScRDGrVfhX6iH2va7jKFYWeU9a6eXUDXN3iYtjHA
SZ4meTQA6crCcp/ZAJace/f8izVQ1cxEj6UywAdBp+K2i0Dvu58ftcwCHA1v
l6+yApKlbvmsXQaY8zhn7YCZUWpNVP34JkBl5K/QpYOMojcW3rynAYw5b+de
28wo+ZFxjV0+YFj2h9UOTYwyHXtbk24AanRjvapGRmU9y+1DSgDB7aPU3MCo
N+1p2pgBUJYVb3asZbSnrkW9tQ2wVOtYxJ4wOqIwT3vaDtioCvbwMTBKynO+
L/8CWLXynG1MJaP8A5Hd5n7Ak+vbklv0jJpheEeJBTB83uJHGTpGn9UOo7LJ
HI/1fax/rGU01OqtSWQcUyoM460ljKSlZ18HyzheSFAmudxh5Dy9+FD1DI5x
FFqtLGTkHl8L8tkcRxThP3bfYrSm89vdTCeOO4b19gk3GAX42/qa53NkfgVO
muuM/gEyAjse
"]], LineBox[CompressedData["
1:eJwVlHk81dkbx7n3e873mkhRKjUalUayTYiQW1mKcI2akiZkXEuUG6KxTCoZ
YVwVpWxRslSWZOkW3wfZQrIbW1rI0iBbuVl+398f53X+OOc5n9f783zOI+/o
ac1liIiI3KHX/3cz0+WG3+XFYGDnSXePayToFQROsOTEoLLTZbdiNAnK8qOJ
MrJikPv0w5YPUSQsn6v+9ouUGOg8q916JJKEpoyL2a4MMeglNsaq/02CDWtm
Xft7Fiz19fXnB5LgUt09kXuPBU2BwtftzvS5hllSaTILYhTs4k5zSTiQXHyw
Pp4FGxyCRhhOJCj5xT74dIMF702jV247ScLYVs7vciEs+JLwc5njcRJ8Q8ur
I5xY4KY/NHGFQ0KIcWaSswILAoJqPI7sJAEi/K02yrPA8a+YnkhNEuabDjI6
N7CATJNxgh0keNuNcU1XsWC+azxYTo0ERz8NFWUmC8SMy8wLfiaBnVny/Es/
CdISMwUGa0nwH486ndVDQp6s2Qa2DAmFWg4b/+ikffKfuqK3igTVcsbl1kYS
ErcYSyiuIGFD936zolISRJ8e5zWQJAjFmzsDE0hwDYk8rzGLYeehe+FacfS9
PeppQ1MYvG776I/dICH43TaVuC8YhhXW3LWPIOGHLYPVA58xdBr87rLPn9bZ
46Or/QFDIW9whrQhodhlUw7/NYbJwqJ0OESC5dIJs846DKoLYcf+pH3Z/PjR
jbW1GNKubi8ZMSEhTajOuVCBISaVF9KgRYJ7w97GviIMZ1uFUjekSQhalueT
mYThpjXHNkaShKJz1Rp/x2N48eZeSuwyEjJ1ubnH4jDgBnP1OCbd15TkO93X
MCRUJlkk0hzXQ6ojTEMwlBtOxSSNYxhTakmtCMbwqWx/T/IohgCt0mi1IAw7
SidOpdIc9t4qWt3nMNQU7gtLb8FQFLXOz8eZrtO81ZjRiGF6lYk2zxGDdP6o
TBbNlb3aq+2YHQa7nJi0RzSHdoK64fRvGKbSB8vz8jHwBQ5BGkYY5G6HLzy/
gWH8YfJ55Q0YjNa8NSrhY9jjqTT1SAaDW6xGZGkErdOzOkViJYan13pkyy7T
PjIKL51CGEzDVXWqvDGEenztO/UZwRnW5QvVnvT7pG5X5QCCmNCOqhp3DDIV
pQPTfQj6LgX/VvcHBtdH2Q7VbxD4BDR7vTmEoTtThS/yBMGdua2CJg4Gmwfe
x6lMBJRfgGjLQQxKm9o12CkIfvDZEt1miOEW5vuf4yNQmzzf0c7GYJgxvEEh
FMFhXoNcpx6GLq98pbBABMkevo+7NDDEK/OMbVwRvBx5Nd2thsFqRtH/xQkE
I64b9Xu3031zUDtRbo1Ak1vz6u1mDFXuhbWJeghenZAd+iiNIdr5W9sAiYCT
rM8WStK5iWCSM3MEtPbb3ZQUx9Dc/GurzygBvX/cM9IlMJi/04lPqSfAMa0y
niNC6y/es9F9QcDg4KdJp3kE97NrpNY/JGDCTTmFP0375MV7dzGUgHMPLb/d
n0AwIdfYNepFgPAzjyOgfb26/Ygd144Aglew8PEjgjmmQjV7BwFX8zoOC/sR
OCykvD0vS8DyqbmHkr0IDm0nIs+KErDG18BWtw3B+35BkF4dExKKHPI4TQi+
OHm5tjxmwk9zl1jcBgR582fTi6KYkKabZu9fi6AX80VfnmaCUmB1Ib8SQYak
6IM3ZkzIKRmWSCtDsM8yvylPgQmaS8u4ghIEpMCCp73IAINLVtIDBQiWZat0
HU1nQEWF1ylhHoJZzaHwP/0YcADFlklmI6gd6+kwN2JAvUnRWoUsBG9vpQyH
SjDg17B/PXUfIMjSFjI6mkWhvfZ7FScVQeRlt5LBG6Jgu0xOjpuEIKRjxz1b
K1FwinKs599EIN43vW38mQgMN4ZsTrtO82fydls4i8CZlen+gigEComSw9Hi
IjBpXdvUGI6AoVus99pmifKLGVUcoHMTXD22TDZikZpvkwgWXkIwI1Vnl5uz
QAWvUe+QvEDXc8LGL1TNU+F3fK7o+iEYMEi/w6kXUpI9N3s43ghy2CnPlEvm
qJgfn2lwPemchoV9lL7/jUq8u/CO74Lg8fUvOrInZqncB6EGA0cRJK6vCMir
+0JpDWXECg8h+CxX3b61aIISbKv7LGmFwGC1uN1723HKwP0/QwVzBNpvDm1b
tP+PqngkGa97AEEXQ3D7woVR6sDYL5McIwRPar2k7QTDVL3aYVPuHppnNEa0
YcUQ1f4k7itfB4FleLuh8fIB6vi0wDJNk86XVfB6w/wPVL9Wb5pAHUG6Z8iL
757vKa7f0nyjMu1vsxsxue8dNVIsf3hAEYGEKtMkWamf8hQaPhRuQWB/Sj1V
5GAfNa3nzFghj0Cj+783+k491EJpVq7uOgSdf5okKs12UvVdauLfVyHQWQyJ
djPtoBJmn7o8X4Fg757heWtBG+UhpVsRII7gbhDzpJlJK6WvSsnpsxDsTpNw
nBhppsTNjPznmQh2qk20yT1oorq5tW0vlgiIOJZ49OO5N5R/Ykvk7q8EnLEX
USmKq6fMntkMLUwSgFWNxVbdfUWta+s1LB0jwK3e5vjFlzVUsfiQ0GCQgJfp
t7d7H6+kwhRPH1l6R8Av1Tm3p2orqKNGk3lULwHS532ZZy3KqdmAedc9bQRk
IxXWZVmKqrp18aVIEwGm7+3WF//zgorNxz+V0f9c2tEheHajgNIcXdGx9yUB
WQnaX3+6XUAR5M0dDCAgw83T+EBgPtW6aX1U+XMC7E9Svw365lHetgrGhvkE
jP61PEf92ENqn2/WXWYOAe4lVy+LvkqnpK6rzVdkERD0avfv2P4+9aR2V75R
KgGSbdqr93fFUxM6lvIm1wg4X9dS7qLhTZW8/tXW9R8CjGLVDI2O9JXyeTad
VWEErKTsw+R2+bJVMxz5+4MJ4PHM9x70j2KfXuO3cMCTAClrh0UPj7tsfY15
GzF3Aiy4W8N2+6ayxTkXn9Y6E/D4b7Pkr7z77EehEe5m9BwKcjxz9bNeBntk
JvnfgxYEFBzdWsjYlcMWrFTQEjclYDEuAzeezmVfVcmKrjciwPnMyU89u/LY
ityn+y30CdhrulZDpuIJ27m1psByOwEN3eu8TiUXsHdOWKyU/JmApX8c+/XN
C9lIvMWjcRMBSdiaMpopZN837N1sRc+5H9XdTK7oF7OD/qoUSC0n4N+6wtET
88/Y/wPlyLo9
"]], LineBox[CompressedData["
1:eJwVins41Gkfh83vGfrNjrqWvDSxy2TXbrFRa1vn7zMVjeqqHGqkNDns5BhK
VnmlkaLS2l3aRVLUeFcp2l7Hwe9xKscMCeudKG1Cuw7lbNJr/7ivz3Vfn5vv
E+r6HaWmpuazxD+bqTEf56nOJr4bNxZ3msphnycV4ctikznj7KDO+3JoYK3c
d2IOkXiT60ZMcQXc22O5On0YkWn/PXZOuVUQNRqZ87wRkZDn6+0Xdaphw+bY
lOF6RJ5ZVOVHOlfDyJWE+LfViEjkQ5NmMdXgZZ8mYZcj8vnuIKH/QDVsTipb
9+VtRBJDzCxu3aoBzXWq30MvIuKU1XpwgVsH9THsW1HnEWmUuAWe2FQHp9s1
r0jjEGGpZ/WHHq6D+BTdXO4pRA62z4cLS+ug65GDbUkgItsO79eNca2H7O+k
tYvbEbnZ5RUwFPMQNmVrdF7mInLdNb8vLrUBbNrCGwXLELGkF0/yixvA4b2y
aopCZEY/7dmJ7gZw2v8gz2ueIpnSy+udeI0g0joc+9UwRT5ae/XD3fRGiJKW
rWt9SJEsafOF2h+aoMw7WKoZR5HtOSV9xS4tUJXcHUliKKLi2aYEBrdATeXm
4IiTFMnpuOQ2fb4Fmng8kTKMItp/6DJp8hboba83yxdTpK/giVLfuBXmBYY9
O+wp0rFaS3zrTSvYrelYf2mWRWQyr50Bkja4sZJr6TbJIuJFdurjf7cBUt9q
rT/OInf8z+p/ktIGTa+Lttx5zSKRiQGO3zNt4HE33aO5i0UKz+39JZWrgAgr
7zhuEYs4cO36fb0VkL9r/OmlcBbxlicf89doB4PoFTFJI2pER+FiMmHbAYrM
DR6SZx8Y1b7nLht3d0Kb/4PX38jfM7bXFg4YhnRB7tp2qz2VC0w1daYs5WAP
WGQowo41zjE8Ht+A79gLfAv7Xt07M0yeh3JsyE8JI+WO6hr3p5hfLIqG5Sf6
QD9I7v190jvGrW6b8GXAc1hs0paNh08wHZYp3JV5L+A4aLoHzIwy6vEXDYRv
BqD061YyVTTCeD6TRRkF/gkFWg8L9SoHmbaJjZVrpl/BcZt05BA6wLQmrg2J
6H4NNufEwqFAJWNi/GS3aGIYxvxVybGJnYzJVhafY/4XnK27agOzj5gdrp/m
LKhGoUtuV3cbFTLmywW9TWgCxGZO791/ugfOxJ00VLyFF0KPvUezGkBSuYIY
vJqEyQIoWGXbCWm1ZYFR5jPQK/v06IRCCfYd2/5MH52D6KNy4d7rA6DS2yva
1aICjxvmoszLg1Cg+9uvu05/AElCI9WvGIG8ErdD4s9ZuDbe1DSpeRSMPLnf
JosofEjrVY3O/glInc4tUfNF+J3NI4cXke/A2qR8uEjMxgWhy+d8wqZAEJQx
p2aqjtmxg1Kf1BnQqqFFFUPqWJpu7KyXOgcxN485NUg18A527Nq6Cwsg+ZfJ
VWv1Zfiec0jPaN57ePxxxBFF3DJslcHUP5Z9gIbTAUjnzdLv04F1B9Vwqbji
SNkeGkuM3E0Fkyxstkkx4XGHxsfHTgk2zLLwisLb0Zn5ND5TlS3iq1g4NGg0
p/8uja8eHDvLYlNY86lxs18hjRVpF5XVKyk8JTyQFFhEY2ut6subv6awMOGj
NX4MjTWp9eNbji31HZVa2k9ozGt317CMpHD92f/OOnfS2ORGtMFnpyj8o6FL
ovQpjQUOjUJ2HIU90B85Y900jjzll133E4UXOLmPHylp3P82w9XxPoWddbp2
+g/S+PeXy4qcximc1ddlHTxH44ie5oeXJilMbzF0ejBP429bk3sUsxS+lrN1
dm6BxhXFeipPFsJJfm27zy/SuP7CF1tCtRFu4cV0ZyAO7jHf1v6rJcKVr36W
PVjOwRmfcV8qrRD2LLYSTK3gYC9e2yTfHuEYmcuFTR9z8AAlWpXviLDsymmj
Um0OfvNUIib7EC5PMI0s0ePge03rwtUPIGzV1H5mehUHhzGjcdvFCO8UmNl+
s5qDJ3+LzO08gvDW2kxy34CDi6/ZlPKCEc7ixSWNfcLBJ39ebDwUhnC9XTnn
K0MOtkuo+d/NCITbZflWgUYcvBh9/u+hKITNKgnvP3wO9jvX6atacs+cgMr8
Jf8/4nzerg==
"]], LineBox[CompressedData["
1:eJwBgQF+/iFib1JlAgAAABcAAAACAAAAwHBK/y19A0Dvfir7EiEJwINlDhFz
gQNADRGdtNYfCcBGWtIiuIUDQPkYuLqaHgnACE+WNP2JA0BMahQNXx0JwMpD
WkZCjgNA8eVKqyMcCcCMOB5Yh5IDQHp59JToGgnAESKmexGbA0CN2QRJcxgJ
wNQWao1WnwNAQLqdEjkXCcCWCy6fm6MDQMLZDSb/FQnAGvW1wiWsA0BnZ9go
jBMJwCTIxQk6vQNAqW/YkqkOCcDmvIkbf8EDQOmz7J9xDQnAqbFNLcTFA0Dh
D6LzOQwJwC6b1VBOzgNAfVpVbcsJCcA3buWXYt8DQDG4o57xBAnAShQFJosB
BEBh+gmfSvsIwCPRRaKHBQRA3LooUyv6CMCwY5hI+AkEQDIx/prr+AjArQcI
F30RBEAkUyW+zvYIwBD7CghvIQRANO9fTVbyCMDV4RDqUkEEQOO4TGhv6QjA
YK8crhqBBEBJn8ZBx9cIwPCzJOhNiQRAg8zuioXVCMABC5uw
"]], LineBox[CompressedData["
1:eJwBsQFO/iFib1JlAgAAABoAAAACAAAAdko0NqoABUDNddjI+rQIwMmo8lLF
+gVAot6UyiZyCMBodPpWHwoHQACs6qs6KgjAm0U0+lMHCED636RmU+UHwKqt
V5e+CwhAhAOtxRbkB8C4FXs0KRAIQGcMebjZ4gfA1uXBbv4YCEDYNjpSXuAH
wBCGT+OoKghA6ncfI2LbB8CGxmrM/U0IQJTOj4tS0QfAckehnqeUCEAFvh6I
x7wHwIGvxDsSmQhArMyOanm7B8CQF+jYfJ0IQEil/pkqugfArecuE1KmCEA+
VrbUircHwOiHvIf8twhANoIEaUKyB8BeyNdwUdsIQIntImKLpwfAbTD7Dbzf
CECJpL6tMKYHwHyYHqsm5AhAYDVHFNWkB8CZaGXl++wIQDVBwiMbogfA1Ajz
Wab+CEDdpHDkm5wHwONwFvcQAwlAJ1/pnzmbB8Dy2DmUewcJQDXM9lfWmQfA
D6mAzlAQCUBy52GtDJcHwB4RpGu7FAlAAuNaQqaVB8AseccIJhkJQDBWHsM+
lAfAO+HqpZAdCUA7Okor1pIHwEpJDkP7IQlAT6todmyRB8DNBLvR
"]]}, {}}, {{}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesLabel->{None, None},
AxesOrigin->{0, 0},
Method->{},
PlotRange->
NCache[{{-Pi, Pi}, {-3.141442791943694,
3.141442756448936}}, {{-3.141592653589793,
3.141592653589793}, {-3.141442791943694, 3.141442756448936}}],
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output",
CellChangeTimes->{3.6805882192654424`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"Arg", "[",
RowBox[{"vel", "[",
RowBox[{
"826", ",", "x", ",", "945", ",", "8.1717302*^12", ",", "142570", ",",
"14175000"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "\[Pi]"}], ",", "\[Pi]"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.6805883670282173`*^9, 3.680588387031139*^9}}],
Cell[BoxData[
GraphicsBox[{{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVUH8w1HkAdd/9fj+f76rJpCadjJ2VuCk1SlSXOfkRh7QNQzRCkZa68yNR
cboiCf24s2vtlUUS20+dIwk9W35MdEqKQz90JZ1md53T1Za72/vjzZs3b+a9
N0+8PSFwB2NiYuJpxP/s5zvL86NYiARFfPSWBoq1dRl63lqIIdHJ8Nn1FA7i
8ZJ5lkI0/1unbqulmPWh4/1ycyFinsmt7a5SPKg+dEXKCLFwKOJQWyVFKD/1
+eMXPG4U3HJM/ZFiZ8eQvqaCR4HP36mb44y+k5+qpZTH0Pkoej6W4uvSBv/u
0zyOZq8p0kdTLE6Tn39dyCPysGDf/ggKrZ0k3DqbR1+IV1NoEEVqjqYjP4bH
0pTiKTdXiuz1alXsIh7C9uVxq0wpkH9gk0jMw7uBX6giFJ8e+DMDVjyeYFHa
NEOxJ0K7w3cuD1dlclLlR4LtaU5LHQQ8smp/EyneEripm29OPKdo3Z1w6kU3
wQHdiW8uDFN01Iz0azoJ6p2jRNEDFB7JZknKOwTLNExWXw9FRbBm7IsmAqsh
H7/rLRRXem686rpAYJjZO5BxhuKUtsbp3RECl6CKPOdiCrN1ykbbQwTJyhRX
bSHFQ9+oTd4ZBG8WWZRF5lOkT9SqtyUTDHwVvtPjAIW4Z7NmYquxP3F0ioZS
KEUl7qtWEPxZf70Kxl+2lE357XUw7pnODdsvocicX6MotyOoPLak+Q9vigXa
0szblgSys4nZ95wpKje6HBMzBEl9BvPCORRTnb2xr7o5FAVKtsjMKNbva+oP
aOfQdL+iXD6DIrA543DZLQ7k3gbHYgHF8EFTNfMzhzNtqoCSSYIZTOR3l+Qc
NJ6TMpWOwLVGPjfsBIfXrT7DpeMEE6NrH43ncFjRoo8/+zuBV+PKI61pHDrr
PXKrHhJ8uhxYbBXKQbtS0VPdQzBe3Xr8jYTDnNrxeRe6jHkBqnKFD4eIq7LK
S7cJwmxivGUuHCarRjXXagkunnvqYzDnYK3Mm75ZSPBr/j/zm+6w8LJ45tV8
0pgXYpi9ppFFnNypoCWfwOap1pB7lcUvPwxbtmYRdDZcW1CiZOGbt2x1+x7j
Xkm8y95dLL7lsw52JBBIw2KypJEsZDn97Z27CN5ZNA/aBrF4evj74K5oAuej
UsHkGhYp6b3J94MI1k2P7bnJsvjpg13jAwnB4p7RE/feCXArLf2zh/4Eg+qC
4DNjApim2J565EmwYUbciLRLgNLdqZcHnQjCc0875x0X4O5Wy7GXcwhW7Gsr
s2cEkJS6uhnMCMYslMn7dQz6nkcUmc0k6M10CIkfZvAkusLrS5bAro/oRHUM
9HEO5Sf/Mv4ce1dftI3B3osb35/Tc3jsvsRW5c/A8DZR0viWg9zGXefozIBN
rJt++ZLDyJV53jrCwN7eszvJqBVS22WrWQb/AYx8Cn8=
"]], LineBox[CompressedData["
1:eJwBcQKO/SFib1JlAgAAACYAAAACAAAAF3Y+WjrVBcCFC2ZciaQBwJy3YLFP
0QXAHsZhz6eKAcAg+YIIZc0FwAWlXVY+cAHAKnzHto/FBcB6tqXHyjkBwD6C
UBPltQXA582O9CjGAMBkjmLMj5YFwBIc9wl4hf+/6c+EI6WSBcAU1cTUPT//
v24Rp3q6jgXAE7Lh29T3/r93lOso5YYFwKkPtaCCZf6/ipp0hTp3BcCs4Vji
YDP9v7Cmhj7lVwXAnacvblyg+r/9vqqwOhkFwOiZ7W+RJvW/zNovDPsUBcBh
YfYnH8n0v5z2tGe7EAXAlYHWe2Vs9L88Lr8ePAgFwBTupEFttfO/ep3TjD33
BMAim0fQXVPyv0q5WOj98gTAjm9Mcaz98b8Z1d1Dvu4EwLuyER81qfG/uAzo
+j7mBMD2xKrZGwTxv/Z7/GhA1QTAxEHNhBGU77/Gl4HEANEEwLLo0Jb3/e6/
lbMGIMHMBMBcZ9kormruvzTrENdBxATA2z+gkJFM7b9zWiVFQ7MEwI4Yncr2
Meu/8DhOIUaRBMACWY9HX33nv8BU03wGjQTAqytWPBoS57+PcFjYxogEwOpS
4BQxqea/Lqhij0eABMCZ3m8aQN7lv20Xd/1IbwTAyz3HdH5i5L88M/xYCWsE
wIJedQm1COS/DE+BtMlmBMCgbeD14bDjv6yGi2tKXgTA/tEWT+0G47/q9Z/Z
S00EwLGsug55yOG/uhElNQxJBMALHpuGEH3hv4ktqpDMRATAO5vaOUAz4b8o
ZbRHTTwEwPmn83I9pOC/ZtTItU4rBMBqDsmcIS/fvxZ1HBAHDgTA7F9l3Qby
2785ez5T
"]], LineBox[CompressedData["
1:eJwVjn881Hccx+/7/Z67NfpBlFR+1LII+Vnh+FxR6Zey/LiZrbLS6EJRNKsL
TYlWTVpKkd0lzq+Ucoz3JxzdSe7sCCPdbdHdtnJJKLrZH8/H6/H84/V4PK3C
Y77YR9JotPBp/t9cxvuUUD06ZoVE2Gwz/AOCQ8n4bwk6Xm9huzGptBceEXOD
j0xQOPuJbbXx3z1QtsPVLEdNYYNw8Y+zv++G/XyWfsELCmfa8NR3rLrBctx3
slhJ4c88V4e1S59CVl7Qs996KNy43WRRiNVTSHx1tOC5hMKHLnjlFj3rBKd1
vCy1mMJ5X3FWOl/sBE326VNvHlL4WmxJxIhvJ3ztdSWCXkPhVMGBJrJCAesy
RbbLiylsVHN3a3bW72BgO1kZc5bCfx+LZllmyEF8nM5PTJvuCbse7eAvhxNy
g+zkFAo7KANEJkZy0CYuTsj6nsLaQNPRsFwZdLV4ez6IorBZ4Sel/Ph2uLkv
ufHjZgorDgZsL3d8DKGi9HvMjRSOEufP+ORJKxjP/Fkwx4fCa5wzpi5zWyHt
XsHpJZ4U9i7w1UelUuBSTVs2rKDwCqvVMcYsCay6yVCc06fwpqWKk3YZzeDR
fkiylklhv10q91TnZvCe6qsfJSnsuu3F+796xbDhy7tFX78ncXrd3rP3HMQQ
YribZ68mcYnDzs4xZSMkJots25pJ3ObrIjSMewg/lH9mmdxA4gUxwdeVyx/C
yf7zJm71JH4+Eao6MoDhjHsE7XoViXuMyhz7/TFc0Rp1cfkkDqaWDZ/n1oNo
DzfZIIXEG/OMooo31UD9+adH8XESpxWU1MikImioW8eNPzb9T7WmWWwTgXTB
gpC+WBInzNt72yKoGnrlYruSXSQ2JxZN2R2+D+/XWnRv8SLxOV7kS35fJbCW
dDhkjBN4ZYNgZuad25A/V99151sCv83Z4p6vLQRKz9d94TCBn+mnrp3pUgjS
oSof4RCBBa4/hT0TCYBTmsNp7SIwQ7u0r7urAOLX7EnRryKwRYvx2aJLV6HE
f7gz4xCBySsOmseqBJjDtvlj50ECM9tsK+wmD0OcU/jzhZEErh54N2HqzgUP
Y4VGuJvAuguyKElLAHh5a35tCCPw6X9p6bNkqvpHPQ90rdsJnHDLeudc64Po
r708GwMHAs95fieL+SAdLUqadTxTQ8OEtYo5y+MGiuRL518fpOEPZpvjHF3y
0P22tMpSFQ3LLm0U8k3yUYClbuhJLw0XiSr7XA1uojTx6y8MW2k4qyK0iDNS
gF7Pli+/LJz21TYdsuJbqIGfpcjj0nCJybaG5uZSJMt14kT060CvtO3GeHwV
ovFjLvcrdMAwHStwFlYhR2GpIvCxDhZb2Ln5KKvQBZFNgE+tDvIkqiL11vto
R5fVZsscHZzKOLkwZMkDJJttxOoN1EHHeuWZT5urkSzljYX/44/w+h87l+0j
taj9u7tDbrVTcMJ7hpvlB4xu2cjX7Kj7AN82nhpN+rwZOV6VxR6WTMB4e6O6
XylFVo5evfOEYxCw3H62Wt6ONDXr9Rh3RsFeomuoWdmBFh6o3ZOQOQL11Qwu
s0yBPkqNBMOHtGCeoZ4SuzxFccggMHLsFbwZ9bGvK+pB1S5teLRKA6FHVgXv
4/WhcsPmivl1g/DOPstqkDOA4jxyKO8YFWjebU5wDlMijx93+b2M6gP6aMf+
Juc/0evvJs/zzijgGmtZ/7vFL1Bq0zUPNN4C39xmab86NoS6allNxVQFrOrQ
4pyParTLbsNU4MUy9EuG7MZxxT9I6ccJir7xCJneE+y1cx5Gb8tRuamnAhnN
L1ycbv4G9QrMo7WyPpQdlM5MTHmLkqJr/YLyVIjhsZuKsR5DnPyVIbnnBpE0
WF3W1T6BIk5LyAGZBnWcMXN0kk6ixlMrVmS2vkL79bFfVaoOfWP4osH4Sy1q
k7XEupIEe8SjxVt5dARxfrf0P+BJsstjZk6Ex46ivIthqbYsik3nDSaHXxpD
L0uSDifZ09nJOUs3zb80gbzb5BKLYTp7C51n05T+AX0+GSmsvKzHLtt0sPtV
0RRqueYzYO7GYK+5CuInAh2KZrGU1tUMdll4B3veII2NzXT7a1cz2TOW1jnd
5kx3FC46sT6Syf4PMu8cLA==
"]], LineBox[CompressedData["
1:eJwVkH0s1HEAh+/0+3bfH8pJsUN5KW+pKC+xRR/JGVoZt7XkLW+lbEzHmJdG
LmGsvC3L6dLJyMvSTGKL6AW3mdcUU9wWkd1pGpuc9Mez5+/nsYhKDIrV4nA4
djv8d5y5yN5rjYsI5zLhQAMPbUpeu1CtheuV435juhT++4l1EZfBrbZC0Wwy
xeERn9p8wmDGSRYqElP4haikv1gG0bZ3c9UpFM08apSlz2CXxK2iPY2Cv6lN
Ey0YVAfWPSrMppiY4y8FeTF4L1/KP1FIEdlq2mKcw+Cve8GdqKcU7glhR40k
DAZP5etM1FLss6t5blDAYMvl2j1fOUVfrdkT3YcMMqxCAo7VU1hXWj7QyBjE
Gous/jRRLGfaJs/3MPggVbflvaZI9Xd2beQQyEO0LGXDFCrv+tpShqCiQS7V
HqG44WGsl0EJek1lh1JGKa46chYD+AQ+U5lO/hMUMFRUrZgRDA2udq99pdCZ
j95yPEvg5lbyO+AHRd70ZLzAm0A/iGvSsUChGfeb5PoSaNwrRJY/KdQfHVpH
LxKYP3PQbCxTTDZvRorDCYR146hfpZCll/V3ZBH4CuBMtygEt3eflOUQKMLr
FGINRWlCuvS+hCBOMJs7t00hiYhIvVJMcPy0yYUuLRY3fextN6sJCgvOxCdR
FkrPmnKljOCysIU3y7IIddPnKuQECYuvBgN0WFyyX/9S3URQ9Hi112YvC2f9
viLPboKYF13FcwYsmrVdN6x7CGwEfvGBB1hYMw0xev0ERzrTkt4ashCsl3h8
HyLY3omoEbAoVXMbPw0TYHGhfI8JC90lseHLMYLO82+Cs0xZSJQLuVWfCcTN
786tHGShmQlR5UwTbLWqosLMWAS78KfCvu38zo4QDJiz+AeQnwNZ
"]], LineBox[CompressedData["
1:eJwVjn0sFHAYgO24+70/xJHPFWVxdcIsttaYXl/dapOopTQ2mg0JU63koqyw
tZlK5Tq1CIlwTR/iqhf5uE1a0WUSydeF1hx9zIn0x7Nnz3+PS3xaZILAyMjI
fZX/Vti7Fw9NCHFUq/jZvZFTsqbR129KiFW63sey1R5vk5RtmhFijUVZp9aZ
U4xascb8hxBdBEXlCU6ctE/Mzy7ohTgQXSFbWMdJU62P7FgU4ofJqG9iR05/
ugQtghUhusc3ZOTYc5JMWW9CYxG2ytVPv9tyynXzmW82E6HKcyhJY83Jv/zE
NdV6EdYlzkWdM+ekUi70KQNE2H5IKdWuAJUU/jqQkSvCwIOSnsBJoMSTiwc3
mDJUNlwOq34EVHckX2pjwbDFM+KSnQpIH2SzBNYMDYNb3C7WA2VZeZXpHRk2
3ja/EVMLVFgfN9MhZfhupdGUVwI90XXnpOxmmH3doRxLgAQxN6ufFzA8utnw
0DUbSBnqu7ifAYbOaZN2yoBq6qJuvjQDHL71piAwFKjJNstXKgZc3lorDQoG
0k60pi47AJZS5nDATiBxXvhYtTtgo09podt2oLzO5J6/YYB9cqy8JwHKkJXd
uV8MmPt+tt/PBOhCw2s/awWgPPXRmLEAqMheNyC/DXgj+65as8KobsprbWQV
oBZ3qfYaGOny1flLzwB/jTW/9dcziu3+mB7xCdA7iTveG2F03NuwRj0CeOWq
5baQz4zkJc61knHAWAvX618HGSkSEyYMs6s/Z0xs7bWM+mH+cNUyYAamLMT0
MBpLt/stNuYojLsbpdMw0g/sKM5iHGenNzeldTGyfHC+d594tcO/yE61MXK2
qjjWYsNx0HI5bv4VI4/MLpA4chSFyMJTXzDyG52uLHLimNTWbjLZzGjPbotg
gwtHj3Nhp6ObGP0DfiMuRg==