-
Notifications
You must be signed in to change notification settings - Fork 0
/
WELCHExactPDF.nb
2472 lines (2439 loc) · 118 KB
/
WELCHExactPDF.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/mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 7.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 145, 7]
NotebookDataLength[ 118451, 2463]
NotebookOptionsPosition[ 117312, 2421]
NotebookOutlinePosition[ 117693, 2437]
CellTagsIndexPosition[ 117650, 2434]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"(*", " ",
RowBox[{
RowBox[{
RowBox[{
"Exact", " ", "density", " ", "of", " ", "the", " ", "Welch", " ", "t"}],
"-",
RowBox[{
"statistic", " ", "for", " ", "odd", " ", "sample", " ", "sizes", " ",
SubscriptBox["n", "1"]}]}], "=",
RowBox[{
RowBox[{
RowBox[{"2", "m"}], "+",
RowBox[{"1", " ", "and", " ",
SubscriptBox["n", "2"]}]}], "=",
RowBox[{
RowBox[{"2", "k"}], "+",
RowBox[{"1.", "\[IndentingNewLine]",
RowBox[{
RowBox[{"Sample", " ", "variances", " ", "are", " ",
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"], " ", "and", " ",
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"]}], ",", " ",
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"], "\[NotEqual]",
RowBox[{
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"], ".", " ", "Ray"}], " ",
"and", " ", "Pitman", " ",
RowBox[{"(", "1961", ")"}]}]}]}]}]}]}]}], " ", "*)"}],
"\[IndentingNewLine]", " ", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"ExactWelchPDF", "[",
RowBox[{
"m_", ",", "k_", ",", "\[Sigma]1_", ",", "\[Sigma]2_", ",", "u_"}],
"]"}], ":=",
RowBox[{"Module", "[", " ",
RowBox[{
RowBox[{"{",
RowBox[{
"a", ",", "b", ",", "\[Sigma]", ",", "\[Alpha]", ",", "\[Beta]"}],
"}"}], ",", "\[IndentingNewLine]", " ", "\[IndentingNewLine]",
" ",
RowBox[{
RowBox[{
SubscriptBox["\[Sigma]", "1"], ":=", "\[Sigma]1"}], ";",
" \
",
RowBox[{
SubscriptBox["\[Sigma]", "2"], ":=", "\[Sigma]2"}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"a", ":=",
FractionBox[
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"], "/",
RowBox[{"(",
RowBox[{"2", "m",
RowBox[{"(",
RowBox[{
RowBox[{"2", "m"}], "+", "1"}], ")"}]}], ")"}]}],
RowBox[{
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"], "/",
RowBox[{"(",
RowBox[{
RowBox[{"2", "m"}], "+", "1"}], ")"}]}], "+",
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"], "/",
RowBox[{"(",
RowBox[{
RowBox[{"2", "k"}], "+", "1"}], ")"}]}]}]]}], ";",
" ",
RowBox[{"b", ":=",
FractionBox[
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"], "/",
RowBox[{"(",
RowBox[{"2", "k",
RowBox[{"(",
RowBox[{
RowBox[{"2", "k"}], "+", "1"}], ")"}]}], ")"}]}],
RowBox[{
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"], "/",
RowBox[{"(",
RowBox[{
RowBox[{"2", "m"}], "+", "1"}], ")"}]}], "+",
RowBox[{
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"], "/",
RowBox[{"(",
RowBox[{
RowBox[{"2", "k"}], "+", "1"}], ")"}]}]}]]}], ";",
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[Alpha]", "[", "r_", "]"}], ":=",
FractionBox[
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"-", "1"}], ")"}],
RowBox[{"m", "-", "r", "-", "1"}]],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "b"}]], "-",
FractionBox["1",
RowBox[{"2", "a"}]]}], ")"}], "r"],
FormBox[
TemplateBox[{
RowBox[{"m", "+", "n", "-", "r", "-", "2"}],
RowBox[{"m", "-", "r", "-", "1"}]},
"Binomial"],
TraditionalForm]}],
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "a"}], ")"}], "m"],
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "b"}], ")"}], "k"],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "b"}]], "-",
FractionBox["1",
RowBox[{"2", "a"}]]}], ")"}],
RowBox[{"m", "+", "k", "-", "1"}]],
RowBox[{"r", "!"}]}]]}], ";", " ",
RowBox[{
RowBox[{"\[Beta]", "[", "r_", "]"}], ":=",
FractionBox[
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"-", "1"}], ")"}],
RowBox[{"k", "-", "r", "-", "1"}]],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "a"}]], "-",
FractionBox["1",
RowBox[{"2", "b"}]]}], ")"}], "r"],
FormBox[
TemplateBox[{
RowBox[{"m", "+", "n", "-", "r", "-", "2"}],
RowBox[{"n", "-", "r", "-", "1"}]},
"Binomial"],
TraditionalForm]}],
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "a"}], ")"}], "m"],
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "b"}], ")"}], "k"],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "a"}]], "-",
FractionBox["1",
RowBox[{"2", "b"}]]}], ")"}],
RowBox[{"m", "+", "k", "-", "1"}]],
RowBox[{"r", "!"}]}]]}], ";", "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "\[Pi]"}], ")"}],
RowBox[{
RowBox[{"-", "1"}], "/", "2"}]],
RowBox[{"(",
RowBox[{
RowBox[{
UnderoverscriptBox["\[Sum]",
RowBox[{"r", "=", "0"}],
RowBox[{"m", "-", "1"}]],
RowBox[{"(",
RowBox[{
RowBox[{"\[Alpha]", "[", "r", "]"}],
RowBox[{"Gamma", "[",
RowBox[{"r", "+",
FractionBox["3", "2"]}], "]"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "a"}]], "+",
FractionBox[
SuperscriptBox["u", "2"], "2"]}], ")"}],
RowBox[{"-",
RowBox[{"(",
RowBox[{"r", "+",
RowBox[{"3", "/", "2"}]}], ")"}]}]]}], ")"}]}], "+",
RowBox[{
UnderoverscriptBox["\[Sum]",
RowBox[{"r", "=", "0"}],
RowBox[{"k", "-", "1"}]],
RowBox[{"(",
RowBox[{
RowBox[{"\[Beta]", "[", "r", "]"}],
RowBox[{"Gamma", "[",
RowBox[{"r", "+",
FractionBox["3", "2"]}], "]"}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"2", "b"}]], "+",
FractionBox[
SuperscriptBox["u", "2"], "2"]}], ")"}],
RowBox[{"-",
RowBox[{"(",
RowBox[{"r", "+",
RowBox[{"3", "/", "2"}]}], ")"}]}]]}], ")"}]}]}],
")"}]}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}],
"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{"Plot", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"ExactWelchPDF", "[",
RowBox[{"1", ",", "2", ",", "1", " ", ",", "1", " ", ",", "x"}],
"]"}], ",", "\[IndentingNewLine]", " ",
RowBox[{"ExactWelchPDF", "[",
RowBox[{"1", ",", "2", ",", "1", " ", ",", "0.5", ",", "x"}], "]"}],
",", "\[IndentingNewLine]", " ",
RowBox[{"ExactWelchPDF", "[",
RowBox[{"1", ",", "2", ",", "0.5", ",", "1", " ", ",", "x"}],
"]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "5"}], ",", "5"}], "}"}]}], "]"}],
"\[IndentingNewLine]"}]}]], "Input",
CellChangeTimes->{{3.5121484034023438`*^9, 3.512148436973633*^9},
3.5121486408759766`*^9, {3.5121491386435547`*^9, 3.5121491420302734`*^9}, {
3.5121492432070312`*^9, 3.512149263265625*^9}, 3.512149310192383*^9, {
3.5121494035810547`*^9, 3.5121494270029297`*^9}, {3.5121607160908203`*^9,
3.512160756665039*^9}, {3.5121608011445312`*^9, 3.512160829330078*^9}, {
3.512161284375*^9, 3.5121613161845703`*^9}, {3.5122345791640625`*^9,
3.5122349122001953`*^9}, {3.512235226727539*^9, 3.5122352755126953`*^9}, {
3.5122357226367188`*^9, 3.5122359452822266`*^9}, {3.512236010004883*^9,
3.512236038942383*^9}, {3.512236490794922*^9, 3.5122364945722656`*^9}, {
3.515074211022461*^9, 3.5150742258935547`*^9}}],
Cell[BoxData[
GraphicsBox[{{}, {},
{Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw1mnk0VV/Yx83jNdx7ydBgTMhQKVF4tlRoFFIZMwvJUAgpERIypZSQeQpR
SHK2MSKZ71WJIoREmX6E17vW+/5zzvqss8+z93mG77PPWlvC5oqhPRMDA4P6
xuV/7zfGTv5ZX6di8tLXP5pqZjDqv2n+8z8qfrsOddKnzeAU/+BixX9UvOi+
o0fcwQzE1d1Xr/yl4hmhrYIXEsygLiKOY2iUip27Sh3bZs2AQ4m2DX+gYvtU
74DOEnNIuHrx+K0nVCyv/PnlX7CEVQ65U+aPqLgjjcRedMESHJ7OGqg9oGI1
jtVrt7wsYX/TbZOZaCqOS5jldsuxhE9CWdYXb1MxQ07nezd+KxB/89MHXaJi
fY60oC1jVlDA6JXJsH9jvUkPZqe7L4KaUT5bkgoVe1hJ16YMX4SmzG9Ou3dR
8bGMU7mWfy/CkK6Boo0cFXOsd7HwUa1hU7Rief3mjflE6qZVjKzhtuj4u1AG
KpYVCI5K6bOG8yoWk6T3FOw0XFHmP2EDoyEJJ7OaKHhWo/KZ4T8b8OprLdas
p+DM0e0kDV5biL5+4KpbNQXzys7mHt1jC/WE0OrHIgo2Vj54TNPfFpROdPHG
x1Mwj9bn6jp+O2C1190jYkHBCttKFZuP2UOX5e5zkxco2GpQXJ+wsIeU85sD
3ppQsIcYTDW724Pqid+N1gYUTF0PFBJ4aA8OKg8v5OtQcJNDgGzkiD00M47d
OihPwVmaQiUCwQ4QmRL60WKRjANsd/cJfXSE84/c55TnyPgac9PpsBFHkI4z
FWGaJeOUzvpJ7mVHeHtHyTZ7goxn+ENHjLc7wbRr38L0ABnLqLmcGQ9wAoOD
MttuNZDxwZAdSnK7LgGV3uiaHkvGBy7aLVXmOAPe+Sw8JJqMk29Fr94knMH1
ZkCmwz0y3u5A7LGkOUOTjMoX+TtkLPhSPvoiuwv4XUs7XuZLxjzjxa9uOrnA
d4q/fIMVGWtj2xBNZVcoPblr/IciGUuWtb52/ngZrJ5xszTLk7GW2PVqvYnL
QJofFcvfQcY6TdU7tVndwCH56Tk3CTJ+0r+XCDroBqKTXO8WBMi448fX94fz
3OB2+I8s9n/8ePBcm83k3StgWP/YVu49P+62zNU9XOIO0j0jX9818WO5mrpI
jN1hYUTJ1KGeH+saBESf6nSHJLZ6g4xqfsyXCCJpf9xhSH9Sc2sxP7685da9
zaoecLnjoDA5gR9L2JB3hWIPCB/40rZkwY+5qdZidcOeYDoto/fQlB8feJSp
m7rgCQrr7vX7zvFjf1U8HsPpBR/FWas8DfhxXxdp7oWyF2yyU8yZOsSPLVx5
Wn8FeEHGxI2goR38GM92v7befBX8f2u6xgzz4bfNu/InnK5BsvX+cTTAh2Mj
qfMCN65Bdfcuu9k+PrzwW5xmFHsN1iqkzA3f82FcOPmar+oa3L7FeULgBR9m
aeV0v03yhghy785HgXy4I38o1fGVNzxWcZ1MFeHDgV1WHn1CvvAmy97RgMqH
r/jWV2ElX/giZDXMwMOHTdY9EuuO+MK2f2c+WzPw4aFIVReeq76Q3rC/VWqM
F7epmVjWdfpC3lnmgpyXvLjw7GFnHHsdKn0eOxef5sWn2vhUW6X9QdLOI4uu
x4ttPfT9uLX9IcpAb4jxEC/etP+Zt4uFP9jILZw13suLlxYJxXuJ/sD9+Yz2
f8K8mIEfVZ3kCgBLLU5hne88uOveGjvjUgAws/o29l3lwYzBvFS1X4HgNnuK
gcGNB+tc6b/9hucm0Ae2H5Rz5MEqK9OvbJVuQmF5d4n/BR6s/M9VQMr9Jpg4
KSeLa/Hg/VEXg23mb0Ju66inMzsPPrbpv6I7wbfgVNxZ8bUkEqYWpP4nnRkE
fnv49HXjSbjmzJawsPIgyO5q9oiJJGHzppaUteYgWKNo1EveJGF71fg2+ekg
eB4v6aBrR8Jxd6jvbqjfBtKD3wUxSiQ8sepKSPXchvcP7+6XqufGt657skiK
hsDCfp2LrtXcWOxyj8tR5RCQpP8Lf/WKGzvKWibcPhwCfsLun3RzuXG/Gv3h
oSshIJdkcsM1ihs/9OTyHWkMgfDHUvWvznHjsF9XQpt97oDu07en9Ka4cM+a
Lv/QdCjMRH35xzLKhcMbNA1jWcPgceBKfu0gF5Y0KbtuuiUMpqzU2Q92c+Hi
2G3f9h4LgzjJV4RiFRc27ds6Q8oJg685hbuo4Vy46cUBC2v7cPAufUL5KsWF
Fay5Bs3m74J4RhV+vJULlwW/NdDii4CW+H63c0Jc2E7wopiGXARsuSbU2sHF
hcfUTx6LsIiAuv3xwfWznLhaM3fr1XcRwPs2Yi6X4MSR6rdNH6Xdg+x31/u8
zDjxub87RtecoiBpj7dOnjEnFo+KPspxNwoin3qWfD3JiZ0STQJV86LA08sl
Qh9x4h9bVsWZJqJAU8xCS2w7J95VzLxkejkaerxR9vtpDswbpPpG3f8+MMmw
e0vc5sC1AZMa3p9iYD6GedjEnwNn5SRJzE3FwPjK+unIqxzYM0rbJGw9Bto7
l+QWHTjwlsI5yrR0LCTdmBhoO86BG4aUHhHusbCr98MRX0EOfD/jq9hdUhxY
hiQIfsxlx0VObW5uZ+OhHRbYfdPZMQ5MSjG/FA+wcm5ZPJkdt4akH754Ix7E
PTcPet5nx5Z/txrkZsXDd6v03E3e7Jg/68MfzcV4cDhYctBShx2jezzbE1IS
wO1Pq/WvATa8VyVT5OHaA/j6XNE4kcaGu9VSq2wEE+H0pftHoZMNC1mjH0cV
EmHXkOHO2AY2rLJuii6YJsKfD/1zewvYsLFJ52xuRSJ4542FBfiy4c8KR19x
+zyEwIvMRdxUNmxyklaRypkEJR91o3tIbLit+8FrRakk+K4V6faUjQ075Qpm
92kkwdGtgsrKy6w4oySz6Lp7EvB+kikx/MaKp52eRbP3J0GK0bHSpGJWfCJM
ztas+DEQR2PLd5xgxbkVfL7FHsnw51Vv4swRVrztpdWzruhkkN4u6vMaWDHp
rJK1QGEyhDNn7D+uwoqLH+aKz40mgwF+WekmyoqlxZ9yy1o+hSF1etWrcRZ8
vqP41A3DFGBQFCN07rBg3WCbO0NyaRCwSdvyTiAL3sTNWzFyMA0W121Wm3xY
8I/dY0arJ9Nguiv74DFnFvzd2XKHq2cafPZVqjQ4zYJz32plHKpOg/IGzRcW
IiyY6+WBp6NGz+CyuXmGz3NmPHzIiqMwPh3GjwQeep3NjBkcet9fyk0HO+W0
b8upzFjn0alBrbfpYMo0IhYYx4xlWNUvKo+lg16uS3KIDzNmzyrxUtDIAOk5
vwex2sx4ZE9IUeFEBny6lxRW0MuEH5ZGtvqYZkF3o3jhxY9M+LjmL8MK9yxo
W8/pEGxhwrJHbWZ4wrKgxqtc5FY1E2bb29jA8DIL0s26C4wzmHD8+03q67zZ
4LyTp2PVnQkfSkH5ni3ZsPL+lrABiQmTPUZ1Og1yYY6FQ5OVjQkr3jrwcZtz
Lkxr3beuWmfc6F8Sh0ODc2Go9Gm+9F9GzMZwTS+jPBcaHldp/NfPiPP23Tol
vC0PIp3nLqbnbIy3NiEz/8mDLVxOeXOHGLFQbc+H4KICsKxpi5TUYMTKKp8i
cj4UwDPP3e4G+xjxtqEq3V9TBSDzeVm1cAcjRqEuMLizEJQLoxpsSYwY67Y+
vVhQCIdOlQ129TJgF4+lhbmS53Apbk3whRMD/p4UrH2EVgyFR23/+2rNgO9u
P7M+ulwM08vvvpDMGLABacj6HF8JeNrGZjidZMDlyRUe51VKIGCv9G7xPQxY
wavtrnZACWx7Hv9o2+I6cZhzcvEt+QW47JNn/bxlnXBkUR/sOl4Knqyvc7z5
1gnrobWCBqtS8OvVPUZhWiek3vY1vvcqhYhr9vf1x9YI/6jsbQLJpZD/Kk2k
8sUawfP6tuDTX6Xwc5+Q0oOjG8+L12MfJJSBoyrLudPuq8Tf7Yn37/15Cbb7
v+bW168QHquKW5tnK8Ahd7/J0YoVYsb5fGcEZyU4i8Qyt+SvELvVdp2zkagE
9xUdy/bYFWLzvyJGgzOVEEDkUfutVgj/C6x/ZUorIUH3WuD0yjKxr8dkc9tG
JTWdIxmL7F0mblBC76hS34C874F1t8wlInxC6JplZA3ggVPyw4lLxO/6NZlN
hTVgomN79tzdJQL9nswfb62B27yRBeC2RJzd+Hf4yk1Af+bAOX61JWKFGqRI
iiIgvPNWyYvWRUImggg/eAzDqPw767k/C8T1/2i43awWng0YNVzXnifubenm
CMmph7EgQQ1d4T9Ek0/DzdN+zVD1QPbIuMkv4ujaPd8nxu2w62f/nPLTUaKQ
eatwvUwXXJ1RlBMKGSVeXNG9WbW7CyoXgyzWnEeJkqerfUUaXaDNJt/UpjZK
HD91kTXasAsMpfweXur5QZRv3SUkHdgF1yxED2Rw/yD6VlqeJvd0QVWn6c1N
fsOEMbeRsWxgN6zRi8pWLw4TtzS4biaEd8OhIabxEd1hguLV+GM5rhtap/MN
ygSHiWLn/TzlOd3whXtFwuDFd2KHoN/x4Y5uWD/ypOHu+Dfim89sJqtkDxyu
+sy5ajJEKLVE+rXU9EBRfpb6u4NDRBl3IdO2lh4QfnLlUqz4EFFTrm7n1N0D
fYtsccd7BonKJDGtr2M9sN3fhP/p8lfCvc5U6Dh/L9QFzZOQ7gAhP6nRwGfV
C/+iVdhCh/oJked/5fz+9kLHlXf2ovn9RN5X1/Az/3ohw8CsscirnxjNjJPZ
zNoHxyjBITS2fqLnNqksZFMfJD3oYpJTpBPSgWe4itT6QDXZY731eh+x+sRr
SduvD9zzipfI1G5Ct+RYUMVsH+zidQs+NNRFFOm7GSgt9cGspwKv1/MuomLY
wShxrQ88NfKkenS7iMCfP3PVuWlw7WP6qcSATqJ0IDi0VYoGfyd26u9r+0is
n2OSJRnS4OrLXMtEh49E6bOsx1HnaLBwY/vVRYaPBPmykeS6OQ2W+MVSX+9r
J3jm7iqVOtLgnyplXjO1jWg6SfcyCaABW/BS2lHPFkLlzz+eiHQahJ24Vp5D
aiFUrXIcWXNowLHpTytHTjNRftv425UCGnDlTS68//yOCGt4v0p+SQO+j19P
nD7SRFhH2dnfbaBBzCMzm5LBRoJPjKoX0EwDig3dh+zXSKyyTX+zaKOBwHxn
endRA0Em+XL96qaB8ObGpXPC9YSE4Zkf6t9pkDSizfu6tI5oPXtV9dkPGogW
1UiJnqwjGrTyvi+M02CrduWpL0G1RP8dGcnLv2mg026iX9tHEGnapyVDVmhw
1NrgajA3QewLuaxssUYDvTn91COohnhxakBPmpEOJ0Q151vyqgmDHIcbAWx0
MHGQetYV8JpQO3GOyZifDhEeXuFxryoJvn8CarEUOrwNqL9iOF1BzCrnTlcL
0EE63lary6qceLVeWj8oTIfzKWXb4x69IlQTrHv7ROlwL4+Zx7DzJXFvzuhU
9RY6/CEyPnceKiOeMxOWhuJ0kGn9WxfrX0rId2yeZJKkw4U+nfwzL18Qlr9D
mp9J0QFPDft2bi8h9Nce2uTL0GH8tyOWfFtEzIzRBAVk6XCajeN+utJz4sx5
A3ZXOTpsU9FTSCfnE+c53BNGd9Lhjv74skRwLtEgZivLqUjf2I+GtzybyyaY
0mUmRZToYOQt+0jCIYsYknXuF1amQ1Vks8MzWgbRHCEww7qLDhIZTvsk9NOJ
yd86Ct82OPw1B8uzqjTi54+T9/J30+Fr6pks7fMphKwnwWWzhw7P22IPaR98
QrhIWuSyq9DB/7/OQbTtEWGwRcrmyQbry1BuIMYHxIj8EdWte+kgZGQoikZi
CclWJamoDW5dl75yMzqKSFYN3jG5wU4XsvICjcOIX65bD+3fRweWsu0jN0SD
NsSQ39Njg5HICc1lfI3gPpxdnrTBM4+S7GaTbIhtjDZ8JRscwJTQeFRNH7pb
Vv1fbrDSY9t3AqpX4Ga3y7/sDf4q8mKamy8ANu/piwnbYCMVenWodghsoaWo
m2xwef+522H+EfA92mqRssE8eFKqwTQG7v77+f7txvogviC0wSceLBN0X5ls
sIeDy8+GhESoTHhdMbDxvenqO080vkgC69XKbsMN7iFNFjW2J0PEnQWu8g1/
sQ3lk5smU2FHjKkV+wYnX9h3wDbkGUxvY+0+uuHffToPXIYpGdBUf9Lp6ob/
2xXmk22fZUJNTL1EzEa8GBherdvW5EAEYzgpbiO+ST8Fdo+cyIMrWZyHfRXo
sKf7qo3d53wQ091WeGIjH+yy9zbaLT2Hm6hxC97In9X7CQsjocUwYKmjaLVj
w/8ay/GlESVg3bUzcGr7xngZ08Y9IqXwNPTy386N/Ny1meur51ApyM85SshL
bLzPV7VQmlMGLB95M9zF6JD4n4isiuorCJlaKm/eTId3bf0RKkaVUJ0p15BF
3Yin76d3Zjtfg0issD8i04FD+jNLCHMVRPjqpDbx0uGY/5fAnrI30D8XafSA
kw5tcoOe1wQIEFOMva65Ub+uvYPFT6cIUEpbz9q7Ud+hyH3BxgZD56mg+0JL
NNgx5a5y70At8D9cbEqbpYGTjufzzxN1ELMWR60bpoFa0p/W3yb1cME70lN4
aEPPfntOMNfXA7V264zpFxrkPfbaofCkAbaFfQvI7KHBxOzVZwHHm6BMI9Tu
+Iaeuab5PNha2AIMC8K35lNpoLGwWLZb+D18eeyRJvCEBqQTvl1HQt7D8PXz
PeKJNHi+6MvnZt4K4kXp1QyRNJg+5Xe3hvQBZsy+hQx506AysKamWLUDvmv/
stylT4O0vLUPr107IP3BhKWcDg3Ce7QG6tM74PTVFxn8mjQ4L0+s0Hg7oS7F
5HT+Lhr810uoMYx1wonTm0NIQjQ4oFhbavCwG2q1RApSh/pgX09fYk1LN3Aw
xHIwftroN/5Tfgr/uoG1TEveuLsPZFqEdDisewByEtk7GvuAau/WTcj3wps+
+1Pi+X0wlbx5Tqlmo48J+kbtde+DVJL3Pp4f/aC2Q1n97Z9eeFx2T8Rf6BN0
seemBE70QqLps9Vx/U/wBhvZ7/3eC5G5bY0NRZ9gU5XXlRudveB7WMokwPcz
LK2ET4QV94JhQIf3JPcAuC4H6cZd6gUmb8EvdiJDME8Xg619PWB6Mb+rcu8Q
+Bn/2P6+tQdKj0ELyWAImlO+vnOt7QEbsUvlL0OHILQ36GV8YQ/UNVfHsMwN
wbsER0+X2z0QvNnuSGb7N9AOPhllrtADLLWlRcPBw7D/5UK1gF83sJMMgmym
R4Gjom7rXuYuCB7q0/HM+AVsqcrTv/nboeRJqEHgn1lQvJ+8MlvRDFfS7fhH
TOahikP4oPrdeoi+npQcbTsPBUxt8pY+9VBk0C6r7j4P4t9PlkXb18P02n7t
6LvzcPmqdgrlUD24mnF7qlXPg5NpjHb2Sh04C5R2R0oswJiaV8i6ex04hK4n
7ptaAGNbM4EIu1qwdH68NSxoCY7zfq84/bUGVvo522ejlsB8xvGgYG0NPNK/
Hmj+eAlM0Y9jsxk10Cl7fnB32RJ8ilhbmbpUA4fHBJ99GVkCUe01yT8Lb0HO
LlZ6j95/MPdJak1J8C3MWYTvHOBdhqONVcl6Nm/g7hlvdZXkFVD7c9LMWLMS
4trkXQvzVmDLdebU21KVkKw7mLK9YgUed1yRb9nYP5do6rEId63Ax5m+4ZS+
CqDLiX78x/4Pqn38XF08KkCWqcauyesf3CUdmA8oKIemUtaYC8dXIT1/T6jY
7lfQoVhV13V+FfYUe60dFH0Fn3Ld5o87bHCF1LIb8yuYTqGZagWtwlnfDhDu
ewmb7uVtl6pYhRvc6pbeAS/B3vbkmynJNTCu+/PI6MOGDgk8GL25vAb7k5/x
y9wsBcp0IC+FfR1M7PS+n3ItBYlmJ9VM6jp8TXuXevdCKWgFaIQ2K6zDpbOU
gr17S+H6yPB2stU6VJ2of2I18QKmX+6xT69fh1yqeO7qhRcwJT/uYXOUAQ2o
x1kwHC2BOaGBOwUnGFDjcnbDyX0lsMLSlTRnyICC2XUDCqRLgGvwTW2oFQNS
rv+W95C5BHbE3ycX+DIgbD6mbBtSDDYrqi/+5jMgf6u/UrUhRUBrvTMTwseI
+LD8/szkQkjQenrk7CZGdNGsWyE3oBAMX7x8vH0rI5KZHiTazQuhPXH4cJM8
I3pPt4fgrYXQZK2dxH6UEa3S7yqHphVA+eI/7YgARtT77tF5emE+JEpdjb//
kxHxHmp+0zOQC8aJ98atZhjRkYq6md+1uUDhzNDctciIBLn+jspl50L0dOdY
BwsTcomwT/3PLRdCXytrUMSZ0MDB4+ssTLngfXpyJMGECWlFtee+UMyBs/42
qo/rmRA/vfv5w8dZkPhUgG2xhQmxyi5z/wrOgj6iqdeogwm1qPQLmV7OAhOW
ndd4BpjQR/vKr34oC85H/S27tcCE2oWnX/WMZYJpashuRzlmJNcmIdl3IBMu
1mcr7L3PjLQqJgif2XRI+3H+X8wDZvSi9rnBo4F0+MbO3fbrCTMKldH5/KEl
HWxOXHHJyWVG+bcs7/mkp4Nt7/580Tpm5PGcRaTeOB0cxpp3MM4zIwONr3Jy
b5+BK/ekZLspCzpJZdfLTk6DyetVLNzWLCjQ19828V4aOI/fHdV1ZEF3arWO
PfJLA8dG2fxaLxZ0X1ZepfN8GtjcdNj9KpIFjWTMhwluSoMLf79Bcg0LGsz2
USzdkQpHP9PNnSVZkaGO+0T+/qfQqJ+rlSPLit4W6716I/IUDlf6iI8osaKd
I83DkyvJoP1g04jFQVbEM5w9lYeTQeO0scsZY1ZU+1P8QNLxZNhT//G6Wigr
ihtJuynn+ATECpoS2SZYEd2mQFjjZRIYp1a+fj/DipibWpibHiZBRHz+l+hF
ViRmUeHs7p8Ec/7RkkKsbOhr6JiojE4SvD9hUrRDnA3tlqwTIPc8Au/pH416
JmyIZEO2a1p5CB92s85H1LGhLyKKx99aJwKTzKLQ6RY2NAU1jFT9RFAT/XmA
2sGGKGx80yG7EiGd6cOtJwNs6Iew8Hj1+gPw7U4gFS6xIfaozRPXUx+A9DVp
6Q9K7GiMu7qbNJIA/pU6RnxP2NHhEgXtulvxsMe8Y8bqGTvaLjdz5KNrPEys
m0eX5LCj6+z2SbPn48FU17vlzEt29FEPdV7bHQ8H+nIh4QM7klau5/YbjoOl
vzwKIgwcaCWOpGR5PA68leksUvYcSOLyIrEqFwtK3XbpXi4cSPbIwc0RwrEw
6j0LDR4cqOCf1iY59lg4W8Plbx/IgSp2oicpIzGw96TmbM5DDhSsR5TZp8XA
X+f0AYX3HKjIWmjBcHMMuGe7lqsqcyIrPWVJfY1o+Ppx8wNeVU5EXqgz27M5
Gk7+995rVIMTtScUNmgvR4H8Cbndicc4kWC9ikdXZRSMzP4omLfnRHU07gvl
+6PgnKbVs5fJnKj20CFZG41I0Oo5E7mHiws9KfT7I3ExAgpXGVy4yFzoNLV3
5IZOBGzeUaL/XYgLcfGHtSzLRMDSdV72uO1cqDA4td1k+i6UirUGzSIuxMKg
GBcfeBe2uxz2KfHhQtVnuEfWMsOBxLTfRukHF3px9GxZMkcYvPZhNhqd5EIN
mkYOh2ZCweHXR52UP1yo6pDCDWZ6KGC603ZeBm7UfcrmHC0nFLyKH4//EuVG
r8n+J+r0QuGT+ZrbcwNuRFAPS0fcvwO5FY3+CtXcqL/iDJl/VwicVYxzHanj
RtSGUqMx0RBgyrC0SG7hRi8NLx3uYw0B8+hFLRKNG7XLT9ctfw4GPns5pqlZ
bmQYPaNudjcYvClR4QU7SGiX7OsMm/HbcOSy8QP5eBJaETA+3lceBFfj3jBc
SSKhm1eihAuzgiCzQtK1LJWEjrXJV8YnBAET0+whjUISCo3ylHjkGQQ4MWrm
ZCMJsT2fD/FSDoIDuOm4xyIJfRDsELErvAVKVHWmSjMexBzyYuiC/U2wUEu7
/M+aB1HKTG336N6ESAv2fuTEg8rr2Njk5W7CZE5vcctVHqSl/0fa+1cg5Gp4
WnyO4kE7m35LBfoEgpRDQeUawYPuCfhf3RV7A4Rfb71yRJoXZfmOkdEnf1iz
zdlM3smLhpV5dYPf+sMI7+7mL7t5N/RN3H00zR9K7A5LXAVe1Dj7joPk5A+6
/C7d6aa8aEpYSb9j0Q+8HSv3M8TyomSKutbWrX7QK2jE8GaVFyUtv/fsvuEL
VfhLYSgLH6rPDOPqsPOFNBeHC4bcfChc+Kbm+HFfcK29XvpTmA/pFWbPOYv6
AsvlNDvhvXzoRN+cQnmlD+xt+NV8zZkPyZadddZb9oY4z/DYXTQ+ZByS1RAe
cw2W2worwwf4kGUT5lIOuAY2OzoHh4b50JHiCLEJx2ug8kVEKXaGD/UryJil
wDWYSD/7Y4mFH23e+Tnv3++rEEkaC7qwjR8pmzt73DC+Cu2DnG9Ez/CjGxKr
y8s7vYCJJ7n9qDE/eln0WZ8s6gWqB5S+e57jR2nMqFGdwwtSEgw5Wy340bsu
0794xBOu6D855+/Mj3JuDvNkp3gCf5nC3OdgfrTmsvVelaAnGIaeVnxazo/c
OE443uP3gOwLM/qvX/MjXqXp/RRGD/hPIdaht5of+WWTTubOukNad1cqTz0/
YgjKZGHqdodf4iaUwI/8KDvafOpSojuEvzFftBznRz5Vjr9XxdyhZuZSrbgo
GXl84Lid5O8GlAaurxpbyajesahA46Ib2D8sWD4vTkY9MbXJy4fdgKT1SyVW
hozuh9jbN/G6wYV77tmMKmRkrFsT/jjjMsxt97n3/TgZzQ8sHNja5QryZiFn
M2+QkdKzEbtubRe44FFOYQkio/Hj/rNRO1wgPGz8o20IGcWVCMvY8bjAj7IT
x6TukRGH3iyXbb8zpJEEIf0RGe18NWib4+kMgjVZsmllZCTpcfCjcv4lYJB4
t/LkJxkJ2Rqneio4gfL+/ypXpjbsrWbedhJ0AsuTO73NZsgIeuQue685wpvr
92dEF8no83P3qf4OR/DuOvcjiYWCui+oFVR7O8Jk8Hj7QzEKKky+p/+pyQH6
RjnT489S0Jzmxb2NPvbQZHw3+fV5CmrotJgYtrOH8jqOh4NmFFSm19UtYGgP
iSnskTttKOjwvqb5QgV7OGvC6l3vRkErF4adX363g55GhmN/wyjIpyCFhd/Q
DroyF2cNqyjIq+6f8Q1NW6ij+Ez5vqUgAz3kRSjYQumthdEUTEFyHrX3yVts
Ic58/vNEEwXphtcG/Vi2gTMCfxtvd1MQ/fStoGOvbaAjeDqpdIqCXGV6x5bV
bKDd5oc2RYyKqGdIMeE61vCNJWfpgSQVuZ9nOftnjzXMZTsVi8hQUZ+Zeupl
SWsQmZzcIqFIRQdSH3s9YrQGO68/i8oHqSiY6eDfanwRloPXn588R0XHB92N
dh26CDJZIqJ371PRbVaeKb90S1DX/dxBit+wTxTN2t+xhBM/k8NiEqnIsFbH
4oqTJXgois8/fEpFPwY3D3xRsoTqVzId2QVUNJal37ztrQUYNqmENryjItMY
eSGJQXOwd5rX0Gulohzam+rv9ebgy1Xxt7Wdii6rtr3Cuebw9NQB6+7eDfvj
slKdnuYw3oc0vg9Tkcs2I1/EYQ7Lvkx/7ceoqH3J4Oi3X2bAs7kh7+cEFXlS
xC4/7TaDPVa6QrOzVFSCHh5zTDWDI4yc7VfnqYjDyLDh6h0zOJ/xPmRpiYoU
jqxveepiBs5HIg8G/KOiu7JBqkNnzOD/zt+g/z9/8z8YCuHr
"]]},
{Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJw1mXc41W/4x4/jmB3HOMrKl4hQFFJZ3R8kGZWSURmZSULHKkRGSBQZZWVm
b2UUPh87e59TEYokStkyf37X9fv99Vyv6xn39dzvezzX9Rywcr5si8fhcHV0
ONz/jg9+nF/Y2SFjifxt78MdnGDKe9/y8CYZW4l6yC3s5wQXOMZWK/+RsVPE
h8UN0U4grOiy5bxIxkIs7nxRqHWChrDnzONTZKzPcZ4yxOEMzDK0/7AuMta9
fNuT8N4ZYtxu6D5MJGOPtUXGaiNcIJ/ONRN3kozp/Nk39ZWJAqcM8hjj5cmY
fGzYPoO9FGjJ/Gove4yMBW5dV+sRocC4lr60lSQZq7yxXjCtSoF9T6UrGgXI
mCWf2VyqGwUC+Kdbg3FkjD7gfPvRSQqYyJvNEtu5MFNnBVtSpysw2GrJ8Zlx
YVlB1VeO1LtDv7ms8exVLmxBS9YortcdXpkI+NQacWFUXSE94rg7nND702yp
z4Xx+0Sp7t9xBzv5F1fzNLiwXI3jl7tUPeAD3Y+HylJc2EnmceUHtR4Q/iq4
x2yVEyNxpHWYtXgC+WOzY3oUJ5Y7mi8VifMC7HBaaNBTTkz41WdkjssLHP18
Mu2ecGIv7B1Dr4l5QYu4/IjUI04swF5j2EzHC7zcU3XL73Fi4lVjoxyxXvCN
y1uqyYIT+8vgvi4n7Q1l549Nf5fmxPCLkcV19j5wuTHBWrKdA+tQIRYNsPvB
wcHJ0dYWDmzsa0kpk4QfrEzKXLNr5MBe2xcY6CN+EM/YqJ9Rw4F5F7XVCd31
g3HtWVXBYg4MEy3QqR7wgzu9yrycMRzYNScVYbcLDyH0y0jnmhkHVrFtvDQu
5A/ef1QdIyfYsXVPS0mGqQBIsjw5jXxhxyIMU5LplwKgZuCYzTyVHVvjNbZh
wwfCdqWo6eV2dixHZpJD479ACHjIosddyo7ZPL64fdI4EMI4hw6/9GXHKr4Z
M1/qCIQEecfZFD52rHFFtCW0JgiqPBMcii+SMFy5fb1nczCI2Nx9/fEcCZuJ
4/dh/xgMEfrnxunUSVjx/cH2yplgsJJcMbxynISlhVnZKHGEwJ7hS2r/eEmY
PF28LotpCJifZuHV+MaGLdeLw8hKCNAz3GumurFhdzb2z3YqPoYLzw2Ft+OJ
GH+4enLL2hPwkmPX1oomYnp8ga9S2MMhq//D3chwInaIw3wrXDwctrlUGkX8
iBjOyFbk9ZVwKIwWsdOyIWLlbG/lu0rDgRj7Jz9Shoh9S3rISHKJgPYXj0+K
Nu7BEkxWjSY2n4JWcu2Fc79YsTYmMDWSiYK/ESObhClWzGF7r1yBchQk+G7k
1Y+xYh1XSvTZtKPgl4Uik/IAK0ZfaRU1bx0Fz0XeotLvWDHbdiMNjYQoGM0u
OEYOZcWcfjqlcLM8B4+yRK5RUVZMQuED0/s/zyGr9T7V9ToL1trTXftyIgbi
5Tw0cq+wYEOBY94PlmMgPJlSMnqeBROQPDNDYYoFiuvtMG2EBet/u8YfeTgW
VIXMTguJsWAakp6ZsW6xMOiBZLXPMWNqJpWGIyxxgBdn8jgQwIy9cLg7MaD2
AsyDYvb25DBhFqHSj0Jn46EbVpjupTNhHIz7TbJxCQAbxuvCSUwY++jmyMe9
CSBMERijPGPCZMP0NNyQBPhmkZ6zz4MJUwqadXsalwB2yiXK5hpMWMNGlfKN
s4ngtNBh+fsLIxbp2fdg400S+N6gL9pDZsTktnvBoj8FSnq0ng4SGbFKdFFA
5VcKfDsd7pTMyIhJ3hwZ4aFPhbOCe48eXWfAHGrcW6UFU4H0Wbzk8tddbmxf
n7iUCq8MdMriixmw2vP2lwJrUwE9G1VxSI8Ba+idjhpOToOFt0NxfzUZsMFq
2Ru6BWlwUIzfsxoYsPXAluXWd2kQSp9xUleeAVs7aJv1g5YG+tibKid+Bixh
fMCxgSsdxhU/vns7TcA0c0Ok/4alA05aCNV4RMAE+U80Bj7KAJ99auaPfAlY
bJbdEv2LDFjdsdpq8SRgpT9F9sTmZMBcf5ayjgMBy+/Xk57pyIDhezJV+hcJ
mC6v7UErciZUNKmWmvERsE6vw9wvX2fCHVPTDM9CeuyRfu/BeNpr+PwkPiR/
CI/dWYl24PLIgYFm4YIbPXiMu5EiOxqWA5072b172/DYgai31JqUHKhzreB7
WIPHflkz4MvbciD9+kD+lQw81sai3BL0Xy44HGbr3XLZPe8mKaSrIxc22h/y
6hPxWKHc+DRFIR+WCMyqDIx4TCyAS71WLx/mTj+zfLdDh9n29Nbvt8mH8bLk
vIOLdJghd5W/QHQ+NCW8U/n3iQ7Tft4sxL+QD+EOSzfSs+mwkV+Z00feFMB+
VvvcJXU67OXEXeXj2kVgXtcZLqJCh/nrPLLyti2CNIqsi74CHXZkyMx41L8I
xIfXTxQcosNW1dX2rr0rgqMFEU3WRDpMwP9k/wOZYlC/UD7WP4TD9DaUWNzZ
SuDW8+29pfY4jPWmQAs9QykUnLX+N2qJw4z2pbVvc5fC3HrrCPE6Dusmmbiy
iJUCxToqw/48DmMiXJs9r1kKPscPygrL4bCvfEFsscGl8F9h9Mv/VndQS0bc
kgSxDG4rSDEM799BFYyyinpEy4HCUJ3twb6DllpFaYQolIPXkJYOF34HFTzS
WnVRqxzC3G2faf/YRu1THH3Jt8sh720qX1XpNnrsxBkO+/Jy+KnAIxN7dhst
Jsocuaz1Bm6eIBhfdNlCnZIIws993oL1ydGcxsYNlENyTZdFuArsck4ana3c
QHtulsnNHq8CB74o+ra8DfTItTjGMe0qcNnQMO+O2kD9LiNZm65V4IPmkj9Z
bKDRoft4ONuqIEbL3XduYx21mmJQy/WohhZj4hW+4+uozJJH2c+xdyB1T2nH
KXMNdXNgkZH8XQvYlwtSE3FrqL9IbTY7Sx0YaVgbGj9eQ63b6QJIYnUQQArP
B6c1VELjd6ihWR18yvxizHFqDd13R/coa08dhPY9LCntWEWxKdo/nkoUpqRa
LZcWVlDuIcV1pi4M0r4YNN1XW0b9hj6dyKA2gNyk3Prg8WWU87M+LX2qAZpm
OI8dk1hGeVN/TxeuNMDUak/iFGkZxbsrWk7tawQpLj1Xg5El9NeTz/njxo1Q
dvaMiLTnEkpHf9tvfbgRmoqPPxwvWEStX16ef/+rCX7471XR4l1AWRR+5vcd
bYV3sRKa00a/UenbvoOpKZ1w7OenpaPJU+iPwCUPYbk+cPsrLckTNIW2yFof
XVDpg6pVf7Nthyn0quO+hAatPlBjlGrpPDWFNgk+l7Uw64PLol4vbg1+R5NP
lXA8Cu0DdzN+pYw931Ez5LuNy2gfvOu75rfPawI1ETatFwvth+2PReVbNyZQ
TzbcrfTn/aA+jp+e1JpA5yoTv/En90PHXJ5++d4J9Fn39ent0n4Y2bNxQL/0
GyrPwucfOdwPO5qJTY+nv6JYgf+bRukBOPNumGXLaBx1lfuj/6ZrAIryXiu2
Ko+jReXiDbofB4A30flWlPA4Si+mRPn8bQCoq4zPdQfH0OS/KOHL6gCIeRtx
JK+PosIrcxKuBwahwX+ZiGh9QSvtH641UQZh86k8Y/D4J1Qis5GOlW0Iep1b
bfnzPqFEb4EuzX1DkKF/vbnI9ROa/+m3t5vQEOhwBQbRGD+h5ywn/uXKDkF8
bD9eUvojymC8WX7XcAhOJN3d6bhPRR9+atI1ShgCl9ziNU7yAFpnXrIWxkuF
YySnQPXxfjTZemjJVogK85QjJNfCfvQ5BMzJilOBopIrOqjVj5555uaZJE8F
9570C3E+fegHOUv56vNUWJw5rK3Q2YMiO5/vb/pSwe1NjnmcXQ/qVsHCOhJE
hZUHYm6ruB70fMv8f/lhVFjjEEqpVuhG3Xs5/HjiqLB5gmtZNaUTTYmsrREv
oAJj4FrqWUobWkKI2SEOUSFEz70im9iGOvW4bH3+RAXmfQsdzNkf0FBlHsUX
o1RgzZ1daR9uRVMyZ4U+/6ACe8+o3kXNFnT5qO+bkH9UiHx53apkrBmVC+zM
bt6iApfVR09Or2Y0BVJ15+lowL3clz5Q1IQG3sFbc7LSgFegec2YtxEdLXjM
5cpHg/hJNVJ1WQO6YszZvl+QBvxFdaL85xvQkYuGb8uFaSCoVnVhxL8exSv0
/c45RAONbiPteiqKyhAptsXHaXDWUt8tcA+KGld3nKw/SYNzS9opmkgd+iKq
uB9VooEev+pyW24NajXxSzYQoYGRnWhav081+p4UaXhAlwZhd11Dn7+tQv+Z
2q88OU+DWp9G58tzlehGI1/7xEUaHIy2Pt1vUYE+Vn/5/doVGpi8Khd7/vIt
Ol70VcDLiAZPcunZLve9QU19tiyDTWiwgGYM96mXo9PsTVu2pjQQ71hsiPIu
Q59+kj2jaE6Dq1SNvEtvSlGtuDXvfxY0wH5N3OsTK0HvX1TPRKxpMP3nJiZS
W4RCHi2k1YYGFxmZn6XLFKLfpGbVle1o8J/8uSPpnHlos+QPiZ/2NHikPb1+
IDAHVWo8pS3qQNt9r4W2pS1loZt8w4d1btPAwEPi5QG712i4HXSaOdLgXfgH
uzRaBorLFhE3v0ODAxn2Cge001GxlggFXScahFYzE9LepaLr6NDKQWcajKZc
eq1m8gq1aUy0m93lws4odTXlRFT6Upf/KxcaeP/rG0P+e4lamGRqnr5LA21x
rgcIXSxaejG7uH2XeQwu8yOTUehdvdxaDQoNOnYOOvs9jUCd9rXfyd1l+6uv
c32vhKDmih51G7tMKBebfMDvj5qQbxYpudIA4dNTXcfc0Zfj6Sp2u/z3ZbzN
fLwVer/Jx953l33wMc1nT2nD+54Y2cBdlkmwbuU+4QxL9B9fuO3yKF/p3B52
H5C/8jDx8i4byH+sCVYLgq10blXBXa74ZBwQ4h0GlrX2DwZ37bNhs6JN1yLB
wiDZ7N4uQ3R+cJNnNPj8Fhpj3uW7drd/NsXEweuxPPrg3fulKx7Way6NByVF
za4/u/4YJM4WNXcnwdUy9dNnd5lxPI+zZTYF7sVpXH2867+kqwpK1kFpIFdw
/ED1rr8VNGJvT3BlgAX/xfihXT26jywnWadlgoUOb/Porl443Nsd67psEBLU
1qrc1Tf+J7fspF4uWBWdzAq+RQO5ATcrm+E8EIj81qa+Gw82WcebbdYKYYXw
3tXflgZbz2JWJoOLITW9g0DYjacOlfXosrAS2FEzd3G12l0vfq1Zjq8MZhSY
/vLuxuMxAdZRyngZiLgbcV40293P/m6lLLscorniD1Gu0yDuH5+E/Im3IIO9
0PczpkFr56cweYMqwP8Zca24sKvnvc+t1w9XQ7LE+VBhPRowHxwmBNG/gz90
+QFu2jTQ8R7xHSx/D1K3Wq991qBBp+QYxZ0bhcQlTpaU3fx1HBorTv6FQlVE
87uru/kdjLisWFlhsCzkgW4fo8GhXy7yT5TqoTXfeYdJateeBqVweKYBZpwj
BPkEaHAqfqHjj1EjaAmbSEjw7Nr/Q5mhb2yEuO5GVhEyDXITXA8dSWwCuXch
KdTdejQz75bmo9sCp3ckBDR365ljqmesYEEbbKpuErT7qaCyslouy9sOs/VS
7aQuKhD17vVrBrXDI8yuFGulQuHqPXYn0w7IzonW+FlDhbkLXo/riF2wc71j
r34WFap86+qKT/TCWTUpgzB3KqTmbndVO/ZC0X7FJkFnKoQOnv7SmN4Lp/Tu
/n5lTwUTKXSDRuqDOOdeqsN1KvwbQk/hfvRBvZ0CzxGECkrS9WX6LwZg9eXw
8zuMVFAYpMbVtQ0A8YHOvMPOEBzz/uV1ZHMArtS9Wbu6NgTibTwazJaDoMWf
Ysc2MwRkW6cBVGoIzjypDOjsHIJfSQJLMnVUkHmfotzwbAhSiB4KbN8/gZaS
V+0Z4hDgPfaO2PCNg9UDy7C07wNw7UZef9XxcfBc/H2fZXgAynSgjag/Dh8q
vjyx7x0AK6FbFW+Cx2Hg76t04vsBaPhQE0lYGof/DrM7/H02AIECNpqZ3V9B
weJ07JlTA0CoLyuaCJwAhBWRIAb1AxNR399qbgoM1B+487D3QeA4VYOS8Rs6
2ewSsy06oSQxWN93YR7s9+9kC9q3Qgc5+Ac7bgFIS/32wwatMPXkkW8a2wJM
Zs3ovoZW2O8TVNAksQCmOtbpZjytEGoawLzHYgFC02NfW7a2wA1BX+xFxwII
v9WfrZNsAY4U92OlmYswGLFzAL/WBM7pNhyTRssgFPrCNPN9Azy9H5/01HoZ
7kczGtTkN0CRfreEossyWGTENowlNsDc9km1p4+XwTz5Nou2TwM4Xt9DOVWz
DFqnD/G8Vm0AB+6ygfADK9DXT573b6gHu+CdOIVfK/CTrCaHDWJg7pAgGOK/
Bjl/dGxBoQ42PrF0z0esgX2EvkaoQB281L7va5qwBsTCmxnf6eqgT8JkTLZ8
DbS1SE1fumvhzI+9aSOTa7C2R8W74VYtSNpEHZQ79w/2uZhJGWbUwJJZ6OEv
pHX4YDrWaiz2Hh5f8lCUT9oAjdRLsXesquB5p5RjQe4GPLiWaFCrWwVJWmOv
xCo3QGz659YBhSooUT1H4O3fgA3Ozk+SzFXwUZK/Z5NpE/aqKxLkiitBAl9n
0+K6CdlCgs376CqhpYwh8qruFjxD8mJayt8CgTt2ym99Gy694pG0vVQOXHO+
JC6mHbBw+c88DymHAx/sT2SSd0DB9M/97aPlcNpHJfjDkR34OvxeYYBUDvcn
J8Q4LXYgUinyjnRXGcy9kbNNb9yBtxw4hwS9MvglNX3X6iwOUQjs3KQ3KoUl
ni+P8vVwyMfkCjZ2rVLYIPTHL13GIf1ihhGHTpUC69j7+mALHHJOXDj/IX8p
HIp+xpl/D4ccNzMS7R8rAauNE6WLeThElqFbetqxBGgdj/4GsdMhYzF2t8UP
FEPM6WRNw310CHUznHiUUAyXS98kiAnSIVxIAt2VH0XQHTdxpkWKDjF/8Np7
obAIWizV4pnO0iFnGc7X+6oUQcXqplqYDx1yHTt2o9a8EOJE3aKf/aRDlEZP
6gtW5oOht9WJhEY8sqibF7JRmgNxydyMq214xKE+tLE3KQeoaMuQQS8eMb9k
erU6JAeMCIfd2b7gEc8pdt0msxwwiVgsf7iCR9ZGhQxQ1hy4lhIke1OSHtms
HXusZJ8NNxqzjhx/Ro98Syw82nUsCxz3zIp0XyMgeEvXUz1LGTB7/x1hjyUB
6ZqaptedzACH6cdTWjcJSN7shaChgQy42SyRV+9KQAYkPO/sK88AKz872bfh
BCRNNsfuLCUDri5+haQ6AqLbaallsJgOZ4c/mjqIMCCM3O5NzTtpIJTfEsc4
w4CQ/vCeuA6pcCWlqrr9LwPiLZaLx8umQlh03sjTVQak/qlUUIVIKix5PxXh
YWBERt/UOp1jTIV2PaOiQ8KMSFpPRdP+pBTwmPvefM6IEeko652Uxl5BlyzD
clgDI6JMuTf5QCAZvKs0DNgTmRBmtc/37vAlgJxp71+LNCZk1fGozgVCAszs
mD4tyWZC/HJf5Kn9iYdrWh5tl94wIQKN/AMWzbt9npoDMV1MiKLX9JgnJR7W
FtmO8OGYkZ6+34zMvS/B4+hHgqgtM5Izhu89/vIFuGQ5Vpw4yoKYlKsaGF6N
hdEegVjSCRbk4nIQ59+zsXD+X7vrlAoLItrRQEo6HgtSepKycTosiHhx18ND
HLEwOf89f9mWBTlf0GXJ+iEGjFUt0t4ksSDSPLNa80oxcHrwUrgcKytSgxOp
7z0SDUT8SSuZ76wInoFl2101Cqo96Q2mZlkRM4XIPPGjUWD3u0fj1QIrMsZ9
/tykcBRgH+3FSLg9yHcCg1gwIQpcixOmf/PvQUTWnzFYd0bCZ9Ntp0L9PYhe
xsFXHuaRkFPZ7H2kZg/yTbNC0prnGWjeuRIrFU1E7nKLetBxhIPb8/c453gi
YinaP1qw/gQyK0Ucy1OISO4l75Vb358AHj+vrlJARDxSfRZ53z8BLC7i7/lm
IrL/RLvjuZtPQAlr0b27SkQUDg/1jTWEgQxZEV91nQ1JKfrCEhL8GHirBZ01
D5IQ1YCh6niVENi2zhbgPExCOhJOardKhsAkSfbDiCwJkfkp78PAEwIlNmcO
uAEJQSXcIsrng0GL4/ZA+jUS4qn/O2MwOxg8bladxEWREPtj9ERBnmAY2muA
e79FQhbW5a+Y0T2C55TQqGM0dsTC5UnGQ7pAWO8sqAr9wo5It7AJ2ywGgNWh
vrHxCXbky/N8NoPvASA/wicT9ZcdOaYjzG3QFgAz6Ybf1wgciCKXXuXnqAAI
J/7wv/ofB7Ilx4vtHAyA7jGW9/yXOBAT1h28uYE/XA6+KJ1cwYF0jKj2fE3y
g6yrf7WrqzmQr79CYtSD/eDfkSi7oRoOZKqQkfWtsx+kDvSnsDVyIDOVg23d
6n7wW9iIy7eHA4lbZN0a+OkLoe9NV82nOZDKzIgbvMq+UPf3Vr0wPydivNsj
xWZ8QOp6kGHmA04k9oqX51VPL7h6t4KL4M+JRJBUZA5ae0FoyHSPdRAncvCv
Z+POBS/4Xq6nI/qEE1GyvKf3U9wLUol7If0lJ2LRGy4bRLsPe+teS6SWcyIf
iVrQqXwfcAdaNxJ/ciI6NQQiwnEPqFMs6dGGXMiYHENv5Xd3aLnyOKnahAsJ
X+9qsKa6Q0UD84ux61wI03AbRaDVHeJeMYUftuJCuss8Jgpz3MHQiMGj0YkL
2Tyu+8/M0R0Gm3E6iyFcSEEQT+q3ZTfoz1ydv/yOC4n9yvlZhdMNuq2+q3EJ
kZGM+jRxrZsU+ErIXosVISNsrCnnRa5SYCnLvphPnIzM9EpfZ9GlAN/s7P4D
0mREPKBwe12GAjauC6tHlcnImqim0u3Vu7AeuFN43piMdI9GnvYIuwvir/n4
Hz8jI3lH8m/OVbmAotZwLzGajFQFRhySK3ABvZ9JIZFxZORz0OfghykucFda
ePlFMhlxkzPbUgl2gZq34r1Z+WSELygXlb/iApdb5IObWslIoQV/m/CkM9ja
L6uc6yAj9/a3Wix2OMM91srFjm4ycouJyvq53BmSLyhZDgyREcz2j+FEoDNM
UxGVbxNkBEdf2j5x0BnW7+EXbX+Qka81MSRDNmdgE2jK/TlDRnizzn78tOwE
chZaPPPzZMTXkKNSsNUJNOlYut2Wycji4sXckWInMMloD1pb252/WDhc9NIJ
HDTDlX02ycjgiSmVGH8n+L//POT///P+By9Zb/4=
"]]},
{Hue[0.1421359549995791, 0.6, 0.6], LineBox[CompressedData["
1:eJw1mnk0V1+0wM3z+P1WQoWEVIZSUco+yFREMlOZlWQeIxSZlVlkHjKVuQzF
vRlD5ukrQyqhxC+zjO/71nrvn3vXZ9297z5n733O3metI2DuoG1FRUFBIUd+
/O/70YzG0t4eEW8/a/LvsYAOTHsfWB3dJuIc8J75tpQOXOf4ul71j4hPRg+U
XFPSAf4LjjsOy0T83lbwhvE9HWgIi2GYnCbiZm2/07XLdYBBfPgI3knE35YE
iReo6EKcq+k1/5dEPFkvXDzYXw92GESvm7wg4nhMWNL7OD2wTl3Ukokn4pH4
f7w0BXog3fJE7+8zIp6+/aWut1cPvnDlmpk+IeIRfTBmdkwf+N//8kD3yPrV
rvT3OvWhiNIlh0KaiDs6b6ZVShiCzM1CuiQpIi4a1zj1VNUQWnK+3T0tScRF
xqzPO5sZwqSKlpi5KBHvjFF/lBJrCAeeib1r5CXinNNTR2X/GcITntnWIAoi
fqf5d1pdmxEYSN2aY2kn4Ci3Nr3F3wSmA+M0clsI+EMuHv+LqSbgMtRRcrmR
gPNVDtd11pjAM6+LrvYfCLhA9q0Z7WUTaMS4drqLCbi0edKBMZtbIK7exxYb
S8AzZq/lMuvfBlorlTPctwj4m6fcU94iptB3+7T+nCEBd80ROfThvCmkGfD6
1OkR8OekSikWZVM4r/5fs5kWAd8Qa1cdszQFa6lEw0JFAp7Htkstl20Knyhn
/GVPEPDHGR+7uATNICItqPvWOie+dXa7q0DCHAxeOK5IrHDiVgbDMkzy5nAs
xoibapETjx1onvXRNoe6p+IWr35z4kzzIw+S3MxhwW5obWGcE5/LRfplH8xB
S1b4iH8TJx7foKnhcN0CiKRmu6xoTny4c16n3t8S8JOZIYHPOPGa93oksThL
sPPzybEO58TzmFKTyvIsoUVYauzEU078IV2i5FqXJTx0y7hW4cmJF4dunHY/
YgXfCd4nmu5w4r/Gpn5ZfbSCcg3J2Z9inDidtbDkGy4buJPJTPPpBCfeGsvm
3SdmAyyr03yFIpy4cBcfB/sVG7BOSdW3FyB/H6N90uVoAzxzTK1r+zhxWc53
dvc6bOBJyM9c+m0OfNSkikbk6V3Qbky2EG3nwNksNyWeMtvCsYGpidYWDtwp
Nit6QdgW1qbEjawbOfBvT4ODPBRsIYmuUSv7Awf+Y2K7edbLFibV5i4fLuHA
4+UMfqv/toUHPbIHOeM48CRBnxm53vsQMj72eeMWBz4fWhI6X/0AjBaEVRON
yPJp36pfkB7AqT3HxnP6HPjtAhr87sYD6OanrXXW4sDfUBizesrYwwFLsbw/
Chx4LmxXGNTaQ/bvR48nRThwM4OLbteaHcD7v8t2UT/YcUt1qbFX246QYiY9
i8bZ8Z7rWqv1rE7woV/ScnGIHe+8ntL+64gT7FYJmmi3s+M3H8tSeso7wRN/
RvV9Zex4XecHBq9gJwjjHDz5wpcdP3bT9oQAlzMkS9nNpXOz47Kiz9+9VnKB
97lWNlpEdvyu00SNqqELjHHd+UHByo5TUrDabtq5wJHtG6NmFOy4ymH8fHqc
C2Q1SXcIzrDh1QUqWotTLlCgS12UV8mGt3CumHaFukK1R7JtiSYbrnlszVf0
pxsctXTKJamy4UYCXJExO24QqaU6SanAhhOPTR3gPuAO5qJrujpn2XCPNsGz
firuwDx6Q/7fQTY8/Zd7fUWRO9yWYzyo+J0Vr/6rSTnp5QHUtJ7NQ66suPMt
NqM74l5gv3idgsKeFSc2hCc1XfUC0riQrKgNK/5wIdJO3cYLXr/rL/U2ZMWF
KW19mzK8QO+uRAq/HCs+RX1W2nj/Q8jvmHa2pWfFeU3Y/lyl9YbrMbr8u0ks
uJpl0A0/ykfw8Ay7mkosmVeTwtYEHsGrvk9OUREsONa4NRan8Ah2CZcaj/qx
4PM7CzEqgY/gTexRaxVLFvzL7uV3YYy+wBL/X1GUOAseNVIkX8nrB+2JodKC
jcy45r6AWKpSf1iTVjS1+8CMi7pKLX5r8IejpO2Qt2+ZcTVPAwnSoD88POj4
RSWfGTetDMYZt/1BNEnvkV0kM46Fi9sEqj2GkGTBxrf6zHhfEV3z4uxjUEmt
u676hwkX0r0cG3gxAP5Gjm3TTDPhOZ/rmCg0AyDZd6vw41cmvKnIsSPGIgD+
3LlAL9vPhHdP0CTsRgRAzNG3mFgtE75mzqC3NBkAE3mvJYkhTDi7wfAhi8hA
cC9/SZgQZMIfR/o5nfv3FPiza/Hkw0x4QO1SVC9HELTFjtjrczHhfM+yckOO
B8EhN66OHiYmPF1PQ0zVIAgapGMDGhcZcXnNL1IcNUHAVhe2ko8x4iNjR4In
/YLhVavXkIsxIy4Z+UXzgWAoJJ1xVyzQYcRdqymyV1EoRKQ6l05oMOK6Zz7n
Zd0OBWeX+2FqiBFvlnbwv58cCpf5bsnxCTHigokC0XXEMBhwR6/aFxjwUAF/
+kjWcKASpncXeMKA6+nuHjvLFwmrUdQ/9LwZcPq084Jl8pEwu7WnGeHKgBsU
FewztIyErt4N0XVrBjyjKa7gRGEkJD36Pf75GgNeLrcQ+0X6GUgOdip57mfA
/4h60KQZP4fbgXH7u/Ppce5khbmjP6KgC9boPbPocQ7XcKa/q1EAW/qb/Cn0
ePH5OYZ+hmjgd+b96vycHhelqWMYFY+G73ey8g+40+MFPKo/Vr2jwVq2VPa2
Ij0es9re/4AnBuyXOszmx+nwvXnG30fvxcLEGzGdhGE6PG54+dpX31jQvPdc
GXrpcCu6w4uVcbEgOal9MrqJDj/Mdt3rNR4LS50jK2eL6HDtn9mV3dxx4F4w
E+zjSYerUhkPUfTGga8pdTEzkQ6PmqLhmLqZAKXdKs8GWOhwu0m7Z4F2CfBd
LsI+lY4OdywpY5B/mgDKh/dLSGzS4tQOP+3YqhKA7YtwqfY3Wty2hOnOLm8i
pN28Wp5UQovnVjXkXp9LBEw5+p2IOi0eyRE5HJeRBEtvBxP+KtHi+3ulRZY/
JMExIR6PGqDF7w3diXQaSYIQ6mzpa1K0eI114vFuQjJo4ZXV9jy0eEzUCyfO
4GSYvECqfTtLg+e8uf2g1+0lUIjxYYpPaXBkfzFK7kEq+ByQv/3Ulwbvy/7V
+yYqFdb3zHdaPGjw1r9PlFUrU2Gh75XsVVsa/PQTwvO9rVQY9RSv1tKkwePV
W/c/ikyDd02Xy25x0+BphyI0JqvT4YGJSbbHG2q8eHAryZA7E2aVfBVqXlHj
fU/iT6yJZYKlRMa3zXRqvMrwQ1KRQiYYUU3x+cZQ4xITQk/M7DJBNf9+SqAH
Nb6BWTzdxjPh2MrD+Gh5alwUaOfpHbLgS3hScNEgFX6twM3m8tds6G/mf23a
TYUv1xlsRq9mw+e9vJ79bVS4Wkm4/h5zDtS7vOP2/0CFu7pL9h+7kANZxv1F
OtlU+OmKxSepcTlge5K1Z8eRCleSvhRtrJULW+3+B7VYqHD26Buf3o++ghUa
hsu0dFS45vEjhkarr2BB7rlZ7R4lftjp6vd97HkwWZ5aeGyZEpfuZ1eZUMiD
puTaS/9GKHF7S3655qI8iLBdMc3Ko8S/1vMd33qSD4eY7hasKFDiZuLUFo0q
hXC7/nPE0UuUuCYp3zPCphAynU87ap2jxOdLDlk9DC4E4dHN869FKPHq86Ka
za2FIPE6ssmChRK3uX1BkfdqEShcr/jaN0iB0ycP5U7rvoZ7Mbv7y+5S4Hl/
SPhKWDG8Vrb4N2FGgfu5CTxKLSmGhc3WMRZjCnxNtzfQYaAYnC2is+9qUOAF
OQfooo6UgM/ZY6f5z1DgJxn/1bVUlsCRN7EvjqzvYd+19p+Pdy2F++dO0I4e
2sOm/uNJMK8pA2famjx39j1s7aNLJnNXGTwcVLlKoNrDCPll1e3fyyDMzeq5
2swuxqcgrRHOUg6FbzO4q8t2sahIzXEus3L4dY5LPF55F6Pw9lizY6sAm/M0
+pqOO1jZSPGSx8NKsJCeyG9s3MJuGhnS2r2sAut8aT3lqi3ssMy378qVVWDL
HU3dVriFRR0SlpDprALHLcXbXdFbWLOy9KXbe1XggxUQR+5sYcf0Zy/XWFZD
nIqb78LWJnYu6QWqkKqBFn0WHe6zm9he83NY/V4LJzwv7tnnbGCCEYXLCiH1
gI9fP/EjYQOTZpx8FlNQD3qKFrr6oRtY0f7DT3bb6+EJW0QR2JP5k6L3STYM
RnLG9TlkNjAIplkdiccgpNe/tKxjHaMk/WNT0sVh+kSr2crSGvbsVvSpJNeP
kDl+s8lLfhXbfJ8+qDjeCDOP919SObiEzXbr+WxOf4La+ONKs3rzGHPIqf7M
pS6Q/DWyIpE6jR3b8zus/aQPXP+KiXIFTmNPx0Pv3Irsg+r1x7d2bacxh8aw
T5Yv+kCe7kTLZ5lpzEdyvs6muA+0BR8m3hv4ie2j+bUh+6UP3G7xXMxm/okF
iH5LZz7dD7W9Rn4HHv7AzlbmSrz90g+7pOKKHdMfmMmLtyyHf/aDwiTV7JTK
D8yUY9jJ979+6Fgo1KrY/wNrJ2X7itEMwBjzloBW2Xfsr+mu2kWxAdhTetkU
OvsNM4qcHtt6NABXakcZd/QmsWY5FBHCPQjFhbkXWmUnsfr7y58bBAfh4EuH
e9H8k5iAQ9jVJbFBGFqni7k28BVr5RGPklQYBCFvPY7UzQmM69x+bRHbQWh4
vMqCVMYx4QgRuemqQdh+JkUXNDmCFRRj6otqQ9Dj0GrFUziCLb6UbkrXHoJs
LePmYpcRbCF/SwoZD8FVQkDgMN0IJvlHP0nPbgiS4vuoRMVIGG/zTsdQ5BCc
T3Ha6/AawiIrLN+xdw2BY0HJBiexH6vzyfT/qDAMkmz2AQqTfVjTzxcX2NWG
YdH5FJvLmz7MTZvu0nXNYXC+VCA4oNKHyfJRjaYbD4Nbd9b1BJ9ebOUmS81z
l2FY/n1S7dznbuzNg4eCcRnD4FqZfzvBmsykm+sLucOw9kjIdZ2iG2PeseK9
WDQMGxx86TXnurDZqimezLfDsH2esHo5/TNWJRNvX9o2DHQBGxnKzm2YU7Zw
7qeFYQhWd3uXx9KGRSr0HDi6PAwMB5Y6GPI+YR1rvrT314eBqWBurX20FVNW
bD7RsTcM7N0T6ppKLZhRnaFOIzsJol4Ym5d+bca0PWzsEokkIJiTPDgfNmOD
BxLzjLhIsG+1N6u/uAnTUCu9VXyEBAd5mzf0DzZifVIrxaGnSJA0Jc9WU96A
vfpH3TwoQQKe4npBHo0GzPm5ZwmHFAkOy1dfH3v8ERuoSRrXvkACxS49tY9D
GHZ8SOYNvRIJlM20XAOYMay7JN+lTYUEqitq6UqoHruYpNvhdZUE6jyXV9sK
PmBejw04XmmSQM9aMLPPpwZTmhcVPWVIgjAnl5CYt9UYPbEL1zMmQZ1Po4P2
QhUWrqVh7niLBMdiLeT67rzDtt3kXjuakcAgrUIo5sVbzHuwT03fggThBdSs
2r2V2LxH8OgpKxIsYdmjvQoV2N+wyNacuyQQ7lhuiPYuxxq/OxxUsSWB4ZBi
4Y3KMozPP+4q6T4J8D8/PHuFSrEexgb5VnsSzP5ngx+tK8ZUozt3jzmSQJOO
4XmW+BtsA+ePdHIiwREp1VNZnIXY1bks9mEXEjxVm90UCMjHEjdVd/5zJZH7
0ZC2zJVXGA+FUMm6Gwluuh9/IWCdi809tjm64E6C2ohP1pnD2ZjKgoNOvwcJ
BLLvnhNQy8KeFX5WyPMkQUgNA01mbQZ2jl927p4XCSbSb+TKG6RhVu1513gf
kuDN52gFedmXmBESN39PZu9/vV/RkRdYCuIWu+pNAjVhwiNEGY/F1H4vbCUz
101tHjQVjZUMSX2R8iFBx94xB79nkVj/zlTNMzLfNcwt8NUJxuInd1RJZKap
EJp6xPMYs/R+G8TxiASIW/3yJu6GVfM6O0uT+e+LJMvFJHOs0KKNWYPMPlRx
zcoyarBz1V/7BpnFky1a9513ANWj4cpXyDzBXbbAzO4DjtGss0JkvilF+hAk
HwhKL4yk18j23o3oPwn2DgOv8JWzb8nMis8JNhlFQRQh8ps5mSG2KKjJIxYo
j9HK7ZHn42R9/1dTXAK4Ujy4FkbmrAsn1ZvLksD5JA0zDZkHWOaKm7tSoMH6
gq8d2T90k4WcLXPpYJJf+KqR7M8Uw3MXLQIzIZ0gFMRI5nOK8fd/ELIh6lfO
ETmy/7tOraZYZObAoxBdZzNyfCgo3u5Z1OeBb8U3HVdyPJN+7Ts9pV4Asmu7
sxbkeJ/pdzW3HC0Erxs9CgrkfLB8dbbZcuMNEPgP6HSQ82fnedzaVFAJNLco
ENwdyP6/tBlbHlYKLK2x2ezkfLMUNmo+w10OH58QNdnJ+SjJyzThPFkOqi/s
ndzukfXZa9fK8yrA8JmyR7sNCRL+cR+XOv8WqjQaQN6SBK2fR8KkblbDXOC1
6+Ym5Hh6fmk1PlkDMj1tKxxGJGA4NkoTSF0LWnnsBUX6JLjqPeY7UPEeoKjZ
MFebBJ9Fvzq77cNgJVzpYyF5/doNfi1J/YMBPSqvqb1CgiDkuGZujgOj8ih1
lTwJRP44SoVf/AgSHn8qXGXJ9hSd34z+boCVJtY0WnESyCQtdfyn1wj0qeWX
vU6Q7f/n/Ju6sRFOr1x5MiJMgoJkF5FTL5tALe/scVM+EvxedM30udYCORWE
88Hk/cwuwyP+8Os2eOfgGqHyZxgura1XnD7YDqeCzfXCZ4aBRd2zTymwHSY0
jid++D4Mb9Y92e1NOqDdU7lmhjQMC9cfhtazdMIr93IJ86ZhqPatry853wPK
Wrxsa4nDkFGw21lj1wOrnIG1B2OGIWRAbrwxqweIZxV+ikYMg8EJbGuYrRcI
Esa7XP7D8G8Qk6GY6YXPHA5KQjbDcFHsY7lWYj/IBEou7DszDOcGhhLq2/ph
+c959fST5Hrj/efhqe1+8LGp9uQVGgbhNi5FBrMBeK+34P+LaxiIVvb92IlB
eOu0ID6xPQR/UnhXxOuHwNWf3+pC0xCks7ifY/05AnNGvzFpcn1Mrgjn9ub6
ApqSBk1H5IcgwShzZ1btCyxJZc38kxmCiPzPzU3FX2BPuO5LxPEh8LwiqOfj
OQqeuVmvROiHQNunx32OeRwKNwQKFBoGgcp9/5gl9yT8u1m8ki45CEamhX3V
ZyfhiMlf2rvCg1B+FdpYtCah5T7POaFDg2DOd+9dZdAkdHauGPnRk+v5pw9R
NCuTUKcm3/p2fAACeC2Vcrq+gbCgi7piyADQfCwv/hHwA14mO+7PHuoHehat
x+YL0/A+5rjvo9t9EDA5pOicPQ/UM3bWF990QenLIC3fpUWgrzlayy7SBg5Z
lhxTeqtghj8dfFLXCM+8klKeWayCswztu/bCRijW6jp+wXEVsKpzlkdeNMLC
rrT8s9BVSGXqVfjl3Ah2xszOMh9WQXiFtoxVtBFs95X3RwisAblAHEpIbADr
oL2Ec3/W4OeZJLemRx/htm3y4eDHG8B855eN9kQ9bI0wdi1GbkD/08rGdKwe
Xqh5+Zokb8Af2seLNJn10Hvc4Ovpig0YkbSpPWBZD1dm9meOTW0AyUnpvf5c
HYhaRh87o/oPuDwPNarvfoCVWyEnx9k24TK9hOnPc+8hqnN9VIV3E9aflo4K
8bwHscs2EeUim2D8WTP50W4tWB9Sng9GmyBISdHi3FoLI1+oS844b0J4qvg+
F6NaqNf3OxMyuAlaTKvE5sAaCL3hfkEqZQsWJM5ec/tdBTGfT9i9LtgCybup
6U7dVZCi8jVNqGoLZL23eUPI/XjpZVWag31bEO/baUjtVwUkUZ7ubfpt6DWY
9Tu9vwqOU9Vbtrhsw2iacrep8jtoKaeNMry2Q14HdsoF7yuhR6y2oc9gBwZm
rHl6sivhS7796jXrHdAflLzBHlEJC2nDRnKPd8D/1MWCIZNKOBBeICRYtQOL
pzs1gvcqwMpC4/2fo7vg4KIbaa5WATT74qf9NnfhdLuN2fOFMiAs+LIR6PdA
f348hnG0DAQ+3T2fQ9yDUD3/i4mtZSDncyno06k9CKzUffsnowy8pn4Icd7Z
A54032FfnTJYqDxjldW4B2KZF8uksFL4c2LWyVyZAikN7lPfvlICK1zjT4vU
KVDQUtj9F0dLYIumL2lFmwIhQ/YtfYoSYPr6/mPQHQpUu9P7D30oBpHY55xF
nhSoWfCkRdW5YjDfOl+2XEiBIrQCNA9LvIHhjqd/A9kpUcmutSmjbBHEyaUq
6R6gRKpnuA11eYpAu6wyWegwJerTtSe1/CuEroQfV1pOUKLfeuwT3DWF0GIm
n0SvTJaPwk75XCiEd+vb8mE+lIjgf0KGKF8ACYKusc9/USIbtoH4F7Z5oJMQ
PnvnLyVqQuv/blzPAwJj9mXJdUo075QZd+pMHjxb6J3poaFCF1zvzFzZfAVB
NRKXCPxU6Mceq/Fu+Ctw15ybitOjQhnrKl4eb3NB19v8fHIjFVq6La+RxpcD
Can76NbbqNCa5ssToww5MIS1DN7soUKtzXzup5eyQY/mpBvrOBX6+Ft8SbQ5
Gwwilyv816hQQEbiexW7bDBKDzxtI0qNTNMelPNjWWDa+OrU2efUyEVDLTLF
PRMyfhpsR8VToxzmkzBsngnf6Jk/z7+kRjoKq9rHNTPBXN3hfl4+NSoliSuy
Hs8Ei0HpQp4GaiT+17P+65cMsJ75JEK5So3st3/yxl7JADvmuaNdRjTIXSxv
QEwzDea8ammYzWjQU76crudiaWA7GzqtYkODUoQMvbhY0sCm+XjhRxcaFHWY
VyGxLRXM/axPv40g67N5aukop4Lh8jdIqadBh6VIKlnKKaA8SjKxPUqL5rWn
IwPMkqFZLV8u7zgtctHM5r1xJRmuVHvwT4nTIh2pnkwQSQb5+ANTt2RpUVPe
HIPffBJc0tS5f0OHFmn/U7EI8E6CM43dXjJBtOjKFfHagNQXwFfUkkD3mxYF
WE9suG4kgE56dU37X1rE60A1sDmWAGGxhWPP1mlRsDHJL+tjAqx4PzvKRUuH
6lKwPL3wBGhX1ysW4adDDxiqBjT4E8B94Wezqh4dYu8SZtPViofO07SrYQ10
aCe3oKW2NRaohNe5NNvokNp3PjOJ0liQ4fl1kdhDhzTp/zJ/fBELWVSd/i/H
yfJTHLSytrHg2R/H8nqDDnH6rRY/ZYuFY27HjnWK06NMaafX+iYx4F2teJP9
JT1aT/xqScMSDWdMev7eyaRHMyYiXyT/RcHvPZNnpXn06JhoGJ3LdBQYqbi3
3aikR08nEvhlPkbBxaF8iOukRyfummQ3eETBxjLrKW4KBkQ1d581uPc5uEuQ
aAStGNBSyAdy6x8J4v2WWS73GdBtq9Qt+thImHZfhCYnBqToNCCe6BwJuvVM
3la+DKiSLiNR6EwknNW4vJiXyIC0Hua/2ymPgGXbrPFT7QxonxTdtNv7cHB8
ZffuvAQjqr7KcN9qOhQmunnj2c4zoqsN/WiwNRQ0/rW7TF9iRJIPf4nfLwiF
E+qipxOuMiKBrM4YhgehMLX4s2jVihHN9V+Xo1kLAf3LdzIrUxiRoWMw4xBb
CMgN3Ig4w8SEFG+kyLwyCoLXOxT3mTiZUMs3syY/FAS8IqVq37mYkIUO+bwk
HAQbXmz0MUJMyCBBO7dr6SmU83U8XkRMSEIGZ1uPeApC9694lHowoQPHCoL4
WgOBhUraXPwnE5r52JmD3QiAGg/qm9NzTOjdMqdRk2wAWM93K6YtMSFYHKz7
IRQAOOmuEBsFM9qfyc3svPkEXEqSZ+d5mNFkmsrhtpwn8MVk1/6NFjPyz+N6
fWX3MeRXNXuf+sCMKNnP/Itv9AddsRi7qQZm1JVe97q81B+osm/fSmljRjKC
4aemU/3B5Nm6HMswM3qzEnUk2tMf2K1Eqf4sMiOv6MyKZXF/cCdEhhSJsKA2
KsbMV+AHSg904k/EsqDSi+b5j48+AteY9xQOSSxoqDv0+2maR5BTddSuIp0F
XaKMo6OZ9gEqqkWFS69ZkEAd9oC10AfwhMi/Gs0siH6v1q5Uygcu4i3XnNZZ
kPHDUuv9Gt4gTrxAVW3Miu7FXlDVS/SCWzIZD7bNWNGE8slrmz5eEHGLfgTd
ZUW2OsSk9+ZeMJc3WNLmyoqYqQ5Hxot7Qf4l51ujkazIGNM6mt3qCYLWRdW7
GCv6JEI3Tb/nAQdrDjsoHWNDp2W/Lxx44g67Fnm8nCfZ0H88exTS99xhiu30
p7HTbEjT3r/PVcsdSi2vCLgCG3Iae39Yjc8dVDju92cZsaGGJ+tKWvVu4G5T
LU0RzYauh72kqKF0g8H9Nyne77Ch8/J6r3dTXKAWH3sdRMOOJBm6faKDXSDj
vrWhNjM7EjpLeHvZ2QXsPnqV/zrIjvaNe41/U3EBmgcZlgfPsiPbmFvJuivO
cLZp/pObLTtq+NrxpvuGM8Q4h0RLDrOjII0HMZePOMHm59fVIePsKHjjF/02
ixOYi/R+nfzBjmKX6ew6txxBaoxbPPovO9LZ+JGXP+IIv7N0f27QcKCyt1fD
5BIcIYJl5rHhEQ6kUdIiOUhwhK6vjO95bnCgar/1mQtn7IGKNaVLWYcDkdJ4
u/j32cP5i+LfnfU5kO5ks/3xtQeQFqfN2HGLA5VoUKkm1T4AB7WX+t62HEju
mr3czpUHwFFxamU0gAMphlxefHPHDrSDNMVS33Ggqpzj+2qKbOGV4V+1mhoO
pKMmzrUabQv/TkVbD37gQGeLT95Q97SFjP6+dNZGDjQ/fua/a0q2MM+vR/Dt
5kB3Ukz9IibvQch7k/XbsxzoxvWkpg3ee1D/995Hfh5OFMXyM60i1wYITUwT
lw5zok8KAgmHn9uAVWLRpgE/JzIOi5NI97QBFrl5qWhhTmTIJ/p19poNGIY7
vqKU4kTe5vu1Y5etYUXII/z7NU50tG5sWUzFGk4YB+rmPOJE0i1J+R27lmDo
9I5A85gTpZDW8q7+soSQ4Nlui0BOdGqru/9bvyX8rFC/KhhO/p+K/CWDAkvI
YNkPWS840VK9ooi5riXsr889nlHBiTxe6iX0llkAhUDr1stfnOhAYu6LO97m
ICH9r3rrDyf6qKx9pfiuOdzWOOlu/JcTpekfr+XUM4f3Xs//8qxzIl1ia8aB
0+bg3qf/M4mGgCSv/c7bP2sGcwGzXYl8BPTIaMul39AMhqYZs2J1Ccj90dJw
4TVTaNEJTakxIKAV0ibp4SVTeNfAkPjVmIBi6zqcjMRMISGNPuKkOQHlfK0Q
UmU3BV09WvdGewLKPUnrUJV3BwaaKa4uBxOQHa1MYgF2G/py1he1awnIkp7u
+wXaW9BA8PjjWUdAdT9ye/x+m0C5/9p0Gk5ATBkeM5PdJhBjsjr6u4WAGARs
Dqwkm8CNfcvNT/oJiFqdWc74jAn0BCwklf8hj3fikhrR2hi6zH/KE/iISKb5
U5PVpCF8o8nbiD9KRDME56+HWwxh5dXdEm5hIgrK8/JYLjIE7rm5QwJiRGTT
ZRK15mEIli5L6xKyRLRxOz3QgNMQNgP23mjoE1F+VAjtAzUDEM7l5gl9TkSn
XKm/zXfqwQWV0R6WWCL6+/lEJLFKD9R/pQRHJRBRVb1TlHaGHjiJ8a8mphJR
w3fvo4dc9ODDW+GeV0VEtP5rWOQLtx5ot0gFNbUSkZPW5Z6Z+7pgdXf1kmoH
Ea16Rq7G6emCJ1PVckcXEb3cbn5pKK8LqdcvmvUPEpHHGY/N41y6MDuELn3/
QUSbAxarY406sOlJtWw1Q0SmFo/3K5boACtvU8Gv30TUsSE1giXrwJk7KlyL
i+T5i179uu2kA0qUjF2uq0S0Pb499+GWDhhktwdubBDR86ffJWPVdMBWKULW
Z5uIXjN7xvqd04H/u8+D/v8+z/8A+6cm+A==
"]]}},
AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
Axes->True,
AxesOrigin->{0, 0},
PlotRange->{{-5, 5}, {0., 0.3826487832713353}},
PlotRangeClipping->True,
PlotRangePadding->{
Scaled[0.02],
Scaled[0.02]}]], "Output",
CellChangeTimes->{3.512237368859375*^9, 3.5690067890478086`*^9}]
}, Open ]],
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{"(*", " ",
RowBox[{
RowBox[{"Constant", " ",
SubscriptBox["K", "g"], " ", "for", " ", "the", " ", "asymptotic", " ",
"approximation", " ", "to", " ", "the", " ", "Welch", " ", "t"}], "-",
RowBox[{
RowBox[{"statistic", ".", " ", "\[IndentingNewLine]", "Zholud"}], " ",
RowBox[{"D", ".", "S", ".", " ",
RowBox[{"(", "2013", ")"}]}]}]}],
" ", "*)"}],
"\[IndentingNewLine]", "\[IndentingNewLine]",
RowBox[{
RowBox[{"KgWelch", "[",
RowBox[{"n1_", ",", "n2_", ",", "\[Sigma]1_", ",", "\[Sigma]2_"}], "]"}],
":=",
RowBox[{"Module", "[", " ",
RowBox[{
RowBox[{"{",
RowBox[{
"n", ",", "a", ",", "b", ",", "\[Sigma]", ",", "\[Alpha]", ",",
"\[Beta]", ",", "w"}], "}"}], ",", "\[IndentingNewLine]",
"\[IndentingNewLine]", "\t\t",
RowBox[{
RowBox[{
SubscriptBox["n", "1"], "=", "n1"}], ";",
" ",
RowBox[{
SubscriptBox["n", "2"], "=", "n2"}], ";", "\[IndentingNewLine]",
"\t\t",
RowBox[{
SubscriptBox["\[Sigma]", "1"], "=", "\[Sigma]1"}], ";",
" ",
RowBox[{
SubscriptBox["\[Sigma]", "2"], "=", "\[Sigma]2"}], ";",
"\[IndentingNewLine]", "\t\t",
RowBox[{"\[Alpha]", "=", " ",
FractionBox["1",
SubscriptBox["n", "1"]]}], ";", " ",
RowBox[{"\[Beta]", "=", " ",
FractionBox["1",
SubscriptBox["n", "2"]]}], ";", "\[IndentingNewLine]", "\t\t ",
RowBox[{
SubscriptBox["w", "0"], "=",
RowBox[{"ArcCos", "[",
SqrtBox[
FractionBox[
SubscriptBox["n", "2"],
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"]}]]], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{
FractionBox[
RowBox[{
SuperscriptBox[
RowBox[{"(",
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "-", "1"}], "\[Alpha]"], ")"}],
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "-", "1"}], "2"]],
SuperscriptBox[
RowBox[{"(",
FractionBox[
RowBox[{
SubscriptBox["n", "2"], "-", "1"}], "\[Beta]"], ")"}],
FractionBox[
RowBox[{
SubscriptBox["n", "2"], "-", "1"}], "2"]],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox["1",
SubscriptBox["n", "1"]], "+",
FractionBox["1",
SubscriptBox["n", "2"]]}], ")"}],
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"], "-", "2"}], "2"]]}],
SuperscriptBox[
RowBox[{"(",
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"], "-", "2"}], ")"}],
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"], "-", "2"}], "2"]]], "\[Times]",
FractionBox[
RowBox[{"Gamma", "[",
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"]}], "2"], "]"}],
RowBox[{
RowBox[{"Gamma", "[",
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"], "-", "1"}], "2"], "]"}],
SqrtBox["\[Pi]"],
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"],
SubscriptBox["n", "1"]],
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"],
SubscriptBox["n", "2"]]}]], "\[IndentingNewLine]", "\t\t",
"\[Times]",
RowBox[{"NIntegrate", "[",
RowBox[{
FractionBox[
SuperscriptBox[
RowBox[{"Cos", "[", "x", "]"}],
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"], "-", "2"}]],
SuperscriptBox[
RowBox[{"(",
RowBox[{
FractionBox[
SuperscriptBox[
RowBox[{"Cos", "[",
RowBox[{"x", "-",
SubscriptBox["w", "0"]}], "]"}], "2"],
SuperscriptBox[
SubscriptBox["\[Sigma]", "1"], "2"]], "+",
FractionBox[
SuperscriptBox[
RowBox[{"Sin", "[",
RowBox[{"x", "-",
SubscriptBox["w", "0"]}], "]"}], "2"],
SuperscriptBox[
SubscriptBox["\[Sigma]", "2"], "2"]]}], ")"}],
FractionBox[
RowBox[{
SubscriptBox["n", "1"], "+",
SubscriptBox["n", "2"]}], "2"]]], ",",
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{
RowBox[{"-", "\[Pi]"}], "/", "2"}], ",",
RowBox[{"\[Pi]", "/", "2"}]}], "}"}]}], "]"}]}]}]}],
"\[IndentingNewLine]", "]"}]}]}]], "Input",
CellChangeTimes->{{3.5121608509335938`*^9, 3.512160851333008*^9},
3.512161090296875*^9, {3.5121612454804688`*^9, 3.512161270104492*^9}, {
3.5121613972890625`*^9, 3.5121614619189453`*^9}, {3.5121615038466797`*^9,