-
Notifications
You must be signed in to change notification settings - Fork 0
/
orth10d.nb
5383 lines (5279 loc) · 266 KB
/
orth10d.nb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 12.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 272368, 5375]
NotebookOptionsPosition[ 268887, 5313]
NotebookOutlinePosition[ 269229, 5328]
CellTagsIndexPosition[ 269186, 5325]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
RowBox[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SeedRandom", "[", "0", "]"}], ";"}]}], "Input",
CellChangeTimes->{{3.808249437390602*^9, 3.808249468789117*^9}, {
3.808249632175363*^9, 3.808249635636956*^9}, {3.808249765110931*^9,
3.808249766532976*^9}, 3.808254098303404*^9, {3.8082563792955236`*^9,
3.808256381597042*^9}, 3.808256704271912*^9, 3.8082567676558867`*^9, {
3.808257135602641*^9, 3.808257140940946*^9}, {3.80825720411176*^9,
3.808257204717432*^9}, 3.808257874664433*^9, 3.808293009163722*^9, {
3.808293569355691*^9, 3.80829356980971*^9}, 3.8083195815073442`*^9, {
3.816823365459618*^9, 3.816823365783063*^9}, {3.8172867370346537`*^9,
3.8172867391812983`*^9}, {3.8408121208094063`*^9, 3.8408121210905905`*^9},
3.840813073673602*^9},
CellLabel->"In[1]:=",ExpressionUUID->"45e836a1-e000-4571-9880-9cea737766a6"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "\"\<Sampler3`\>\""}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"STEPS", "=", "3"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CHAINS", "=", "3"}], ";"}]}], "Input",
CellChangeTimes->{{3.8168124021590147`*^9, 3.816812402159195*^9}, {
3.816820205778687*^9, 3.81682020614793*^9}, 3.8172501919117117`*^9,
3.817285432667527*^9, 3.817435500804758*^9, 3.817435604876099*^9, {
3.81744024760004*^9, 3.817440253709037*^9}, {3.820740856282202*^9,
3.8207409008730907`*^9}, 3.8207410047465763`*^9, 3.820741052776198*^9, {
3.8207410968162527`*^9, 3.820741100943769*^9}, {3.820980551767469*^9,
3.820980563997992*^9}, 3.820980617916573*^9, {3.823756418630746*^9,
3.823756424005458*^9}, {3.830112106490539*^9, 3.830112110185746*^9}, {
3.830112225007798*^9, 3.8301122717751417`*^9}, 3.8301123103879013`*^9, {
3.8318681355310917`*^9, 3.831868136241496*^9}, {3.832296299563047*^9,
3.832296325010687*^9}, {3.83852739275846*^9, 3.838527396171302*^9}, {
3.840779563262048*^9, 3.840779566979927*^9}, {3.8408103676709385`*^9,
3.840810368155196*^9}, {3.84081042060872*^9, 3.840810426607334*^9}, {
3.8408104704156914`*^9, 3.840810476554878*^9}, 3.8408111788743734`*^9, {
3.8418050430196953`*^9, 3.8418050474717274`*^9}, {3.8418063356485415`*^9,
3.8418063869332466`*^9}, {3.841806494209369*^9, 3.841806494678008*^9}, {
3.8418065939406333`*^9, 3.841806619512684*^9}, {3.8418066622869616`*^9,
3.841806662849296*^9}, {3.841806769417779*^9, 3.8418067701207323`*^9}, {
3.841807017057907*^9, 3.841807018729353*^9}, {3.84180751865545*^9,
3.8418075197334757`*^9}, {3.8418076141376753`*^9,
3.8418076177149467`*^9}, {3.8418078428852053`*^9,
3.8418078455094395`*^9}, {3.842840340000411*^9, 3.8428403420727468`*^9}, {
3.843020901590493*^9, 3.843020904152385*^9}, {3.843021520371374*^9,
3.843021520855631*^9}, 3.8430218072117414`*^9, {3.8450885857512116`*^9,
3.8450885859977894`*^9}, {3.8489418598391466`*^9, 3.848941860022662*^9}, {
3.848942696371207*^9, 3.848942698337946*^9}, {3.848970730658123*^9,
3.8489707319836636`*^9}},
CellLabel->"In[3]:=",ExpressionUUID->"aac9ca2d-c078-432d-afca-9e864ccd42b9"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"SIGMA", "=",
RowBox[{"Table", "[",
RowBox[{"0", ",", "10", ",", "10"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"For", "[",
RowBox[{
RowBox[{"i", "=", "1"}], ",",
RowBox[{"i", "\[LessEqual]", "10"}], ",",
RowBox[{"i", "++"}], ",",
RowBox[{
RowBox[{"SIGMA", "[",
RowBox[{"[",
RowBox[{"i", ",", "i"}], "]"}], "]"}], "=",
FractionBox["1",
SuperscriptBox["11",
RowBox[{"2",
RowBox[{"(",
RowBox[{"i", "-", "1"}], ")"}]}]]]}]}], "]"}],
";"}], "\[IndentingNewLine]",
SqrtBox[
RowBox[{"Diagonal", "[", "SIGMA", "]"}]], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"U", "[",
RowBox[{
"x1_", ",", "x2_", ",", "x3_", ",", "x4_", ",", "x5_", ",", "x6_", ",",
"x7_", ",", "x8_", ",", "x9_", ",", "x10_"}], "]"}], "=",
RowBox[{
FractionBox["1", "2"],
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"{",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "}"}], ".",
RowBox[{"LinearSolve", "[",
RowBox[{"SIGMA", ",",
RowBox[{"{",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",",
"x7", ",", "x8", ",", "x9", ",", "x10"}], "}"}]}], "]"}]}],
"]"}]}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"dU", "[",
RowBox[{
"x1_", ",", "x2_", ",", "x3_", ",", "x4_", ",", "x5_", ",", "x6_", ",",
"x7_", ",", "x8_", ",", "x9_", ",", "x10_"}], "]"}], "=",
RowBox[{"GradientG", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "}"}]}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"ddU", "[",
RowBox[{
"x1_", ",", "x2_", ",", "x3_", ",", "x4_", ",", "x5_", ",", "x6_", ",",
"x7_", ",", "x8_", ",", "x9_", ",", "x10_"}], "]"}], "=",
RowBox[{"HessianH", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"dddU", "[",
RowBox[{
"x1_", ",", "x2_", ",", "x3_", ",", "x4_", ",", "x5_", ",", "x6_", ",",
"x7_", ",", "x8_", ",", "x9_", ",", "x10_"}], "]"}], "=",
RowBox[{"D3", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6", ",", "x7",
",", "x8", ",", "x9", ",", "x10"}], "}"}]}], "]"}]}], ";"}]}], "Input",\
CellChangeTimes->{{3.808319417169717*^9, 3.808319426953169*^9}, {
3.8083195116541348`*^9, 3.808319512297307*^9}, 3.8083195999310637`*^9, {
3.814330433516584*^9, 3.814330446640849*^9}, {3.817440287245225*^9,
3.8174403199989433`*^9}, 3.817440673136303*^9, 3.81744086413822*^9, {
3.817503864474145*^9, 3.817503908025481*^9}, {3.81750395539356*^9,
3.817504043249189*^9}, {3.8175040797242107`*^9, 3.817504129818478*^9}, {
3.817504253290216*^9, 3.817504314357682*^9}, 3.817504422272829*^9, {
3.817506453495331*^9, 3.817506515588607*^9}, 3.817506560648061*^9, {
3.8175066227361107`*^9, 3.817506625250626*^9}, 3.817507141139452*^9,
3.832029041813363*^9, 3.838527436371852*^9, 3.838527548753686*^9, {
3.840779432029445*^9, 3.8407794484632454`*^9}, 3.840811172360276*^9, {
3.8418057845835667`*^9, 3.841805785083457*^9}, {3.8418058289986124`*^9,
3.841806018820402*^9}, 3.841806086823037*^9, {3.8418061963348684`*^9,
3.8418061968973303`*^9}, 3.8418062502751436`*^9, {3.841807381066459*^9,
3.841807389895855*^9}, {3.8418074440283194`*^9, 3.841807450276975*^9}, {
3.8418074813475847`*^9, 3.841807482940955*^9}, {3.8418078598546057`*^9,
3.841807876444446*^9}, {3.8418079575345454`*^9, 3.8418079644080524`*^9}, {
3.8418080033236895`*^9, 3.8418080855540257`*^9}, {3.8418082011085997`*^9,
3.8418082158760004`*^9}, {3.841808285468857*^9, 3.8418082867654266`*^9}, {
3.8418089683433685`*^9, 3.8418089966177692`*^9}, {3.8418100852801003`*^9,
3.841810146301073*^9}, {3.8418101934939985`*^9, 3.8418101942438173`*^9}, {
3.841810266563077*^9, 3.841810309755947*^9}, {3.8418196078691325`*^9,
3.84181963551235*^9}, {3.841821911584039*^9, 3.8418219120059543`*^9}, {
3.8423375976092153`*^9, 3.8423376166986237`*^9}, {3.8424136414611907`*^9,
3.842413651641822*^9}, {3.8424139426879463`*^9, 3.842413945718592*^9}, {
3.8425951986123714`*^9, 3.8425951990964775`*^9}, {3.8425953206715794`*^9,
3.8425953374958878`*^9}, {3.842720034675085*^9, 3.842720064587289*^9},
3.842828066582609*^9, {3.8428397129333453`*^9, 3.8428397131347866`*^9}, {
3.843018969009003*^9, 3.8430189695869617`*^9}, 3.8430816652221837`*^9,
3.8430818430314736`*^9, {3.843081885705332*^9, 3.84308188591581*^9},
3.843081921109686*^9, {3.8430841212989445`*^9, 3.8430841215372944`*^9},
3.8430841539179535`*^9, 3.8430850588133893`*^9, {3.8433528912335787`*^9,
3.8433528935235367`*^9}, {3.8433529477864714`*^9,
3.8433529772364206`*^9}, {3.843353522617175*^9, 3.84335352434713*^9}, {
3.8448354502928543`*^9, 3.8448356268010197`*^9}, {3.844836002317313*^9,
3.8448360850203357`*^9}, {3.844836181864067*^9, 3.8448362109895277`*^9}, {
3.844838833014307*^9, 3.844838836931977*^9}, 3.844841112731256*^9, {
3.8448411913625307`*^9, 3.8448411940032415`*^9}, {3.844841233740271*^9,
3.8448412362254496`*^9}, {3.8448412830525*^9, 3.8448412838061943`*^9}, {
3.845088568583091*^9, 3.8450885689366117`*^9}, {3.845088621836038*^9,
3.8450886338368154`*^9}, 3.84508866856703*^9, {3.8459327853940105`*^9,
3.8459327865460496`*^9}, {3.8459329340827837`*^9, 3.845932936293355*^9}, {
3.845932981880821*^9, 3.845932982155089*^9}, {3.845933363836052*^9,
3.8459333659054255`*^9}, {3.8459336371326237`*^9, 3.845933642607249*^9},
3.8459336860459023`*^9, {3.8465289453341227`*^9, 3.84652894590823*^9}, {
3.846532000416685*^9, 3.8465320305682936`*^9}, {3.8465327423202133`*^9,
3.8465327438780174`*^9}, 3.8465335272570133`*^9, {3.8465344007792497`*^9,
3.8465344014717445`*^9}, {3.8489418884233685`*^9,
3.8489419172019787`*^9}, {3.848942095255129*^9, 3.8489421276782703`*^9}, {
3.848942224315063*^9, 3.848942245347438*^9}, {3.8489423242551203`*^9,
3.8489423247763386`*^9}, {3.848942365713991*^9, 3.848942366506889*^9}, {
3.848942464628497*^9, 3.8489424694956675`*^9}, {3.848942575869368*^9,
3.8489425997552223`*^9}, {3.848942666015644*^9, 3.8489426677119293`*^9}, {
3.848942706861102*^9, 3.848942720757658*^9}, {3.8489430034373465`*^9,
3.8489430036071815`*^9}, {3.848956977971016*^9, 3.848956995771182*^9}, {
3.848957080704894*^9, 3.8489570820710635`*^9}, {3.8489706656565638`*^9,
3.848970680525428*^9}, {3.8489707915245113`*^9, 3.848970792337356*^9}, {
3.8491533559902735`*^9, 3.8491533766443987`*^9}},
CellLabel->"In[7]:=",ExpressionUUID->"b9b3fe6e-f312-47f1-8786-232c79fe6d00"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"1", ",",
FractionBox["1", "11"], ",",
FractionBox["1", "121"], ",",
FractionBox["1", "1331"], ",",
FractionBox["1", "14641"], ",",
FractionBox["1", "161051"], ",",
FractionBox["1", "1771561"], ",",
FractionBox["1", "19487171"], ",",
FractionBox["1", "214358881"], ",",
FractionBox["1", "2357947691"]}], "}"}]], "Output",
CellChangeTimes->{
3.8489426025933847`*^9, 3.848942671546382*^9, {3.848942703396655*^9,
3.848942722907362*^9}, 3.848942829010537*^9, 3.848943015662877*^9,
3.848957017168104*^9, 3.848957085102908*^9, 3.848970707353242*^9,
3.8489707661035967`*^9, {3.848970835830801*^9, 3.8489708505490403`*^9}, {
3.848970914456676*^9, 3.848970941565317*^9}, 3.8489709886754885`*^9,
3.8491534659283504`*^9, 3.849162466176684*^9, 3.8491625960180187`*^9},
CellLabel->"Out[9]=",ExpressionUUID->"20a2933f-fc32-4464-a726-b353373ed440"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
SqrtBox[
RowBox[{"Diagonal", "[", "SIGMA", "]"}]], "//", "N"}]], "Input",
CellChangeTimes->{{3.840773518367429*^9, 3.8407735279617853`*^9}, {
3.8418062077540145`*^9, 3.8418062232660537`*^9}},
CellLabel->"In[14]:=",ExpressionUUID->"da99e213-762b-4c24-af15-1c498fa89576"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
"1.`", ",", "0.09090909090909091`", ",", "0.008264462809917356`", ",",
"0.0007513148009015778`", ",", "0.00006830134553650706`", ",",
"6.209213230591551`*^-6", ",", "5.644739300537774`*^-7", ",",
"5.1315811823070673`*^-8", ",", "4.665073802097334`*^-9", ",",
"4.2409761837248493`*^-10"}], "}"}]], "Output",
CellChangeTimes->{
3.840774008254961*^9, 3.840778842348008*^9, {3.8407794508271666`*^9,
3.840779455966586*^9}, 3.8407795252472177`*^9, 3.840779570327211*^9,
3.8408103051194973`*^9, 3.8408103714576917`*^9, 3.8408104330338297`*^9,
3.840810479366717*^9, 3.840811280006689*^9, 3.8408113387167835`*^9,
3.8408113799883833`*^9, 3.840811478700859*^9, 3.840812029815527*^9,
3.8408121243242073`*^9, 3.840812606932948*^9, 3.840812679654438*^9,
3.840813075657529*^9, 3.8418050508459473`*^9, 3.8418051674832*^9,
3.8418052122382016`*^9, 3.8418052809439445`*^9, 3.8418053214499507`*^9,
3.8418053737079887`*^9, 3.8418056259213033`*^9, 3.84180571483436*^9,
3.8418060471729097`*^9, 3.841806118114992*^9, {3.8418062020522137`*^9,
3.8418062402150297`*^9}, {3.841806345458686*^9, 3.8418063962904015`*^9},
3.841806430224533*^9, {3.8418064750576344`*^9, 3.8418065001454287`*^9}, {
3.841806600517084*^9, 3.8418066222150574`*^9}, 3.8418066843910575`*^9,
3.8418067878387947`*^9, 3.8418070312732763`*^9, {3.841807092618102*^9,
3.8418071124415135`*^9}, 3.841807268620556*^9, 3.841807395207095*^9,
3.8418074864869943`*^9, 3.841807522732607*^9, 3.8418076213390837`*^9, {
3.841807862807028*^9, 3.841807881911867*^9}, 3.8418079696097975`*^9,
3.8418080121651993`*^9, {3.8418080532804337`*^9, 3.841808088600175*^9},
3.841808156447324*^9, 3.8418082517581024`*^9, 3.8418082889055347`*^9,
3.8418086674493237`*^9, {3.841808971545592*^9, 3.8418090177221303`*^9}, {
3.8418100978082485`*^9, 3.841810149597165*^9}, 3.8418101975086718`*^9, {
3.841810276435729*^9, 3.841810312708369*^9}, 3.841818270769696*^9,
3.8418195939505663`*^9, 3.8418196394801545`*^9, {3.8418197037617702`*^9,
3.841819715462118*^9}, 3.8418197456620398`*^9, 3.841821275498319*^9,
3.841821368742087*^9, {3.841821444680805*^9, 3.841821471330825*^9}, {
3.8418215123680325`*^9, 3.8418215490514903`*^9}, 3.8418216129427347`*^9,
3.8418216832434225`*^9, 3.841821713658154*^9, 3.8418217988869047`*^9,
3.8418218510309353`*^9, 3.841821914926998*^9, 3.8420693354810886`*^9, {
3.842069376309746*^9, 3.8420694032253027`*^9}, {3.84206952794975*^9,
3.8420695516004853`*^9}, 3.8423372647103786`*^9, 3.842337359231823*^9,
3.8423374528757267`*^9, {3.8423376015926533`*^9, 3.8423376498938236`*^9}, {
3.842413647166402*^9, 3.8424136574529257`*^9}, 3.842413706209156*^9, {
3.8424138593771796`*^9, 3.8424139081819973`*^9}, 3.842413949998724*^9,
3.8424140002066183`*^9, 3.8424140485118556`*^9, 3.842414090411417*^9,
3.842414136264843*^9, 3.842414213938777*^9, 3.842595171646274*^9,
3.8425952022363605`*^9, {3.842595327716798*^9, 3.842595340354454*^9},
3.842600116907176*^9, 3.8427200481909*^9, 3.842720081846768*^9,
3.8428135668055596`*^9, 3.8428136465217147`*^9, 3.842813762342836*^9,
3.842813845342249*^9, 3.8428279983338184`*^9, 3.842828069863003*^9,
3.8428281497849874`*^9, {3.8428282437316456`*^9, 3.8428282715250616`*^9},
3.8428283203551407`*^9, 3.8428284022155914`*^9, 3.842828459045972*^9, {
3.8428285036792793`*^9, 3.8428285302355323`*^9}, 3.842828683062028*^9,
3.8428288488333883`*^9, 3.8428289611144104`*^9, 3.842828992741126*^9,
3.8428290294878607`*^9, 3.842829077318731*^9, {3.8428291264997463`*^9,
3.8428291548212323`*^9}, 3.8428293054487257`*^9, 3.842829410209633*^9,
3.8428296231210556`*^9, 3.8428297448612413`*^9, {3.842830007935007*^9,
3.842830094903678*^9}, {3.8428301653714757`*^9, 3.842830192774476*^9},
3.842830235045782*^9, 3.842830269709503*^9, {3.84283040291708*^9,
3.842830429103552*^9}, {3.842830482059876*^9, 3.8428305275022793`*^9},
3.842830594505308*^9, 3.842830630012581*^9, 3.8428306644463005`*^9,
3.842830709793781*^9, 3.8428307418174973`*^9, 3.842830852060021*^9,
3.842830940870228*^9, 3.842831003214943*^9, 3.842831048789923*^9,
3.842831115852272*^9, {3.842831176928621*^9, 3.8428312056573477`*^9},
3.8428312840646544`*^9, 3.8428313586451163`*^9, {3.8428314336151295`*^9,
3.8428314604213047`*^9}, 3.8428315195824146`*^9, 3.8428316031252594`*^9,
3.8428316950001316`*^9, 3.84283179668251*^9, 3.842833732272644*^9,
3.84283384958066*^9, 3.842833890448251*^9, 3.84283398640561*^9,
3.842834339477416*^9, 3.8428344023881283`*^9, 3.8428345037752476`*^9,
3.842834539141957*^9, 3.8428346038298283`*^9, 3.8428352707776675`*^9,
3.8428394751739664`*^9, 3.8428395757408295`*^9, 3.8428396347622547`*^9,
3.8428397160785275`*^9, 3.8428397599213896`*^9, 3.842839835635334*^9,
3.8428398749538765`*^9, 3.842839967427287*^9, 3.842840048071313*^9,
3.8428401651908555`*^9, 3.842840242638155*^9, 3.842840344880452*^9,
3.843016899243758*^9, {3.8430172862515883`*^9, 3.843017303266927*^9},
3.8430180049478693`*^9, 3.843018116301386*^9, 3.843018631662936*^9,
3.8430189722269926`*^9, 3.8430191810931945`*^9, 3.843019241613795*^9,
3.8430203388755984`*^9, 3.8430207019826326`*^9, 3.843020863533249*^9,
3.843020909604186*^9, {3.8430209477982597`*^9, 3.8430209675435805`*^9},
3.843021006506836*^9, 3.8430210532613797`*^9, 3.8430211304669533`*^9,
3.8430212677930665`*^9, 3.8430213481646223`*^9, 3.843021549630072*^9,
3.843021682542788*^9, 3.8430218096329036`*^9, 3.843021903881338*^9,
3.843080797624956*^9, 3.843081531566362*^9, {3.843081638750016*^9,
3.843081689235669*^9}, 3.8430817516935434`*^9, {3.8430818459375906`*^9,
3.8430818888052225`*^9}, 3.843081923517377*^9, {3.8430836263455296`*^9,
3.843083641554844*^9}, 3.843083692121623*^9, {3.8430841246038847`*^9,
3.8430841563554354`*^9}, 3.843085085474568*^9, 3.8430851214783196`*^9,
3.843085199232632*^9, {3.8430868147759504`*^9, 3.8430868444015083`*^9},
3.8430870017613826`*^9, 3.8430871394127283`*^9, 3.8430871825039234`*^9,
3.843196974075267*^9, 3.8433263364298196`*^9, 3.843326375979635*^9,
3.8433276352616863`*^9, 3.843327676816004*^9, 3.8433279167983875`*^9,
3.8433281223505545`*^9, 3.843328226260437*^9, 3.843328273062622*^9,
3.843329171102326*^9, {3.843329310057688*^9, 3.843329336306964*^9},
3.8433294295942764`*^9, 3.843329544053123*^9, 3.843329577650732*^9,
3.8433529143435583`*^9, {3.843352951746463*^9, 3.843352979626416*^9},
3.8433532483276587`*^9, 3.843353458097499*^9, 3.8433535275071545`*^9,
3.843353597227128*^9, 3.843353710679614*^9, 3.843353777319584*^9,
3.8433539188919487`*^9, 3.84335403748454*^9, 3.844568627377505*^9,
3.844568714301239*^9, 3.8445688020642405`*^9, 3.844569119185997*^9,
3.844569160467536*^9, 3.8445692643386436`*^9, {3.844569411741331*^9,
3.844569433239775*^9}, 3.844569494561277*^9, 3.8448356409257917`*^9,
3.8448356991278286`*^9, {3.8448358108255167`*^9, 3.8448358331668386`*^9},
3.8448360886621256`*^9, 3.8448361341623325`*^9, {3.844836188099186*^9,
3.844836213656756*^9}, 3.8448363407480164`*^9, 3.8448364914633107`*^9,
3.844838713584502*^9, 3.8448388441450634`*^9, 3.844838974911293*^9,
3.844839168825938*^9, 3.844839434824155*^9, 3.844839695696764*^9,
3.844839924284713*^9, 3.844840183186515*^9, 3.844840553186381*^9,
3.8448408358358293`*^9, 3.8448410674597635`*^9, {3.8448411154785614`*^9,
3.8448411553013787`*^9}, {3.844841196796488*^9, 3.844841238747655*^9},
3.844841295071714*^9, {3.8450885776905317`*^9, 3.84508858842013*^9},
3.8450886389476233`*^9, 3.845088670790506*^9, 3.8450888453535852`*^9,
3.8459326065123963`*^9, 3.8459327528423986`*^9, 3.84593278880884*^9,
3.845932938586872*^9, 3.8459329849424534`*^9, 3.8459333269172335`*^9,
3.84593336858327*^9, 3.845933557281503*^9, 3.8459336450171556`*^9,
3.8459336883450565`*^9, 3.8459339174719014`*^9, 3.8465288128036237`*^9,
3.8465288719465733`*^9, 3.846528914742838*^9, 3.8465289488116903`*^9,
3.846529022529578*^9, 3.846529067616542*^9, 3.846529121740415*^9,
3.846529169751582*^9, {3.846529207943014*^9, 3.8465292316077023`*^9},
3.846529318553253*^9, 3.8465293740326056`*^9, 3.8465294613488245`*^9,
3.8465295306370015`*^9, 3.8465296308313627`*^9, 3.8465297623262215`*^9,
3.846529871850564*^9, 3.8465300878491774`*^9, 3.846530246081091*^9,
3.8465304655912333`*^9, 3.846530583100012*^9, 3.8465308235927486`*^9,
3.846531042792621*^9, 3.8465310926701713`*^9, 3.8465311462078094`*^9,
3.846531204684164*^9, 3.8465314574981203`*^9, 3.8465315089984365`*^9,
3.8465316077964993`*^9, {3.846531693222903*^9, 3.8465317179157715`*^9},
3.8465318353520465`*^9, 3.846531913219776*^9, 3.846531967240836*^9,
3.846532004372154*^9, 3.8465320351199284`*^9, 3.8465320947836795`*^9,
3.846532216957611*^9, 3.8465322581551094`*^9, 3.8465323442031436`*^9,
3.8465325073247995`*^9, 3.8465326805200987`*^9, 3.8465327504174733`*^9,
3.846532923542731*^9, 3.8465329732766857`*^9, 3.8465330193272486`*^9,
3.8465331167241197`*^9, 3.8465331544325953`*^9, 3.8465331903211746`*^9,
3.84653323652709*^9, 3.846533271048812*^9, 3.846533529365591*^9, {
3.846533684342716*^9, 3.8465337122300673`*^9}, 3.8465338354972973`*^9,
3.8465339497054615`*^9, 3.846533983108399*^9, 3.846534104198867*^9,
3.846534202058828*^9, 3.846534325742277*^9, 3.8465344045639544`*^9,
3.846534471339347*^9, 3.846534546853032*^9, 3.84653467617525*^9,
3.8465347948998203`*^9, 3.8465349122266755`*^9, 3.8465349654403887`*^9,
3.8465350101075335`*^9, 3.8465351519028444`*^9, 3.8465351823597794`*^9,
3.846535309246544*^9, 3.8465353631503496`*^9, 3.8465354124494514`*^9,
3.8489421491428685`*^9, 3.848942186373189*^9, 3.8489422605984945`*^9,
3.8489423276330504`*^9, 3.848942369414426*^9, 3.8489424722667475`*^9,
3.848942602757333*^9, 3.8489426717130013`*^9, {3.848942703527319*^9,
3.8489427230212936`*^9}, 3.848942829138376*^9, 3.8489430158168397`*^9,
3.848957017436861*^9, 3.8489570852500916`*^9, 3.84897070773215*^9,
3.84897076632122*^9, {3.848970835918686*^9, 3.848970850655794*^9}, {
3.8489709145708027`*^9, 3.848970941653328*^9}, 3.8489709887595916`*^9,
3.849153466156425*^9, 3.8491624662995987`*^9, 3.8491625962988963`*^9},
CellLabel->"Out[14]=",ExpressionUUID->"ea6f8627-88eb-4cc5-b388-c72882c49e10"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"dKdq", "[",
RowBox[{"p_", ",", "q_", ",", "i_"}], "]"}], ":=", "0"}], ";"}]], "Input",\
CellChangeTimes->{{3.848970906306204*^9, 3.848970908270934*^9}},
CellLabel->"In[15]:=",ExpressionUUID->"b29abeb7-f628-4238-9d5f-4ba966b94d47"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"QS", "=",
RowBox[{"hmc", "[",
RowBox[{
"U", ",", "dU", ",", "ddU", ",", "dddU", ",", "10", ",", "5000", ",",
"10000", ",",
RowBox[{"{", "}"}]}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.8083300207373047`*^9, 3.808330021545154*^9}, {
3.808330057606822*^9, 3.8083300651679897`*^9}, {3.8084239541699953`*^9,
3.80842395728685*^9}, {3.8084239901973343`*^9, 3.808423995188442*^9},
3.808424069693726*^9, {3.808424279402*^9, 3.808424280309854*^9}, {
3.808424475946063*^9, 3.808424476285922*^9}, {3.808424556851679*^9,
3.808424559741782*^9}, {3.814329296202931*^9, 3.8143292994796143`*^9}, {
3.8143295041114073`*^9, 3.8143295054599323`*^9}, {3.814329810538707*^9,
3.814329811379229*^9}, {3.814329988519002*^9, 3.8143299916986837`*^9}, {
3.8143305561948767`*^9, 3.814330559434862*^9}, {3.814330647247871*^9,
3.814330648698496*^9}, {3.814330761246181*^9, 3.814330762762741*^9}, {
3.814331829179998*^9, 3.8143318412578373`*^9}, {3.816820029273882*^9,
3.816820032992787*^9}, {3.816820354397133*^9, 3.816820365714497*^9}, {
3.8168232128438177`*^9, 3.8168232191226807`*^9}, 3.816823302483911*^9, {
3.817250728188867*^9, 3.817250735153274*^9}, 3.817251209164393*^9, {
3.8172515707942343`*^9, 3.817251578534383*^9}, {3.817253057791007*^9,
3.817253059916018*^9}, {3.817440652184346*^9, 3.8174406526768913`*^9}, {
3.817504347067464*^9, 3.8175044007622213`*^9}, {3.8175069155531683`*^9,
3.817506931370562*^9}, {3.817507457551218*^9, 3.8175074637568502`*^9}, {
3.817507501262887*^9, 3.8175075272461367`*^9}, {3.817507561779567*^9,
3.8175075886919003`*^9}, {3.81750764435394*^9, 3.8175076452708817`*^9}, {
3.8237564121436577`*^9, 3.8237564131344557`*^9}, 3.830112406886807*^9, {
3.832023346481739*^9, 3.832023356466961*^9}, {3.832029094230179*^9,
3.8320290964168653`*^9}, {3.832029462117546*^9, 3.832029470132633*^9},
3.832296187659807*^9, {3.832296239252947*^9, 3.832296255554605*^9}, {
3.8322963299030237`*^9, 3.83229633284181*^9}, 3.83229641034589*^9,
3.832296454608801*^9, {3.832296498322076*^9, 3.832296500159094*^9}, {
3.838527121487178*^9, 3.838527129028344*^9}, 3.8385272264975977`*^9,
3.838527366065442*^9, {3.838527452828074*^9, 3.838527456556349*^9},
3.8385275422634487`*^9, {3.838528228323004*^9, 3.838528266447521*^9}, {
3.8407739680721984`*^9, 3.8407739696187105`*^9}, 3.8408112452120385`*^9, {
3.8408114061384697`*^9, 3.8408114237126036`*^9}, {3.840812488175477*^9,
3.840812624663289*^9}, {3.841806038393789*^9, 3.841806043376937*^9}, {
3.841807260185065*^9, 3.8418072622002077`*^9}, {3.841807632242778*^9,
3.841807635570116*^9}, {3.8418181137901273`*^9, 3.8418181142431493`*^9}, {
3.8418195861555343`*^9, 3.841819588983119*^9}, {3.8418197373477807`*^9,
3.8418197425221295`*^9}, {3.8418212579242935`*^9,
3.8418212636885962`*^9}, {3.841821527708317*^9, 3.8418215299890404`*^9}, {
3.8418216083345733`*^9, 3.841821610396589*^9}, {3.841821695115665*^9,
3.8418217045353284`*^9}, {3.841821794762909*^9, 3.8418217953410025`*^9}, {
3.8418218481723547`*^9, 3.841821871541855*^9}, {3.841822046261695*^9,
3.841822047948801*^9}, {3.842069330029236*^9, 3.8420693304978495`*^9}, {
3.8420693647630935`*^9, 3.8420694006790714`*^9}, {3.84206952373199*^9,
3.842069548648167*^9}, {3.842337253681719*^9, 3.84233726122696*^9},
3.8423373547641525`*^9, {3.842337449204733*^9, 3.84233745056382*^9}, {
3.842337639786823*^9, 3.8423376468634386`*^9}, 3.842413703616016*^9, {
3.8424138543158813`*^9, 3.842413905010866*^9}, {3.842413996082621*^9,
3.842413997707345*^9}, {3.842414044184737*^9, 3.8424140450127797`*^9}, {
3.8424140874120865`*^9, 3.842414087646433*^9}, {3.842414122278944*^9,
3.842414133781046*^9}, {3.842414209580457*^9, 3.8424142105645742`*^9}, {
3.842595158930546*^9, 3.84259516739741*^9}, {3.8426001132674484`*^9,
3.842600113439242*^9}, {3.8427200436816893`*^9, 3.8427200442318497`*^9}, {
3.842720077336485*^9, 3.842720078664167*^9}, {3.8428135627909136`*^9,
3.842813563228426*^9}, {3.842813643428687*^9, 3.842813643491151*^9}, {
3.8428137475182285`*^9, 3.842813749611457*^9}, {3.842813842514822*^9,
3.8428138430459137`*^9}, {3.842828146785694*^9, 3.842828147035638*^9}, {
3.842828232484284*^9, 3.8428282634957304`*^9}, {3.8428282946134057`*^9,
3.8428283114353633`*^9}, {3.8428283936863475`*^9,
3.8428283994037547`*^9}, {3.8428284444244156`*^9,
3.8428284562030425`*^9}, {3.8428284865707455`*^9,
3.8428285248617983`*^9}, {3.842828679376136*^9, 3.8428286803151703`*^9}, {
3.842829406554243*^9, 3.8428294069605427`*^9}, {3.842829738659577*^9,
3.8428297421120415`*^9}, {3.842830246746155*^9, 3.8428302665853777`*^9}, {
3.842830421308514*^9, 3.8428304266510077`*^9}, {3.842830496728301*^9,
3.842830524502989*^9}, {3.8428306542768197`*^9, 3.842830662274938*^9}, {
3.842830706826956*^9, 3.842830739349332*^9}, {3.842830806086615*^9,
3.8428308064771323`*^9}, {3.84283084318711*^9, 3.8428308489826226`*^9}, {
3.8428309251553307`*^9, 3.842830938074049*^9}, {3.8428309973725643`*^9,
3.842831000356234*^9}, {3.8428310410105033`*^9, 3.842831056475647*^9}, {
3.8428311053234987`*^9, 3.8428311127281446`*^9}, {3.8428315142086773`*^9,
3.842831517207971*^9}, {3.8428315967831397`*^9, 3.8428316004383936`*^9}, {
3.8428316828311186`*^9, 3.842831687236478*^9}, {3.8428317868098326`*^9,
3.8428317932928233`*^9}, {3.842833721618902*^9, 3.842833728789093*^9}, {
3.842833863218045*^9, 3.8428338702476416`*^9}, {3.84283397122165*^9,
3.8428339794228745`*^9}, 3.8428340778893127`*^9, {3.842834327058565*^9,
3.8428343363376436`*^9}, {3.8428345277541084`*^9, 3.842834534393077*^9}, {
3.8428401604295535`*^9, 3.8428401624672117`*^9}, {3.842840238557994*^9,
3.842840240154443*^9}, {3.8430172895809207`*^9, 3.8430172930332537`*^9}, {
3.8430192247271967`*^9, 3.843019229929088*^9}, {3.843020935707376*^9,
3.8430209381599216`*^9}, {3.84302112412471*^9, 3.8430211263116894`*^9}, {
3.843021539648193*^9, 3.843021546224619*^9}, {3.8430216642971544`*^9,
3.8430216724828143`*^9}, 3.8430218813711886`*^9, {3.843080617375188*^9,
3.843080636947776*^9}, {3.8430816333072166`*^9, 3.843081635096697*^9}, {
3.843081740993554*^9, 3.843081748290124*^9}, {3.8430818576657033`*^9,
3.843081864829067*^9}, {3.843082136077016*^9, 3.843082136410294*^9},
3.8430841388108144`*^9, {3.8430851925519524`*^9, 3.8430851962499743`*^9}, {
3.843086837994474*^9, 3.843086841177266*^9}, 3.8433282229017887`*^9, {
3.843328267557907*^9, 3.8433282694318466`*^9}, {3.8433293041468153`*^9,
3.8433293338943624`*^9}, 3.843329534147857*^9, {3.8433295710625615`*^9,
3.843329575134762*^9}, {3.8433535932371016`*^9, 3.8433535946772313`*^9}, {
3.8433537727296047`*^9, 3.843353774699587*^9}, 3.844569408518029*^9, {
3.8448358042248645`*^9, 3.844835806706053*^9}, 3.844838546981146*^9, {
3.844838595393461*^9, 3.844838605454103*^9}, {3.844841061727195*^9,
3.8448410637249355`*^9}, {3.84484112564672*^9, 3.8448411287572856`*^9},
3.8448412094393167`*^9, {3.845088551900528*^9, 3.845088555512825*^9}, {
3.845088809546076*^9, 3.845088809854854*^9}, {3.8465287251967087`*^9,
3.846528739728691*^9}, {3.8465288056427803`*^9, 3.8465288070340586`*^9}, {
3.846529114279646*^9, 3.846529117941084*^9}, {3.8465291637977996`*^9,
3.846529166878932*^9}, {3.8465316893206377`*^9, 3.846531689703659*^9}, {
3.846531736915019*^9, 3.8465317371298976`*^9}, {3.8465330108299246`*^9,
3.8465330128993626`*^9}, {3.8465332645715876`*^9, 3.846533266228161*^9}, {
3.846533692815788*^9, 3.8465337061072035`*^9}, {3.8465339717249355`*^9,
3.846533975322338*^9}, {3.8465343572367697`*^9, 3.846534360659613*^9}, {
3.846534538921303*^9, 3.8465345425468135`*^9}, {3.8465347881334195`*^9,
3.846534790958355*^9}, {3.846535338134053*^9, 3.846535341372508*^9}, {
3.8465354014840703`*^9, 3.846535406008998*^9}, {3.8489421581385536`*^9,
3.8489421807544866`*^9}, {3.848942252495022*^9, 3.848942252669842*^9}, {
3.848942422741745*^9, 3.848942427043209*^9}, {3.848942478946973*^9,
3.848942479514456*^9}, {3.8489428213290415`*^9, 3.848942825388419*^9}, {
3.8489430100961237`*^9, 3.848943012639429*^9}, {3.848970693043109*^9,
3.848970698169097*^9}, 3.848970762243975*^9, {3.848970830875198*^9,
3.848970846094433*^9}, {3.8489709350050545`*^9, 3.8489709390698247`*^9}},
CellLabel->"In[16]:=",ExpressionUUID->"8ecbaa4f-6d9e-4355-98b5-0089406492fe"],
Cell[CellGroupData[{
Cell[BoxData[
TemplateBox[{
" ","\" \"","1001","4.0488505997285295`*^15",
"9.577911076445036`*^-20","0.000045399929762484854`","0.`",
"4.412126545338385`*^10","1",RowBox[{"{", "1", "}"}],
RowBox[{"{", "4", "}"}],RowBox[{"{",
RowBox[{
"0.002430442096739797`", ",", "0.004305676431342455`", ",",
"0.000851855127950063`", ",", "0.0009370406407450693`", ",",
"0.00010464697253675069`", ",", "0.00003334375043029813`", ",",
"0.000014141005145095954`", ",", "0.000010624346465135951`", ",",
"7.2565715901481806`*^-6", ",", "0.000010624346465135951`"}], "}"}],
RowBox[{"{",
RowBox[{
"3.6808173937459345`*^15", ",", "2.7288813118496245`*^15", ",",
"1.9084545048565625`*^15", ",", "2.792996672634261`*^15", ",",
"1.9084524327152372`*^15", ",", "2.309198223798012`*^15", ",",
"2.3091989915384965`*^15", ",", "1.908452432104913`*^15", ",",
"1.908452432104913`*^15", ",", "1.9084524320323812`*^15"}], "}"}],
"0.002430442096739797`","False","True"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491625986928444`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"fc175277-3898-4fd7-a9e1-7dcda0194954"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","2002","2.899449042491806`*^11",
"1.7334860205774068`*^-17","0.000045399929762484854`","0.`",
"2806.8213972578756`","2",RowBox[{"{",
RowBox[{"1", ",", "2"}], "}"}],RowBox[{"{", "4", "}"}],RowBox[{"{",
RowBox[{
"0.00007862281933640175`", ",", "0.00010464697253675088`", ",",
"0.00003667812547332795`", ",", "0.00001711061622556611`", ",",
"4.505760032628285`*^-6", ",", "7.367277700816571`*^-7", ",",
"7.367277700816577`*^-7", ",", "8.104005470898228`*^-7", ",",
"4.1586361975774865`*^-7", ",", "2.8404044789136573`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"4.732088122786947`*^11", ",", "2.635862793969856`*^11", ",",
"2.0277971876299185`*^11", ",", "2.967564387212573`*^11", ",",
"1.843453378759692`*^11", ",", "8.470318813461284`*^11", ",",
"9.317256628042821`*^11", ",", "1.2401633340631514`*^12", ",",
"3.5383872048705586`*^12", ",", "5.698606988876715`*^12"}], "}"}],
"0.00010464697253675088`","False","True"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491626000232887`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"bfbbd5af-47e5-4ee1-82d5-99c207c2db1a"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","3003","7.215734416963543`*^9",
"1.1147073539419691`*^-15","0.07870635316840176`","0.06465083806036002`",
"14.077616533842242`","3",RowBox[{"{", "1", "}"}],RowBox[{"{", "4", "}"}],
RowBox[{"{",
RowBox[{
"0.000014141005145095959`", ",", "0.000011686781111649564`", ",",
"8.780451624079308`*^-6", ",", "6.596883263771078`*^-6", ",",
"4.0961454842075365`*^-6", ",", "7.367277700816571`*^-7", ",",
"4.1586361975774913`*^-7", ",", "6.697525182560524`*^-7", ",",
"3.124444926805026`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"4.031053561991353`*^9", ",", "2.0412484232184708`*^9", ",",
"6.559758574953565`*^9", ",", "2.263589413655894`*^10", ",",
"5.339827827372975`*^10", ",", "1.3641533148731873`*^12", ",",
"3.2165696857058774`*^12", ",", "2.416727496672653`*^12", ",",
"3.8922259252630664`*^12", ",", "1.2215470173936764`*^13"}], "}"}],
"8.780451624079308`*^-6","False","True"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.849162601415991*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"8cf455f7-6958-49bd-a8c9-f2edc4e5e3c0"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","4004","1.5460620259130095`*^10",
"2.0601944429362775`*^-12","0.4049675704025664`","0.4371648667768654`",
"16.44983612532279`","4",RowBox[{"{", "1", "}"}],RowBox[{"{", "4", "}"}],
RowBox[{"{",
RowBox[{
"0.00001168678111164956`", ",", "9.658496786487249`*^-6", ",",
"3.7237686220068497`*^-6", ",", "4.9563360358911175`*^-6", ",",
"5.997166603428262`*^-6", ",", "1.305158185093633`*^-6", ",",
"3.436889419485531`*^-7", ",", "6.088659256873205`*^-7", ",",
"3.436889419485531`*^-7", ",", "2.134037925554964`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"2.7532638289157796`*^9", ",", "3.9778156949058533`*^9", ",",
"3.3156104725732246`*^10", ",", "1.5460620275579931`*^10", ",",
"2.4910690913115685`*^10", ",", "4.346612862106719`*^11", ",",
"5.180317644523594`*^12", ",", "1.500597634713506`*^12", ",",
"5.180552706519197`*^12", ",", "1.3437017191328205`*^13"}], "}"}],
"4.9563360358911175`*^-6","False","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491626027065363`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"9fbef24a-5ce8-44e4-9e39-2fc79f7cd41e"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","5005","7.937313809758326`*^9",
"1.7011279223098866`*^-10","0.38173239858256974`","0.5354487290216836`",
"15.34359302452774`","5",RowBox[{"{",
RowBox[{"1", ",", "2", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{
"4.9563360358911175`*^-6", ",", "5.997166603428265`*^-6", ",",
"7.982228749163016`*^-6", ",", "0.00001414100514509597`", ",",
"7.98222874916302`*^-6", ",", "8.914406017988065`*^-7", ",",
"5.535144778975644`*^-7", ",", "4.5744998173352423`*^-7", ",",
"4.158636197577493`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"1.6838711087238682`*^10", ",", "8.526801199637009`*^9", ",",
"2.4910672225234646`*^10", ",", "6.556812245278325`*^9", ",",
"7.937313825101918`*^9", ",", "8.470318811373715`*^11", ",",
"2.6583220543087056`*^12", ",", "2.6584002463489307`*^12", ",",
"3.8922259252640493`*^12", ",", "2.380449561247762`*^13"}], "}"}],
"7.98222874916302`*^-6","True","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491626039532022`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"913858ae-0be2-4998-a23e-a5513c400535"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","6006","8.470318811215796`*^11",
"3.0601468029385494`*^-9","0.667910235907535`","0.5751963440817119`",
"15.79182268444168`","6",RowBox[{"{",
RowBox[{"1", ",", "2", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{
"4.9563360358911175`*^-6", ",", "5.997166603428265`*^-6", ",",
"7.982228749163016`*^-6", ",", "0.00001414100514509597`", ",",
"7.98222874916302`*^-6", ",", "8.914406017988065`*^-7", ",",
"5.535144778975644`*^-7", ",", "4.5744998173352423`*^-7", ",",
"4.158636197577493`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"1.6838711087238682`*^10", ",", "8.526801199637009`*^9", ",",
"2.4910672225234646`*^10", ",", "6.556812245278325`*^9", ",",
"7.937313825101918`*^9", ",", "8.470318811373715`*^11", ",",
"2.6583220543087056`*^12", ",", "2.6584002463489307`*^12", ",",
"3.8922259252640493`*^12", ",", "2.380449561247762`*^13"}], "}"}],
"8.914406017988065`*^-7","True","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491626050063853`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"6847d478-e081-4450-b33c-6ddd3f5ee3ff"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","7007","2.6583220542969897`*^12",
"1.9146491144164694`*^-7","0.6598382633469563`","0.5714845948133904`",
"11.715606492264893`","7",RowBox[{"{",
RowBox[{"1", ",", "2", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{
"4.9563360358911175`*^-6", ",", "5.997166603428265`*^-6", ",",
"7.982228749163016`*^-6", ",", "0.00001414100514509597`", ",",
"7.98222874916302`*^-6", ",", "8.914406017988065`*^-7", ",",
"5.535144778975644`*^-7", ",", "4.5744998173352423`*^-7", ",",
"4.158636197577493`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"1.6838711087238682`*^10", ",", "8.526801199637009`*^9", ",",
"2.4910672225234646`*^10", ",", "6.556812245278325`*^9", ",",
"7.937313825101918`*^9", ",", "8.470318811373715`*^11", ",",
"2.6583220543087056`*^12", ",", "2.6584002463489307`*^12", ",",
"3.8922259252640493`*^12", ",", "2.380449561247762`*^13"}], "}"}],
"5.535144778975644`*^-7","True","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.849162606044608*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"8fb2dbbf-b573-4336-a52b-098e8ab8b9a3"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","8008","2.6584002463323945`*^12",
"0.0002822542818737737`","0.4603409569900581`","0.5009519417739828`",
"16.535967709876026`","8",RowBox[{"{",
RowBox[{"1", ",", "2", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{
"4.9563360358911175`*^-6", ",", "5.997166603428265`*^-6", ",",
"7.982228749163016`*^-6", ",", "0.00001414100514509597`", ",",
"7.98222874916302`*^-6", ",", "8.914406017988065`*^-7", ",",
"5.535144778975644`*^-7", ",", "4.5744998173352423`*^-7", ",",
"4.158636197577493`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"1.6838711087238682`*^10", ",", "8.526801199637009`*^9", ",",
"2.4910672225234646`*^10", ",", "6.556812245278325`*^9", ",",
"7.937313825101918`*^9", ",", "8.470318811373715`*^11", ",",
"2.6583220543087056`*^12", ",", "2.6584002463489307`*^12", ",",
"3.8922259252640493`*^12", ",", "2.380449561247762`*^13"}], "}"}],
"4.5744998173352423`*^-7","True","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.849162607107765*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"5c10645f-d18b-4603-9745-0da6acb538a5"],
Cell[BoxData[
TemplateBox[{
" ","\" \"","9009","3.8922259252440825`*^12",
"0.009028426783095525`","0.13536917954474287`","0.1778259957944823`",
"19.96703153117168`","9",RowBox[{"{",
RowBox[{"1", ",", "2", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{"1", ",", "4"}], "}"}],RowBox[{"{",
RowBox[{
"4.9563360358911175`*^-6", ",", "5.997166603428265`*^-6", ",",
"7.982228749163016`*^-6", ",", "0.00001414100514509597`", ",",
"7.98222874916302`*^-6", ",", "8.914406017988065`*^-7", ",",
"5.535144778975644`*^-7", ",", "4.5744998173352423`*^-7", ",",
"4.158636197577493`*^-7", ",", "1.6033342791547435`*^-7"}], "}"}],
RowBox[{"{",
RowBox[{
"1.6838711087238682`*^10", ",", "8.526801199637009`*^9", ",",
"2.4910672225234646`*^10", ",", "6.556812245278325`*^9", ",",
"7.937313825101918`*^9", ",", "8.470318811373715`*^11", ",",
"2.6583220543087056`*^12", ",", "2.6584002463489307`*^12", ",",
"3.8922259252640493`*^12", ",", "2.380449561247762`*^13"}], "}"}],
"4.158636197577493`*^-7","True","False"},
"RowWithSeparators"]], "Print",
CellChangeTimes->{3.8491626082267714`*^9},
CellLabel->
"During evaluation of \
In[16]:=",ExpressionUUID->"e36c0224-fefb-457d-9fa2-4f8ff7dac1b0"]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"QS1", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QS", "[",
RowBox[{"[", "i", "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",",
RowBox[{"Length", "[", "QS", "]"}], ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS1", "=",
RowBox[{"QS1", ".",
RowBox[{"MatrixPower", "[",
RowBox[{"SIGMA", ",",
RowBox[{"-", ".5"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS2", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QS", "[",
RowBox[{"[", "i", "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "2", ",",
RowBox[{"Length", "[", "QS", "]"}], ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS2", "=",
RowBox[{"QS2", ".",
RowBox[{"MatrixPower", "[",
RowBox[{"SIGMA", ",",
RowBox[{"-", ".5"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS3", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QS", "[",
RowBox[{"[", "i", "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "3", ",",
RowBox[{"Length", "[", "QS", "]"}], ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"QS3", "=",
RowBox[{"QS3", ".",
RowBox[{"MatrixPower", "[",
RowBox[{"SIGMA", ",",
RowBox[{"-", ".5"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{"BarChart", "[",
RowBox[{
RowBox[{"StandardDeviation", "[", "QS1", "]"}], ",",
RowBox[{"ChartLabels", "\[Rule]",
RowBox[{"{",
RowBox[{
"1", ",", "2", ",", "3", ",", "4", ",", "5", ",", "6", ",", "7", ",",
"8", ",", "9", ",", "10"}], "}"}]}]}], "]"}]}], "Input",
CellChangeTimes->{{3.84484107463823*^9, 3.844841097001574*^9}, {
3.8459330714662647`*^9, 3.845933075878519*^9}},
CellLabel->"In[17]:=",ExpressionUUID->"2e810517-a198-4c46-8ad8-c903b244c9f2"],
Cell[BoxData[
GraphicsBox[{
{Opacity[0],
PointBox[{{0.3479427549194991, 0.}, {10.453901197192788`, 0.}}]}, {{},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{0.5460988028072107, 0.}, {1.4539011971927893`,
1.0598404157767793`}, "RoundingRadius" -> 0]},
ImageSizeCache->{{24.17139208624627,
54.08939458136197}, {-103.38733618625487`, 95.49959016371413}}],
"DelayedMouseEffectStyle"],
StatusArea[#, 1.0598404157767793`]& ,
TagBoxNote->"1.0598404157767793"],
StyleBox["1.0598404157767793`", {}, StripOnInput -> False]],
Annotation[#,
Style[1.0598404157767793`, {}], "Tooltip"]& ]},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{1.5460988028072107`, 0.}, {2.4539011971927893`,
1.0361116790303313`}, "RoundingRadius" -> 0]},
ImageSizeCache->{{56.577125621490445`,
86.49512811660614}, {-98.94565696811976, 95.49959016371413}}],
"DelayedMouseEffectStyle"],
StatusArea[#, 1.0361116790303313`]& ,
TagBoxNote->"1.0361116790303313"],
StyleBox["1.0361116790303313`", {}, StripOnInput -> False]],
Annotation[#,
Style[1.0361116790303313`, {}], "Tooltip"]& ]},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{2.5460988028072107`, 0.}, {3.4539011971927893`,
1.0684422163339689`}, "RoundingRadius" -> 0]},
ImageSizeCache->{{88.9828591567346,
118.90086165185032`}, {-104.99746988960639`, 95.49959016371413}}],
"DelayedMouseEffectStyle"],
StatusArea[#, 1.0684422163339689`]& ,
TagBoxNote->"1.0684422163339689"],
StyleBox["1.0684422163339689`", {}, StripOnInput -> False]],
Annotation[#,
Style[1.0684422163339689`, {}], "Tooltip"]& ]},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{3.5460988028072107`, 0.}, {4.453901197192789,
1.010891561338539}, "RoundingRadius" -> 0]},
ImageSizeCache->{{121.38859269197879`,
151.3065951870945}, {-94.22481262534946, 95.49959016371413}}],
"DelayedMouseEffectStyle"],
StatusArea[#, 1.010891561338539]& ,
TagBoxNote->"1.010891561338539"],
StyleBox["1.010891561338539`", {}, StripOnInput -> False]],
Annotation[#,
Style[1.010891561338539, {}], "Tooltip"]& ]},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{4.546098802807211, 0.}, {5.453901197192789,
1.0088968218579641`}, "RoundingRadius" -> 0]},
ImageSizeCache->{{153.79432622722297`,
183.71232872233867`}, {-93.85142600173204, 95.49959016371413}}],
"DelayedMouseEffectStyle"],
StatusArea[#, 1.0088968218579641`]& ,
TagBoxNote->"1.0088968218579641"],
StyleBox["1.0088968218579641`", {}, StripOnInput -> False]],
Annotation[#,
Style[1.0088968218579641`, {}], "Tooltip"]& ]},
{RGBColor[0.982864, 0.7431472, 0.3262672], EdgeForm[{Opacity[0.63],
Thickness[Small]}],
TagBox[
TooltipBox[
TagBox[
TagBox[
DynamicBox[{
FEPrivate`If[
CurrentValue["MouseOver"],
EdgeForm[{
GrayLevel[0.5],
AbsoluteThickness[1.5],
Opacity[0.66]}], {}, {}],
RectangleBox[{5.546098802807211, 0.}, {6.453901197192789,