-
Notifications
You must be signed in to change notification settings - Fork 10
/
Copy pathimportanceSampling.nb
2236 lines (2147 loc) · 109 KB
/
importanceSampling.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 10.1' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 111371, 2227]
NotebookOptionsPosition[ 108066, 2108]
NotebookOutlinePosition[ 108423, 2124]
CellTagsIndexPosition[ 108380, 2121]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Solve", "[",
RowBox[{
RowBox[{
RowBox[{"Exp", "[", "x", "]"}], "\[Equal]", "1"}], ",", "x"}],
"]"}]], "Input",
CellChangeTimes->{{3.7016881162165813`*^9, 3.701688133508973*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"x", "\[Rule]",
RowBox[{"ConditionalExpression", "[",
RowBox[{
RowBox[{"2", " ", "\[ImaginaryI]", " ", "\[Pi]", " ",
RowBox[{"C", "[", "1", "]"}]}], ",",
RowBox[{
RowBox[{"C", "[", "1", "]"}], "\[Element]", "Integers"}]}], "]"}]}],
"}"}], "}"}]], "Output",
CellChangeTimes->{3.701688134065237*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SubsuperscriptBox["\[Integral]", "0", "1"],
RowBox[{
SuperscriptBox["x", "2"],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.702764742919292*^9, 3.70276475459186*^9}}],
Cell[BoxData[
FractionBox["1", "3"]], "Output",
CellChangeTimes->{{3.702764743689941*^9, 3.7027647551185083`*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SubsuperscriptBox["\[Integral]", "0",
SqrtBox["2"]],
RowBox[{
SuperscriptBox["x", "3"],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.702637081224292*^9, 3.7026370932091007`*^9}, {
3.7026670051522837`*^9, 3.702667005456141*^9}}],
Cell[BoxData["1"], "Output",
CellChangeTimes->{3.702667006236586*^9, 3.702761332931199*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{
FractionBox[
RowBox[{"2", " ",
SqrtBox["2"]}], "3"], "//", "N"}]], "Input",
CellChangeTimes->{{3.702637097711726*^9, 3.702637098192266*^9}}],
Cell[CellGroupData[{
Cell[BoxData["0.9428090415820634`"], "Input",
CellChangeTimes->{3.7027613377238913`*^9}],
Cell[BoxData["0.9428090415820634`"], "Output",
CellChangeTimes->{3.702761337813093*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SubsuperscriptBox["\[Integral]",
RowBox[{"-", "\[Infinity]"}], "0"],
RowBox[{
RowBox[{"Exp", "[", "x", "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.7016885482178497`*^9, 3.7016885780336742`*^9}}],
Cell[BoxData["1"], "Output",
CellChangeTimes->{{3.701688563073257*^9, 3.701688579067445*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SubsuperscriptBox["\[Integral]",
RowBox[{"-", "\[Infinity]"}], "0"],
RowBox[{"x", " ",
RowBox[{"Exp", "[", "x", "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.701688595953409*^9, 3.7016885966642857`*^9}}],
Cell[BoxData[
RowBox[{"-", "1"}]], "Output",
CellChangeTimes->{3.7016885971841516`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SubsuperscriptBox["\[Integral]",
RowBox[{"-", "\[Infinity]"}], "0"],
RowBox[{
SuperscriptBox["x", "2"], " ",
RowBox[{"Exp", "[", "x", "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.701688622223023*^9, 3.701688622567279*^9}}],
Cell[BoxData["2"], "Output",
CellChangeTimes->{3.701688623247532*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{
SubsuperscriptBox["\[Integral]", "0", "1"], " ",
RowBox[{
RowBox[{"Exp", "[", "x", "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}], "//", "N"}]], "Input",
CellChangeTimes->{{3.701689018904676*^9, 3.7016890386616173`*^9}}],
Cell[BoxData["1.718281828459045`"], "Output",
CellChangeTimes->{{3.701689024120564*^9, 3.701689039202794*^9}}]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"f", "[", "x_", "]"}], ":=", " ",
RowBox[{"Exp", "[",
RowBox[{
RowBox[{"0.4",
SuperscriptBox[
RowBox[{"(",
RowBox[{"x", "-", "0.4"}], ")"}], "2"]}], "-",
RowBox[{"0.08",
SuperscriptBox["x", "4"]}]}], "]"}]}]], "Input",
CellChangeTimes->{{3.701756629495447*^9, 3.701756631755775*^9}}],
Cell[BoxData[
RowBox[{
RowBox[{"Q1", "[",
RowBox[{"x_", ",", "s_"}], "]"}], ":=",
RowBox[{
SqrtBox[
FractionBox["1",
RowBox[{"2", "\[Pi]", " ",
SuperscriptBox["s", "2"]}]]],
RowBox[{"Exp", "[",
FractionBox[
RowBox[{"-",
SuperscriptBox["x", "2"]}],
RowBox[{"2",
SuperscriptBox["s", "2"]}]], "]"}]}]}]], "Input"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{
SubsuperscriptBox["\[Integral]",
RowBox[{"-", "5"}], "5"],
RowBox[{"10", "*",
RowBox[{"Q1", "[",
RowBox[{"x", ",", "1"}], "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}], "//", "N"}]], "Input",
CellChangeTimes->{{3.701756845535008*^9, 3.701756901499094*^9}}],
Cell[BoxData["9.999994266968562`"], "Output",
CellChangeTimes->{{3.701756890391841*^9, 3.7017569021129503`*^9}}]
}, Open ]],
Cell[BoxData[
RowBox[{"Clear", "[", "x", "]"}]], "Input",
CellChangeTimes->{{3.701813872893661*^9, 3.701813875688485*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"\[Integral]",
RowBox[{"10", "*",
RowBox[{"Q1", "[",
RowBox[{"x", ",", "0.1"}], "]"}],
RowBox[{"\[DifferentialD]", "x"}]}]}]], "Input",
CellChangeTimes->{{3.701756971817506*^9, 3.701756981131122*^9}, {
3.701813862708601*^9, 3.701813863536263*^9}}],
Cell[BoxData[
RowBox[{"5.`", " ",
RowBox[{"Erf", "[",
RowBox[{"7.071067811865475`", " ", "x"}], "]"}]}]], "Output",
CellChangeTimes->{
3.70175698167822*^9, {3.7018138681134663`*^9, 3.701813878193446*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Solve", "[",
RowBox[{
RowBox[{
RowBox[{"5", " ",
RowBox[{"Erf", "[",
FractionBox["x",
SqrtBox["2"]], "]"}]}], "\[Equal]", "y"}], ",", "x"}], "]"}]], "Input",\
CellChangeTimes->{{3.701767532545272*^9, 3.701767545922358*^9}}],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"Solve", "::", "ifun"}], "MessageName"],
RowBox[{
":", " "}], "\<\"Inverse functions are being used by \
\[NoBreak]\\!\\(Solve\\)\[NoBreak], so some solutions may not be found; use \
Reduce for complete solution information. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", ButtonNote -> \
\\\"Solve::ifun\\\"]\\)\"\>"}]], "Message", "MSG",
CellChangeTimes->{3.701767546869196*^9}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"x", "\[Rule]",
RowBox[{
SqrtBox["2"], " ",
RowBox[{"InverseErf", "[",
FractionBox["y", "5"], "]"}]}]}], "}"}], "}"}]], "Output",
CellChangeTimes->{3.701767546880191*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Solve", "[",
RowBox[{
RowBox[{
RowBox[{"5", " ",
RowBox[{"Erf", "[",
FractionBox["x",
SqrtBox["2"]], "]"}]}], "\[Equal]", "0"}], ",", "x"}], "]"}]], "Input",\
CellChangeTimes->{{3.701756997375812*^9, 3.70175702127534*^9}}],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"Solve", "::", "ifun"}], "MessageName"],
RowBox[{
":", " "}], "\<\"Inverse functions are being used by \
\[NoBreak]\\!\\(Solve\\)\[NoBreak], so some solutions may not be found; use \
Reduce for complete solution information. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", ButtonNote -> \
\\\"Solve::ifun\\\"]\\)\"\>"}]], "Message", "MSG",
CellChangeTimes->{{3.7017570107287397`*^9, 3.7017570216063843`*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"x", "\[Rule]", "0"}], "}"}], "}"}]], "Output",
CellChangeTimes->{{3.701757010732788*^9, 3.701757021621243*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Solve", "[",
RowBox[{
RowBox[{
RowBox[{"5", " ",
RowBox[{"Erf", "[",
FractionBox["x",
SqrtBox["2"]], "]"}]}], "\[Equal]", "1"}], ",", "x"}], "]"}]], "Input",\
CellChangeTimes->{
3.701757025373885*^9, {3.701767339978753*^9, 3.701767342242855*^9}, {
3.7017673931072283`*^9, 3.7017673931775503`*^9}}],
Cell[BoxData[
RowBox[{
StyleBox[
RowBox[{"Solve", "::", "ifun"}], "MessageName"],
RowBox[{
":", " "}], "\<\"Inverse functions are being used by \
\[NoBreak]\\!\\(Solve\\)\[NoBreak], so some solutions may not be found; use \
Reduce for complete solution information. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/Solve/ifun\\\", ButtonNote -> \
\\\"Solve::ifun\\\"]\\)\"\>"}]], "Message", "MSG",
CellChangeTimes->{3.701767343504314*^9, 3.70176739358664*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"x", "\[Rule]",
RowBox[{
SqrtBox["2"], " ",
RowBox[{"InverseErf", "[",
FractionBox["1", "5"], "]"}]}]}], "}"}], "}"}], "//", "N"}]], "Input",
CellChangeTimes->{{3.701767418400836*^9, 3.701767418864509*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{"{",
RowBox[{"x", "\[Rule]", "0.25334710313579983`"}], "}"}], "}"}]], "Output",
CellChangeTimes->{3.701767419290539*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"u", "=",
RowBox[{
RowBox[{"RandomReal", "[",
RowBox[{"1", ",", "100000"}], "]"}], "//", "N"}]}], ";"}]], "Input",
CellChangeTimes->{{3.701767957445096*^9, 3.701767995501946*^9},
3.701768199994409*^9, {3.701769332327073*^9, 3.701769332926605*^9}}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Histogram", "[",
RowBox[{"u", ",", "100"}], "]"}]], "Input",
CellChangeTimes->{{3.701768093357284*^9, 3.70176817200307*^9}}],
Cell[BoxData[
GraphicsBox[{
{RGBColor[0.987148, 0.8073604000000001, 0.49470040000000004`], EdgeForm[
Opacity[0.]], {},
{RGBColor[0.987148, 0.8073604000000001, 0.49470040000000004`], EdgeForm[
Opacity[0.]], RectangleBox[{0., 0}, {0.01, 977},
RoundingRadius->0], RectangleBox[{0.01, 0}, {0.02, 1011},
RoundingRadius->0], RectangleBox[{0.02, 0}, {0.03, 991},
RoundingRadius->0], RectangleBox[{0.03, 0}, {0.04, 973},
RoundingRadius->0], RectangleBox[{0.04, 0}, {0.05, 1003},
RoundingRadius->0], RectangleBox[{0.05, 0}, {0.06, 1036},
RoundingRadius->0], RectangleBox[{0.06, 0}, {0.07, 972},
RoundingRadius->0], RectangleBox[{0.07, 0}, {0.08, 970},
RoundingRadius->0], RectangleBox[{0.08, 0}, {0.09, 999},
RoundingRadius->0], RectangleBox[{0.09, 0}, {0.1, 975},
RoundingRadius->0], RectangleBox[{0.1, 0}, {0.11, 1005},
RoundingRadius->0], RectangleBox[{0.11, 0}, {0.12, 956},
RoundingRadius->0], RectangleBox[{0.12, 0}, {0.13, 1020},
RoundingRadius->0], RectangleBox[{0.13, 0}, {0.14, 1004},
RoundingRadius->0], RectangleBox[{0.14, 0}, {0.15, 963},
RoundingRadius->0], RectangleBox[{0.15, 0}, {0.16, 1023},
RoundingRadius->0], RectangleBox[{0.16, 0}, {0.17, 977},
RoundingRadius->0], RectangleBox[{0.17, 0}, {0.18, 1005},
RoundingRadius->0], RectangleBox[{0.18, 0}, {0.19, 992},
RoundingRadius->0], RectangleBox[{0.19, 0}, {0.2, 1013},
RoundingRadius->0], RectangleBox[{0.2, 0}, {0.21, 988},
RoundingRadius->0], RectangleBox[{0.21, 0}, {0.22, 934},
RoundingRadius->0], RectangleBox[{0.22, 0}, {0.23, 979},
RoundingRadius->0], RectangleBox[{0.23, 0}, {0.24, 1013},
RoundingRadius->0], RectangleBox[{0.24, 0}, {0.25, 1096},
RoundingRadius->0], RectangleBox[{0.25, 0}, {0.26, 993},
RoundingRadius->0], RectangleBox[{0.26, 0}, {0.27, 995},
RoundingRadius->0], RectangleBox[{0.27, 0}, {0.28, 980},
RoundingRadius->0], RectangleBox[{0.28, 0}, {0.29, 1000},
RoundingRadius->0], RectangleBox[{0.29, 0}, {0.3, 996},
RoundingRadius->0], RectangleBox[{0.3, 0}, {0.31, 1042},
RoundingRadius->0], RectangleBox[{0.31, 0}, {0.32, 1053},
RoundingRadius->0], RectangleBox[{0.32, 0}, {0.33, 1016},
RoundingRadius->0], RectangleBox[{0.33, 0}, {0.34, 980},
RoundingRadius->0], RectangleBox[{0.34, 0}, {0.35, 1018},
RoundingRadius->0], RectangleBox[{0.35, 0}, {0.36, 955},
RoundingRadius->0], RectangleBox[{0.36, 0}, {0.37, 1032},
RoundingRadius->0], RectangleBox[{0.37, 0}, {0.38, 1027},
RoundingRadius->0], RectangleBox[{0.38, 0}, {0.39, 1028},
RoundingRadius->0], RectangleBox[{0.39, 0}, {0.4, 1003},
RoundingRadius->0], RectangleBox[{0.4, 0}, {0.41, 1024},
RoundingRadius->0], RectangleBox[{0.41, 0}, {0.42, 1056},
RoundingRadius->0], RectangleBox[{0.42, 0}, {0.43, 999},
RoundingRadius->0], RectangleBox[{0.43, 0}, {0.44, 971},
RoundingRadius->0], RectangleBox[{0.44, 0}, {0.45, 1054},
RoundingRadius->0], RectangleBox[{0.45, 0}, {0.46, 966},
RoundingRadius->0], RectangleBox[{0.46, 0}, {0.47, 981},
RoundingRadius->0], RectangleBox[{0.47, 0}, {0.48, 960},
RoundingRadius->0], RectangleBox[{0.48, 0}, {0.49, 981},
RoundingRadius->0], RectangleBox[{0.49, 0}, {0.5, 995},
RoundingRadius->0], RectangleBox[{0.5, 0}, {0.51, 989},
RoundingRadius->0], RectangleBox[{0.51, 0}, {0.52, 1037},
RoundingRadius->0], RectangleBox[{0.52, 0}, {0.53, 1047},
RoundingRadius->0], RectangleBox[{0.53, 0}, {0.54, 1015},
RoundingRadius->0], RectangleBox[{0.54, 0}, {0.55, 996},
RoundingRadius->0], RectangleBox[{0.55, 0}, {0.56, 992},
RoundingRadius->0], RectangleBox[{0.56, 0}, {0.57, 955},
RoundingRadius->0], RectangleBox[{0.57, 0}, {0.58, 982},
RoundingRadius->0], RectangleBox[{0.58, 0}, {0.59, 1005},
RoundingRadius->0], RectangleBox[{0.59, 0}, {0.6, 939},
RoundingRadius->0], RectangleBox[{0.6, 0}, {0.61, 975},
RoundingRadius->0], RectangleBox[{0.61, 0}, {0.62, 996},
RoundingRadius->0], RectangleBox[{0.62, 0}, {0.63, 985},
RoundingRadius->0], RectangleBox[{0.63, 0}, {0.64, 1034},
RoundingRadius->0], RectangleBox[{0.64, 0}, {0.65, 998},
RoundingRadius->0], RectangleBox[{0.65, 0}, {0.66, 1025},
RoundingRadius->0], RectangleBox[{0.66, 0}, {0.67, 1048},
RoundingRadius->0], RectangleBox[{0.67, 0}, {0.68, 970},
RoundingRadius->0], RectangleBox[{0.68, 0}, {0.69, 946},
RoundingRadius->0], RectangleBox[{0.69, 0}, {0.7, 999},
RoundingRadius->0], RectangleBox[{0.7, 0}, {0.71, 1026},
RoundingRadius->0], RectangleBox[{0.71, 0}, {0.72, 1052},
RoundingRadius->0], RectangleBox[{0.72, 0}, {0.73, 1019},
RoundingRadius->0], RectangleBox[{0.73, 0}, {0.74, 1000},
RoundingRadius->0], RectangleBox[{0.74, 0}, {0.75, 1030},
RoundingRadius->0], RectangleBox[{0.75, 0}, {0.76, 1006},
RoundingRadius->0], RectangleBox[{0.76, 0}, {0.77, 985},
RoundingRadius->0], RectangleBox[{0.77, 0}, {0.78, 974},
RoundingRadius->0], RectangleBox[{0.78, 0}, {0.79, 988},
RoundingRadius->0], RectangleBox[{0.79, 0}, {0.8, 1035},
RoundingRadius->0], RectangleBox[{0.8, 0}, {0.81, 1026},
RoundingRadius->0], RectangleBox[{0.81, 0}, {0.82, 1033},
RoundingRadius->0], RectangleBox[{0.82, 0}, {0.83, 1019},
RoundingRadius->0], RectangleBox[{0.83, 0}, {0.84, 955},
RoundingRadius->0], RectangleBox[{0.84, 0}, {0.85, 961},
RoundingRadius->0], RectangleBox[{0.85, 0}, {0.86, 998},
RoundingRadius->0], RectangleBox[{0.86, 0}, {0.87, 1022},
RoundingRadius->0], RectangleBox[{0.87, 0}, {0.88, 987},
RoundingRadius->0], RectangleBox[{0.88, 0}, {0.89, 1008},
RoundingRadius->0], RectangleBox[{0.89, 0}, {0.9, 1015},
RoundingRadius->0], RectangleBox[{0.9, 0}, {0.91, 983},
RoundingRadius->0], RectangleBox[{0.91, 0}, {0.92, 1001},
RoundingRadius->0], RectangleBox[{0.92, 0}, {0.93, 937},
RoundingRadius->0], RectangleBox[{0.93, 0}, {0.94, 1009},
RoundingRadius->0], RectangleBox[{0.94, 0}, {0.95, 1049},
RoundingRadius->0], RectangleBox[{0.95, 0}, {0.96, 987},
RoundingRadius->0], RectangleBox[{0.96, 0}, {0.97, 998},
RoundingRadius->0], RectangleBox[{0.97, 0}, {0.98, 1007},
RoundingRadius->0], RectangleBox[{0.98, 0}, {0.99, 954},
RoundingRadius->0], RectangleBox[{0.99, 0}, {1., 1023},
RoundingRadius->
0]}, {}, {}}, {{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{-0.02, 0},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
PlotRange->{{0., 1.}, {All, All}},
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{
3.7017681034740677`*^9, {3.701768137844781*^9, 3.701768202343725*^9},
3.701768896605385*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"u", "[",
RowBox[{"[", "100000", "]"}], "]"}]], "Input",
CellChangeTimes->{{3.701768852262238*^9, 3.7017688687633667`*^9}, {
3.701768901619648*^9, 3.701768902325993*^9}}],
Cell[BoxData["0.9629800337389016`"], "Output",
CellChangeTimes->{{3.701768853393807*^9, 3.701768869865025*^9},
3.701768902726317*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Histogram", "[",
RowBox[{
SqrtBox[
RowBox[{"2", "u"}]], ",", "100"}], "]"}]], "Input",
CellChangeTimes->{{3.701768183226972*^9, 3.7017681909410477`*^9},
3.701768746932829*^9, {3.701768946751513*^9, 3.701768957738257*^9}}],
Cell[BoxData[
GraphicsBox[{
{RGBColor[0.987148, 0.8073604000000001, 0.49470040000000004`], EdgeForm[
Opacity[0.]], {},
{RGBColor[0.987148, 0.8073604000000001, 0.49470040000000004`], EdgeForm[
Opacity[0.]], RectangleBox[{0., 0}, {0.01, 4},
RoundingRadius->0], RectangleBox[{0.01, 0}, {0.02, 16},
RoundingRadius->0], RectangleBox[{0.02, 0}, {0.03, 25},
RoundingRadius->0], RectangleBox[{0.03, 0}, {0.04, 27},
RoundingRadius->0], RectangleBox[{0.04, 0}, {0.05, 42},
RoundingRadius->0], RectangleBox[{0.05, 0}, {0.06, 51},
RoundingRadius->0], RectangleBox[{0.06, 0}, {0.07, 73},
RoundingRadius->0], RectangleBox[{0.07, 0}, {0.08, 66},
RoundingRadius->0], RectangleBox[{0.08, 0}, {0.09, 99},
RoundingRadius->0], RectangleBox[{0.09, 0}, {0.1, 78},
RoundingRadius->0], RectangleBox[{0.1, 0}, {0.11, 92},
RoundingRadius->0], RectangleBox[{0.11, 0}, {0.12, 119},
RoundingRadius->0], RectangleBox[{0.12, 0}, {0.13, 128},
RoundingRadius->0], RectangleBox[{0.13, 0}, {0.14, 130},
RoundingRadius->0], RectangleBox[{0.14, 0}, {0.15, 157},
RoundingRadius->0], RectangleBox[{0.15, 0}, {0.16, 172},
RoundingRadius->0], RectangleBox[{0.16, 0}, {0.17, 157},
RoundingRadius->0], RectangleBox[{0.17, 0}, {0.18, 165},
RoundingRadius->0], RectangleBox[{0.18, 0}, {0.19, 203},
RoundingRadius->0], RectangleBox[{0.19, 0}, {0.2, 184},
RoundingRadius->0], RectangleBox[{0.2, 0}, {0.21, 197},
RoundingRadius->0], RectangleBox[{0.21, 0}, {0.22, 208},
RoundingRadius->0], RectangleBox[{0.22, 0}, {0.23, 221},
RoundingRadius->0], RectangleBox[{0.23, 0}, {0.24, 258},
RoundingRadius->0], RectangleBox[{0.24, 0}, {0.25, 211},
RoundingRadius->0], RectangleBox[{0.25, 0}, {0.26, 255},
RoundingRadius->0], RectangleBox[{0.26, 0}, {0.27, 244},
RoundingRadius->0], RectangleBox[{0.27, 0}, {0.28, 279},
RoundingRadius->0], RectangleBox[{0.28, 0}, {0.29, 293},
RoundingRadius->0], RectangleBox[{0.29, 0}, {0.3, 290},
RoundingRadius->0], RectangleBox[{0.3, 0}, {0.31, 328},
RoundingRadius->0], RectangleBox[{0.31, 0}, {0.32, 289},
RoundingRadius->0], RectangleBox[{0.32, 0}, {0.33, 380},
RoundingRadius->0], RectangleBox[{0.33, 0}, {0.34, 339},
RoundingRadius->0], RectangleBox[{0.34, 0}, {0.35, 339},
RoundingRadius->0], RectangleBox[{0.35, 0}, {0.36, 330},
RoundingRadius->0], RectangleBox[{0.36, 0}, {0.37, 378},
RoundingRadius->0], RectangleBox[{0.37, 0}, {0.38, 343},
RoundingRadius->0], RectangleBox[{0.38, 0}, {0.39, 375},
RoundingRadius->0], RectangleBox[{0.39, 0}, {0.4, 388},
RoundingRadius->0], RectangleBox[{0.4, 0}, {0.41, 416},
RoundingRadius->0], RectangleBox[{0.41, 0}, {0.42, 396},
RoundingRadius->0], RectangleBox[{0.42, 0}, {0.43, 435},
RoundingRadius->0], RectangleBox[{0.43, 0}, {0.44, 421},
RoundingRadius->0], RectangleBox[{0.44, 0}, {0.45, 418},
RoundingRadius->0], RectangleBox[{0.45, 0}, {0.46, 482},
RoundingRadius->0], RectangleBox[{0.46, 0}, {0.47, 467},
RoundingRadius->0], RectangleBox[{0.47, 0}, {0.48, 429},
RoundingRadius->0], RectangleBox[{0.48, 0}, {0.49, 475},
RoundingRadius->0], RectangleBox[{0.49, 0}, {0.5, 486},
RoundingRadius->0], RectangleBox[{0.5, 0}, {0.51, 535},
RoundingRadius->0], RectangleBox[{0.51, 0}, {0.52, 549},
RoundingRadius->0], RectangleBox[{0.52, 0}, {0.53, 482},
RoundingRadius->0], RectangleBox[{0.53, 0}, {0.54, 531},
RoundingRadius->0], RectangleBox[{0.54, 0}, {0.55, 519},
RoundingRadius->0], RectangleBox[{0.55, 0}, {0.56, 591},
RoundingRadius->0], RectangleBox[{0.56, 0}, {0.57, 528},
RoundingRadius->0], RectangleBox[{0.57, 0}, {0.58, 588},
RoundingRadius->0], RectangleBox[{0.58, 0}, {0.59, 580},
RoundingRadius->0], RectangleBox[{0.59, 0}, {0.6, 599},
RoundingRadius->0], RectangleBox[{0.6, 0}, {0.61, 597},
RoundingRadius->0], RectangleBox[{0.61, 0}, {0.62, 618},
RoundingRadius->0], RectangleBox[{0.62, 0}, {0.63, 616},
RoundingRadius->0], RectangleBox[{0.63, 0}, {0.64, 665},
RoundingRadius->0], RectangleBox[{0.64, 0}, {0.65, 618},
RoundingRadius->0], RectangleBox[{0.65, 0}, {0.66, 618},
RoundingRadius->0], RectangleBox[{0.66, 0}, {0.67, 620},
RoundingRadius->0], RectangleBox[{0.67, 0}, {0.68, 679},
RoundingRadius->0], RectangleBox[{0.68, 0}, {0.69, 690},
RoundingRadius->0], RectangleBox[{0.69, 0}, {0.7, 767},
RoundingRadius->0], RectangleBox[{0.7, 0}, {0.71, 735},
RoundingRadius->0], RectangleBox[{0.71, 0}, {0.72, 702},
RoundingRadius->0], RectangleBox[{0.72, 0}, {0.73, 722},
RoundingRadius->0], RectangleBox[{0.73, 0}, {0.74, 727},
RoundingRadius->0], RectangleBox[{0.74, 0}, {0.75, 735},
RoundingRadius->0], RectangleBox[{0.75, 0}, {0.76, 751},
RoundingRadius->0], RectangleBox[{0.76, 0}, {0.77, 778},
RoundingRadius->0], RectangleBox[{0.77, 0}, {0.78, 776},
RoundingRadius->0], RectangleBox[{0.78, 0}, {0.79, 818},
RoundingRadius->0], RectangleBox[{0.79, 0}, {0.8, 842},
RoundingRadius->0], RectangleBox[{0.8, 0}, {0.81, 826},
RoundingRadius->0], RectangleBox[{0.81, 0}, {0.82, 819},
RoundingRadius->0], RectangleBox[{0.82, 0}, {0.83, 792},
RoundingRadius->0], RectangleBox[{0.83, 0}, {0.84, 843},
RoundingRadius->0], RectangleBox[{0.84, 0}, {0.85, 826},
RoundingRadius->0], RectangleBox[{0.85, 0}, {0.86, 878},
RoundingRadius->0], RectangleBox[{0.86, 0}, {0.87, 866},
RoundingRadius->0], RectangleBox[{0.87, 0}, {0.88, 935},
RoundingRadius->0], RectangleBox[{0.88, 0}, {0.89, 845},
RoundingRadius->0], RectangleBox[{0.89, 0}, {0.9, 948},
RoundingRadius->0], RectangleBox[{0.9, 0}, {0.91, 927},
RoundingRadius->0], RectangleBox[{0.91, 0}, {0.92, 954},
RoundingRadius->0], RectangleBox[{0.92, 0}, {0.93, 925},
RoundingRadius->0], RectangleBox[{0.93, 0}, {0.94, 927},
RoundingRadius->0], RectangleBox[{0.94, 0}, {0.95, 967},
RoundingRadius->0], RectangleBox[{0.95, 0}, {0.96, 936},
RoundingRadius->0], RectangleBox[{0.96, 0}, {0.97, 942},
RoundingRadius->0], RectangleBox[{0.97, 0}, {0.98, 928},
RoundingRadius->0], RectangleBox[{0.98, 0}, {0.99, 970},
RoundingRadius->0], RectangleBox[{0.99, 0}, {1., 992},
RoundingRadius->0], RectangleBox[{1., 0}, {1.01, 998},
RoundingRadius->0], RectangleBox[{1.01, 0}, {1.02, 1042},
RoundingRadius->0], RectangleBox[{1.02, 0}, {1.03, 1074},
RoundingRadius->0], RectangleBox[{1.03, 0}, {1.04, 1059},
RoundingRadius->0], RectangleBox[{1.04, 0}, {1.05, 1035},
RoundingRadius->0], RectangleBox[{1.05, 0}, {1.06, 1041},
RoundingRadius->0], RectangleBox[{1.06, 0}, {1.07, 1029},
RoundingRadius->0], RectangleBox[{1.07, 0}, {1.08, 1028},
RoundingRadius->0], RectangleBox[{1.08, 0}, {1.09, 1125},
RoundingRadius->0], RectangleBox[{1.09, 0}, {1.1, 1001},
RoundingRadius->0], RectangleBox[{1.1, 0}, {1.11, 1097},
RoundingRadius->0], RectangleBox[{1.11, 0}, {1.12, 1104},
RoundingRadius->0], RectangleBox[{1.12, 0}, {1.13, 1149},
RoundingRadius->0], RectangleBox[{1.13, 0}, {1.14, 1144},
RoundingRadius->0], RectangleBox[{1.14, 0}, {1.15, 1163},
RoundingRadius->0], RectangleBox[{1.15, 0}, {1.16, 1228},
RoundingRadius->0], RectangleBox[{1.16, 0}, {1.17, 1105},
RoundingRadius->0], RectangleBox[{1.17, 0}, {1.18, 1134},
RoundingRadius->0], RectangleBox[{1.18, 0}, {1.19, 1229},
RoundingRadius->0], RectangleBox[{1.19, 0}, {1.2, 1226},
RoundingRadius->0], RectangleBox[{1.2, 0}, {1.21, 1239},
RoundingRadius->0], RectangleBox[{1.21, 0}, {1.22, 1213},
RoundingRadius->0], RectangleBox[{1.22, 0}, {1.23, 1231},
RoundingRadius->0], RectangleBox[{1.23, 0}, {1.24, 1241},
RoundingRadius->0], RectangleBox[{1.24, 0}, {1.25, 1205},
RoundingRadius->0], RectangleBox[{1.25, 0}, {1.26, 1239},
RoundingRadius->0], RectangleBox[{1.26, 0}, {1.27, 1344},
RoundingRadius->0], RectangleBox[{1.27, 0}, {1.28, 1295},
RoundingRadius->0], RectangleBox[{1.28, 0}, {1.29, 1281},
RoundingRadius->0], RectangleBox[{1.29, 0}, {1.3, 1246},
RoundingRadius->0], RectangleBox[{1.3, 0}, {1.31, 1285},
RoundingRadius->0], RectangleBox[{1.31, 0}, {1.32, 1344},
RoundingRadius->0], RectangleBox[{1.32, 0}, {1.33, 1328},
RoundingRadius->0], RectangleBox[{1.33, 0}, {1.34, 1349},
RoundingRadius->0], RectangleBox[{1.34, 0}, {1.35, 1301},
RoundingRadius->0], RectangleBox[{1.35, 0}, {1.36, 1358},
RoundingRadius->0], RectangleBox[{1.36, 0}, {1.37, 1326},
RoundingRadius->0], RectangleBox[{1.37, 0}, {1.38, 1419},
RoundingRadius->0], RectangleBox[{1.38, 0}, {1.39, 1353},
RoundingRadius->0], RectangleBox[{1.39, 0}, {1.4, 1435},
RoundingRadius->0], RectangleBox[{1.4, 0}, {1.41, 1356},
RoundingRadius->0], RectangleBox[{1.41, 0}, {1.42, 621},
RoundingRadius->
0]}, {}, {}}, {{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, \
{}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->{True, True},
AxesLabel->{None, None},
AxesOrigin->{-0.028399999999999998`, 0},
FrameLabel->{{None, None}, {None, None}},
FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
GridLines->{None, None},
GridLinesStyle->Directive[
GrayLevel[0.5, 0.4]],
PlotRange->{{0., 1.42}, {All, All}},
PlotRangePadding->{{
Scaled[0.02],
Scaled[0.02]}, {
Scaled[0.02],
Scaled[0.05]}},
Ticks->{Automatic, Automatic}]], "Output",
CellChangeTimes->{{3.701768191767829*^9, 3.701768204867261*^9},
3.7017688291979113`*^9, 3.701768905660675*^9, {3.7017689473898706`*^9,
3.701768958499846*^9}},ImageCache->GraphicsData["CompressedBitmap", "\<\
eJzFfXdwnMl1J5ZZG6SVLcnlKpe0ts91Z8lXdXb5SvfP1bku+K7qSv9eKFln
r2WvdsVlQBykQRjkTASSIIjACBIECBAgQIAgkUnkMMgAkYmccxhgZn7Xr7u/
b2aIATEgqbupwgDT0/3e65dfd3+NX3zp9vU/n/vS7Ztff/nFf3X58uLX3/za
9Yu/veDCmo5/5OT00ddOTk6DXzjR32B/yjf++obe5IdD/v5v/NdxGI0myJfT
Tx1sOynbzCZL258o/Uz7+xmNxn3waKxZtn0pSPpb/uu72FyaQFpiJAJ1Oly9
+RBru6zr1iyuJl3BwOQkslOvIjwiEhHh4QgLDUZYTCpmVgw/4cPPYG1uBNfi
wvnw1MwCbBnN2F0dR1LCFYwubh5TiFqfREraPSxvGp2+4G2fwbC5hPsZl+Ef
EIj4azcxvWYAjFu4dy0eNT3TgHkXL4qzEaQLRHBwBF52DKljt1dmcOPaJQQE
BuByxn0sbe3BvLOCjMvxaBtbxt7WCrIzktlYHULDotHQ81oy7RNUPMpBRUOP
5JEBD+9cR1ljH/vbjObKAj6VkLAI1HWPyTGfYqCxFJnFL8UQ8xYyM66hoXsQ
5fn3ERYWjsiICMadEOhC4tA1vPDHfNgpdL0oRM7zZkUeqHx4A/dLm45Lrrx8
nImKtleyMxfLf+G/vsfeTHh4NRxBcTfQ/6oLkb5uuFFYzxi0htKiIsysrqH4
Zix++U9n8eh5DZob61HfqMfmjtHpxwoM8x7uxukQfvke+vvaEaxxQVZ5O2PU
EooLn2BudR3Li4tYmJ/HeF8tXNwC0DM8C5PZbIHBmFKRlwZNYBx6+vpwLdIf
cTcfsz5GVJc8Qd/UOoabn8LF2Qv17b1MYPfh6q7F+Mq2FQwT8m/Ewy8qBX19
PYgL9MT13CqmlQaUsbm8XtxAadZVaHSX0DMwgMrHd+Hi4YfxVYPTH0sBZCZF
4/4TKQDsoaGyFJ0jU1ga74DzBTeUN3WimsZ5BmNqbVcZ1lKUgdCUfCk3A16U
lWBgah76sgf41a/+AekPS9Da3Ija2iYsrGwpw2pyriL2Vqkqt87aSrzQD2Jz
fRXzC/OYn59DRqQfrj8qYxzftZHdaZi3ZxHo5YOW4QU+uqMqD8FRyUzV4CQM
4gRT/kWkRwdAGxaPscWtLySr9tbGofXwQc/0Oh/aUHIPEYk3oJi0YWkYAVoN
tMFhCPb3wtdf/Qa6sBQsb+xKdp9m3NlCcogXCuuG+ZjpriomvlisGcwS/Snk
3YjGpVsPJV82Ead1Re3wgiIx0/Yqov01qOma4T2G6wrgxdBsq+5mD9diApH9
VC/pWuDa1T69JiV2HA+S45BdUic/CbQN+WkIvfxAjNldQWSAN6q7ppQhrcW3
EJn22GrICabDOyi8kQg3Lz+0Dc38qWI6uSmIv/vcBnpxehLcPTSIiAjF+XNn
cc5Fg+dtg29Ix7g2Cl/fUAxJFo+3ViIkirHHZEFq3FlDTmosvmXa1DOx8hPJ
lu3FV/DyCcfk8g4fOlhbjLBLV7AtmWLaXcXVMB8k3n+O2eEWaDyD0DdCBgXV
GMyGNVwO1ODFkFCO9ZlOaHyjsLChyNiM/NsxcNfFYnZ5Da976+FykbF2ck2F
YdxYQLSfBvpJMYWFkXpo/BOxbjCrMGqYMXj5RqC2qQW5N5mX1cVgbnP3reKp
zLyMq9lPJYxtXI4Lw9O6AXviOaG8M6Oqyr+Fby64o6pjTAlM9sTTUZnHXEMA
hqZmkRkXhPSCCmxuv2k8xrUxJp4QDE6vcSpet1QgOEKIR0VqNCA7JRza6OvY
3IOi+CQeT+8wTCwJ8Qy8fILQOC4eqfgnmb/ohr+vHxr0DfALiMHs6p6TIl4w
0SQFeuDF4DwfvzbVAQ+fSFvR3E3ABTdvFpCC4O7mjlt5z2FUxXuaiybKz4OJ
RpC/MFQLD78EK+M7ic2pPvh++w38g0Nw4ZvfIj23Rga2E9hcXcLaFgtKKQkk
nhNWkrnyoESSsYWkmBCSzE/kKMP2OpaW19Hx7C5J6LTk+4uCG3D2CsLU6q6T
0nV7cw3LK5toKEgjCZ2QNDHdRdaVKKTkPEdBSjSyKlucTltLhmnu7hKCvTxR
3TvFqWh+eg9BLBswQHmZ8Sz7Jp42DvBPKwONcL7ozdy7wcKerRn4MRNtHF3i
fWoepSM86Q75N9X49BV5cHa+iPOaCCYegyod0w5SI33woLybDx1tLoUmIBbr
uxalL81KQWxaHhbmpvF6bBgtrXps7BpV9ObddRZ9NHjaPMpHdJdnwzv8OnZM
UF1jc34G83Wp/Pvh2kdw8b2EjT0o/BusL4K/fwB8PD2R+7xREVBbyS34xdyQ
8plHsLcGL3tnlFEL453w1/ojLMgLMTefqKPK85FZVCtGrU1De/ECanvmlFFT
A83QslGhjOSkrIqTigiXJxCl88KFb3+L7Oo2GzGRNZlQejcJXrpoVFeXwdvd
BXnVvYoVvixIZZwNRHVtA25ejoB3RDIp54/V0XvIT49hcSkB1VXP4OHigpLG
YYu5MzW5Gowvv9ayFGT3J5KvxPum0ge46KFDZXUVwrxdkZFXTTmnytfFkXZo
XJ2R96wahZlXcd49GItbRgtiM8pz0uCmjWRkVyDA7QJyylsVsoebSpgyeaGs
5gXSmPPQxWeCZaoKp/Z2NlH3LBe/+fJLZD9TpbLCLM2TOc/7T8rwIC0O3oGX
sGwwK6PMTOl7m6vh8fU/ISy1UBk11FyMsxdcUVLxEgX3UnHRS4expR1llMm4
g/baZ3D+53/EpczyU9LYRtur8M9f/hL3KlsVy/lbaVy7W6t4kn0bUdGxyC+r
h0ENNSeZK9tEVXEuYiKjkXzjPl7Pr0lncBKGjSVGQAaiYuLwpLoFe2ZLmr+9
OIbbN7OxtLmn9jfvbaOqhMGKZslS/jOm+ibVZ472NCMpLhbxSdfRPTL9hYqe
5V4siY2OjkFuSTV2jGaVNPIJzZVPEMvApd7Kwczyplp5tFUVo6a1X1EL/j7S
04TL8bG4nHKTRddlpet4dz0elTbKTwroPcbD54iPiUHC1TR0Dk+roAfbalBY
qVd6ylFmVBfnoLp9UGHuH8qv1tbWsLGxwT9xO/qF+P77kipiYlNtDV42tGBl
c/eEHNba2oqenp59wz6XSmBi1URKZCC0ulCWafkg5FIG1UZK6q7X69Hb23vg
8FF9FTw9dZhlwce8MQGtuzvqeieU4V1dXRgYGDhweHNZFkISZKLEXndiA5H3
slMpp2jODx8+RGFhITlMvZ7eyzKYF4pDa0k8Wp8mnmZ/JB5D29PE76h/fWyn
zdFv3xEKsao0kQjin8442PY+Yz80vA9GS4k1vJPWbVxobeU3qVWv/4iU8yMM
DQ19ZKsVTeVZCE14qGpFTryOBQC9ohVkBY8fP0ZpaSnTh76+PvZewD4/fV6B
vvYK9ORfPIWeAjf6ssgduVfOoyT9Av192qr54eXzeCqaz1g157DmZxkX0Vdo
0zuHASm7ua85K+kcym9Rs5s1kKzEc6i47UzN1r3vJ5xDJWvutW2m3jV39zVn
xsvmxzaw77LmF5ku1Gzd+86lb/Hy3r7mzPhvUX/fBT22QKi3bLbunRp5Fg1Z
rm82X2fNjQ/2NadHf4umbFcL7OOid8RZNOe4Otm2/RYth7UVMgKizqLtoRvT
HUu/lPDfQp9rO5ba2vP2t3WwNnVsoRufUOcjNyscbrjG+nU+soxlHJNtln49
sq0rX6HlFG8jeN35btTC1MwZPc+iT5IGfkRO8iPcvn37DVVurXqIoNjbclXK
iPRwP3Jwin/s7OzEq1evbGJGR0cHBoZY4mZgP60aQK+l+NDuTwGoyx9DJV6Y
qvBm1br/aavmgWK1+Yx18xMvzFT6sOrc0rs7AI05Hhgtpd4B1kD6nnhivpr3
tgbSV+SJhRpfGyCsubfQE4svfN+EXXvfHdMVPm8S2PPYE8svfd+E3c2aV2v3
we4q0GC1bh/sqkw3LL3YB6ST9d6o3wek5aEGmw3aN5ubWfNWo/ZN2MUZrliv
11pgH+e9m3I02G7SOtm2eWCn2aqtU7QZrNu6AlCU5sLG+jEdsuqX7YG9Fj/b
sazN2Gq/zXpsYaqLbb8OfzSyfuY2e23++9qg95fwTvE2gsfaqIX9ZiIbyuB5
3kdobGx8Q5Wn+hvg6qJB3wyrrMbb4e7qifaRRXuZgmj7DC0tLaytj9c05o1x
mGdewtikgbHFF8Zmb/bjxaajpSnRFGBqo+nSBI8rnyUL9red5u8n+DReZnlg
kOk/m8qZtw3XW9rO2GmzBll11x3jz73pe3td3wVk2S03zFb50Hc2w1tth+OQ
NmuQBddduU1ZgaSuRjvD39Z2WmljavAoxYWswOnj94Bn3ZZ3zQU7zHzehz7R
dopzleAx8xHftXjD9CrtjK2ikrYV3roMD40Xy001SMl6hj0TDlbUTygHwcTE
uJpxwGQQCrs5AdNYAUyDt2FsDyVldfpjyX8zY1Yvc4zkdJmt/blVc1eBcKOM
Xuvm9kcaRVrWzeSLye+ozaQX/mhg9sq8lNPPrDS6/oHwMT+1Ygi12WeIUhpQ
4JkYG8Ho+DQrmnHCDiM+lW2m3svMQpnjafGRE7U1pZ/a0fG3tVkr7k/tSPZn
dtpEv5PWbcrEnM7YFkU0BaW++FTGT/NEMZeXqe8qczCefCogIH8lcUlYvPsb
RVZ7ezvy8vIwODjo9F3JPfMW49rATQbvGndcBOUvHQP19OlTplYTqmaaVwek
z6MZaVX/Rw3/+nCQNFslW/hcdjAv98DY6CbcansYTEOZqqZyuH/hGFxFEQTc
j5nmjzMu3uI/5pkaYRS9ycRNSb4vnwJ9+POjofih5Iap96oAx0CZV3pgnq3j
HCFo/8oxkExOViBP24CEYVkwnKuzr+DGv3RMcLm5uaRXEu53YOq5IoTFgpZ5
a4aBYgFtc5L9PQVTd7wF/p85Bp+VLpiampLwmUFuzzN4U/wHpl3FVanTYTA5
+H9xNCX5IwVOd4Jgg6IUDsJRJPZHiilMlEgFC1H19k+PBurH1o6GUgBG0jEV
2BfvAox5/J0Fbqbqz+wLVZN+cjRNUlZnSCG591BY/2PHJEvLFYz9Es4ZmJc6
mH262qQ6Iv3x4qKlhj90DPSzZ88wPT2tbvCaN14zadyFaThLSIQ0pSuGf+bO
ilyCQCum8KOjac+fqcrPFLzhIkwjWdwNkJLa4FXg//Bo0vszhdVcq+7APFfn
pMQmY0+Sreb/4GDYYvHsc8xPDKKhsQH1dXWora3Fi5cvMT6/wmke7GpGSXEJ
2nqGlBMC2FqZRXXZUzyvfIH51W17YfEvpOIbexItnk/q6u8fRtH30fL0Dv7X
//if0OqCERsTjZCQUNT2j6O/8Smcz7vhaspVuJx3RmnLIMyGVSSFaKELj0FU
iD8CIq9hZcfoZC9t+Qup95wsJmUKd7Qba6333z+cvN2NBVyO0MI/6joW1nYl
WwxIDvbF/Wet/FND8R3oYtLR3/4CXt4hWKRu2zPw17jjZfe4DXmKDdmS58EU
J0coznIXI09zGHlipflTGNZmkRwTBI1vBMYW14G9BegCgtAm961m+hsRHh6N
pyXZCE/KUZaH7sYGIvdFB/+kbEMrPr2+vh79/f0/480fWxGYLebO3DqIjzPV
IqhaxefPDib4tGKQe+vITAiHpy4eC8szCNGFQD8giJ1/1YKIiEgUM2IjEh46
KdTKta191JaXl2N4eBgzMzM/k12xs8iUdhbYXVUXC3iYHr7HZ2Bji58eLv6V
+Ql0941Iqc/Dz8sDz+obEBUchLo+sfk10fUCQaHRKC97hODYmxAb40akhWnx
qLbLRvyK3/i3/P3MG0ZDP/7A+jDT0peqlp46nMyhtmc4d06DwblVrC+NwMv1
Airb+5GTHIHo9Dxsbm0gOyUS0WmFmB1th5uLOzrGFzE31MTqVS90ji3bNSKF
TOwsMaHPwdgdJ8llZDJJmpe7eQRwkEzTzioepCXBk1UdXp4eSEi5izWDiZXK
vYjwZ5WIlxe8dTHon1rhHHyWlQoPDw2rUjxx81E13121R+bPVWNKsJC3u8oC
S6cj5J2SOkXbZfMz05iamSVUivrtbK1hcnIS69u7qkrSJsT87DRm5hfVjQlr
8/65aj32KbIJcdJ6jjkW4rKzs/mC7s+lRe1DsdjODNNFhLbOaJhGc2zV/iPH
8JB1zc7OSjwnbfFQ3rgxztjBAhwlPQstQukpiZL5vwNoiGPM1Vhx7ASda+Gp
HuUrMBuViEMMN1pnkQ6CV1TkP3xQOCcknARRohBHjgTlbyRPiY88MWPVgyx6
KNNr9iJ+78+EnLaOlqj9R6ki5qV2Cd5feMLRXEs21H2Jfzb1p3KmbDqmHQ8e
PKDZSBTHRVggke0sqtu8nM8UOl6LU1fmxQ5SfgcxVFZWYm5uTu6csg5mMyeW
QxzLExAXWgmiI0Rb69oHBKlI1AZkZ5RIBkdzPhRI02Qpi18PVENTQKperS2A
S5c+LRwNzX+XGi1BCp1r8z9GMGlFiukIj/0sBPCYr9TmvNTz5G5l5mhqaYOR
vBQjHusjx8j2mc/kMXxvS6pMi8Tiz+dPUjLqgwnrtGNqlJWVxRXVLtLdlZNy
nHmhWTYHqsiNXbHcJBzEVFVVhYWFBfuYDEt2MRko0Jl5SeIgJms9/t1gOqMM
p5PCZqUSIGe8y3ebjSZLFN7eWGfJxY6k4y3KNc9IaXDm6GBYFBSe4urMFUkf
ZKGwO+4t+iSI+y5K71xFSVMfJ2x1fhQJof7w8vFCYHgSxpfoCIAJL55kwYel
GR6evsgtraeTRfayB5XG9TGmYw+ZuypgKrApaPyYr3PyoES+2ixPPbI4S96O
ryKJpcCJw9Ofpdkpvoo6NDiAK8E+SH7wBMsbm8hPiUJIYibmFudw65IOsTeL
sDTRB81FV9R2j2JIXwlnZw16J22zNMWefsnfzwgnMfZIceVvmunI4RQ+uhyJ
8wxrdHQkLpw7iwsevqhpbUWUToeX8qzXWEcVQsIvoaaiAEHRGTLd3UNymB+e
NPQ5KWe3lChRU1NDpyv+TlXKNl7gcItmDtI0lq9YtNPwYRR+xstCjacO4wss
lYwPxr2yFuwZFhAUGIL2Qds64klxNiKTctWkLfuAOoIonJiYwOLiokrlfJPF
dIziUKexU5iOA1SS9hXcSET8zTxkXw5Hdg2rtUzLCOKl2RwHN9vXhLAwWe0k
ZqtU3mO1WXZNuxWVtkb/K4XEuQZpUj78AC3pJuV7RKKaSrZq1dyh/7AE+HOm
9AtICPHGV199hbzaLi7WBJ0PHlZ3cpL15dnQRV9HZ2sFvPwisEqy31tEqLcn
KjtG7RnXryRHzOujIs14XSgYaqFWBC7GaPIC5OCJOprFCetZyIXv9sNKzO9h
5lUTzn7198gSooa+LBuubt7Iyc2Bxs0Vj170wbS1iBg/d1xKuYX0K1HQBF3G
0qb9lQTxlMR3yAEKkilbGH0oHBY5AnJb3JF5cmJFPJZTYjG04zBv9hn0tWVo
G5gQvsVkRNvL50hPz0BFQyedGOav1dkxPMy8jbvZ+ZhY3LDrbr+UWs1J7Yrj
sjeNPzkuZ2Ceb7CEeOkceN7DZmKeruSJ51sYbC+W/6Nsw+6aYIhx2+lgdJ5E
jWDYHP/OEXRkovfu3UN3d7fT128FXSRB1zsKmmZC1r+0tPQW0Jp3Ba2YrF3Q
PGPzhWmi2Aq0xlHQisQF6I8laI1gxFiuUEzpuMzzjRwVQ3najknVvgtGSjqa
BFQxCeV0Jg+XVIZOlYsyVAQf1ZbbdHLq74yYziLy6YkpCp/7ifAwLJhQ+CZP
Ypki+63X8bRSUXIHMStK7sbfz4jEXgmsHGwQNzWpzA5AVXS5q6vLyVVqBc3H
2Bm5HzJrVxTOQYJfvnyJ5eXlI4E+Zkclqo6m3wLdp46gs+i+Pui4HcxlR9MJ
VykZgTnCZgFaxTz78hipIu19kbPbXYexI1xW7n5CV9oCncrfAzMti5G3t8bM
jI+p5DHSR9rjWKFhLLvtiBDs4RnvMtPICtJIB6dtq5DM2PY2eG5qbA8XUyc9
39s4oYxntmhBGMyDrnmujovdAYykrJmZmVxZ3XnbSUtkoekyLjJsTgeiM+0x
dLWOoqMJ0p7EysqKk4di63ubfKGRTxTmD4ipr48XD05+Hw6kohqOgLRndSWH
5TbHsb66gqXlVavT9Ky7YQuLC4vY2bXkL2vLS1hZXVe4eACJH7+FRI0gjVtq
cOnhpQM9DlBTdA+eGg+WZbkj+c4jbLOKa2VmEFEB3vBkOaKvLgbD8/Q0jBHl
ebehcXeDm7sn7hVWs0zT/tqun+pbDMJU+R6tn7BsCnpiuULEUaYlxo4wovwt
vFQonn7VCLcLLmh+NYXFqVfwd72Iqq4hXo5FJudgdX0F9y6HITq9AAtjPawc
c0fL4DQme+tYOeaBrvElu+WYn8J1ZiomFgVN48WsNHt2jEg+RjTyv4g2Zq0y
ZvNVP8+iw4luq7wPr5A4Yhd/ZSeGIKv0GaKCg1HXN83bXnfV8AqtqjzfZkMi
JdwPRQ29NhUaBSNaaevs7AyVcyEdME2VEUdpgXeihLVRAsRDN6Pv8eH1T3vN
Q/z6m/OoaO7FUE8D3C46o6y+HuEh4fuqNKp/IhNz37rbQ+ylfan5+Xmsrq5G
KZTSYls7yVuUPFoeUSjV0D44nMi26hy4eAbiUkw4NBpvZBZUYJclCcGBQWiV
RdpMXyMv0kpo/yzxgUpkZkwAcuwUaYpXiVYoNO6IulEfTBTqQ6ioEXyUwSdY
baYiWKnWW3wylQlQU5MbTH3JyskeRRdG2spZGRaJ2TUDDDurKC8txfjsFK6E
+iGroo1PoLmUnsxKRw89hOAbiiVSHMMcgujpra7Xdk1OJZ7cjHGLnDBNkr2f
EGldkydfPmaRRzSf5Om/aeo5aQuM/BCQj9MdPoMT3Lfw/faJJ5SKvjEJ8/YK
UmMDEBxzBbdTE3Huog/zEZvofVkAV1cN7mTegburG4rqB2HeWWb1qCdTl+u4
eikU3mEpWNky2Z1EvFQd88wLefAnlKShbNaQgitnF26rrPbmBwYgcieolfH3
sLu5gudFj3Dvfi4GJuak8ZnQ3fyChch7qNP3q3vmG0tTKMrLxsPHTzG7um3X
/1qoq7GiblulbrLUPnVrQ9bU2UsNPixkchF3794l9/A7oJlMem1tzQKZ0nNa
0qB1zPeCrNihCnl9hJee5I+Zf3sPyIoEUxQQLO7wOodp/DuBOAnz6isBpj+F
+ytaCrt5NKmnKtTMVImY3R7iCAgSb1paGpqbmyUI8TwWrxxZdqrsG6qZuvBN
/NM1xyhsbGzE+vq6U7rSgYFWzkZZ58u0t2B6lc756CBoRcgC9MfcV+0DzdIB
bqrMOXGlbdUmH01Kh0M3w0RZPGXE7aE8Dl15FxSfvoFCK8rrjlCxuWa2Egst
4reH8Owx8V1QnbZFRWkVrWSzMoahclKfzCO0hI4Wveh7vc7JQXSKXgp0pyzo
mAio5OGgTTvcvR6TSE3jhUreyez2SLqbLs3cBs3Mi48Vq6CgRckWm4ADkGkC
TU1N9BDVG5BZhG70IFX9WGnmUtkR+VGz56V3UdwTwiZos4mpEKuT1YPGFEyl
d4o7mpgzZBtfjeHaY7ScXyZaG50taWfMwbCV+DfW04obqddxPTUDzd0j3KTM
ph3UlOTiWvI15Je+oEdZ+Wt2tAe3067h+o1M9L+etxv/cqUakl8my6fzLvxH
MSSmA1EHk6UkD3NDrXB3dsGN7AIUPbwNZxd3tI0uoOXpPbi4B+JJ6RP4uTnj
flkby1BmEeLphuRb2biflgh33wjMrBvsJg8FKvt21c0JdemMPLQ8Mx14OInN
ZZnwC82A8robE4gH5TW4HBqIwnq+foKumkfQRSajvbEM2sA4bHL2riJc64lK
/bDdKkMlkWoeSaLYDT4p9myHMjlHTa9SOakBh5H6GQb1ZXBz9UfH0AQmhtoR
6OGO8rpahLGsXS/vqJh71cyy9igUFz9AhNXeStYl3RsJsagvMjIy6NFUstL8
t1HMGErLlPSSau8IxbRTmx7ig6/PXsCFs1/DNyQV29sLCAm0OlXWz+qM8Eg8
KWEUJ1pOlWVfsl9nMEppCQKbm5v5in0SxW1Bcl1LrmWeEr6TlaNUMHPChRM4
lO7vY228B1pW+DYPzGGB6bCnhz/6R18jNsgflZ08A+XZvS40Di+qixAYkSIr
PgOuBPlYVXG2DqVAJXhTOBRa59IHicWG74glL1np85i70GxNOO91jHppHZjC
JD1PfwFlLQMYbC3HhYveGJpdQ8ntS/CPuY7XE69xJUKLy/crsTr7Cp7O1LcP
XXVPcN7VD0MzG3ZNT9wp8rlY5WXey0Kbv1ipYE6Y6n7u2eR27L4+HeGHzoAO
8RvRUl4If18vePsEoLCihefwW8vTyEiMgsbTEzHJdyiH51jaq4rg78P6+geh
tL7XnlvjuwR0bJ+VOvxwJCtO+RKWIFbM97Sgl9YxOiNFTN9ZEvGJWYTPwYQr
i0B7u6ze292z0lxRq+3s7PCFoeNW/Xb3jIpOHEDtCVELEfqeJEtyfJof2zHW
XxClHnF6d0Ol0mZrjjjP7NXtaKlBkUK52SiFvSvDjRIWauRxFzMvLxXUxxXr
I2VmWYl5tlYUDnxXQeviWOaQmppKMd6pWqHCLgZfXrzzF4VoXrg7goHmSQ/2
bW1tSQyfWmOwjiESS7Rk8doRsSiWb4OFOEH7QMP3+bK7qT9NLCO8MxZFYyxY
1lm2GMsdCZ3qEkZYLJX6fbHwJ2mcPhGVCGmGYuavCw/C4OQgCkX56lQdfy+Q
1pr0gUAqqrO9vf0hQSp68gFBKuKqk0rxFpC2Ci+KpxPUcNoKlRKB/v5d0H/G
T74IdxVqi542uUeyeXrBTxaQJbC6jVkG856slfYOp6uUXUmLS1Of4vLmJxf4
skeb/7uR9ynfnqS9F/OiXrq0VUEelSryVJ8aF7pilacueSSnhJi8BFFNe7l9
ySK6NDjDvKQXWWpbgKOUKQ8DNikddlcE+PqLnEJhagUCRWcUrx2OKeQYViQ5
PnwbiF7mjTExM+ZqsDMvZ+YQOdb2o5JDz5XONwtw27OCnMnnoq7rT7Ei5+Qb
5MTy79QTw1S5yp3m/+MYZ+giGLI7vQKC1IjCk1nZWZEoKc9QHep7oVTUpEcB
T1pLR1JpG7X/uhV4cTTaYki+vzoaBv7Ow+1JBRdxmqWAdNhJnBnzdvq7g4EK
pv+A32v2pLgYc/JxGBNT0eeP7iEhIQFZBc+xtSfc9eSAHqlXE5CUnI7OQXF4
zMRyiGd5mbzvg8fldEeovSTwtWzjh8PkQ3m/PJyyndlB+PsF4NWMuPWtNv8G
XD1DUFFdjiCNC24VN8KwOgOdhyvSsx4jPzMFrl7BmF03oC4/A25eYbyvzsMZ
t4sbDqbspFg5FNm1EsWd/vdhFP4QO3NDCA4KweCcQVy2GOiLYnlJWl99EXTh
l9HW8Ax+uniIQ6gbiA7wRWWLHonBAXjawsGi52UB/EMSsG223X9QbFshU9Fg
XlNRZtwZ7fR7b/pcRvp/Poz0T6xIZ2Lfm0dQYDDalBprQO7lsKowXK2xziA7
IRQ5z54iLDgcHcOLvO9sfxOCQ8Iwb4B1oZiens5vrGlsbBxX9J0On5L5sYqO
mdXnanXrzZ09qa6pL4VXt//pYPr5cTwzs+sd7C2OISRYTsE4j+AANoVXYgoL
r1pZmRjBysRshCfkqFPISQhBNp9CBNplEUxH94KDQzEnp3CM37t0nEpGnoTT
yAnFJ1BuL33C96yr3SVxdx0/i9/s9dYZ/AGo8LuTkoSU68nwDQjB8DwlzOuI
CtCiTD/CIY20PmcF4yXUviiGf9hVCOs0IClIi5K6elZcBqBCFpdDzaXQ6mKw
brLVH0XNp+xQz2874Ns5omwX01hsF9MYyaFpvG0WP+IdXw/qEeHtil+f98XI
vIHL5umdRPiEJWFg6BXiQ7xxNbsG6/PD8GY1Y9ELVpJX5uG8WyDGFzdRejcB
vhFXeN9LQV64+qDKnp3am8BJVYG0VkdU+HEUIrkrjtZ4aK3T/tzeJiLFA20v
vUbylVS8XhQeaGdtHvdSEuDj44vE9GwsrosrEnvryxDCZOcXFI7KVmH/O6tz
yLwWz2pRXyRl5GBxw3DozMTzIbIcOy3XpCOUU97CxClMi5hpNlENyWLqO8xQ
uR3Sch5cvEwmkxr/zHS+RHRQi1CT1b3U1tNRPJU6HYqr9MBIV4xlOpxEvQiK
XSwodjHhDNDNisoNmqd5xvEO07HOQK5fv05ux2laIYWcJvc5/vKchtjE5Gs9
i63C4omZnwDmN/jxHtQoXsRgUK/RlNw5zQ+v8TOlyukpOnHEUgTzco9cFnhv
TljrGL/ql8RKaJUzxF2xYtZyY1hSePJ3hV2mWj3xQpsHMrjYv/Ph0fEUk0/q
4w8PfOvDgVQsRoA8Yx/k0H3xKH2LNy3k+f7N0SzgbaB5AbPQylDcO0bK8B0r
0Mqza3/t2EwI5e7urkR34k101s/8SW+l3ujDgz8dJKAHZFu1DiK0lYZESFlb
kzs5GPVqSXrEVR7GOqUgJJpWX4klj3dHSKk7PS7IPJh5pe9ghKcEwpUe4d8p
2+d1qPavjoZ4x5qH0lRPHIj1tAWrYVXYHD3tqRe7v395NNQC2SmBRlngUg+2
+FrvJdPjsNRZSUH4dQlSl/6NY1iVJ6qhGtqiIJ55SNPkM/7so6lfZInHPhxq
a6vhqJVHmugpNVmyWHR6Sb3N7h1mSJfdkbXYoqGEns5hKnu4743GRnynpcLw
2UjZHR284O13+ZNjdHTLYKQDRbuYm5niT1SPj49j/PVrTM/OiUNqxh0Mv+pF
/+AIdvbMimiW5ybR3d2DmYVVu+vYH5ziH8DAss9gVikMLewxxzcFnfs3+Oq3
30Lrp4XG3R1hSdexZthCdkosNF7e8NJ4ID4tB3QB9URvPbxdnBEYHAgXVx+0
DM0ceFXsh6TaqrqcJfM3ofrxLXz9rSvqu8awY9jBtmEXkz310Lix3HtpB1vz
A/B09YB+aAJ3Y3VIeVjJHcbTzCQExd/mV3FbVwc2tva7JB18baI07wbOnnVH
TbtIkFsrshF8KVN6NDNuxQajoKIc4WyccmBvsqcWOlaqLe/Z7kPSSYGCggJ+
HIVepxyn/hM7ge5PDpvRcVZsjigzUhePmp/l4tw359D4ahL6mlyExlsKzYeJ
VGiWIpxR3zEod1t5rWwpNOUWJV2JSVni3t7eUSejbtzwu6l4ZHE6bh1byWXO
1lKN9JY58op0F7VVVehtb+C1sBCcEfNzc9jeFTn/8+xr8I9OR131I4TEZQrB
MsHdjtYht+w5m2oIsw4huKmeOgToQrG4Z78iParOHVdLPq248Yduy5HHmdSk
noavDYt1uxbfwyZsxJPsNHz7zW9w1iMQIwu7fAnnks4bOWXiJpiWygfw0MWi
p7Me7h6+eL22xyL/CHzcPZmGvkZGdACu54ubwp4/uIqgS3f3GZn9CQ/fF1GF
EpalDjHh71Mupjx/KsKr+nC2J7/1jBvS0N39AYkSmu15Xju+bdY/kiIzoaOy
ABdc/DA0tyNtMR9uLh6IiY2BxtUVec+bYTJuIzMpDN5+OgRovRCZkoMdNr3X
HTXwcnNFeGQ4XNy1qO+bOtwhTj7lt4aZui/RDWWi/+/zsGeiyM6qQLoqx0Sn
OylJIr4sdwli6XI7qiStQ6MyYXmlFGnBF4fbMCnryNAwNnaMatv06CC/Pql/
eEL1MCbDJjrbmtHS3oMtS+CaHR9GfX0DxqYXDw5cvydhmFm1Y+pNEjNeG1Ie
ojMN3RYZfaOLeIhkZ5npcqK4HEzq80k7Rv0H75A+/aFCysZrnh/xHGl3VTUv
utbMrnnF2pDzo8NRk0tOSUmxpE9fvA0H+2xe6ZfafIfXTA7gUDyl0Wh0FEff
u+Cwkeb/LxzHlGxeWafgzxupeN4dzafWaKBukLHUk7SQNuo5msEbsu70/eE7
ozplhUruwTGQbCbKAwtcL5e76Sw86fkPHMOkPKW2HxNdC0bWxnhFYFmJpx4M
Jqf5ZshUNwZ9nT51XMUbGhrsiI2eIe67Jk/0s7bBO/vVhU5xsVKMyOMXdbIK
5hPHpkz/U4AW5g7Am3wIXpYG0NaTqj7ejuI9WE1t8H5qjdfKa8nN1v4UqVO3
Rc7Fr0Zq8TnzQWn4nhUNWrHySHdxyb3q/TRoRf0u+jid/l3RQi6CDr+otbMd
WqgP7bluz3CdcJAWO2ZwnF9dyq+q6Umy0DJwQ+z9sh/z6iC/LVI8puTB9XAf
LUTv1pSjtBxoF3SRMCUt5vVRSctJ0UbmSe5sb5P4Qs8k9V4hHZI39tDFs2pu
RHenCVosvorHfK3lei0f4SHXBp0+3s+m7u5ui+n8tSSDPzrQe4VHOQJv3prk
B+V86NYNcoMsOlkuItXyBQjqz8i3g8JGK/6d5BH5NovE0+TMbllxeVie4D0q
rFQLlxRYbOoSlpKq0P9Wo5Kcpe6yWjnD/5POFGuT+7msOt9AT0cr9J292DQY
laFzE8Nobm7B2JT9s77vSxqVhjpdEAZmlasFvovl4Q5eny8RaXvruJUQDi9f
f/h5eyLqyl36L04Y66iGx0VnhEWG8etT6vomDkw735k8S+U6NL+LnZU5lD0p
RGFRIe6mJ+OCsysKqhow0FHHz3iOrxmxuzIKbzcPtAyM4VZ0IDIKxOnR8uxk
6GJv8Ach7T1S9qHoHF4wYn2iFxf/8e8Rn5yGtKvxOHf+Iu4/qUDN82xWk94X
4ma57524YOSXl7EKOxRtQ2KjlAo1exU2nePIz88nu+b9fv6upLJYszgqS2cD
f6606O5lxKcXYX60F2GhYaDnWvVV2e9aPtP/+jFZXW70DpTKLdnK0qfQN9Zw
VVRWLvY25xEX6o+0tHSEhTM+MVVsr85BUOwdy8pFlA55jK8RbFzToDjjMtld
C39d2NsLYIv8x1SfR9nQ20ilnMaE5/m3ce7rr/Bb9wCMzu9IFJ9jdqgFzr/5
B7joBKkzfXX8tpShxW1szvXD080T7cOTuBWnw9WcMpjMJhTfTURI4n0qXQ+l
lD8Hz6ov/rMx7ojB99eXiv8OOGcx+NqiO3D1C8EC7Y0xkA9TYuHprYWvtwbx
NwvoH4BhqrcBWg9Xpp06uHkFonVo9kCD//dSH/gjR43uItitDUhKlGGTExPY
Es81qIpm3t3G2NgYdpV9TeMOBnq70fNqmNqsFy3b2zswvfiWRUuFCn5P9uqA
KFP3Nu34eRsv8Avb7/dF1F/sH086z/T9wO/t/QOr/1ffH0Z/W1ub3f/PpXy/
vr7Oj168+b3TR/8Xh8xoqA==\
\>"]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"x", "=",
RowBox[{
RowBox[{"RandomReal", "[",
RowBox[{"1", ",", "10"}], "]"}], "//", "N"}]}]], "Input",
CellChangeTimes->{{3.70176856649321*^9, 3.701768569875935*^9}, {
3.7017693185764103`*^9, 3.701769319383185*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
"0.9557201322384856`", ",", "0.6176147953859905`", ",",
"0.8716354877689405`", ",", "0.8406792527120477`", ",",
"0.8767124694838984`", ",", "0.2886033560411032`", ",",
"0.43531125109089763`", ",", "0.8501332657200067`", ",",
"0.15170260286550774`", ",", "0.8035678273637223`"}], "}"}]], "Output",
CellChangeTimes->{
3.701768575451021*^9, {3.701769319998618*^9, 3.70176932421996*^9}}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
SqrtBox[
RowBox[{"2", "u"}]], "[",
RowBox[{"[", "i", "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "1000"}], "}"}]}], "]"}], "//",
"N"}]], "Input",
CellChangeTimes->{{3.7017692548769417`*^9, 3.701769263994534*^9}, {
3.701769344403982*^9, 3.701769350211391*^9}}],
Cell[BoxData["933.0263605467495`"], "Output",
CellChangeTimes->{
3.701769220928855*^9, {3.701769255365707*^9, 3.701769265609683*^9},
3.701769352291404*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"x", "=",
RowBox[{"RandomReal", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", "100000"}], "]"}]}],
";"}]], "Input",