-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathhierarchical_weibull2.nb
13065 lines (12538 loc) · 541 KB
/
hierarchical_weibull2.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[ 553795, 13057]
NotebookOptionsPosition[ 546213, 12932]
NotebookOutlinePosition[ 546620, 12950]
CellTagsIndexPosition[ 546577, 12947]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell["ref: \[OpenCurlyDoubleQuote]Practical Applications of Bayesian \
Reliability\[CloseCurlyDoubleQuote]", "Text",
CellChangeTimes->{{3.830148608074943*^9,
3.830148616129074*^9}},ExpressionUUID->"8d909b17-d9be-478f-84f2-\
5054b734427a"],
Cell[CellGroupData[{
Cell["sample using reliability function", "Subsubsection",
CellChangeTimes->{{3.8235354491013203`*^9, 3.8235354630056562`*^9}, {
3.8235354969016848`*^9,
3.8235355290697107`*^9}},ExpressionUUID->"4ffea42d-64c0-4274-82ce-\
6eca70b40ce4"],
Cell[BoxData[{
RowBox[{
RowBox[{"Clear", "[", "\"\<Global`*\>\"", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SeedRandom", "[", "0", "]"}], ";"}]}], "Input",
CellChangeTimes->{{3.808249437390602*^9, 3.808249468789117*^9}, {
3.808249632175363*^9, 3.808249635636956*^9}, {3.808249765110931*^9,
3.808249766532976*^9}, 3.808254098303404*^9, {3.8082563792955236`*^9,
3.808256381597042*^9}, 3.808256704271912*^9, 3.8082567676558867`*^9, {
3.808257135602641*^9, 3.808257140940946*^9}, {3.80825720411176*^9,
3.808257204717432*^9}, 3.808257874664433*^9, 3.808293009163722*^9, {
3.808293569355691*^9, 3.80829356980971*^9}, 3.8083195815073442`*^9, {
3.816823365459618*^9, 3.816823365783063*^9}, {3.8172867370346537`*^9,
3.8172867391812983`*^9}},
CellLabel->
"(Debug) In[99]:=",ExpressionUUID->"af0fe3cb-4374-408d-afe0-f6d381632569"],
Cell[BoxData[{
RowBox[{
RowBox[{"SetDirectory", "[",
RowBox[{"NotebookDirectory", "[", "]"}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"<<", "\"\<Sampler1`\>\""}], ";"}]}], "Input",
CellChangeTimes->{{3.8168124021590147`*^9, 3.816812402159195*^9}, {
3.816820205778687*^9, 3.81682020614793*^9}, 3.8172501919117117`*^9,
3.817285432667527*^9, 3.817435500804758*^9, 3.817435604876099*^9, {
3.81744024760004*^9, 3.817440253709037*^9}, {3.820740856282202*^9,
3.8207409008730907`*^9}, 3.8207410047465763`*^9, 3.820741052776198*^9, {
3.8207410968162527`*^9, 3.820741100943769*^9}, {3.820980551767469*^9,
3.820980563997992*^9}, 3.820980617916573*^9, {3.8227306941465607`*^9,
3.822730694659491*^9}, {3.822914077673979*^9, 3.822914078055475*^9},
3.822914843886628*^9, 3.82291516448693*^9, 3.822915747812892*^9, {
3.8229158419818974`*^9, 3.822915843668749*^9}, {3.8229171391032887`*^9,
3.8229171405172873`*^9}, {3.823535536950664*^9, 3.8235355377496967`*^9}, {
3.824373488903586*^9, 3.824373489917871*^9}},
CellLabel->
"(Debug) In[1]:=",ExpressionUUID->"3975a13b-a9aa-4697-8b52-d582c128fc8c"],
Cell[BoxData[{
RowBox[{
RowBox[{"data", "=",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
"\"\<~/Downloads/Computer_scripts_and_data_files/Example8.2Data_ttf_\
Table8.4.txt\>\"", ",", "\"\<TSV\>\""}], "]"}], "[",
RowBox[{"[",
RowBox[{
RowBox[{"2", ";;"}], ",",
RowBox[{"2", ";;"}]}], "]"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"cenlim", "=",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
"\"\<~/Downloads/Computer_scripts_and_data_files/Example8.2Data_CenLimit_\
Table8.5.txt\>\"", ",", "\"\<TSV\>\""}], "]"}], "[",
RowBox[{"[",
RowBox[{
RowBox[{"2", ";;"}], ",",
RowBox[{"2", ";;"}]}], "]"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"censor", "=",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
"\"\<~/Downloads/Computer_scripts_and_data_files/Example8.2Data_CensorLG_\
Table8.6.txt\>\"", ",", "\"\<TSV\>\""}], "]"}], "[",
RowBox[{"[",
RowBox[{
RowBox[{"2", ";;"}], ",",
RowBox[{"2", ";;"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"censor", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"If", "[",
RowBox[{
RowBox[{
RowBox[{"censor", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[Equal]",
"\"\<FALSE\>\""}], ",", "False", ",", "True"}], "]"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "10"}], "}"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "100"}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"data1", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"If", "[",
RowBox[{
RowBox[{
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[Equal]", "\"\<NA\>\""}],
",",
RowBox[{"cenlim", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], ",",
RowBox[{"data", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "10"}], "}"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "100"}], "}"}]}], "]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.822824264333557*^9, 3.822824283825055*^9}, {
3.822824331997476*^9, 3.8228243716842823`*^9}, {3.822824407625411*^9,
3.822824408384859*^9}, {3.822883752252137*^9, 3.8228837626286097`*^9}, {
3.8228840320844316`*^9, 3.822884073395755*^9}, {3.822886790846822*^9,
3.822886791361621*^9}, {3.822889740477742*^9, 3.822889753431324*^9}, {
3.823540069197192*^9, 3.823540069412766*^9}, 3.823915544105348*^9,
3.823917270864052*^9},
CellLabel->
"(Debug) In[103]:=",ExpressionUUID->"c3c2a907-b40f-42ec-8efb-61331ab6fcbd"],
Cell[BoxData[
RowBox[{
RowBox[{"(*",
RowBox[{
RowBox[{"trueShape", "=",
RowBox[{"{",
RowBox[{
"1.05271", ",", "1.30417", ",", "1.25654", ",", "1.53451", ",",
"1.58226", ",", "1.33131", ",", "1.38997", ",", "2.15500", ",",
"1.47402", ",", "1.28420"}], "}"}]}], ";", "\n",
RowBox[{"trueScale", "=",
RowBox[{"{",
RowBox[{
"681.011", ",", "610.300", ",", "489.721", ",", "439.992", ",",
"570.550", ",", "549.467", ",", "674.867", ",", "510.555", ",",
"610.915", ",", "667.784"}], "}"}]}], ";", "\n",
RowBox[{"cenlim", "=",
RowBox[{
RowBox[{"Import", "[",
RowBox[{
"\"\<~/Downloads/Computer_scripts_and_data_files/Example8.2Data_\
CenLimit_Table8.5.txt\>\"", ",", "\"\<TSV\>\""}], "]"}], "[",
RowBox[{"[",
RowBox[{
RowBox[{"2", ";;"}], ",",
RowBox[{"2", ";;"}]}], "]"}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{"censor", "=",
RowBox[{"Table", "[",
RowBox[{"False", ",", "100", ",", "10"}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"data1", "=",
RowBox[{"Table", "[",
RowBox[{"0", ",", "100", ",", "10"}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"Do", "[",
RowBox[{
RowBox[{"Do", "[",
RowBox[{
RowBox[{
RowBox[{"v", "=",
RowBox[{"RandomVariate", "[",
RowBox[{"WeibullDistribution", "[",
RowBox[{
RowBox[{"trueShape", "[",
RowBox[{"[", "j", "]"}], "]"}], ",",
RowBox[{"trueScale", "[",
RowBox[{"[", "j", "]"}], "]"}]}], "]"}], "]"}]}], ";",
RowBox[{"chmc", "-",
RowBox[{"mcmcIf", "[",
RowBox[{
RowBox[{"v", ">",
RowBox[{"cenlim", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",",
RowBox[{
RowBox[{
RowBox[{"censor", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "True"}], ";",
RowBox[{
RowBox[{"data1", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], "=",
RowBox[{"cenlim", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}], ",",
RowBox[{
RowBox[{"data1", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", "v"}]}], "]"}]}]}],
",",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "10"}], "}"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "100"}], "}"}]}], "]"}]}], "*)"}],
"\[IndentingNewLine]"}]], "Input",
CellChangeTimes->{{3.823915527278084*^9, 3.823915575855986*^9},
3.823917265985948*^9,
3.830148868936532*^9},ExpressionUUID->"35f1697c-3cf1-44ac-9fa0-\
4fdcfe3148df"],
Cell[BoxData[{
RowBox[{
RowBox[{"\[CapitalAlpha]", "=",
RowBox[{"{",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8", ",",
"\[Alpha]9", ",", "\[Alpha]10"}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[CapitalBeta]", "=",
RowBox[{"{",
RowBox[{
"\[Beta]1", ",", "\[Beta]2", ",", "\[Beta]3", ",", "\[Beta]4", ",",
"\[Beta]5", ",", "\[Beta]6", ",", "\[Beta]7", ",", "\[Beta]8", ",",
"\[Beta]9", ",", "\[Beta]10"}], "}"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.8228863881708517`*^9, 3.8228864600215187`*^9}},
CellLabel->
"(Debug) In[109]:=",ExpressionUUID->"c63b7b28-8467-47e7-857d-7692d5b2181c"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"U0", "[",
RowBox[{
"\[Alpha]1_", ",", "\[Alpha]2_", ",", "\[Alpha]3_", ",", "\[Alpha]4_",
",", "\[Alpha]5_", ",", "\[Alpha]6_", ",", "\[Alpha]7_", ",",
"\[Alpha]8_", ",", "\[Alpha]9_", ",", "\[Alpha]10_", ",", "\[Beta]1_",
",", "\[Beta]2_", ",", "\[Beta]3_", ",", "\[Beta]4_", ",", "\[Beta]5_",
",", "\[Beta]6_", ",", "\[Beta]7_", ",", "\[Beta]8_", ",", "\[Beta]9_",
",", "\[Beta]10_", ",", "a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}],
"=",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"PowerExpand", "[",
RowBox[{"-",
RowBox[{"Log", "[",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"If", "[",
RowBox[{
RowBox[{"censor", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}], ",",
RowBox[{"1", "-",
RowBox[{"CDF", "[",
RowBox[{
RowBox[{"WeibullDistribution", "[",
RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], ",", "y"}],
"]"}]}], ",",
RowBox[{"PDF", "[",
RowBox[{
RowBox[{"WeibullDistribution", "[",
RowBox[{"\[Alpha]", ",", "\[Beta]"}], "]"}], ",", "y"}],
"]"}]}], "]"}], ",",
RowBox[{"y", ">", "0"}]}], "]"}], "]"}]}], "]"}], "/.",
RowBox[{"{",
RowBox[{
RowBox[{"y", "\[Rule]",
RowBox[{"data1", "[",
RowBox[{"[",
RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",",
RowBox[{"\[Alpha]", "\[Rule]",
RowBox[{"\[CapitalAlpha]", "[",
RowBox[{"[", "j", "]"}], "]"}]}], ",",
RowBox[{"\[Beta]", "\[Rule]",
RowBox[{"\[CapitalBeta]", "[",
RowBox[{"[", "j", "]"}], "]"}]}]}], "}"}]}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "100"}], "}"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "10"}], "}"}]}], "]"}], "+",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"PowerExpand", "[",
RowBox[{"-",
RowBox[{"Log", "[",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"PDF", "[",
RowBox[{
RowBox[{"GammaDistribution", "[",
RowBox[{"a", ",", "b1"}], "]"}], ",", "\[Alpha]"}], "]"}],
",",
RowBox[{"\[Alpha]", ">", "0"}]}], "]"}], "]"}]}], "]"}], "+",
RowBox[{"PowerExpand", "[",
RowBox[{"-",
RowBox[{"Log", "[",
RowBox[{"Simplify", "[",
RowBox[{
RowBox[{"PDF", "[",
RowBox[{
RowBox[{"GammaDistribution", "[",
RowBox[{"c", ",", "d1"}], "]"}], ",", "\[Beta]"}], "]"}],
",",
RowBox[{"\[Beta]", ">", "0"}]}], "]"}], "]"}]}], "]"}]}],
")"}], "/.",
RowBox[{"{",
RowBox[{
RowBox[{"b1", "\[Rule]",
RowBox[{"1", "/", "b"}]}], ",",
RowBox[{"d1", "\[Rule]",
RowBox[{"1", "/", "d"}]}], ",",
RowBox[{"\[Alpha]", "\[Rule]",
RowBox[{"\[CapitalAlpha]", "[",
RowBox[{"[", "j", "]"}], "]"}]}], ",",
RowBox[{"\[Beta]", "\[Rule]",
RowBox[{"\[CapitalBeta]", "[",
RowBox[{"[", "j", "]"}], "]"}]}]}], "}"}]}], ",",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "10"}], "}"}]}], "]"}]}], "]"}]}],
";"}]], "Input",
CellChangeTimes->{
3.823786041680241*^9, {3.823788807014958*^9, 3.8237888095093*^9}, {
3.8237888558039503`*^9, 3.8237888755486*^9}},
CellLabel->
"(Debug) In[111]:=",ExpressionUUID->"442ea949-53a0-4a23-b80b-4634aa541f2b"],
Cell[BoxData[""], "Input",
CellChangeTimes->{3.824373176063591*^9, 3.824378270252618*^9},
CellLabel->
"(Debug) In[112]:=",ExpressionUUID->"0653d4d1-80d5-4458-b653-abfa6cd949e9"],
Cell[BoxData[
RowBox[{
RowBox[{"mle", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"Last", "[",
RowBox[{"FindMinimum", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"U0", "[",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4",
",", "\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",",
"\[Alpha]8", ",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1",
",", "\[Beta]2", ",", "\[Beta]3", ",", "\[Beta]4", ",",
"\[Beta]5", ",", "\[Beta]6", ",", "\[Beta]7", ",", "\[Beta]8",
",", "\[Beta]9", ",", "\[Beta]10", ",", "a", ",", "b", ",", "c",
",", "d"}], "]"}], ",",
RowBox[{
RowBox[{"\[Alpha]1", ">", "0"}], "&&",
RowBox[{"\[Alpha]2", ">", "0"}], "&&",
RowBox[{"\[Alpha]3", ">", "0"}], "&&",
RowBox[{"\[Alpha]4", ">", "0"}], "&&",
RowBox[{"\[Alpha]5", ">", "0"}], "&&",
RowBox[{"\[Alpha]6", ">", "0"}], "&&",
RowBox[{"\[Alpha]7", ">", "0"}], "&&",
RowBox[{"\[Alpha]8", ">", "0"}], "&&",
RowBox[{"\[Alpha]9", ">", "0"}], "&&",
RowBox[{"\[Alpha]10", ">", "0"}], "&&",
RowBox[{"\[Alpha]1", "<", "3"}], "&&",
RowBox[{"\[Alpha]2", "<", "3"}], "&&",
RowBox[{"\[Alpha]3", "<", "3"}], "&&",
RowBox[{"\[Alpha]4", "<", "3"}], "&&",
RowBox[{"\[Alpha]5", "<", "3"}], "&&",
RowBox[{"\[Alpha]6", "<", "3"}], "&&",
RowBox[{"\[Alpha]7", "<", "3"}], "&&",
RowBox[{"\[Alpha]8", "<", "3"}], "&&",
RowBox[{"\[Alpha]9", "<", "3"}], "&&",
RowBox[{"\[Alpha]10", "<", "3"}], "&&",
RowBox[{"\[Beta]1", ">", "0"}], "&&",
RowBox[{"\[Beta]2", ">", "0"}], "&&",
RowBox[{"\[Beta]3", ">", "0"}], "&&",
RowBox[{"\[Beta]4", ">", "0"}], "&&",
RowBox[{"\[Beta]5", ">", "0"}], "&&",
RowBox[{"\[Beta]6", ">", "0"}], "&&",
RowBox[{"\[Beta]7", ">", "0"}], "&&",
RowBox[{"\[Beta]8", ">", "0"}], "&&",
RowBox[{"\[Beta]9", ">", "0"}], "&&",
RowBox[{"\[Beta]10", ">", "0"}], "&&",
RowBox[{"\[Beta]1", "<", "1000"}], "&&",
RowBox[{"\[Beta]2", "<", "1000"}], "&&",
RowBox[{"\[Beta]3", "<", "1000"}], "&&",
RowBox[{"\[Beta]4", "<", "1000"}], "&&",
RowBox[{"\[Beta]5", "<", "1000"}], "&&",
RowBox[{"\[Beta]6", "<", "1000"}], "&&",
RowBox[{"\[Beta]7", "<", "1000"}], "&&",
RowBox[{"\[Beta]8", "<", "1000"}], "&&",
RowBox[{"\[Beta]9", "<", "1000"}], "&&",
RowBox[{"\[Beta]10", "<", "1000"}], "&&",
RowBox[{"a", ">", "0"}], "&&",
RowBox[{"b", ">", "0"}], "&&",
RowBox[{"c", ">", "0"}], "&&",
RowBox[{"d", ">", "0"}], "&&",
RowBox[{"a", "<",
RowBox[{"90", "-",
RowBox[{"10", "ii"}]}]}], "&&",
RowBox[{"b", "<",
RowBox[{"90", "-",
RowBox[{"10", "ii"}]}]}], "&&",
RowBox[{"c", "<",
RowBox[{"90", "-", " ",
RowBox[{"10", "ii"}]}]}], "&&",
RowBox[{"d", "<",
RowBox[{"90", "-",
RowBox[{"10", "ii"}]}]}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4",
",", "\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",",
"\[Alpha]8", ",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1",
",", "\[Beta]2", ",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5",
",", "\[Beta]6", ",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9",
",", "\[Beta]10", ",", "a", ",", "b", ",", "c", ",", "d"}], "}"}]}],
"]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"ii", ",", "1", ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}]], "Input",
CellChangeTimes->{{3.8236544957929697`*^9, 3.823654513328127*^9}, {
3.823654666691367*^9, 3.823654679594281*^9}, {3.82377130758873*^9,
3.823771310034842*^9}, 3.8237860476160107`*^9, {3.8238158741916523`*^9,
3.823815892102949*^9}, {3.8238159719348593`*^9, 3.823816007174426*^9}, {
3.8238242987974997`*^9, 3.8238244733954363`*^9}, 3.823824514838274*^9, {
3.823824545828907*^9, 3.823824554284561*^9}, 3.8239017021846943`*^9, {
3.824373045562697*^9, 3.824373050858698*^9}},
CellLabel->
"(Debug) In[113]:=",ExpressionUUID->"3d2d32fa-4b51-4786-bc85-bbc0d62faaa3"],
Cell[BoxData[
RowBox[{
RowBox[{"qinit", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8",
",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2",
",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6",
",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10",
",", "a", ",", "b", ",", "c", ",", "d"}], "}"}], "/.",
RowBox[{"mle", "[",
RowBox[{"[", "ii", "]"}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"ii", ",", "1", ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}]], "Input",
CellChangeTimes->{{3.823654536318873*^9, 3.82365460462508*^9}, {
3.8236547245932293`*^9, 3.82365475322642*^9}, {3.8237915465482693`*^9,
3.8237915573204927`*^9}, {3.823791593449441*^9, 3.823791623601123*^9},
3.823792386723074*^9, {3.82379374700254*^9, 3.823793747598248*^9}, {
3.8238160193823442`*^9, 3.8238160226450443`*^9}, {3.823824486601076*^9,
3.823824497415742*^9}, {3.8238245825056562`*^9, 3.823824583256405*^9}, {
3.823917223318974*^9, 3.8239172493364964`*^9}},
CellLabel->
"(Debug) In[114]:=",ExpressionUUID->"301b8c47-ea19-4332-9d0d-0b0c28f25f05"],
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.8239168526241417`*^9, 3.823916891867565*^9},
3.823917245069109*^9},
CellLabel->
"(Debug) In[115]:=",ExpressionUUID->"88f77ce9-4313-4a60-999e-6c387e740874"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"outbnd", "[", "q_", "]"}], ":=",
RowBox[{
RowBox[{"AnyTrue", "[",
RowBox[{
RowBox[{"q", "[",
RowBox[{"[",
RowBox[{"1", ";;", "20"}], "]"}], "]"}], ",",
RowBox[{
RowBox[{"#", "<=", "0.0"}], "&"}]}], "]"}], "||",
RowBox[{"AnyTrue", "[",
RowBox[{
RowBox[{"q", "[",
RowBox[{"[",
RowBox[{"21", ";;", "24"}], "]"}], "]"}], ",",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"#", "<=", "0.0"}], "||",
RowBox[{"#", "\[GreaterEqual]", "100"}]}], ")"}], "&"}]}], "]"}]}]}],
";"}]], "Input",
CellChangeTimes->{
3.823558739081317*^9, {3.8235755686792316`*^9, 3.8235755729419727`*^9},
3.823575923497136*^9, {3.823575985728202*^9, 3.8235760060088053`*^9},
3.823576657101342*^9, {3.8235812970246277`*^9, 3.8235813570124207`*^9},
3.823582013177043*^9, {3.8235825046157827`*^9, 3.8235825057342443`*^9}, {
3.8236484158028316`*^9, 3.823648418412325*^9}, 3.823648448829712*^9, {
3.823648725693712*^9, 3.823648726220564*^9}, {3.8236488434030857`*^9,
3.823648864306568*^9}, {3.823648985790164*^9, 3.823648986532823*^9},
3.823649278942437*^9, {3.823649580126604*^9, 3.8236495806633587`*^9}, {
3.823652378503459*^9, 3.823652409676248*^9}, {3.823652624559615*^9,
3.823652639855153*^9}, 3.823653535056409*^9, {3.82365413822464*^9,
3.8236542114634657`*^9}, 3.823654775656323*^9, {3.8236564571115*^9,
3.823656457742106*^9}, {3.823657043207876*^9, 3.823657046622532*^9}, {
3.823657767295886*^9, 3.823657776342586*^9}, {3.823657808326536*^9,
3.8236578260144777`*^9}, {3.823660701009714*^9, 3.823660711032131*^9}, {
3.823742415822846*^9, 3.823742424630402*^9}, {3.823766646979649*^9,
3.82376666039491*^9}, {3.823766921255192*^9, 3.8237669393590717`*^9}, {
3.823766987046273*^9, 3.8237670017488194`*^9}, {3.823767110316333*^9,
3.8237671112842903`*^9}, {3.823767271891533*^9, 3.823767272089736*^9}, {
3.823767442048998*^9, 3.8237674628711576`*^9}, {3.8237681028837957`*^9,
3.823768114219171*^9}, {3.823768147636348*^9, 3.823768153276569*^9}, {
3.823768591775681*^9, 3.823768597189554*^9}, {3.823770640359736*^9,
3.8237706433599167`*^9}, {3.823770782043687*^9, 3.823770799018283*^9}, {
3.823770959016947*^9, 3.823770982353352*^9}, {3.823771060910932*^9,
3.8237710620867863`*^9}, {3.8237718749017067`*^9, 3.823771940155128*^9},
3.8237722251363783`*^9, {3.823772680097311*^9, 3.8237726806734457`*^9}, {
3.8237733654849777`*^9, 3.823773397060663*^9}, {3.8237734655238743`*^9,
3.8237734860836763`*^9}, 3.8237736040432663`*^9, {3.8237736556825542`*^9,
3.823773657088943*^9}, {3.823774024226273*^9, 3.823774073336605*^9},
3.823774248381708*^9, 3.823774736730171*^9, 3.82377527653802*^9, {
3.823775910575478*^9, 3.823775915206647*^9}, {3.823775991997881*^9,
3.8237760003953047`*^9}, {3.823776455160084*^9, 3.823776460870358*^9}, {
3.8237765276680098`*^9, 3.82377654944508*^9}, {3.8237770785244102`*^9,
3.823777080141191*^9}, {3.823777523407855*^9, 3.823777536326314*^9},
3.823780166247161*^9, 3.823780628398589*^9, 3.82378110651502*^9,
3.823781884532795*^9, {3.823784611612419*^9, 3.823784618130858*^9}, {
3.823785307713908*^9, 3.8237853088802757`*^9}, {3.823786832685403*^9,
3.823786838191695*^9}, 3.8237915409716263`*^9, {3.823792392206279*^9,
3.823792398836319*^9}, 3.823792719033441*^9, {3.823793057052104*^9,
3.823793064448214*^9}, 3.823813620317466*^9, 3.823814037373522*^9, {
3.823816044373618*^9, 3.823816047020776*^9}, 3.823920410190888*^9,
3.823927376223477*^9, 3.823928797790786*^9, {3.823936236533553*^9,
3.8239362369396057`*^9}, 3.8239446157242107`*^9, 3.823946657265383*^9},
CellLabel->
"(Debug) In[116]:=",ExpressionUUID->"fb599c1a-2a95-4211-86d9-17ccdc4acb9f"],
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"U", "[",
RowBox[{
"\[Alpha]1_", ",", "\[Alpha]2_", ",", "\[Alpha]3_", ",", "\[Alpha]4_",
",", "\[Alpha]5_", ",", "\[Alpha]6_", ",", "\[Alpha]7_", ",",
"\[Alpha]8_", ",", "\[Alpha]9_", ",", "\[Alpha]10_", ",", "\[Beta]1_",
",", "\[Beta]2_", ",", "\[Beta]3_", ",", "\[Beta]4_", ",", "\[Beta]5_",
",", "\[Beta]6_", ",", "\[Beta]7_", ",", "\[Beta]8_", ",", "\[Beta]9_",
",", "\[Beta]10_", ",", "a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}],
"=",
RowBox[{"U0", "[",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8", ",",
"\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2", ",",
"\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6", ",",
"\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10", ",", "a",
",", "b", ",", "c", ",", "d"}], "]"}]}], ";"}]], "Input",
CellChangeTimes->{{3.82378609179562*^9, 3.8237861200783157`*^9}, {
3.823786810491288*^9, 3.8237868107215967`*^9}},
CellLabel->
"(Debug) In[117]:=",ExpressionUUID->"1b8fe7a8-cab7-4b3c-8d8c-d6c6d625fdbe"],
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"dU", "[",
RowBox[{
"\[Alpha]1_", ",", "\[Alpha]2_", ",", "\[Alpha]3_", ",", "\[Alpha]4_",
",", "\[Alpha]5_", ",", "\[Alpha]6_", ",", "\[Alpha]7_", ",",
"\[Alpha]8_", ",", "\[Alpha]9_", ",", "\[Alpha]10_", ",", "\[Beta]1_",
",", "\[Beta]2_", ",", "\[Beta]3_", ",", "\[Beta]4_", ",", "\[Beta]5_",
",", "\[Beta]6_", ",", "\[Beta]7_", ",", "\[Beta]8_", ",", "\[Beta]9_",
",", "\[Beta]10_", ",", "a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}],
"=",
RowBox[{"Simplify", "[",
RowBox[{"GradientG", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8",
",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2",
",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6",
",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10",
",", "a", ",", "b", ",", "c", ",", "d"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8",
",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2",
",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6",
",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10",
",", "a", ",", "b", ",", "c", ",", "d"}], "}"}]}], "]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"ddU", "[",
RowBox[{
"\[Alpha]1_", ",", "\[Alpha]2_", ",", "\[Alpha]3_", ",", "\[Alpha]4_",
",", "\[Alpha]5_", ",", "\[Alpha]6_", ",", "\[Alpha]7_", ",",
"\[Alpha]8_", ",", "\[Alpha]9_", ",", "\[Alpha]10_", ",", "\[Beta]1_",
",", "\[Beta]2_", ",", "\[Beta]3_", ",", "\[Beta]4_", ",", "\[Beta]5_",
",", "\[Beta]6_", ",", "\[Beta]7_", ",", "\[Beta]8_", ",", "\[Beta]9_",
",", "\[Beta]10_", ",", "a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}],
"=",
RowBox[{"Simplify", "[",
RowBox[{"HessianH", "[",
RowBox[{
RowBox[{"U", "[",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8",
",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2",
",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6",
",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10",
",", "a", ",", "b", ",", "c", ",", "d"}], "]"}], ",",
RowBox[{"{",
RowBox[{
"\[Alpha]1", ",", "\[Alpha]2", ",", "\[Alpha]3", ",", "\[Alpha]4", ",",
"\[Alpha]5", ",", "\[Alpha]6", ",", "\[Alpha]7", ",", "\[Alpha]8",
",", "\[Alpha]9", ",", "\[Alpha]10", ",", "\[Beta]1", ",", "\[Beta]2",
",", "\[Beta]3", ",", "\[Beta]4", ",", "\[Beta]5", ",", "\[Beta]6",
",", "\[Beta]7", ",", "\[Beta]8", ",", "\[Beta]9", ",", "\[Beta]10",
",", "a", ",", "b", ",", "c", ",", "d"}], "}"}]}], "]"}], "]"}]}],
";"}]}], "Input",
CellChangeTimes->{
3.822888516444146*^9, 3.822890197374745*^9, {3.822890540861368*^9,
3.8228905530694637`*^9}, 3.82289085284464*^9, {3.822891010253562*^9,
3.822891017926022*^9}, 3.822891845463046*^9, {3.8228962998120213`*^9,
3.822896329580885*^9}, {3.822896367987722*^9, 3.822896411299499*^9}, {
3.822896540836049*^9, 3.8228965616585417`*^9}, 3.822896785279277*^9, {
3.822896837557857*^9, 3.822896839981758*^9}, {3.822896995804059*^9,
3.822897000612232*^9}, {3.822897135894601*^9, 3.822897137058505*^9}, {
3.822897244163577*^9, 3.822897245493267*^9}, {3.822897595157508*^9,
3.822897597060042*^9}, {3.822908034852561*^9, 3.8229080985231*^9}, {
3.822908457816498*^9, 3.822908464269712*^9}, {3.822908505655816*^9,
3.822908508526498*^9}, {3.8229087159197683`*^9, 3.8229087498649893`*^9}, {
3.822908789853404*^9, 3.822908793218178*^9}, {3.822909009192171*^9,
3.822909011189533*^9}, {3.822913323991028*^9, 3.822913329471593*^9}, {
3.822913480286107*^9, 3.822913481811852*^9}, 3.822913647367165*^9, {
3.822913809821418*^9, 3.8229138125160313`*^9}, {3.822915138260703*^9,
3.822915140731805*^9}, {3.8229154751890993`*^9, 3.822915488596508*^9}, {
3.822915682527038*^9, 3.822915688776388*^9}, {3.822915828711012*^9,
3.822915829748595*^9}, {3.8229367594980497`*^9, 3.822936785847254*^9}, {
3.822940103880726*^9, 3.822940105644102*^9}, {3.823039808023817*^9,
3.823039833202338*^9}, {3.8230422562757597`*^9, 3.8230422614942007`*^9}, {
3.8230520470599823`*^9, 3.823052052362673*^9}, {3.8230597263122997`*^9,
3.823059726874763*^9}, {3.823062124551766*^9, 3.823062153058711*^9}, {
3.823062614585812*^9, 3.8230626279585037`*^9}, {3.823062839903404*^9,
3.823062841973922*^9}, {3.8230705413402863`*^9, 3.823070549034774*^9}, {
3.823113827587571*^9, 3.823113847203849*^9}, {3.823113992684887*^9,
3.8231139955155783`*^9}, {3.823114148921176*^9, 3.823114151436599*^9}, {
3.8231142746215143`*^9, 3.8231142748514557`*^9}, {3.82311558566929*^9,
3.823115598691505*^9}, {3.8231157946333857`*^9, 3.823115807839656*^9}, {
3.823133111919525*^9, 3.8231331339235067`*^9}, {3.823133398210232*^9,
3.8231333992119837`*^9}, {3.8235358797714*^9, 3.823535897405941*^9}, {
3.8235360003521023`*^9, 3.823536006092115*^9}, {3.823536089431979*^9,
3.8235361402566843`*^9}, {3.823536401458234*^9, 3.823536444977562*^9}, {
3.823536538572878*^9, 3.823536539284709*^9}, {3.8235400077077923`*^9,
3.823540023860713*^9}, {3.823557450862138*^9, 3.823557451824428*^9}, {
3.8235575009821873`*^9, 3.823557502013237*^9}, {3.823557674687855*^9,
3.823557677150516*^9}, {3.823557898280394*^9, 3.823557898575472*^9}, {
3.823558655196516*^9, 3.823558658513027*^9}, {3.823558698901227*^9,
3.8235587016011257`*^9}, {3.823573832601261*^9, 3.8235738418717117`*^9}, {
3.823575224195706*^9, 3.82357524428261*^9}, {3.823582035404727*^9,
3.8235820368366327`*^9}, {3.823582437394333*^9, 3.823582442166894*^9}, {
3.823582763988064*^9, 3.823582765945578*^9}, {3.823646679520197*^9,
3.823646679695684*^9}, 3.8236469621539373`*^9, 3.8236474901953907`*^9},
CellLabel->
"(Debug) In[118]:=",ExpressionUUID->"0678a4c5-b7f3-4748-9226-b21e382aec6f"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"INTERVAL", "=", "1001"}], ";"}], "\[IndentingNewLine]",
"\[IndentingNewLine]",
RowBox[{
RowBox[{"QS", "=",
RowBox[{"hmc", "[",
RowBox[{
"U", ",", "dU", ",", "ddU", ",", "24", ",", "20000", ",", "50000", ",",
"True", ",", "True", ",", "qinit"}], "]"}]}], ";"}]}]}]], "Input",
CellChangeTimes->{{3.823653542656541*^9, 3.8236535480786953`*^9}, {
3.823657052150523*^9, 3.823657054958152*^9}, {3.823657182657817*^9,
3.823657185254086*^9}, {3.8236611176527576`*^9, 3.8236611275884943`*^9}, {
3.823662239647232*^9, 3.8236622403978987`*^9}, {3.823742300775391*^9,
3.823742302015888*^9}, {3.823757747144318*^9, 3.823757754684988*^9}, {
3.823758543613056*^9, 3.823758562251978*^9}, 3.82375955835699*^9, {
3.8237607849891787`*^9, 3.82376079867185*^9}, {3.8237610896923523`*^9,
3.823761092429503*^9}, {3.823761240067227*^9, 3.82376124194127*^9}, {
3.823765773897637*^9, 3.8237657750242434`*^9}, 3.823765997536989*^9, {
3.8237665269989758`*^9, 3.823766537844503*^9}, {3.8237665708931913`*^9,
3.823766597832539*^9}, {3.823766722350765*^9, 3.82376675337549*^9}, {
3.8237676213318768`*^9, 3.823767623620388*^9}, {3.823768003319765*^9,
3.8237680284460993`*^9}, {3.82377067186095*^9, 3.823770679359745*^9}, {
3.8237708524580793`*^9, 3.823770853072996*^9}, {3.823771089739222*^9,
3.823771113037861*^9}, {3.82377138586124*^9, 3.823771390331168*^9},
3.823771963874886*^9, {3.8237723978886147`*^9, 3.8237724142596493`*^9}, {
3.823772693493894*^9, 3.8237726974072733`*^9}, {3.8237738728191767`*^9,
3.823773876403441*^9}, 3.823780645897386*^9, {3.823781127881515*^9,
3.8237811281662617`*^9}, {3.823781498584598*^9, 3.8237814991373987`*^9}, {
3.8237818990126543`*^9, 3.823781899179173*^9}, {3.823782372007275*^9,
3.823782376172927*^9}, {3.8237824625948772`*^9, 3.8237824716668158`*^9}, {
3.823782958280595*^9, 3.823782967670479*^9}, 3.8237837145709553`*^9, {
3.823786152850913*^9, 3.823786159384758*^9}, {3.8237863546853533`*^9,
3.823786360271864*^9}, 3.823786796816478*^9, {3.8237868741275997`*^9,
3.823786877918337*^9}, {3.823787141748961*^9, 3.823787155030878*^9},
3.823787356930394*^9, {3.823787999696252*^9, 3.823788024528076*^9}, {
3.823788059137882*^9, 3.8237880819180326`*^9}, {3.823788580725322*^9,
3.8237885878532352`*^9}, {3.823789581947536*^9, 3.82378961161458*^9}, {
3.823789773849523*^9, 3.823789774108107*^9}, {3.823790353115671*^9,
3.823790371254569*^9}, {3.823790827351779*^9, 3.8237908571908503`*^9}, {
3.8237909258157797`*^9, 3.82379092661047*^9}, {3.8237909896067333`*^9,
3.82379100087175*^9}, {3.8237910760052156`*^9, 3.8237910778091*^9}, {
3.823791158179996*^9, 3.823791164374681*^9}, 3.823791884988154*^9, {
3.8237923077675447`*^9, 3.823792313389432*^9}, 3.823792729718523*^9,
3.8237931256492567`*^9, 3.823793208591569*^9, {3.8237936927787523`*^9,
3.8237936948243113`*^9}, {3.8238132948136997`*^9, 3.823813300796671*^9}, {
3.823813637785597*^9, 3.823813647172495*^9}, {3.823813685909141*^9,
3.823813785716745*^9}, {3.823813887415691*^9, 3.823813907932822*^9}, {
3.823814042738749*^9, 3.823814046884782*^9}, 3.823814087716757*^9,
3.82381412755676*^9, {3.823816091327276*^9, 3.8238161096041*^9},
3.8238174849674253`*^9, {3.8238199279435663`*^9, 3.823819928294438*^9}, {
3.823901629872361*^9, 3.8239016308926067`*^9}, {3.823916763656814*^9,
3.823916773213578*^9}, {3.823917258690255*^9, 3.823917279613681*^9}, {
3.8239180553483057`*^9, 3.823918060878334*^9}, 3.823920403175095*^9,
3.823920437069236*^9, {3.823940402255073*^9, 3.823940416229658*^9},
3.823940510669094*^9, {3.8239407118775463`*^9, 3.823940742653534*^9}, {
3.8239407811814127`*^9, 3.823940781429578*^9}, {3.8239408652771482`*^9,
3.8239408654292583`*^9}, {3.8239446223945093`*^9, 3.8239446230748*^9},
3.824372794320005*^9, {3.8243728574213467`*^9, 3.824372857716372*^9}, {
3.824373022822405*^9, 3.8243730229662037`*^9}, {3.824373184409905*^9,
3.824373185064501*^9}, {3.824373216776923*^9, 3.824373237027176*^9},
3.824378023860491*^9, {3.824378180679962*^9, 3.824378180846188*^9},
3.824378234422249*^9, {3.8243782740153513`*^9, 3.8243782759801083`*^9}, {
3.824378825843894*^9, 3.82437882703001*^9},
3.83014888429713*^9},ExpressionUUID->"4dcce1fe-7175-44ec-8938-\
4a221682106a"],
Cell[BoxData[
TemplateBox[{
"LinearSolve","luc",
"\"Result for \\!\\(\\*RowBox[{\\\"LinearSolve\\\"}]\\) of badly \
conditioned matrix \\!\\(\\*RowBox[{\\\"{\\\", RowBox[{RowBox[{\\\"{\\\", \
RowBox[{\\\"79.48059005384266`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.050046415386488796`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", RowBox[{\\\"-\\\", \
\\\"0.7589749758595121`\\\"}], \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"\
\[LeftSkeleton]\\\", \\\"22\\\", \\\"\[RightSkeleton]\\\"}], \\\",\\\", \
RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\
\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\
\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\
\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", RowBox[{\\\"-\\\", \\\"67.65900844889494`\\\"}], \
\\\",\\\", \\\"31802.25288080832`\\\"}], \\\"}\\\"}]}], \\\"}\\\"}]\\) may \
contain significant numerical errors.\"",2,121,27,19644157943612182906,
"Local"},
"MessageTemplate"]], "Message", "MSG",
CellChangeTimes->{
3.824373063410585*^9, 3.8243731904754057`*^9, {3.824373220478573*^9,
3.8243732427551327`*^9}, 3.8243732985113487`*^9, 3.824373396352665*^9,
3.824373497381762*^9, 3.824373551960793*^9, 3.82437803190342*^9,
3.824378140017164*^9, 3.8243781846906967`*^9, 3.824378238949347*^9,
3.824378326153111*^9, 3.824378879150226*^9},
CellLabel->
"(Debug) During evaluation of \
In[120]:=",ExpressionUUID->"3825a398-a346-4089-9974-f27b301b75c4"],
Cell[BoxData[
TemplateBox[{
"LinearSolve","luc",
"\"Result for \\!\\(\\*RowBox[{\\\"LinearSolve\\\"}]\\) of badly \
conditioned matrix \\!\\(\\*RowBox[{\\\"{\\\", RowBox[{RowBox[{\\\"{\\\", \
RowBox[{\\\"79.48059005384266`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.050046415386488796`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", RowBox[{\\\"-\\\", \
\\\"0.7589749758595121`\\\"}], \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"\
\[LeftSkeleton]\\\", \\\"22\\\", \\\"\[RightSkeleton]\\\"}], \\\",\\\", \
RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\
\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\
\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\
\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", RowBox[{\\\"-\\\", \\\"67.65900844889494`\\\"}], \
\\\",\\\", \\\"31802.25288080832`\\\"}], \\\"}\\\"}]}], \\\"}\\\"}]\\) may \
contain significant numerical errors.\"",2,121,28,19644157943612182906,
"Local"},
"MessageTemplate"]], "Message", "MSG",
CellChangeTimes->{
3.824373063410585*^9, 3.8243731904754057`*^9, {3.824373220478573*^9,
3.8243732427551327`*^9}, 3.8243732985113487`*^9, 3.824373396352665*^9,
3.824373497381762*^9, 3.824373551960793*^9, 3.82437803190342*^9,
3.824378140017164*^9, 3.8243781846906967`*^9, 3.824378238949347*^9,
3.824378326153111*^9, 3.824378879278075*^9},
CellLabel->
"(Debug) During evaluation of \
In[120]:=",ExpressionUUID->"e012b6de-a7cd-4224-980c-a2e4bb1dd2fe"],
Cell[BoxData[
TemplateBox[{
"LinearSolve","luc",
"\"Result for \\!\\(\\*RowBox[{\\\"LinearSolve\\\"}]\\) of badly \
conditioned matrix \\!\\(\\*RowBox[{\\\"{\\\", RowBox[{RowBox[{\\\"{\\\", \
RowBox[{\\\"79.38536663105015`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.05000372164225661`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
RowBox[{\\\"-\\\", \\\"0.7589150188034967`\\\"}], \\\",\\\", \\\"1.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}], \\\",\\\", \
RowBox[{\\\"\[LeftSkeleton]\\\", \\\"22\\\", \\\"\[RightSkeleton]\\\"}], \
\\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\
\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\
\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \\\"1.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", RowBox[{\\\"-\\\", \
\\\"67.41301385617216`\\\"}], \\\",\\\", \\\"31671.154502691166`\\\"}], \\\"}\
\\\"}]}], \\\"}\\\"}]\\) may contain significant numerical errors.\"",2,121,
29,19644157943612182906,"Local"},
"MessageTemplate"]], "Message", "MSG",
CellChangeTimes->{
3.824373063410585*^9, 3.8243731904754057`*^9, {3.824373220478573*^9,
3.8243732427551327`*^9}, 3.8243732985113487`*^9, 3.824373396352665*^9,
3.824373497381762*^9, 3.824373551960793*^9, 3.82437803190342*^9,
3.824378140017164*^9, 3.8243781846906967`*^9, 3.824378238949347*^9,
3.824378326153111*^9, 3.8243788792845173`*^9},
CellLabel->
"(Debug) During evaluation of \
In[120]:=",ExpressionUUID->"16701f0f-ba8f-4fc8-a89f-9cbfe002255f"],
Cell[BoxData[
TemplateBox[{
"General","stop",
"\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"LinearSolve\\\", \\\"::\
\\\", \\\"luc\\\"}], \\\"MessageName\\\"]\\) will be suppressed during this \
calculation.\"",2,121,30,19644157943612182906,"Local"},
"MessageTemplate"]], "Message", "MSG",
CellChangeTimes->{
3.824373063410585*^9, 3.8243731904754057`*^9, {3.824373220478573*^9,
3.8243732427551327`*^9}, 3.8243732985113487`*^9, 3.824373396352665*^9,
3.824373497381762*^9, 3.824373551960793*^9, 3.82437803190342*^9,
3.824378140017164*^9, 3.8243781846906967`*^9, 3.824378238949347*^9,
3.824378326153111*^9, 3.824378879290165*^9},
CellLabel->
"(Debug) During evaluation of \
In[120]:=",ExpressionUUID->"95838c5d-6996-4c0f-a33f-8eea924f9877"]
}, Open ]],
Cell[BoxData[
RowBox[{
RowBox[{"QQ", "=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QS", "[",
RowBox[{"[",
RowBox[{
RowBox[{"i", ";;",
RowBox[{"Length", "[", "QS", "]"}], ";;", "CHAINS"}], ",", ";;"}],
"]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "CHAINS"}], "}"}]}], "]"}]}],
";"}]], "Input",
CellChangeTimes->{{3.8238154676797733`*^9, 3.823815601564576*^9}},
CellLabel->
"(Debug) In[122]:=",ExpressionUUID->"eb80f3ab-288b-46a0-b080-d0d61556d445"],
Cell[BoxData[
RowBox[{"Table", "[",
RowBox[{
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QQ", "[",
RowBox[{"[",
RowBox[{"i", ",", ";;", ",", "j"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "CHAINS"}], "}"}]}], "]"}], ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"Opacity", "[", ".1", "]"}]}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "24"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{3.8238247427428417`*^9},
CellLabel->
"(Debug) In[123]:=",ExpressionUUID->"80b62389-6cc7-493d-89d8-2b847c9afa0b"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"Table", "[",
RowBox[{"j", ",",
RowBox[{"{",
RowBox[{"j", ",", "21", ",", "24"}], "}"}]}], "]"}]], "Input",
CellChangeTimes->{{3.823592199609214*^9, 3.823592231189793*^9}, {
3.823592333711832*^9, 3.823592341166078*^9}},
CellLabel->
"(Debug) In[124]:=",ExpressionUUID->"e935c722-ab87-4bda-8aa5-397224b12e8d"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"21", ",", "22", ",", "23", ",", "24"}], "}"}]], "Output",
CellChangeTimes->{
3.823592161409912*^9, {3.8235922005952682`*^9, 3.8235922317322702`*^9},
3.823592341690515*^9, 3.823596202897415*^9, 3.823638083486721*^9,
3.8236388000063972`*^9, 3.823640379974907*^9, 3.823641264782216*^9,
3.82364610487952*^9, 3.8236478660347652`*^9, 3.823655386555579*^9,
3.823656969333384*^9, 3.823657150935753*^9, 3.823658319381926*^9,
3.823754097495673*^9, 3.8237558077057953`*^9, 3.8237562898709803`*^9,
3.82375695854307*^9, 3.823757520245529*^9, 3.823758333752824*^9,
3.823758861936489*^9, 3.823759971316832*^9, 3.823766490092043*^9,
3.8237746452835417`*^9, 3.823775148011742*^9, 3.823775705654636*^9,
3.823789534928873*^9, 3.8237939792198143`*^9, 3.823813372960672*^9,
3.823821650593072*^9, 3.823826226450574*^9, 3.823902677513239*^9,
3.823916467915444*^9, 3.8239181962768497`*^9, 3.82392139820231*^9,
3.823928360895007*^9, 3.8239406016833067`*^9, 3.8239463967625504`*^9,
3.823948456297625*^9, 3.8243729876895123`*^9, 3.824373678058221*^9,
3.8243786603329077`*^9, 3.8243826097306957`*^9},
CellLabel->
"(Debug) Out[124]=",ExpressionUUID->"815a4b94-bb2b-4982-8467-98a9668e464e"]
}, Open ]],
Cell[BoxData[
RowBox[{"GraphicsGrid", "[", "\[IndentingNewLine]",
RowBox[{"ArrayReshape", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"ListPlot", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"QQ", "[",
RowBox[{"[",
RowBox[{"i", ",", ";;", ",", "j"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "CHAINS"}], "}"}]}], "]"}], ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"Opacity", "[", ".1", "]"}]}], ",",
RowBox[{"Ticks", "\[Rule]",