-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathJelinski–Moranda.nb
4193 lines (4146 loc) · 236 KB
/
Jelinski–Moranda.nb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 241100, 4185]
NotebookOptionsPosition[ 237516, 4122]
NotebookOutlinePosition[ 237880, 4138]
CellTagsIndexPosition[ 237837, 4135]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
RowBox[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SeedRandom", "[", "0", "]"}], ";"}]}], "Input",
CellChangeTimes->{{3.808249437390602*^9, 3.808249468789117*^9}, {
3.808249632175363*^9, 3.808249635636956*^9}, {3.808249765110931*^9,
3.808249766532976*^9}, 3.808254098303404*^9, {3.8082563792955236`*^9,
3.808256381597042*^9}, 3.808256704271912*^9, 3.8082567676558867`*^9, {
3.808257135602641*^9, 3.808257140940946*^9}, {3.80825720411176*^9,
3.808257204717432*^9}, 3.808257874664433*^9, 3.808293009163722*^9, {
3.808293569355691*^9, 3.80829356980971*^9}, 3.8083195815073442`*^9, {
3.816823365459618*^9, 3.816823365783063*^9}, {3.8172867370346537`*^9,
3.8172867391812983`*^9}},
CellLabel->
"(Debug) In[1]:=",ExpressionUUID->"6aa5b155-c8b6-4950-bc31-a07e80a911f7"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "\"\<Sampler1`\>\""}], ";"}]}], "Input",
CellChangeTimes->{{3.8168124021590147`*^9, 3.816812402159195*^9}, {
3.816820205778687*^9, 3.81682020614793*^9}, 3.8172501919117117`*^9,
3.817285432667527*^9, 3.817435500804758*^9, 3.817435604876099*^9, {
3.81744024760004*^9, 3.817440253709037*^9}, {3.820740856282202*^9,
3.8207409008730907`*^9}, 3.8207410047465763`*^9, 3.820741052776198*^9, {
3.8207410968162527`*^9, 3.820741100943769*^9}, {3.820980551767469*^9,
3.820980563997992*^9}, 3.820980617916573*^9, {3.8227306941465607`*^9,
3.822730694659491*^9}, {3.822914077673979*^9, 3.822914078055475*^9},
3.822914843886628*^9, 3.82291516448693*^9, 3.822915747812892*^9, {
3.8229158419818974`*^9, 3.822915843668749*^9}, {3.8229171391032887`*^9,
3.8229171405172873`*^9}, {3.823535536950664*^9, 3.8235355377496967`*^9}},
CellLabel->
"(Debug) In[3]:=",ExpressionUUID->"7eff355e-9bc7-4669-b9b3-db6f7109a16b"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"-",
RowBox[{"PowerExpand", "[",
RowBox[{"Log", "[",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"PDF", "[",
RowBox[{
RowBox[{"ExponentialDistribution", "[", "\[Lambda]", "]"}], ",",
"t"}], "]"}], ",",
RowBox[{"t", "\[GreaterEqual]", "0"}]}], "]"}], "]"}], "]"}]}]], "Input",\
CellChangeTimes->{{3.825717233744521*^9, 3.8257172620865583`*^9}},
CellLabel->
"(Debug) In[5]:=",ExpressionUUID->"212e05e2-2736-4d27-97ed-6652963802c8"],
Cell[BoxData[
RowBox[{
RowBox[{"t", " ", "\[Lambda]"}], "-",
RowBox[{"Log", "[", "\[Lambda]", "]"}]}]], "Output",
CellChangeTimes->{{3.825717243505425*^9, 3.825717249683043*^9},
3.825717867426518*^9, 3.825718205967695*^9, 3.825718399934609*^9,
3.825718487903735*^9, {3.8257186864533367`*^9, 3.825718725696171*^9}, {
3.825718833398089*^9, 3.825718885674385*^9}, 3.82572066558392*^9,
3.825720854712288*^9, 3.8257210674439707`*^9, 3.825721108159718*^9,
3.82572704001877*^9, 3.825727162361059*^9, 3.825727671670089*^9,
3.8257277773379993`*^9, 3.825728321052681*^9, 3.825728960858798*^9,
3.8257296667525*^9, 3.825729736387432*^9, 3.825729823970634*^9,
3.825729945157805*^9, 3.82573039794762*^9, 3.825803141738451*^9,
3.8288261537142076`*^9, 3.828826198339302*^9, 3.828826250520075*^9},
CellLabel->
"(Debug) Out[57]=",ExpressionUUID->"ea3dabcf-7809-47c1-b66b-5dc402cfc799"]
}, Open ]],
Cell[BoxData[
RowBox[{"(*",
RowBox[{
RowBox[{
RowBox[{"N0", "=", "800"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"n0", "=", "60"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Lambda]10", "=", "3"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Lambda]20", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Tau]0", "=", "30"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"NN", "=",
RowBox[{"Length", "[", "data", "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"INTERVAL", "=", "1"}], ";"}]}], "*)"}]], "Input",
CellChangeTimes->{{3.824797877729094*^9, 3.824797890266142*^9}, {
3.824797958946247*^9, 3.824797963474185*^9}, {3.8247980956117268`*^9,
3.8247980985627832`*^9}, {3.824798171626272*^9, 3.824798174250143*^9}, {
3.824871013588922*^9, 3.8248710184668417`*^9}, 3.824873861501131*^9, {
3.8248739741262417`*^9, 3.8248739754441137`*^9}, 3.8248740451185627`*^9, {
3.824876095188408*^9, 3.824876096507303*^9}, 3.824876172034321*^9, {
3.82487641698841*^9, 3.822675627713685*^9}, {3.8226780729608994`*^9,
3.82267808764895*^9}, 3.8226781591609087`*^9, {3.8226782855130587`*^9,
3.822678286848304*^9}, {3.822682935825254*^9, 3.822683056730687*^9},
3.8226832112833223`*^9, {3.825107830097363*^9, 3.825107840400428*^9}, {
3.825376350944564*^9, 3.8253763799736137`*^9}, {3.8253803715752287`*^9,
3.825380375357637*^9}, 3.8253804482205563`*^9, 3.8253847353913813`*^9,
3.825386157364841*^9, 3.8253893753457327`*^9, {3.825390557993737*^9,
3.82539056751958*^9}, {3.825390764526593*^9, 3.825390765284821*^9},
3.825390796614028*^9, {3.825390863269362*^9, 3.825390863851976*^9}, {
3.825390946836618*^9, 3.825390951803205*^9}, {3.825391072864626*^9,
3.825391073032188*^9}, 3.8253926623116693`*^9, {3.825392722693924*^9,
3.825392726060487*^9}, {3.825393313614963*^9, 3.825393313901581*^9}, {
3.8254002350999002`*^9, 3.825400235434561*^9}, 3.825731842816152*^9},
CellLabel->
"(Debug) In[174]:=",ExpressionUUID->"12056f8f-f46b-4ca3-964a-4b9656492af3"],
Cell[BoxData[
RowBox[{"(*",
RowBox[{
RowBox[{"data", "=",
RowBox[{"Flatten", "[",
RowBox[{"{",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"RandomVariate", "[",
RowBox[{"ExponentialDistribution", "[",
RowBox[{"\[Lambda]10",
RowBox[{"(",
RowBox[{"N0", "-", "i", "+", "1"}], ")"}]}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "\[Tau]0"}], "}"}]}], "]"}], ",",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"RandomVariate", "[",
RowBox[{"ExponentialDistribution", "[",
RowBox[{"\[Lambda]20",
RowBox[{"(",
RowBox[{"N0", "-", "j", "+", "1"}], ")"}]}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"\[Tau]0", "+", "1"}], ",", "n0"}], "}"}]}], "]"}]}],
"}"}], "]"}]}], ";"}], "*)"}]], "Input",
CellChangeTimes->{{3.822678115648424*^9, 3.822678152344617*^9}, {
3.8226783151528387`*^9, 3.822678338152616*^9}, 3.825386163605687*^9, {
3.825719779112687*^9, 3.825719779999361*^9}},
CellLabel->
"(Debug) In[181]:=",ExpressionUUID->"9aee1e6a-ca1d-45c6-82fe-5c134f19ba06"],
Cell[BoxData[{
RowBox[{
RowBox[{"data", "=",
RowBox[{
RowBox[{
RowBox[{
"Import", "[",
"\"\<~Downloads/handbooksoftwarereliabilityengineering/DATA/CH4/SYS1.\
DAT\>\"", "]"}], "[",
RowBox[{"[",
RowBox[{"12", ";;",
RowBox[{"-", "2"}]}], "]"}], "]"}], "[",
RowBox[{"[",
RowBox[{";;", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"N0", "=",
RowBox[{"2", " ",
RowBox[{"Length", "[", "data", "]"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"data", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"RandomVariate", "[",
RowBox[{"ExponentialDistribution", "[",
RowBox[{
RowBox[{"1", "/", "3"}],
RowBox[{"(",
RowBox[{"N0", "-", "i", "+", "1"}], ")"}]}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "data", "]"}]}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]"}], "Input",
CellChangeTimes->{{3.8257153966885843`*^9, 3.825715531713738*^9}, {
3.825718452618618*^9, 3.8257184850359707`*^9}, {3.825718706959814*^9,
3.825718722933845*^9}, {3.8257187579092693`*^9, 3.825718780347151*^9}, {
3.825718814489687*^9, 3.825718881831851*^9}, {3.825720216763358*^9,
3.825720217474118*^9}, 3.825729819493807*^9, {3.8257303895413017`*^9,
3.8257303950680723`*^9}, 3.8258031323673697`*^9, {3.830115871759104*^9,
3.830115878273086*^9}},ExpressionUUID->"88eef10d-427c-4169-98b0-\
ea1989c0587f"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"U", "[", "\[Lambda]1_", "]"}], "=",
RowBox[{"Simplify", "[",
RowBox[{"Total", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{
RowBox[{"PowerExpand", "[",
RowBox[{"-",
RowBox[{"Log", "[",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"PDF", "[",
RowBox[{
RowBox[{"ExponentialDistribution", "[", "\[Lambda]", "]"}],
",", "t"}], "]"}], ",",
RowBox[{"t", "\[GreaterEqual]", "0"}]}], "]"}], "]"}]}], "]"}], "/.",
RowBox[{"{",
RowBox[{
RowBox[{"t", "\[Rule]",
RowBox[{"data", "[",
RowBox[{"[", "i", "]"}], "]"}]}], ",",
RowBox[{"\[Lambda]", "->",
RowBox[{"\[Lambda]1",
RowBox[{"(",
RowBox[{"N0", "-", "i", "+", "1"}], ")"}]}]}]}], "}"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "data", "]"}]}], "}"}]}], "]"}], "]"}],
"]"}]}],
RowBox[{"(*",
RowBox[{"-",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"N0", "-",
RowBox[{"Length", "[", "data", "]"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "t"}], " ", "\[Lambda]1"}], ")"}]}], "/.",
RowBox[{"t", "\[Rule]",
RowBox[{"Total", "[", "data", "]"}]}]}], ")"}]}], "*)"}],
";"}]], "Input",
CellChangeTimes->{{3.824798078420321*^9, 3.82479820892268*^9}, {
3.824798351268783*^9, 3.824798444435299*^9}, {3.824798683138339*^9,
3.824798812667768*^9}, {3.824798947051549*^9, 3.8247989617558327`*^9}, {
3.824799102941595*^9, 3.824799115564328*^9}, {3.824799206580085*^9,
3.824799240611725*^9}, 3.824799312381834*^9, 3.824799394557672*^9, {
3.82479945879666*^9, 3.82479946842881*^9}, {3.824799542413204*^9,
3.824799543211771*^9}, 3.8248636868846817`*^9, {3.8248638148680964`*^9,
3.8248638786282597`*^9}, {3.824864711514967*^9, 3.824864750247938*^9}, {
3.824864963489201*^9, 3.824864969247636*^9}, {3.824865052465322*^9,
3.8248650529718037`*^9}, {3.824865100024632*^9, 3.824865137830184*^9}, {
3.824865316002928*^9, 3.824865344850099*^9}, 3.824865387483205*^9, {
3.824865417600205*^9, 3.824865438978271*^9}, {3.8248656387956753`*^9,
3.824865713187004*^9}, {3.824867049592964*^9, 3.824867051138234*^9}, {
3.824868417510271*^9, 3.824868421406719*^9}, {3.8248710089994383`*^9,
3.824871009693295*^9}, 3.824872195274089*^9, 3.824873755578967*^9, {
3.8248739082234287`*^9, 3.824873911142541*^9}, {3.824874127601722*^9,
3.822678388489332*^9}, {3.822678428280992*^9, 3.822678453770197*^9}, {
3.8226784849770603`*^9, 3.822678488089197*^9}, {3.822678600369632*^9,
3.822678619842081*^9}, {3.822678672865734*^9, 3.8226786798892107`*^9}, {
3.822678790074387*^9, 3.8226788177538643`*^9}, 3.822678873721414*^9, {
3.822683099789221*^9, 3.822683116780435*^9}, {3.8251080045428953`*^9,
3.825108022139064*^9}, {3.825108242074274*^9, 3.825108262841107*^9}, {
3.825108370158791*^9, 3.8251083816934547`*^9}, {3.825384542528989*^9,
3.825384648918984*^9}, {3.825384718908483*^9, 3.825384720323079*^9}, {
3.825385623898158*^9, 3.8253856936745567`*^9}, {3.8253863587667103`*^9,
3.82538636514184*^9}, {3.825387549587175*^9, 3.825387550682221*^9}, {
3.825387584747428*^9, 3.825387617124914*^9}, {3.825387669134452*^9,
3.825387686545822*^9}, {3.8253877299135942`*^9, 3.8253877602565737`*^9}, {
3.825387841171155*^9, 3.825387851825273*^9}, 3.8253881279966717`*^9, {
3.825390121694598*^9, 3.825390123430393*^9}, {3.82539093107045*^9,
3.825390931764135*^9}, {3.825715536804729*^9, 3.825715581738703*^9},
3.825715660262227*^9, {3.825717536374125*^9, 3.825717536790167*^9}, {
3.82571759347431*^9, 3.8257176716322536`*^9}, {3.825717706017742*^9,
3.825717794461769*^9}, {3.82571783282381*^9, 3.825717847824052*^9},
3.825717915845983*^9, {3.8257179948592243`*^9, 3.825718056378707*^9}, {
3.825718089488236*^9, 3.825718154587749*^9}, {3.82571819381328*^9,
3.8257181987858133`*^9}, {3.825718235683639*^9, 3.825718239232774*^9}, {
3.8257182795160418`*^9, 3.825718380926961*^9}, {3.8257185886214447`*^9,
3.825718664306615*^9}, 3.8257205156765137`*^9, {3.82572064535575*^9,
3.825720659913652*^9}, 3.8257271528028603`*^9, {3.8257275724050007`*^9,
3.825727578965185*^9}, {3.8257276158318653`*^9, 3.825727663693797*^9}, {
3.825727829352857*^9, 3.8257278396933613`*^9}, {3.825728291163577*^9,
3.825728296519849*^9}, {3.8257296609543333`*^9, 3.825729663640209*^9},
3.825729694568244*^9},
CellLabel->
"(Debug) In[61]:=",ExpressionUUID->"b2c727d7-3fe5-437f-81d5-c3b42c3c16d6"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"dU", "[", "\[Lambda]1_", "]"}], "=",
RowBox[{"GradientG", "[",
RowBox[{
RowBox[{"U", "[", "\[Lambda]1", "]"}], ",",
RowBox[{"{", "\[Lambda]1", "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.8248651442325773`*^9, 3.824865149303952*^9}, {
3.82486705526359*^9, 3.824867056696405*^9}, {3.824867109671906*^9,
3.824867112086319*^9}, 3.824868660580357*^9, 3.8248687042461557`*^9,
3.824868832116352*^9, {3.824871023348235*^9, 3.824871037010594*^9}, {
3.825715620570724*^9, 3.825715642561338*^9}, {3.825715680567915*^9,
3.825715683808449*^9}, {3.825720756290715*^9, 3.825720762225649*^9}, {
3.825727912367008*^9, 3.8257279175888042`*^9}, 3.8257282855898*^9},
CellLabel->
"(Debug) In[62]:=",ExpressionUUID->"6f2d8446-9c47-4ab3-896f-654eb317b18c"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"ddU", "[", "\[Lambda]1_", "]"}], "=",
RowBox[{"HessianH", "[",
RowBox[{
RowBox[{"U", "[", "\[Lambda]1", "]"}], ",",
RowBox[{"{", "\[Lambda]1", "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.824866400470134*^9, 3.824866432251779*^9}, {
3.824867060562065*^9, 3.824867064181905*^9}, {3.8248671162401133`*^9,
3.824867118966241*^9}, 3.824868833228517*^9, {3.824871025588564*^9,
3.824871035651095*^9}, {3.825392193559525*^9, 3.825392194805924*^9}, {
3.825715626122882*^9, 3.8257156501931677`*^9}, {3.825715687913241*^9,
3.8257156898416033`*^9}, {3.825715734010895*^9, 3.8257157387212353`*^9}, {
3.825720766650605*^9, 3.825720807370337*^9}, {3.825727921990048*^9,
3.825727927211775*^9}, 3.825728287298674*^9},
CellLabel->
"(Debug) In[63]:=",ExpressionUUID->"39feadf1-d3cc-4b3a-a1b8-83e1af3f021f"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"INTERVAL", "=", "1001"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"STEPS", "=", "5"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CHAINS", "=", "5"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dt10", "=", "0.0000000001"}], ";"}], "\n",
RowBox[{
RowBox[{"dt20", "=", "0.0000000001"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"outbnd", "[", "q_", "]"}], ":=",
RowBox[{"AnyTrue", "[",
RowBox[{"q", ",",
RowBox[{
RowBox[{"(",
RowBox[{"#", "\[LessEqual]", "0"}], ")"}], "&"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Dim", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"qinit", "=",
RowBox[{"RandomVariate", "[",
RowBox[{
RowBox[{"UniformDistribution", "[", "]"}], ",",
RowBox[{"{",
RowBox[{"CHAINS", ",", "Dim"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS", "=",
RowBox[{"hmc", "[",
RowBox[{
"U", ",", "dU", ",", "ddU", ",", "Dim", ",", "5000", ",", "10000", ",",
"True", ",", "True", ",", "qinit"}], "]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.825716104103396*^9, 3.825716133516217*^9}, {
3.825716177102722*^9, 3.8257162199397593`*^9}, {3.825718677579792*^9,
3.82571868435189*^9}, {3.825718872995072*^9, 3.825718892734666*^9}, {
3.8257207788653393`*^9, 3.825720780273355*^9}, {3.82572083330447*^9,
3.825720869464089*^9}, {3.825720992833641*^9, 3.825721062045289*^9}, {
3.825721102861444*^9, 3.8257211051487513`*^9}, {3.825727933223653*^9,
3.8257279450505667`*^9}, {3.8257283148183107`*^9,
3.8257283151717787`*^9}, {3.825729940414426*^9, 3.8257299422184153`*^9},
3.828826195177326*^9, {3.828826246831999*^9, 3.828826248308856*^9}},
CellLabel->
"(Debug) In[64]:=",ExpressionUUID->"28d96fc8-8011-4139-bd6d-447897055a94"],
Cell[CellGroupData[{
Cell[BoxData[
TemplateBox[{
" ","\" \"","1001",RowBox[{"-", "922054.2685561789`"}],
"2.49023616730644`","0.21672681194251667`","0.2196453793696374`",
RowBox[{"-", "2287.835697645888`"}],RowBox[{"-", "924342.1042538248`"}],
RowBox[{"-", "2.271989817795716`*^8"}],"0.00003200568495929872`",
"0.00006237002557739829`","True",RowBox[{"{",
RowBox[{"1", ",", "5"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.828826258600786*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"d8ea747b-1c3e-444c-a3b2-82a6b5836c31"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","2002",RowBox[{"-", "1.4322344965278095`*^7"}],
"0.003251441248375656`","0.24132678334742869`","0.4309884779075283`",
RowBox[{"-", "2286.1248943405117`"}],RowBox[{"-", "66225.16654740366`"}],
RowBox[{"-", "1.4324631090172434`*^7"}],"0.0001336956889071248`",
"0.00028658858268655423`","False",RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.828826266570957*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"4cee0c23-d39b-48d2-9b87-61803530efba"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","3003",RowBox[{"-", "4031.692964582932`"}],
"2.848596219379267`","0.4535218801489728`","0.4995596809860869`",
RowBox[{"-", "2285.343792372783`"}],RowBox[{"-", "6317.036756955715`"}],
RowBox[{"-", "823079.6879303065`"}],"0.00038144940355580375`",
"0.0006757608868127241`","True",RowBox[{"{",
RowBox[{"2", ",", "5", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.828826274557069*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"0001f45a-c305-4f28-aab3-f7f1843be5f9"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","4004",RowBox[{"-", "35341.655224427195`"}],
"0.002503756878417551`","0.46795696517732666`","0.4443286822185363`",
RowBox[{"-", "2285.511982145788`"}],RowBox[{"-", "3347.703388981167`"}],
RowBox[{"-", "37627.16720657298`"}],"0.0008176706730433963`",
"0.005500879411464487`","False",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.8288262824336348`*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"439a71d7-37bd-4512-8fb1-1c547dd15416"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","5005",RowBox[{"-", "207.78670557760552`"}],
"2.7705153528666155`","0.22995756962877212`","0.43214592118565964`",
RowBox[{"-", "2288.382048949268`"}],RowBox[{"-", "2496.1687545268737`"}],
RowBox[{"-", "7539.616015759533`"}],"0.0019280246755010938`",
"0.009745143431053444`","True",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.828826290479362*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"a004c807-6fcc-4e8c-a249-acd6589a6ddb"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","6006",RowBox[{"-", "5253.091932882098`"}],
"0.0016346153190342004`","0.42601920456267417`","0.5248931883436466`",
RowBox[{"-", "2286.5240828774354`"}],RowBox[{"-", "2496.1687545268737`"}],
RowBox[{"-", "7539.616015759533`"}],"0.0019280246755010938`",
"0.009745143431053444`","False",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.8288262937424507`*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"11dfaf1e-6130-4bba-8dde-6fe2f2b09995"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","7007",RowBox[{"-", "211.10594896505336`"}],
"2.3973911842181614`","0.4741573745751183`","0.4358496437698161`",
RowBox[{"-", "2285.0628055618204`"}],RowBox[{"-", "2496.1687545268737`"}],
RowBox[{"-", "7539.616015759533`"}],"0.0019280246755010938`",
"0.009745143431053444`","True",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.8288262970466948`*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"d1684562-20fd-42aa-9587-b306dfff616b"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","8008",RowBox[{"-", "5252.605037708419`"}],
"0.0032299741578167285`","0.6302603341271352`","0.5090964571206207`",
RowBox[{"-", "2287.0109780511134`"}],RowBox[{"-", "2496.1687545268737`"}],
RowBox[{"-", "7539.616015759533`"}],"0.0019280246755010938`",
"0.009745143431053444`","False",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.8288263004067793`*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"d2e11b76-26b9-46e8-b38f-428141161527"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","9009",RowBox[{"-", "210.0727460096291`"}],
"1.489932811683944`","0.38553874630856666`","0.35214746099851585`",
RowBox[{"-", "2286.0960085172446`"}],RowBox[{"-", "2496.1687545268737`"}],
RowBox[{"-", "7539.616015759533`"}],"0.0019280246755010938`",
"0.009745143431053444`","True",RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "6"}], "}"}]},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.825727783548335*^9, 3.825728963805464*^9,
3.825729668983017*^9, 3.8257297387246447`*^9, 3.825729826318863*^9,
3.825729947252529*^9, 3.825730400301259*^9, 3.8258031465255527`*^9,
3.8288261571300364`*^9, 3.828826201304373*^9, 3.828826303656089*^9},
CellLabel->
"(Debug) During evaluation of \
In[64]:=",ExpressionUUID->"3d3d48e1-9450-4488-bd36-3db65af41138"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"QS", "//", "Dimensions"}]], "Input",
CellChangeTimes->{{3.825716232112521*^9, 3.825716237875525*^9}},
CellLabel->
"(Debug) In[73]:=",ExpressionUUID->"5ebd8265-412b-4c06-a76b-41f8cc62703b"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"25000", ",", "1"}], "}"}]], "Output",
CellChangeTimes->{
3.825716238369419*^9, 3.8257178948033123`*^9, 3.825718426222116*^9,
3.825718514655002*^9, {3.8257186893758717`*^9, 3.8257187289268303`*^9}, {
3.8257188367799253`*^9, 3.825718862437345*^9}, 3.8257189147746677`*^9,
3.825721138417251*^9, 3.825728978196062*^9, 3.825729683386032*^9,
3.825729753340722*^9, 3.825729840214404*^9, 3.825729960518037*^9,
3.825730415006413*^9, 3.825803161276286*^9, 3.82882617055722*^9,
3.82882621619363*^9, 3.8288263069936047`*^9},
CellLabel->
"(Debug) Out[73]=",ExpressionUUID->"c7470434-28b8-4fb2-a594-a7f807ce35fa"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"ListPlot", "[",
RowBox[{"QS", "[",
RowBox[{"[",
RowBox[{";;", ",", "1"}], "]"}], "]"}], "]"}]], "Input",
CellChangeTimes->{{3.8257162396994534`*^9, 3.8257162579599047`*^9}, {
3.8288263251739273`*^9, 3.828826326851267*^9}},
CellLabel->
"(Debug) In[79]:=",ExpressionUUID->"640eac33-c814-46ae-938d-635f6fe81224"],
Cell[BoxData[
GraphicsBox[{{}, {
{RGBColor[0.368417, 0.506779, 0.709798], PointSize[
NCache[
Rational[1, 360], 0.002777777777777778]], AbsoluteThickness[1.6],
PointBox[CompressedData["
1:eJxc3XlYDf3/P/Bk66ab7CEkIcQdQggvxB1CCCE59pCEJEkqSXuntG9OdapT
tiTpJhyEEEIIIXv2EEL49b1+79fnuubpH9c0M48zM+99mfd0X7Ju5nJNDQ0N
1WINjf/7////qx4zInv0NfnKu2PEH2iSvEsfffMHYluLXnueH9jE657Y1qHY
MUFmvxL4+LY0xUw9Su1+X2zrktTTo7nr/EYGD6wU2/r0IN/eunzuHbFtAJ4h
rf7rUmhoU/Z6gWdElzsOvHjH+LbY7ktDqq7LgiPKxbYxeAMoP81r+1lrvn4T
auG1fb3uJT5+IDk62/1r8Zm3B1HHi+rUgrF8/mDwTMEbQuaywbrmU3j/UPCG
gWcG3nDwRoA3krZqfv70PP+m2Dan9VYruz59wM9vFHijwRsDHkk9LwJvrNTz
GgveOKnnNY7aPXg+KLqUw2O81PMaT7vqUpI6ht8S+y2knpcFXN8EemzZ7+OI
auF5TaBRhkVDRkbz85wo9bwmwvX9C9f3r9TTsJR6XpYU8DLcbIInP99JcH2T
4Pomw/VNhuubQinuMw71fMH7p8D1WVHEHKXaIYd9K7i+qXB9U6WexjQqeaM4
vngy+9OknsZ0OtU0s+d/hnz/0+H6rKUeWYNnLfXU1hQUYtZqSh+OHzOkHs2g
pl02Xdk6rkKcP0PqqWfQGB2991Uz+fdnUk2b8MsW3uJ4mgn3O1PqqWdKPY1Z
Uo9m0eYP7S2nzeX7nUUnvgUFJ/wR+9WzpJ6GDZ3bc0zL4oW4P7KhfgNvLsyY
w+FhI/XUNpT0JOq48ijnR7PJ797YvT49RP5Is6We12ypp54t9TTmUF4HR99p
lXz+HLqVlHt15wy+/zlSTz1H6mnMlXo0l6ZMeu78ZTk//7mU3nneKr0r4nrV
c2nQlHlq1/ecn9uS4eukbTe0+P5t4fnZSj21rdTTmEfHdvy+2rGUr38edVv3
dEsTXf79efD85tGPpYuSsmfw+fOlHs2Xel7zpZ56Pjy/BVSUFRR57iOH/wKp
57WApg5Z2Gfqc+GpF1Dm9a9H85Zz/LGD52cn9bzs4PnZST2NhdRufdUat46c
fhbC81so9dQLqVdRr483/Pn67aUe2Us9L3t4fvZST2OR1KNFVDZy+d5xGzj+
LaLXO94bTqzi619Euz98jUku5d+XST19GXgymvy6VcSvPHE9Mhk9nFmUklPC
+a1M6inAU4NXKaMWBgEJG7pz+ltMNya19CkJKxO/v5i0LDxK437z7y+WerLF
Us9rsdRTLKYtyeOmjooT16deTKZxIYdWtBD5VSV4GkvIukXQ6ok7xfH6S6Qe
LSH34iP3OzqJ82VLaNSlBd9eafLzXSL1FEsoOszk5I/DHF7gVYKnsVTq6S+l
lv5zNza5LfbTUipp0Mp6ywKxX7aUckvcIrWncvxcSjlDlqfl/+/3wVODV7mU
lmt+mDstmPOPZaR76MvY3GwO/2VSj5ZJPdkymvhZw7XVE/79ZVJPAZ4avErw
NJZLPf3lNHLwI+1zPfn3l0s92XKp57Vc6inAU4NXCZ7GClpyTHtqw8b8+yuk
Hq2gmq7q1jff8O+voDE5q06Ma8rxb4XUU4CnBq8SPI2VUk9/JX0s6pJzcgHn
hyulnmwlxL+VUk+xkp5kz1of05nzy5X0d9n24AVrOP2tpBNTLlWtaMj5nwPE
Pwfacc9yUsBrfv4OUk/mIPW8HGheg2HD+nYS8UsBntqBfHe+b3X0tkhflQ6U
/Ze209AFnP5XST39VVKPVlHxjL0/D3tz/FtF1tezNshm8f2vknoK8NTgVYKn
sZpaLrdc3stInK+/mlxzS3PiR/Pvr5Z6stU03E7zJEVw+bZa6inAU4NXuZr2
v3/VsvYvvv819LPig1XPc2K//hqpR2uo80Nz8/51D8Xvr6Hg8se//f/i8m0N
nV57utuUDH7+4KnXUILD/PDOp0V4V4Kn4Uix7hn90zbz83eUeuQo9WSOUs/L
UeopHOniXrOIW/e4vgdepSM5XTJop/2En/9aqae/VurRWqknWyv1vNbSL1nX
ihtcH1KsJb+c2Aby9Rz/1lLzUW+7HdHk+L+W0gcVJXUewPmvk9TTdyKZqvvL
RTs4/jtJPZkT6Qxv4XroAz9/J1qYe2HJkGkc/5zoWl7c+UZlHP7gVYKnsY7s
ph/TmHhIePrrpB6tk3qyddTVfOHR9hs5/1tH10tu7nz6D6d/8NTr6K9XOyw3
cPuuEjwNZ/Ju7XHeJF8cr+Ms9fSdpZ6JM93+dOB8TrNH4vqcacjzyevtroj7
sQZP5kwDbTVX53J73Rk8L2dK2bYxfN9vsS0HT+FMD1emLY/bKa43Bzy1Mz0Z
su10ZANxfCl4lc5U0rCJ7bZjYn81eBrrpZ7Oeqmnv17qmayXerSeOp5d3nS0
UtRfrMGTradZ3wv+cY7l+wfPCzw5eArwctZTeAPzdO2+fP/glYJXuZ7GeZ26
1vOaiN/V4GlskHo6G6Se/gapZ7JB6tEGym5qUPr4pjjfegN9uG64+aArx9cN
9M509UqzsSL9OoPntYEGHTvkFbFMHC8HTwFeDnhq8ErBqwSvGjyNjRQwY+XB
13fFts5Gunpd95w1l3/6G6lRo+A3Oev4/jeSavrKKT2Pc31/I5mEkt31IRz/
wZNtpKz4zYk/DMXzdAbPCzz5RrrtH1TttILz8410sneDgOV3xP4c8NTglW6k
ug51WaZWnN+AV72Rvg7/7S3n8kPDRerpuEg9fRe683X75ZE+4n5MXKhooEGr
x4kifpILDV+UGePsx+EPngw8Zxd60DM5Y44llx8u5LNiRcce38T1yl1olEfX
HIdj4nwFeDku1Ml5srn6Iefv4JWCV+lCzmabqry/8/2Dp7GJfr3J8Ljyv/vf
JPX0N0k9k01SjzZJPWvwZJsoKOHKga0f+P430bGgYX/Wcn3EaxPNSLx8cccz
UT+Sb6JWSbOX6HJ/iwK8HPDUm6g636pH6jMRHqWb6PIf45UfXvH9g1e9iXKv
Ke2/Peb7d5V6Oq6Uu/TK7xN9uTx0lXomrlKPXOmNge3eiDjO/1xJu+TfVM9P
XF8Az9mVEju4FX7uyfEfPLkrZYeO+3bnGqd/8HLAU7tS42Tj41P/E8eXglcJ
XjV4GpspIT9y3/n5Irx1Nks9/c1Sz2QzTR5aOLR5IdfPN0s9681Eg8Z5957O
9QXwnDeT2xmtsgEvuP0Mnhw8xWbyT3fve/e7eL454Knrzz80O2lwf87/N9NP
G2+r2CLO/8Cr3kxDG+/b9R/X5zTcpJ6Om9TTd6PMHrcS2hzm+3eTeuRG0Z3u
WlxfxfkfeDI36nn4i0FFAed/4Hm50bg6i2v7Vonj5eApwMsBTw1eqRu5dBr7
X91sru+60RLVnfuffcXzrgZPY4vU09lCKxv7PApew/nfFqlnsoUeb7z7YRiX
57SFRp616LdkEKd/8GTgOW+hoi63lTEWnP63UOcRjbqP1uH8DzwFeDngqbdQ
7q9Nh3s1F/dXuoW+tf26pvV2rv+AVw2ehrvU03Gn040X2+j9xeWfOxV9mGTh
Himu18Rd6pE7TXs4ZtI/3D9nDZ4MPGfwvMCTg6cAL8ed/lo8dtvSIi7/wSsF
rxK8anfqFOiR/0eTxx+2Sj2drVJPfytdnxwW2e4g3/9WqUdb6X7ry/2e9Obw
30pLSgfdSpvH8X8r7dd94d1jOt8/eF7gycFTbCXt3JUL1vmJ83PAU2+lw7Pn
j9OyFtul4FVupV0VUTpDz4v91eBpeEg9HQ+pp+8h9Uw8pB550Gp3i+aRvcS2
NXgyD1KmRvr05/qYM3he4Mk96J+4vjXFC/n+wcsBT+1BK0/227V4E8d/8CrB
qwZPYxvdqTM6OvenOF9nm9TT30ael7okLLMR/bMm26QebZN61tvo6JBhVDZJ
nC/bRi0mz+4/u5bzv2300PhU57HteLwNPPk2yl5l0KSS+3sV20i1vi7y6QIu
/8BTb6N8T9ukQg2+/20UcD1t1+8K0T6rBK96G2V0XFDTrQW3Pz2lno4nrboR
euafdM7/PaWeiafUI0+pZ+1J34e91lsdx/mfJy03DjA9VMH3D56XJxlfuGzs
YMT5P3gKT7LOa3ki5wXnf+CpwSsFr9KTGhzsOO9LHt8/eBrbqavJgcT1V7j8
2y719LdLPZPtUo+2U17d1a2zBojnb72dzHxL/2uSx/cPnvN2uuhy5NvBj1z+
gScHTwFeDnjq7eQ24P2nPx+4/rud/LY57dXewfe/nTru7B3825/LP/A0vOh+
2IRPJly/0/KSejpeNCuvwffRx4Sn6yX19L2knpEX3Wp4aPSMOs4vwTMDj7zq
85uQus7DxPOw9KI9S5cOes/9G9bg2YInA88BPGfw3Lxo6PvGSTvcuf+k/vr3
Wu38Nk/4/uDJwYsFTwGeCrwc8ArAU4NXDF6pF40/4trmYqHI78rBq/Si1sd6
Pyi4L55/FXjV4NV6UefZvg7Nz3L+501rCnJHBAaJ69XyJtdE+ZSkWxz+3lJP
15s+9Y3y/qHL5aM3Tdt7YHOFCYe/N10JbVsrN+D2Anhm4BF4luBZg2frTdpP
OocvD+P6hTctTOzu+HaK2O/gTV/X+E/JHsjhD54beF7g+YMn96abiReCr9uI
648FTwGeCrwcb3qwIL/Zoovi/gq8qf/jNI8TL0V9Xu1Nd50bJsdeF8cXg1cK
Xjl4ld5k1bjtwVNWIj1WgVftTV558+0dMzj8wdPwoaIrB8K+vRX3p+VDL1a2
eP+MxPXp+Eg9XR+pp+8j9YzAMwHPzIeaXHHaovjN453gWYJn7UPtc86uWDWS
0z94Mh/af38+nYgQ4eEAnrMP7fjUoVnzm5z+wfPyoeb9g2YFcPvKHzy5D0WW
RSx9GyTSWyx4Ch+6ZD5lYIduIn2qwMsBr6D++bTq0HsTj1+owSsGr9SHeryz
23xzqNhf7kMeVklNCmK5futDjy5p2p47w+EPXjV4tT40+cLtysUHOP3vIPPP
utt7+3H476DffUt/lF7n8N8h9XR30IJhLd3TXTj8d0g9ox10sXnp1jRLLi93
0KyDBx0cxon7NQOPdtCKyd/0xzTn8AfPegflbRxn9nak8GzBk4HnAJ4zeG7g
ee2goLqPv9JrOP8HTw5eLHiKHTT7WaXN9++c/sHL2UGDK953XeTA6R88NXjF
O6hokqbO8wbcngGvHLxK8KrAq95BEw1nrO7oI+JHLXgavlJPy5dmeSW5FHXn
9qCv1NP1lXr6vnRG1kwZe1mcbwSeCXhmvpTYu/xEUC2X/74U9vvkMD1uv1uC
Z+1L8zY/dFztIeK7rS8piu/GWGpw/yp4DuA5g+cGnhd4/uDJwYv1Jdcx+3/u
OMz5f/35tusednsu9qvAywGvwJc8FT8WtDjL/TfgFfvS0vs1wy8P5fo8eOXg
VYJXBV41eLXgaeykDVvtFH14voXWTurxp5fK3pDDfyeVPO1x//ISkT/o7pR6
+julnhF4JuCZ7aSojzMeaq7l/iPwLHdSsyXZn1LWcfkPni14MvAcwHPeSa5x
t7dozhLbbuB5gee/k/r1Sn1clc398TvJY/DKvr4/Of/fSV8KJu46Hib2K8BT
gZcDXgF4avCKd9LoPps1Fu7i8AevfCfJnddl580U5WcleFU7SVkZMGTYWO7v
A692J3XavjXnUBn3f/lJPS0/qafjJ/V0/Ui2TuP7FUsOfz+pZ+RHZ76aFPzH
/Vkm4JmBR370wa22+FAXzv/96HRMqueUIk7/ftShaHyBxVhO/+DJ/OiTnUa/
BtM4/YPnDJ6bH50q7eK08B8Of/D8/Wh8T49vJr6c/4MX60fTAxp/ccnn/lnw
VODl+NFolzZXlkRw/u9HM8d1kB/azeU/eMXglYJXXn+/5aYWaQHcvwlelR9N
nNbq+nQd7u+qP/9V32adgrj+B57GLppp87HzaRNO/7ukns4uuhS+Kkd7O6f/
XVJPfxdNPG5vOZHzH6NddLB44Opd3L42Ac8MPNpFJnv73c8x5vwfPOtdNLBj
Qm1EY7HfFjwZeA67aMGn2vQ2PN7uvIsuu+Qo3jUWnht4XuD5gycHL3YX5X20
uDQ+ktP/Lvq5f+J/3ZaL+1XtojPp1f3mX+P0D14BeGrwindRkK/Pt8RALv/B
KwevchdZ2b1qoMn9NVXgVe8ig8TBw7ccEefX7qLmqvguUxtx+PuTsuzETdtE
EV+0/KmH24DTE324P8hf6un6Sz19f2o33efulzEi/zDyp1k+VzJtrDj/B88M
PALPEjxrf3pXO61P5gJu/4EnA88BPGfw3Pzp18jAt+obXP8Hz9+fmhzuUuXm
wONR4MX6U0jjhKtttTn9g6eqP791bYEN9z/m+FNy+WRFOce3An9yGmSbFfKO
+y/9KevB8NadW3L570/6AxLGFBzm9h945eBV+tPyxZUbCr4Lrwq8avBq/Wny
ouAHqwo5/APolH7HaBOej6gVIPV0AqSeboDU0w+QekbgmYBnBh4F0IC5+2Jm
jeP8Hzxr8GwDyKJhptrkKoc/eA7gOYPnBp5XQH352/fgYCdRnvkHUNvktJzw
7Vz/Ay8WPAV4qgDyavTn/t8nuf0XQOYvqNUEng9WAJ4avGLwSgNo15txJx3n
ivlQ5QFUmtqx99bnXP8LIDMH2cszCeL4KvCqA6ixwenGTbZy+AfQ9giVa9IX
nr8cKPW0AqnxoXlXli4Xnk6g1NMNlHr6gTQj+/X5HTw+YwSeSSBpDOhoEGDH
+X8gHdyhOXoMt6cJPEvwrMGzBU8GnkMgjdrXbkJTbm84g+cGnlcgTRt619XG
kdt/gTSxaF2l5t9c/gfSTesu87/MFvtjA+nP56zTw3px+IOnAi8HvALw1OAV
g1caSH+foqDUkZz/B1L4YLV5cjMOf/CqwKuuv54WJ58Y6Iv9tYFkc/O5U1wd
9/8EST2tIDIKPBStbMvjQUFSTzeIIicOm5N6kedHBEk9I/BMwDMDj8CzDCLf
abSl4Riu/4NnC54siLx80n53ecrhD55zEP1TkNl/bH9u/4PnBZ5/EA1WOwbu
XcrpH7xY8BRBNGzd6d3be4j9KvBywCsIomRXGvMhl8cvwSsOoguOtx8d8eP6
fxDVHrmTnNxepNdy8CqDaO0rq+r+hzj/B68avNogupuYbbKuNYd/MO31yl4d
rcn5fzB9G/pgzfsaLv+DpZ5uMC3VPO54uhvX/4OlnhF4JsEUmFdjlRTN5T94
BJ4leNbg2YInC6bxgzwNmnfk/h/wnMFzC6anipdR7s5c/wPPP5g2NJ/bq+YC
h38wTV9ygV42FPtjgykgQXWhaxaX/8HUYHbYklbB3P8DXk4wWThPC+l7mvv/
gunUqVP3l2zi8h+8YvBKwSsHrxK8qmAqe//04pauPN4DXm0wvfQJuzhQl8M/
ROpphUg9nRCppxsi9fRDpJ4ReCbgmYWQ453vjwx4vI1C6Ed0odfvNZz+wbMO
oV0Lnz1uG8v1f/Bk4DmA5wyeG3he4PmDJwcvNoROFJxxPr+Rwz+EfN0DE77w
+IYKvJwQqrYePWfoBk7/4KlDqOXgWU/la0V9oTiEjr9fc1N5hMMfvHLwKkNo
wpszPUtvcP0/hGZnXD4SvYzb/yH04FKDA7+/cviH0NlHw/e+uSZ+TyOULBVN
p8cFcPoPpUYVQ+K63+Lx8FCppxsq9fRDpZ4ReCah9ME8KvWSs/DMQmlIxkz9
QHtO/+BZhtJDw8Hlz3m+nTV4tuDJQumXUjtgGs+ncAilwbcLRt4x5vofeG7g
eYXSZc/1CyP3ivDwB08OXix4ivrz4490L70n8mcVeDmhlPvg7MBFh7n9D546
lLq1vPDfYJ4/UAxeaSg5zvQ9MLwX9/+G0v7PRnstBnH4h9LK0ms3U3k+RBV4
1eDVhpLNr3mhb77z+x9hUk8rjDqPu577bysxvqwTRoftbHKv8nxI3TCppx9G
BvrPq67tEPuNwqhp+lKbxdzeNgHPLIz2auoOPtWQ5w+G0ZrxuRtfTBaeZRj5
jh0eu8aR6//g2YInCyNnjcPna3n83AE85zC6HeXUr8MqHv8Lo+j4Cbr7tDn9
g+cfRr1eR72IbMftP/BiwVOE0Y++Pm5XunH+H/Z/82uca9nPCaNB3XYnzezN
7T/w1GGUZWC5f+FQbv+BVwpeOXiVYXTc6HHaqP+N/4VRu6rfvvb3OP2DVxtG
+/ZrkLyA8385eR2RT47vJOJPI7nU05LT31MW9BhXJ65PW04bh/W2HMD9rTpy
upmglef5VhzfFjxdOR28HTzTtb2IH3rg6YNnKKdRwbt8di/g9qScdtGVOZcP
im1jOaV9GNNjZxyPL8rJqPXj5Zt/ivAzlVN2TIOZk89wf6OcNs3reG+QofDN
wSPwLMCzBM8KPGvwbOT09nTqPdtzwrOV06w6mVbnZuJ524EnA2+ZnDZ7954T
puD4DJ4jeM7gucjJIPbde8Oj4vm7yenJm4OzPsaI/R5y8swcO3DBLx7/AM8X
PH855bSfVfZTJY4PBk8OXiR4seAlyilj79lNw8K5/SSnfdOaJswcKe5PCZ5K
Th72a4frzRHH7wMvB7w88Ark5NLbqWDXfyJ+FoKnBq9ITitUGocacfu3WE5/
/Zo7/sla8TxLwCuVU/vNBmV3p4rfK5OT45Ynsvv2XF+TU6s5c8IK+H2+CvAq
wXsGXhV4b+VU0inX2DaF0zd4NeDVymlc4sL+cbNF+q+T06Rjs8pClNz/E04R
O7Pe/u4g+hcahUs9rXCppx0u9XTC6WzNzQke3F/eFjzdcDL7/fJRo1Pi9/TA
0w8n0/92Zl0bJfYbgmcEnjF4JuCZhtNA86BRQ/ZwezOcAj7/Pa/1FnG8OXgE
nkU4rYx6MzSyl3jeluBZgWcNnk04JTxPGhdQyfUV8OzAk4G3LJzONHhWoFUs
6lMO4DmC5wyeC3hu4HmE08voVe+s+X1ur3AK27PsrHUtp3/w/MELDie3l97b
TwTw+wrgRYIXC14ieIpwerozOXftenF9SvBU4TR51LU9uUvF/e8DLyecbi/+
+1X3BiK+5oXT73tvvVMmifmVBeAVhlOS4uzVGRf5/XDwisLp16C0cdFO3B8f
Tqc3pTfosld4JeCVglcGXjl4FeBVhtP16ysW7Dsk7u8ZeFXhtCgwpU+nKeJ5
vgWvGrwa8GrBqwNPI4KCPS0Pp/YV240ipJ5WhNTTjpB6OhFkr2pY7cXz4duC
pwueXgR16KofvE+Pxw8iqMl9tcllfv/TEDwj8IzBM4mg7GHb7NO5fW4Knhl4
5uAReBYRNLp/8nmvJ9wejaARP/Qsa0LF8VbgWUfQogXdJ+V1EOFtA54teHbg
ycBbBp5DBIV1vlRTLhfP2zGCLkw/qvWnit+niaA/e6fNsl0t7t8lgr7fGttz
sQ+Pd4DnEUHOne9c+fWF57+B5wuefwTNzu/1reQfsT8YPDl4kRH0yL59iG0A
t5fAS4ygg18OZGzL5Pn94CnBU4G3D7wc8PLAKwCvEDw1eEUR1OPg8GfpG8Tz
K46gnK9TlQ/4/Z+SCFoR821qY19+3wS8MvDKwauIoGHBn+ObbRfHV4L3LIKO
nwgf12Eoj+eC9xa8avBqImj/ps8/ggLE8bXg1YGnsZt6GnyfePIVl/+76fuz
hztUa4WntVvqae+Wejq7pV5b8HTB09tNVSt9I5re5fmj4BnuprS/5paM3Mrz
iXZTdBtZyEOeb268m7r9fNMr+C3PLwDPFDyz3dRos/bSy7x+gjl4tJvaBYx5
8+mNON8CPEvwrHZTi7U7HWLiuH9qN4VlLL7mxPOpbHZTztVTFz1Hcf0fPDvw
ZLtJJ9XObCj3Jy3bTV3XOXRekMjjWbvpXqbn4RlPxfU7gucMnstu+u95cfja
htzfsZu2uAQU/l4owsMDPC/wfMHzBy8YPDl4keDFgpcIngI85W5yq1hm3zeO
x9PA27eb/k6K+qOpyekfvLzddDene/PVmdzfupuunN5pFOoufr8QPDV4RfXx
o7SjmW1vLv930xzN9G2qLiK+lNSH7za9u3U8X7p0N+U+yvLsOEL4ZeCVg1cB
XiV4z8CrAu/tbvJ+Ocxhrorr/+DVgFdbf79e53RSIzn9g6cRKfUaRUo9rUip
px1J6R0C3087xfPDI6Ve20jSb3HAqM6Ay3/w9CLJcobehhHeYls/kvbp/hdp
0kzk34aR9P7EmThvHm8xAs84kmq62139sp/TP3imkfTQY4jLU34/3Qw8c/AI
PItI6vR3fsL529wfGUk/M+flXP4jPCvwrMGziaS9VoZOp7n9YhtJjQ88bGZw
VGzbgScDbxl4DuA5gucMnkv980zbbrTpMqd/8DzA8wLPN5Jqmzx7YMHlmT94
wZEUYtRf6XuX2/+RlD39c+9eHUX4RUZS3XHPeW1rOf3X/15mq/R2jcR2IniK
SCqze/17AI/fKcFTgbcPvJxIMr2c/ijvpbi/vEjSuDrylfoN97eDVxhJvR/G
bLlqw+1/8IrAK66P34e0rFN3ifApiaTQ8vHBUc95fjZ4ZeCVR9KQ2J7yHTw/
tCKS/pjtPF31v/WwwHsGXhV4b8GrjqSVx5q66oaI/TXg1YJXB55GlNRrFCX1
tKKknnYUrc2Ys3+oD7f/o6Re2yjSvTCtb+xg7v8DTw88/Sjyy+7p0NyM038U
Va7/XnTNgMv/KPryM02rJ49HGINnEkVp9//czt8r2o+mUbQnr2PJv39z+R9F
TQ4d+nbPkct/8Kh+271N9ZFJIv+yAM8SPCvwrKNoRL8Ltw7e4vZ/FO2a/2Nx
9h/ur44iJ7dBI5z/CM8uiny/TnzqXcD9f+Ati6K6SNPqxyvE+qgO4DmC5wye
C3huUXT6rz6yLtOE5wGeF3i+4PnX778yYkbSL27/R9E657FdIm7x+xHgRYIX
G0WXq7xeVPB8+8Qo2unxdkD0EG7/g6cETwXeviiqXqrReK8L96+DlxdFVn2W
fjDg+moBeIVR1HLIraNFPL9HDV4ReMXglYBXCl5ZFDUznWjufJHH68GrqE8f
k/WCh3Xi/nvwnoFXFUVXy09sObFQeG+jSGbauWXfDzy+A15NFPXIuzLRy5n7
+6Mo836vxqXvOP1H0ftLl1utaMH9f9FSr1G01NOKlnra0TT/xozhHT6L83Wi
pV5b8HTB0wNPHzzDaLIbHpIo28flfzS9/PW01fJG4v6No+lX9z2Fmwdw+Q+e
KXhm4JmDR+BZgGcZTUZzBs64xesnWkXTnaFP3DIPcf0fPJtoku97PvGVmYiv
tuDZgScDb1k0/emw5u/fdSJ/coimnf8ezy3k8t0RPOdoOv+Xz5GeA8Xvu4Dn
Bp5HNMUoTB8FzeX+f/B8wfOPpq5fQmoX7uH2P3jyaFp8Z8uMh/9y/z94sdE0
S3O34S3O3xPBU4CnjKbM4Oy4rZf4fQrw9kWTi1X2iXyeX5ADXl40rUyb7XZ/
ibj+AvAKo6mJbKZq5W6u/0dT43tlJvEXuPwHrziaJuevXJ/dg9v/0dRi5OGP
72y4/K/3Nra8XTJRXE8ZeOXR5LP2ZrxNlTi+ArxK8J6BVwXeW/Cqo6ng+/MF
LqtE/K6JpoX33j2anMr1f/DqwNOIoZYRw3N6t+b2f4zU04qRetoxUk8nhkKD
x+aW7+f6fwzp7W81220CzxcFTw88/Ri6sy4yosSKy3/wjMAzjiH5Sdudmvq8
HgF4pjEU8rBt68vm3P8Hnjl4BJ5FDMXnD7rYPIT7/2JoUFrS1vRQHv+Lodrw
9Xe3jBbPyxo8mxhy09u+59oQrv/HUOLk8x3G8XrfduDJYqhRumZJS16vfFkM
LW7607FTvoi/DuA5guccQ38vat/6Jq8P7AKeWwxN9uv0aWU69/+B5wWeL3j+
4AXHULhG7d2Rz3l+PHiR4MWClwieAjxlDBUNiR6dye1fVf223aJzP+zF89gH
Xg54eeAVgFcInjqGypL+szjN62MXgVccQ/1dCo9oVYj4VQJeKXhlMaTVYVrz
bx+EVw5eBXiV4D2LoY4WGa9b8/rpVeC9jaHbq/p++Sef+//BqwGvFrw68DRi
pV6jWFr58vWelx7c/o+VetqxtHDt94nD2nD5H0tNjIKvzOf5BG3B042lvyN/
HEjJ5vY/ePqxNEbr0fpPX7n/P5amjh3oepfXzzACzxg8k1hSapYemHtSeKbg
mYFnDh7F0r3mo36M4/WGLGJpy/odN17x+jyW4FnFUrbHi7YNnHi9K/BswLON
pZtuU1Tjh3P5D54MvGXgOcRSfsLiYWbc3nGMpbzvTR2ucPvZGTyXWPoYXBn+
uAG/71bvHY7p0eM4p//6+9ObZdPuf+v5g+cLnn8sPagrnNfalev/4MljqWqF
+loTBbf/wYuNpbXDPK83juXxP/AU4CljqWGy+4rDCu7/iyUNpdGRe7s4/YOX
A15eLL1c0W/jpDSebwVeIXjqWNoc+H1SVAzX/8ErBq8klho4v5iUsprr/7F0
aV5y9MvunP5jyTa5xlLuy+kfvArwKsF7Bl5VLAV9c1g77ZcIj7fgVcfSm4Co
OOLxjRrwamMpNikoySNU/F5dLI0dnGj9ZiKX/3FSr1Gc1NOKk3racdR+3N9T
OiVx+R8n9drGUd+usZc3fOb+P/D04ujOkz7hhrweqz54huAZxZHMLPjUNF7v
whg8E/BMwTMDzxw8As8ijtK/bOk+l/ufLeNo/aSRpyKmi/RkBZ41eDbg2YJn
B54sjoZcyzBOHsr9f+A5gOcYR+PftTxz9TyX/+C5xNFRh1GZcg9O/+B5xFFt
pwZyr8Y8/g+ebxxlnuv2Vz+uj/mDFwyePI5+xiyz6r+X2//gxYKXGEc3tK+4
mdtw+R9HOz1zP8Y2FPFdGUfh/dU3OvD6CCrw9oGXA15eHOX7PYl4H8Xt/zjK
CF5T5mPO43/gqeuf7/2AeI9s7v8Hrxi8EvBKwSuLo2lhi43W8vth5eBVxJFt
k21VZn48/we8Z3E09ciuBx15fKoqjlrNfdekqRmn/zjq5H9w7GADTv/g1YBX
C14deBrxUq9RvNTTipd62vF0oeeTwubPef5fvNRrC55uPM0Y+ep7Y34/RA88
/XiS3R6y4cUCEf6G4BmBZxxPvX7Yb9vWh8t/8Ezj6dN0g/wb4TzfHDzzeEru
pJh63JL7/8GziKcfmiVXL90T25bgWYFnDZ5NPA0ydbvVoYDTP3h24MnAWxZP
i/Ov5y/9l+f/gOcInjN4LvF00bdZyKRPPP8nno7U3ZoSNlzEXw/wvOJpyYou
Jk2P8/gfeP7gBcfTbNeTfXYmcPkfT3s3dVZXuvP8H/Bi4+n7i/WHVL48/g+e
ov56gmWt3+hy+gdPBd6+eHJO2n4vn9cry4mnqz0fOvRNFtt59feb12hUDs+f
LIgn92VZ+03Wc/kPnjqevl4f/vL8PtGeKwKvOJ7M1WMXVWwS91sSTzrfazqs
esvjf+CVgVcOXkU8Te167YYHz/esBO8ZeFXxlNL9bEohj1+8Ba8avJp4+rX5
aUffq1z+g1cXTw70csy/FTz/J0HqNUqQeloJ9KhRu+7fXER4aCeQeUz2WePH
Irx0EqReW/B0E+jvd2lbvkwV16+XQK3v9hnfkZ+nPniGCVStHm7mOp/n/4Fn
nECWKx+/NmrJ/f8J1Ct2zTY/P67/J5B6exfz6by+txl45gm0aWGvFKdI7v9P
oKrNl3cc4/VoLcCzTKB/ly6Ye6YHt//Bs06gvVuCkvW4/LIBzzaBar3+a596
msv/BJJn+LdY9orH/8FbBp5DAoUYxBxy4/UVHMFzTqA83ZWx5tu4/E+gL4Xv
++7n+bdu9c9j7urwx224/p9AcTkPnp7i8Vov8Hzrz18Zu/K9jN+HTaBVn4NN
H87i9A+ePIFc8yYFLNzF5T94seAlgqeov7+F4+ynPeP6P3gq8PYlkElkzMDe
r3n8D7w88ArAKwRPnUBb/IdWW5ty+z+BXj28J7txgcv/BNJdPaz/aJ6/UpJA
b189zEk5yOU/eGXglYNXAV4leM8SaKDN9PN+93n8P4Ge2pZPt+T1V94mkNeX
7F4PXnL7vz5+jgrv82oI9/+BVwteXQJ935LRbDDXNzUSpV6jRHo7Y/A5azPu
/0+kdY/HDisJ4fZ/Iv07YbRzEa9/ppMo9dom0gvZDgNnNff/g6cHnn4irUm+
vn42v59qmEj6m/K0W4dy+Z9ID51e9z5vzf3/4JmAZwqeGXjmidSvbMXi/fy+
NIFnAZ5lIq1vq3nx9mlxP1aJZH72/XXDXrzeYyJ9abK33w4uX20SKazJ+e/2
bTj9g2cHngy8ZeA5JNKA+PTCJQPE+Y7gOYPnkkhRf5ZcOteL5/8nUonGlj+2
U7j/P5EWXBjndJ7f5/NKJHdltm6hhng+vuD5J5LBHXL59xynf/Dk4EUm0u69
h9ok5YjtWPASwVOAp0yk/AxZ3228/pgKvH3g5SRS0Mbbf9u24/SfSI7DQtxf
D+f5v4nUWbvpqH3FPP83kTRObNP10eD2P3hF4BUnktJeEXWW1xMqAa80kT4b
NpE9ucXpP5FORh5cd4/rq+XgVSTSib7bmycOFvdXmUhOA4YP/8PzMZ6BV1V/
fWNXVR3k9QjfgledSH+NHzuiDX//rga8WvDqEsmCjl9//i+X/0l0fMiRZTN2
8/hfktTTSqJlp/MGpZtz/3+S1NNJot75PZ4alnH5D54ueHpJZD/Xb8zo6fz+
GHiG4BmBZwyeCXimSVQ1NLFhwb88/geeOXiURM8eJR/9wv2XFuBZgmeVRC1j
n6Xs5/f3rMGzSaI9a5xdBvcR59smkSyg1KR3H07/4MnAWwaeA3iOSdRtY5uK
hY35/TbwXMBzS6K7E8rmT+X50R5JNLDXgUtJbtz+T6LWBiPuvuH3ZXyT6OLf
76zde3L9P4k69NT87BTG43/gycGLTKLok0mze/P7e7FJdMSzRbdG8Vz/B08B
njKJ1n8LW9+d3zdUgbcviTxMLiwZw+tj54CXB14BeIXgqZPod1JyfPlPcb9F
9c/bfciLJvx9wmLwSuqf99JxO8bx+2ql4JWBV55EmgVDbtk85vl/4FUm0Vm/
DtGPk8X+Z+BVJdG71cenOW7m+j941eDVJFGk3udWvo15/i94dfW/VzC+l6oD
z/9LlnqNkunQA7MOtYu4/E+WetrJUk8nWeq1BU8XPL1k2t765rYDvN6UfjL1
GNOz7eUTIr4agmeUTIrZNRMN0rj9n0wTm4+adpfX5zEBzxQ8M/DMwSPwLMCz
BM8qmXJbp9mP2sT9f+DZJNOABelzWj/n9crAswNPBt6yZPolszF8yd8rdADP
MZk+TXp05AKPTzsnU7vgk/tbZHH/P3huyXSn44Cbc37y/D/wvMDzTaZ93131
Xa7x/H/wgsGTJ1NXo4WnnE3FdiR4seAlgqcATwmeCrx9yVT8/OP51iNF+zoH
vLxkKoto3mQrt68Kkslw/by+NbyefSF46mRyyz5+MZzHu4rAKwavBLxS8MqS
qWZO8PeH+lz+g1eRTD8GvL+0gt9PqATvGXhV4L0Frxq8GvBqk6liwROftTzf
oC6Z3qS3f7b3GY//7ZF6jfZIPa095Hj76lNfft9Qe4/U09kj9dqCpwue3h7K
uJ7Xb99Qnv8HniF4RuAZg2cCnuke2l1o+rd3Cc//A88cPALPAjxL8Kz2UK8J
Z2V6E3j+zx66d+XBhCmuXP6DZwueHXiyPZSU8Gv77o+c/vfQ2C471g1qyfN/
99DG9jm7nA/x+797yP1G1fPd/H6T8x4yaGPc9vth7v8Dzw08D/C8wPOtv941
jjMf8Pdn/ffQ4XyTuz3O8fj/HhqyWmZ2m9eXl4MXCV4seIngKcBT7qFPNgua
PnvI83/B27eHzljF/Uk8xOtrgZcHXsEespBfi25zmMt/8NR7aHhKUtVgXp+h
CLxi8ErAKwWvbA/dcTKJud1CnF++h6x3e2ff5e+TVYBXCd4z8KrAewteNXg1
4NWCV7eHZLfvLSzm9SY1FFJPU0GHZxS/PbmS2wMKmhgy5NFiHXF8EwUVbOj8
aWklzw9SSL1m4GmD1wI8HfBaK2htl/KoS57i+LYKstLKu5VtJI5vr6BrT9wP
5fJ8H13wOilo0vCniTN5voseeF3B0wfPQEFDSrPymvHxhuD1As8IvL4KGuv/
J/TuF56vDN4A8EzAGwSeKXhDwTNT0OJEg4zAY2L/CPDMFRQSFH9xA7/vPRo8
Am+cgva2O558kNtXFgq6PtdR86euCJ+J4FkqqGvZ+q3rnUT5Phk8KwW9PFZt
FPRTbE8Dzxq8meDZgDdHQa72/112HMvtKwVtczzqvjFVePPBs1PQkjkPdh3y
F8/THjyZgkZ82vxkPH9fYgl4yxTUKsug2pvXK1sBnoOCFr7WfOL2VcT31eA5
gucEnrOCXh15NMJxtrjfDQqy1DRc6tOb23MK+jXIoZtqnNjvCp4beO7geSgo
8u6yX79yxf16Kqh1brQyaR2P/yio6pIi2Jjnm/iA5wueH3j+4AUqSKXdv3/g
GPE8gsELVdD5wODByjncf6ygRvYrbbR5fksEeJEKmvBh691x/L5JNHix4MWD
lwheMngK8FLBUyrIc8LtzfOWifSSoaAFLwIX5Zdw/7SCnlywi9B8J55PNnj7
wDsAXg54ueDlKWjXrPxZA+aK55kPXgF4x8ArVFDU+Cyf7bvE8zgJnhq8M/Xx
YUXajFxeH6hIQRvIMrCOxzPOg1cM3iXwShSkG3/pl0mOOP4qeKXg3QCvTEFe
kcVuX7i9dBu8cgX5fbv+Qsbf57kHXoWCVjT8k/SH23cPwausv97c86+nFYn9
T8B7Bt4LBYW1u9wuIIzn5yroQPTY7XvMhPcavLcK2v3zRcyIxWL/e/Cqwfuk
oK+zS43f8PzGGvC+glcL3g8F7TQ7mX0tRNSf6urL343WV4dz/eA3eBopVHby
sGXyEbFfM0XqNUqRek1SpJ5WCg0L/9C1315x/c3A006hC716dDo6WcTPFuDp
gNcavLbgtU+hjn36DBnYlPsXwOuUQorvVQeD+Hu2eim0XN3lrpzH57uCpw+e
AXiG4PUCzwi8vuAZgzcAPBPwBoFnCt7QFKpu+PNS8Rdu34A3IoUG3Tzvemgo
t3dSaGmDX4t67RHPfzR4BN448CzAm5hCOYa/d0xrwuujpNDx57aOFdy/MzmF
bBMTXz3m9dKtwJsGnjV4M1Mo+/v3hrNfiPhnA96cFJpoNWO2KX+f2ha8+eDZ
gWcPngy8JSm0b+v90C48n3tZCr0+0fWzaTtR31iRQgP+/No79r4ITwfwVoPn
mELheyPCPA+I63NKIeNBlxa68vqCzuBtAM8FPFfw3MBzB88jhdIbJ7keGiLu
3xM8L/B8wPNNobRTrQ8H8/wNP/D8wQsELxi8UPDk4EWAF5lCP18uanTtvHj+
0eDFghefQj6PnK/P3cPzzVPo/ch/zA694PIfPAV4qeApwcsAT5VC+blJk0ck
iP3Z4O1LIYPFEw92/y3u70AKjW3YYW0qz9/JAS8XvDzw8lMoZPTiWZ8+cvsW
vGPgFabQf+NSQh64Ce9kff6+es9TP16fXQ3eGfCKUqhJYffHZ3qL+HQevOIU
2vB4zNPMaJGfXAKvBLyr4JWm0M02b5wWn+LyP4U+635vp6rh/nTwboNXnkLz
tE3Hd+f+uXvgVYD3ELzKFJqjyE3twN8reALesxTqccEi0Iy/B/4CvKoUevN2
VnBj7v95Dd7bFLJe8M/Qq1vF8e/Bq04hmfWNki+J4vl9SiH10zpf9W2ez5NC
c9e0mLQ1WIxnfq3fr8zcYXGa5/enUP9Pb6yv/uLyH7w68H6Dp5FKWTPOTNEO
EtermSr1GqVKvSap9OPcn+TdSVz+p0q9ZqkUNDvc8iHP59EGr0UqtRp3xNHX
j+cPgNc6lYyaGm3r/POGKP/Baw+eLnidUmnT9XFtFzQR+/VSqWj08ambef5q
11RqPvTS3BHfef0x8AzAMwSvV/3293XFx67z+AR4fVPJ/tZwu6/8/okxeAPA
MwFvUCpFFfQYlMrrmZuCNxQ8M/BGgGeeSn/5XjutyevFjk4l76q/ZD8O8PtP
4I0DzyKVtly2m7ZprYhvE1PpwLeB2xeZc/9oKp3eNL+0RzMRnyeDZ5VKsy4M
OvbhP3H901LpzpXwGUVPefwkleKvjE78+UKE70zwbMCbA54tePNTqaBFyPVm
3YRnl0qH+191avmXuH978GT18aNV6ykT+whvCXjLwFsBngN4q1MpdOQS+/38
fopjKvU6PiM8YbPY71R/vxvbrY8czuM34G1IpWLd2w0rrvF8LvBcwXMDzx08
j1RqH5eiWsPf//EEzyuVHvfoujn9hAhfH/B8wfMDzz+V3NcZxm46J84PBC84
lWaOXtTuAa9fHgqePJWc8g70se0h7icCvEjwosGLBS8+lca+71tY3ozb/+Al
p9L0mvBh+xXc/5xKe/OmFn5ozuV/fXq7nxVWxt+3UabS4PtbR3TrL+JHBngq
8LLB25dKbbr+/N6/nMv/VNq8xUxfK477r8HLBS8PvHzwClLp+Lnh0xL5+5nH
UinVztX/f+9fF6bS+9JrA15niPR6Ejw1eGfAKwLvPHjFqbRx3xGZD8/nvQRe
CXhXwStNpVt5O48O/U94N1JJ/2Xv0Vu5Pl8G3m3wysG7B14FeA/BqwTvCXjP
wHuRSquHLLXwn8/z88B7nUrmOi+djrcW4fsWvPfgVYP3Cbwa8L6mUt1bn7sb
eb3/WvB+gFeXSg7dw1acbczlP3gaaVJPM43yq42nK85w/38aRX/pTf5q4TVJ
k3paaZQ4IWHv7gjhNUuj5Yez5209z+uFgtcCPJ00Cpp85/wvnt/QOo2OFHU6
3+iFWH+3bRpNfu095zh/H6Q9eLrgdQJPD7yuaTRUnXvK/AbPPwTPADxD8HqB
Z5RGhXWtB7nxeox90yijs96S7ByerwTegDRqfNP104Q5/P5yGv3dS9e4WaBY
n2NQGjkGu+8O/8DrmYI3NI3Un9Vjh3P4moE3Io3WXbR47xkhts3BGw0egTcu
jVq+vDZ0KvF6KGl0dc3r9vaPuP8/jbScq4x/veDvNYA3uf78Vmv02vL3iazA
mwaeNXgzwbOpj589O33sxfNr54BnC978NNrsPbNnQT+ebwWePXgy8Jak0aaI
g+s0XLj9n0a5rc4t3DJdbK8AzwG81eA5plH29DV+H/l7vU7gOYO3IY0u3Quw
6ndW7HcBzxU8N/Dc02jHoNHHtZ+K+OkBnmd9/K3+5OrP3x/zSqNmqbrqzfw9
Jx/wfMHzS6OzFX0Xx2fx+hDgBaaRs1r7k8ElHi8GLzSNxuc45iRZC0+eRj1/
pLfsx+uDR4AXCV50Gi1e7G/r4M/lfxoFRgXfmjFRnB+fRj/bFpqsmcrzzcFL
TiONExmrXHm+gSKN/sr23j5ER8wnSQVPCV4GeCrwssHbB94B8HLSqEVuxb0J
/P2IXPDywMtPo2Wjy0w/m/D8FvCOgVeYRuWd17x4sJPb/+CpwTsDXhF458Er
Bu8SeCXgXQWvFLwbaTTq7vPrKzfz+7RplNPR2foVrw92O40sPGbpqj15fb00
2hnR7NJe/t73PfAq6vP7F+8nP7EQ3kPwKtPodI+Cvlu3ifzySRopHX67Hpsj
/GfgvUij2UVWU/fw+upV4L1Oo99+XWQmlVz+p5GLtlZx0+bi+PfgVafRxXUb
ultUcv8/eDXgfU0j63en8w24PVebRvPLJt9ydOH2fxpNX+Pv+iaD5/uB9xs8
DaXU01SS4sQZ/WMjufxXSr0mSrLLHO//31oe/1fSiTGH/1noLsK3GXja4LUA
T0dJ/Uv/yVLP5fF/8NqC1x48XfA6KakkrnagywWebwheV/D0ldQitOeaPTy/
ygA8Q/B6KemG2Wmtx7z+s5GSWgd9sfvG8637gmcM3gDwTMAbBJ4peEOVdD6l
9y298Tz+r6TtqwND42Zw+x88c/BGKylG33TbUf7+Jynpjuz3vycecPkPnoWS
iv7RNN67RJw/ETxL8CaDZwXeNPCslbRRr3pkl59i/0wlrcyPXvlxKY//gzcH
PFvw5iupNnDx4F7BPP4Pnn29l3xmTOhhsV8G3hLwloG3AjwHJV361/1dgpLL
fyV1bH3TqrU39/+D5wSeM3gblLQ+vVhzMK8H6aIk79nrI9z4fUhX8NzAcwfP
Q0mPZrp8W8XzlT2VZPs2WL76jciPvJQ0umPHbr4/efwfPF8lDdbQclzXVJzv
p6SKdc8/xL8X+/3BCwQvWEk+t/ffKLskzg8FTw5ehJIeurkUOg0X+XUkeNHg
xYIXr6R4DasrK3l9hUQl/Rjlui5iLJf/Sgrpu//3P+7c/lfSzoxbZ1fqi/ww
FTylkmZOWOzQ5Jo4PgM8FXjZ4O0D7wB4OeDlKskv4LvN/95XyFOS67w52Y1K
xf588AqUFDnt6Y2OM0V4H1PSoLQHLXx4vb1C8E6Cp1ZS5rVBd2e/5v5/Jf16
8nqNbwmP/4N3HrxiJf186N3Y9Kn4/UvglSip9F/ljS+m3P5XUpq952kNfl6l
4N0Arwy82+CVK+nu0dDiEaEivtwDr0JJJzUG3Pp9SngPwatUUpbjX+o5y3n8
X0ndWyz7p86f1+MC7wV4VeC9VtKZc3PmHU3i8X8lHQ0Kf3pzvjj+vZL0dhTW
/uL1c6vB+wReTX3+FDZAL3Ct2P4KXi14P5Q0++Qpt7Y8vlCnpGf/mq8P4PXt
foOnkS71NNOlXqN0qdckXepppUu9ZumUcbCDfAZ/n1UbvBbg6YDXOp1Cb2pE
up/n9wnBaw+ebjods3yxN06Xy/90CtlYot/dj+f/1R9/fclxnX08/w88ffAM
0ikxv1mGLvcHG6bTqqHzKzz4fZxe6VR8Uqeo7gd/7xi8vun0+fmPS6VPRfwx
TqdFs11expHwBqTTrbF9ujTk/NUknd411vse31fUnweBZ5pOreVhPfd1FvuH
gmeWTjtDI9/eM+b2P3jm6fRB++MYDR8RH0enU82wpzq/zvD853Q6o0wy3dVB
bI9Lp8zgHI2uu3j+H3gT00lmOKzPwUni9y3Bm5xOXdvNSnsp5/IfvGnppPiT
/9hvEvf/gzcznaZu+LVxI+eHNuDNSad161sf+esmv18J3vx0ejG+4z+KLdz+
B88+ne4ovzhtu8/lP3hLwFsG3grwHMBbDZ5jOn18Hh3VhceXnMBzBm8DeC7g
uYLnBp47eB7pNGHRQeOrHlz+g+eVTn+C3zg8eyXimw94vun0Zur3na84P/NL
p5lut+694u8t+afTkUFdI4tNRXoIBC8YvFDw5Om094VRRfeP3P+fTsbzqmYZ
DOP+f/Ci02lQ6eo2Ou25/AcvHrxE8JLTaVn5lSw7Wx7/T6co40l7Iq+I55UK
njKdoh2P7mheLsrzDPBU6dTgTkqDIQ9E+ssGbx94B9LJ/Hnw9X9a8Psr6XQh
Ieb1Rk4/ueDlpdNB21t7V/qK9JGfTq5z9mhPJx7/T6eJvYZM1ggW/rH68P42
3DOnj/AK0+n2qBbxXXl9pJPgqcE7k07aTis6NtjA/f/1+Utxu8zu/H298+mk
tfvMoBlXRXuqOJ3KrWK8ll4Uz+cSeCXgXa0/fl+Z1lENnv8H3g3wysC7XZ8/
J0dFNOL19MvBu5dO/c6sWn6V10utAO8heJXplHvRd/o2GZf/4D0D70U6DRj6
dsH7zfx+DnivwXsL3vv68tJ+0dd3/H3IavA+gVeTTqM+JY/LqeL5f+m0ZvHy
fr9+iONrwftRX778c/ft6khu/4P3GzyNDKmnmUHTr5c8mMXfE2mUIfWaZFDD
1eMc7JO5/M+Qes0yqI9rC39D/r6mdgat+x1g8oX3t8igibnfPvzm+bc6GRRQ
psp4WsblP3htM+j7T52LG3n+QnvwdMHrBJ5evXfwq8EOe7G/K3j64Blk0Laa
Uy+vbeD+f/B6gWcEXl/wjMEbkEHeS++bbt8p8icT8AZl0I5XiVNN/7eeYQbd
K6wOmftd+EPBM8ugtECPP4FPeP4feObgjQaPwBsHngV4E8GzzKBxkeZlVbx+
5uQMMpyjFa+9QmxbgTcNPGvwZmbQbJMvOopWvJ4KeHPAswVvPnh2GZS7IPRI
yWORn9qDJwNvCXjLwFuRQdZdXfpf5e9hOWRQ+wcet6/w+7erM+iTfHPInI48
/g+eUwbturP+8ip+390ZvA3guWRQgenn4DarxfW7gucGnjt4HuB5ZpBKdrBB
8gex3ws8H/B8wfMDzx+8wAxK7/dePaqW5/+BF5pBHV/ZbjnxTuRXcvAi6uOT
7HqrUIXYjgQvGrxY8OIzaNi5iDjdtuL4RPCSwVOAl5pBxv5bu997wf3/4GWA
p8qgj3uWejzj771n18enlr2tYyeK+LIvg6Y2Tmyk/Mjz/zNo0OwzWZeOCD8H
vFzw8sDLB68ggzosKr62k+cPHAOvMINORZzQWnqd2/8ZdGiY6lXYCS7/M6hq
4Md/3PqK48+AVwTeefCKwbuUQSH6r5L+mcz9/+BdBa80g84mOHY0txT3cyOD
FHM2jv55ncf/wbsNXnkGuVx3PWLRUKSPe+BVgPcwg4qdskaeaCautxK8J+A9
qy/PrPu8LuXv/b4Ar6r++mZkagRx/eY1eG8zqM3yPJOcljz+n0Ezqjb0qub1
Z6rB+wReTQbVxTx6E7JQnP8VvNoMuntn+qMIL3F9PzLoW9hHA61uXP6D9xs8
jUx6/Cn6ujnnt5qZUq9RptRrkin1tDKlXjPwtDPpaKvasi9GXP6DpwNe60yK
8t7UwY7Xu2mbSaWLG1dH1ojwaA+eLnidwNPLpLTNc5xfDhTbXTOpxVGN49f5
fUV98AzAM8yk2zO+xWfx92R6gWcEXl/wjDPJtMee/kFruP8fPBPwBmXS94Ge
w25Wcf9/Jrk7+A7oy/0ZQ8EzA28EeObgjQaPMum8TkSP41t4/h94Fpl0Mv/o
p5HFwpsIniV4kzMp/eEG4w4RPP8PvGmZ1CbH9lDhAC7/M8kk/txpJ/7+2kzw
bMCbA54tePPBswPPHjwZeEvAWwbeikyapjtq7FQdEb4O4K3OpFO32+28+pv7
/zOp54JKy3e8HrQTeM7gbQDPJZOGDGzfv8NLnv8Pnht47pkU/qLN7cu8npVH
Jp27fVFHwesHe4LnBZ5PJtmtPdL44Xju/wfPDzx/8ALBC86koIMvhoemifZX
KHhy8CLAiwQvOpO+tm/YIuEBv18OXjx4iZnUYWGnyR+3cf8/eArwUjPJu+8d
7emWPP8fvAzwVOBlg7cvk1Zl/Vj9t5lo7x8ALwe8XPDy6uPjLq28aR1E/pcP
XkEmnVg1t/nHblz+g1cI3slM6nT9ye7j4fz9W/DOgFcE3vlMWv9h1slJXF4V
Z1LY9pwmQR7c/1//e1Yrf7c9JcK7BLyr4JVm0hfHw3du9uP5f+CVgXcbvPJM
sk6Nku/l8YR7mXSsdf7Nfz3F/VXU529fdHsvKuX+/0wq3rppa//eonyqBO8J
eM/Ae5FJ1/1G3NTm9aGqwHudSY+anpsh8xLHvwXvPXjVmdRIUZBjc47b/+DV
1P9e93+3nloljv8KXm0mbZ+90HWLvtj/A7y6TDq44YW6h7u4/9/gaaiknqaK
9tvXdBjI88MbqehY8fhLMUVc/qvo6KYhrsVcvmipyP9wlWyBCY//q6i49wjl
mxhu/6uo4fqf31vx955bgKcDXmvw2oLXHjxdFS050j3DkL//3Ak8PRU1su3o
dLY5j/+Dpw+eAXiG4PVS0cNnrVpqLeTxfxW9rvN0muDN5T94xuANAM9ERXty
u4Xu4PfJB4FnqqLK9dbbV2/m8l9FazIU69fw9/rMVPSt632dOP7e0QjwzMEb
DR6paNl5b+NF/cT541Q06di/8tXdRXqwUNGCBgVZ6/h7CBPBswRvsoqmxF4o
/zlPbFup6NSDfM14d5G/TAPPGryZ4NmANwc8W/Dmg2enIsMaj432g7n8B08G
3hLwltXHn+B2+x9yfFmhorueg26/vcDj/+Ctrt/fMMFufDav16aiyx2Pf9Cf
xuW/imrHzx1gyN+rdFbRk3Z3LjzaLY7fAJ6Livrmfbp0LUnEX1cV/ewzfHTN
MR7/B89dRVrDFbZ+07j9D54neF7g+YDnq6JXqy2Cv/H7/H4qKlMljndez+vD
qahN1BnDSQ7CCwQvGLxQ8OTgRYAXCV60inZ1cxvVZL14/rEqMlFqmdls4Pf/
VLQobUULh9Fc/oOXDJ5CRZFnAkyb/uD+fxU5rvEvUnfk/n/wMsBTgZetopmZ
XcucrnL7H7wD4OWAl6uiqgAP5XhevzlPRe0PfvDY0FuEZ76Krret1LzA62UV