-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathSampleOneMatrix.nb
More file actions
9218 lines (9150 loc) · 451 KB
/
Copy pathSampleOneMatrix.nb
File metadata and controls
9218 lines (9150 loc) · 451 KB
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 13.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 461389, 9210]
NotebookOptionsPosition[ 459134, 9160]
NotebookOutlinePosition[ 459527, 9176]
CellTagsIndexPosition[ 459484, 9173]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}]], "Input",
CellLabel->
"6/7/24 08:54:38 \
In[265]:=",ExpressionUUID->"4cb67bbd-0bcc-480b-9211-7b5394994b7e"],
Cell[CellGroupData[{
Cell["Exact solution", "Section",
CellChangeTimes->{{3.926731675482567*^9, 3.926731681141871*^9}, {
3.926731863932665*^9,
3.926731864944477*^9}},ExpressionUUID->"db79803c-091a-42bc-9975-\
e4601f860f5c"],
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"intanalytic", "[",
RowBox[{"n_", ",", "p_", ",", " ", "\[Beta]_"}], "]"}], ":=",
RowBox[{"BesselI", "[",
RowBox[{
RowBox[{"n", "+", "p"}], ",", "\[Beta]"}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"matanalytic", "[",
RowBox[{"p_", ",", "\[Beta]_", ",", "N0_"}], "]"}], ":=",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"intanalytic", "[",
RowBox[{
RowBox[{"i", "-", "j"}], ",", "p", ",",
RowBox[{"2", "N0", "*", "\[Beta]"}]}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "1", ",", "N0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",", "1", ",", "N0"}], "}"}]}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"partitionfunc", "[",
RowBox[{"\[Beta]_", ",", "N0_"}], "]"}], ":=",
RowBox[{"Sum", "[",
RowBox[{
RowBox[{"Det", "[",
RowBox[{"matanalytic", "[",
RowBox[{"p", ",", "\[Beta]", ",", "N0"}], "]"}], "]"}], ",",
RowBox[{"{",
RowBox[{"p", ",",
RowBox[{"-", "35"}], ",", "35"}], "}"}]}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"w1analytic", "[",
RowBox[{"\[Lambda]_", ",", "N0_"}], "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{"partitionfunc", "[",
RowBox[{"\[Beta]", ",", "N0"}], "]"}], ",", "\[Beta]"}], "]"}], "/",
RowBox[{"partitionfunc", "[",
RowBox[{"\[Beta]", ",", "N0"}], "]"}]}], "/",
RowBox[{"(",
RowBox[{"2",
RowBox[{"N0", "^", "2"}]}], ")"}]}], ")"}], "/.",
RowBox[{"\[Beta]", "->",
RowBox[{"N", "[",
RowBox[{"1", "/", "\[Lambda]"}], "]"}]}]}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"partitionfuncuN", "[",
RowBox[{"\[Beta]_", ",", "N0_"}], "]"}], ":=",
RowBox[{"Det", "[",
RowBox[{"matanalytic", "[",
RowBox[{"0", ",", "\[Beta]", ",", "N0"}], "]"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"w1analyticuN", "[",
RowBox[{"\[Lambda]_", ",", "N0_"}], "]"}], ":=",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{"partitionfuncuN", "[",
RowBox[{"\[Beta]", ",", "N0"}], "]"}], ",", "\[Beta]"}], "]"}], "/",
RowBox[{"partitionfuncuN", "[",
RowBox[{"\[Beta]", ",", "N0"}], "]"}]}], "/",
RowBox[{"(",
RowBox[{"2",
RowBox[{"N0", "^", "2"}]}], ")"}]}], ")"}], "/.",
RowBox[{"\[Beta]", "->",
RowBox[{"N", "[",
RowBox[{"1", "/", "\[Lambda]"}], "]"}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"w1analyticthooft", "[", "\[Lambda]_", "]"}], "/;",
RowBox[{"0", "<=", "\[Lambda]", "<=", "2"}]}], ":=",
RowBox[{"1", "-",
RowBox[{"\[Lambda]", "/", "4"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"w1analyticthooft", "[", "\[Lambda]_", "]"}], "/;",
RowBox[{"\[Lambda]", ">", "2"}]}], ":=",
RowBox[{"1", "/", "\[Lambda]"}]}], ";"}]}], "Input",
CellChangeTimes->{{3.926692243371737*^9, 3.9266922436710377`*^9},
3.926731862654697*^9},
CellLabel->
"6/7/24 08:54:39 \
In[266]:=",ExpressionUUID->"aee34aa8-cf24-4815-92c7-0e4ff5065552"]
}, Open ]],
Cell[CellGroupData[{
Cell["Some function", "Section",
CellChangeTimes->{{3.926731997072426*^9,
3.926732007234192*^9}},ExpressionUUID->"1c4470c1-22fe-45ab-803b-\
6bf7cebcadcd"],
Cell[BoxData[{
RowBox[{
RowBox[{"lastsort", " ", "=", " ",
RowBox[{
RowBox[{"#", "~", "Extract", "~",
RowBox[{"Ordering", "[",
RowBox[{"#", ",", " ",
RowBox[{"-", "1"}]}], "]"}]}], " ", "&"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"canonize", "[",
RowBox[{"w", "[",
RowBox[{"m_", ",", "n_"}], "]"}], "]"}], ":=",
RowBox[{"w", "@@",
RowBox[{"(",
RowBox[{"lastsort", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"m", ",", "n"}], "}"}], ",",
RowBox[{"{",
RowBox[{"n", ",", "m"}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "m"}], ",",
RowBox[{"-", "n"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "n"}], ",",
RowBox[{"-", "m"}]}], "}"}]}], "}"}], "]"}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"canonize", "[", "1", "]"}], "=", "1"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"canonize", "[",
RowBox[{"w", "[", "m_", "]"}], "]"}], ":=",
RowBox[{"w", "@@",
RowBox[{"(",
RowBox[{"lastsort", "[",
RowBox[{"{",
RowBox[{
RowBox[{"{", "m", "}"}], ",",
RowBox[{"{",
RowBox[{"-", "m"}], "}"}]}], "}"}], "]"}], ")"}]}]}],
";"}]}], "Input",
CellChangeTimes->{{3.926732015240726*^9, 3.926732028087635*^9}, {
3.92673207487372*^9, 3.926732075093041*^9}},
CellLabel->
"6/7/24 08:54:41 \
In[274]:=",ExpressionUUID->"f0058aa6-3b9e-4f71-b51b-dc23d0ba7562"]
}, Open ]],
Cell[CellGroupData[{
Cell["SU(3)", "Section",
CellChangeTimes->{{3.9267316690633287`*^9,
3.926731670825782*^9}},ExpressionUUID->"cc8696ab-a102-48f5-b368-\
3147c8281fcb"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"N0", "=", "3"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"w", "[",
RowBox[{"i_", ",", "j_", ",", "k_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"-", "2"}], "/",
RowBox[{"N0", "^", "2"}]}],
RowBox[{"canonize", "@", " ",
RowBox[{"w", "[",
RowBox[{
RowBox[{"i", "+", "j", "+", "k"}], ",", "0"}], "]"}]}]}], "+",
RowBox[{
RowBox[{"1", "/", "N0"}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{"i", ",",
RowBox[{"j", "+", "k"}]}], "]"}]}]}], "+",
RowBox[{
RowBox[{"1", "/", "N0"}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{"j", ",",
RowBox[{"i", "+", "k"}]}], "]"}]}]}], "+",
RowBox[{
RowBox[{"1", "/", "N0"}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{"k", ",",
RowBox[{"j", "+", "i"}]}], "]"}]}]}], "+",
RowBox[{
RowBox[{"1", "/", "N0"}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{
RowBox[{"j", "-", "i"}], ",",
RowBox[{"k", "-", "i"}]}], "]"}]}]}], "-",
RowBox[{
RowBox[{"1", "/",
RowBox[{"N0", "^", "2"}]}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{
RowBox[{"j", "+", "k", "-",
RowBox[{"2", "i"}]}], ",", "0"}], "]"}]}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"w0", "[",
RowBox[{"i_", ",", "j_", ",", "k_"}], "]"}], ":=",
RowBox[{"w", "@@",
RowBox[{"SortBy", "[",
RowBox[{
RowBox[{"{",
RowBox[{"i", ",", "j", ",", "k"}], "}"}], ",", "Abs"}], "]"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"w", "[",
RowBox[{"0", ",", "0"}], "]"}], "=", "1"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"eqmn", "[",
RowBox[{"n_", ",", "m_"}], "]"}], "/;",
RowBox[{"n", ">=", "0"}]}], ":=",
RowBox[{
RowBox[{
RowBox[{"Sum", "[",
RowBox[{
RowBox[{"w0", "[",
RowBox[{"i", ",",
RowBox[{"n", "-", "i"}], ",", "m"}], "]"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{"n", "-", "1"}]}], "}"}]}], "]"}], "+",
RowBox[{
RowBox[{"m", " ", "/",
RowBox[{"N0", "^", "2"}]}], " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{
RowBox[{"m", "+", "n"}], ",", "0"}], "]"}]}]}], "+",
RowBox[{"\[Beta]", " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{
RowBox[{"n", "+", "1"}], ",", "m"}], "]"}]}]}], "-",
RowBox[{"\[Beta]", " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{
RowBox[{"n", "-", "1"}], ",", "m"}], "]"}]}]}], "-",
RowBox[{
RowBox[{"1", "/",
RowBox[{"N0", "^", "2"}]}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"n", " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{"n", ",", "m"}], "]"}]}]}], "+",
RowBox[{"m", " ",
RowBox[{"canonize", "@",
RowBox[{"w", "[",
RowBox[{"n", ",", "m"}], "]"}]}]}], "+",
RowBox[{"\[Beta]", " ",
RowBox[{"N0", "^", "2"}],
RowBox[{"(",
RowBox[{
RowBox[{"w0", "[",
RowBox[{"n", ",", "m", " ", ",", "1"}], "]"}], "-",
RowBox[{"w0", "[",
RowBox[{"n", ",", "m", " ", ",",
RowBox[{"-", "1"}]}], "]"}]}], ")"}]}]}], ")"}]}]}], "//",
"Expand"}]}], ";"}], "\n",
RowBox[{
RowBox[{
RowBox[{"updownpair", "[",
RowBox[{"\[Lambda]_", ",", "matrix_", ",", "eq_"}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"var1", "=",
RowBox[{
RowBox[{"Cases", "[",
RowBox[{"eq", ",",
RowBox[{"w", "[", "x__", "]"}], ",", "Infinity"}], "]"}], "//",
"Union"}]}], ",", "rule", ",", "matrixwithrule", ",", "varsdp"}],
"}"}], ",", "\[IndentingNewLine]",
RowBox[{
RowBox[{"rule", "=",
RowBox[{
RowBox[{"Reduce", "[",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"eq", "/.",
RowBox[{"\[Beta]", "->",
RowBox[{"1", "/", "\[Lambda]"}]}]}], ")"}], "==", "0"}], ",",
"var1"}], "]"}], "//", "ToRules"}]}], ";", "\[IndentingNewLine]",
RowBox[{"matrixwithrule", "=",
RowBox[{
RowBox[{"matrix", "/.", "rule"}], "//", "Expand"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"varsdp", "=",
RowBox[{
RowBox[{"Cases", "[",
RowBox[{"matrix", ",",
RowBox[{"w", "[", "__", "]"}], ",", "Infinity"}], "]"}], "//",
"Union"}]}], ";", "\[IndentingNewLine]",
RowBox[{
RowBox[{"w", "[",
RowBox[{"1", ",", "0"}], "]"}], "/.",
RowBox[{"{",
RowBox[{
RowBox[{"ConicOptimization", "[",
RowBox[{
RowBox[{"w", "[",
RowBox[{"1", ",", "0"}], "]"}], ",",
RowBox[{
RowBox[{
RowBox[{"VectorGreaterEqual", "[",
RowBox[{
RowBox[{"{",
RowBox[{"#", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\"\<SemidefiniteCone\>\"", ",", " ",
RowBox[{"Length", "[",
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], "]"}],
"&"}], "[", "matrixwithrule", "]"}], ",", "varsdp"}], "]"}], ",",
RowBox[{"ConicOptimization", "[",
RowBox[{
RowBox[{"-",
RowBox[{"w", "[",
RowBox[{"1", ",", "0"}], "]"}]}], ",",
RowBox[{
RowBox[{
RowBox[{"VectorGreaterEqual", "[",
RowBox[{
RowBox[{"{",
RowBox[{"#", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\"\<SemidefiniteCone\>\"", ",", " ",
RowBox[{"Length", "[",
RowBox[{"#", "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], "]"}],
"&"}], "[", "matrixwithrule", "]"}], ",", "varsdp"}], "]"}]}],
"}"}]}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"updowncurve", "[",
RowBox[{"\[Lambda]list_", ",", "\[CapitalLambda]_"}], "]"}], ":=",
RowBox[{"Block", "[",
RowBox[{
RowBox[{"{",
RowBox[{"eq", ",", "op", ",", "matrix"}], "}"}], ",",
"\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{
RowBox[{"var", "=",
RowBox[{
RowBox[{"canonize", "/@",
RowBox[{"(",
RowBox[{"w", "@@@",
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"i", ",", "j"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",",
RowBox[{"-", "\[CapitalLambda]"}], ",",
"\[CapitalLambda]"}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"-",
RowBox[{"Abs", "[",
RowBox[{"\[CapitalLambda]", "-",
RowBox[{"Abs", "[", "i", "]"}]}], "]"}]}], ",",
RowBox[{"\[CapitalLambda]", "-",
RowBox[{"Abs", "[", "i", "]"}]}]}], "}"}]}], "]"}], ",",
"1"}], "]"}]}], ")"}]}], "//", "Union"}]}], ";"}], "*)"}],
"\[IndentingNewLine]",
RowBox[{
RowBox[{"eq", "=",
RowBox[{"(",
RowBox[{
RowBox[{"eqmn", "@@@",
RowBox[{"(",
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"i", ",", "j"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",", "0", ",",
RowBox[{"\[CapitalLambda]", "+", "3"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"-",
RowBox[{"Abs", "[",
RowBox[{"\[CapitalLambda]", "+", "3", "-",
RowBox[{"Abs", "[", "i", "]"}]}], "]"}]}], ",",
RowBox[{"\[CapitalLambda]", "+", "3", "-",
RowBox[{"Abs", "[", "i", "]"}]}]}], "}"}]}], "]"}], ",",
"1"}], "]"}], ")"}]}], "//", "Union"}], ")"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"op", "=",
RowBox[{"Cases", "[",
RowBox[{
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"i", ",", "j"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",",
RowBox[{"-",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], ",",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"-",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], ",",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], "}"}]}],
"]"}], ",", "1"}], "]"}], ",",
RowBox[{"x_", "/;",
RowBox[{
RowBox[{"Norm", "[",
RowBox[{"x", ",", "1"}], "]"}], "<=",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}]}]}], "]"}]}], ";",
"\[IndentingNewLine]",
RowBox[{"matrix", "=",
RowBox[{"Outer", "[",
RowBox[{
RowBox[{
RowBox[{"canonize", "[",
RowBox[{"w", "@@",
RowBox[{"(",
RowBox[{"#1", "+", "#2"}], ")"}]}], "]"}], "&"}], ",", "op", ",",
RowBox[{"-", "op"}], ",", "1"}], "]"}]}], ";", "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"updownpair", "[",
RowBox[{"#", ",", "matrix", ",", "eq"}], "]"}], "&"}], "/@",
"\[Lambda]list"}], "//", "Transpose"}]}]}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"\[CapitalLambda]", "=", "2"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"op", "=",
RowBox[{"Cases", "[",
RowBox[{
RowBox[{"Flatten", "[",
RowBox[{
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"i", ",", "j"}], "}"}], ",",
RowBox[{"{",
RowBox[{"i", ",",
RowBox[{"-",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], ",",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"j", ",",
RowBox[{"-",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], ",",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}], "}"}]}], "]"}],
",", "1"}], "]"}], ",",
RowBox[{"x_", "/;",
RowBox[{
RowBox[{"Norm", "[",
RowBox[{"x", ",", "1"}], "]"}], "<=",
RowBox[{"Floor", "[",
RowBox[{"\[CapitalLambda]", "/", "2"}], "]"}]}]}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"matrix", "=",
RowBox[{"Outer", "[",
RowBox[{
RowBox[{
RowBox[{"canonize", "[",
RowBox[{"w", "@@",
RowBox[{"(",
RowBox[{"#1", "+", "#2"}], ")"}]}], "]"}], "&"}], ",", "op", ",",
RowBox[{"-", "op"}], ",", "1"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"\[Lambda]list", "=",
RowBox[{"Range", "[",
RowBox[{
RowBox[{"1", "/", "40"}], ",", "12", ",",
RowBox[{"1", "/", "40"}]}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"l2up", "=",
RowBox[{
RowBox[{"Join", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}], ",",
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}]}], "}"}], ",",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#1"}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#2"}], "}"}], "]"}]}], "}"}],
"&"}], "@@",
RowBox[{"updowncurve", "[",
RowBox[{"\[Lambda]list", ",", "2"}], "]"}]}], ")"}], ",", "2"}],
"]"}], "[",
RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"l4updown", "=",
RowBox[{"Join", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}], ",",
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}]}], "}"}], ",",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#1"}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#2"}], "}"}], "]"}]}], "}"}],
"&"}], "@@",
RowBox[{"updowncurve", "[",
RowBox[{"\[Lambda]list", ",", "4"}], "]"}]}], ")"}], ",", "2"}],
"]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"l6updown", "=",
RowBox[{"Join", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}], ",",
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}]}], "}"}], ",",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#1"}], "}"}], "]"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",", "#2"}], "}"}], "]"}]}], "}"}],
"&"}], "@@",
RowBox[{"updowncurve", "[",
RowBox[{"\[Lambda]list", ",", "6"}], "]"}]}], ")"}], ",", "2"}],
"]"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"analist", "=",
RowBox[{"Join", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}], "}"}], ",",
RowBox[{"Transpose", "[",
RowBox[{"{",
RowBox[{"\[Lambda]list", ",",
RowBox[{
RowBox[{
RowBox[{"w1analytic", "[",
RowBox[{"#", ",", "3"}], "]"}], "&"}], "/@", "\[Lambda]list"}]}],
"}"}], "]"}], ",", "1"}], "]"}]}], ";"}], "\n",
RowBox[{
RowBox[{"lowplot", "=",
RowBox[{"ListPlot", "[",
RowBox[{"l2up", ",",
RowBox[{"Joined", "\[Rule]", "True"}],
RowBox[{"(*",
RowBox[{",",
RowBox[{"PlotMarkers", "\[Rule]",
RowBox[{"{",
RowBox[{"cross", ",", "0.023"}], "}"}]}]}], "*)"}], ",",
RowBox[{"Frame", "\[Rule]", "True"}], ",",
RowBox[{"PlotStyle", "\[Rule]", "Gray"}], ",",
RowBox[{"Filling", "\[Rule]", "Axis"}], ",",
RowBox[{"PlotRange", "\[Rule]",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"0", ",", "12"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",",
RowBox[{"RotateLabel", "->", "False"}], ",",
RowBox[{"FrameStyle", "\[Rule]", "Thick"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[CapitalLambda]4pl", "=",
RowBox[{"ListPlot", "[",
RowBox[{"l4updown", ",",
RowBox[{"Joined", "->", "True"}], ",",
RowBox[{"PlotStyle", "->",
RowBox[{"{",
RowBox[{"Blue", ",", "Blue"}], "}"}]}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"\[CapitalLambda]6pl", "=",
RowBox[{"ListPlot", "[",
RowBox[{"l6updown", ",",
RowBox[{"Joined", "->", "True"}], ",",
RowBox[{"PlotStyle", "->",
RowBox[{"{",
RowBox[{"Red", ",", "Red"}], "}"}]}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"anapl", "=",
RowBox[{"ListPlot", "[",
RowBox[{"analist", ",",
RowBox[{"Joined", "->", "True"}], ",",
RowBox[{"PlotStyle", "->",
RowBox[{"{",
RowBox[{"Dashed", ",", "Black"}], "}"}]}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"su3pl", "=",
RowBox[{"Show", "[",
RowBox[{
"lowplot", ",", "\[CapitalLambda]4pl", ",", "\[CapitalLambda]6pl", ",",
"anapl"}], "]"}]}]}], "Input",
CellChangeTimes->{{3.926731740673234*^9, 3.9267317449560137`*^9},
3.9267318119967537`*^9, {3.926732011584798*^9, 3.926732029997505*^9},
3.926732086735091*^9},
CellLabel->
"6/7/24 08:54:51 \
In[278]:=",ExpressionUUID->"8d454947-9214-48ed-836d-9c0f031cc3e1"],
Cell[BoxData[
GraphicsBox[{{{},
InterpretationBox[{
TagBox[GraphicsComplexBox[CompressedData["
1:eJzt2flbzHv8//FpT4vWqaZ1ptnK0iIJp7yfWSpEOEXZkxNlK1k6liTH1iJE
FEkHZStLkoOInGhfpL20TvtepFKfjp6vH77f6/MPfK5r5hfXaOY293lcuC7X
k7Vlz6o/RGk02gcxGu2/X//fRzcVGfHrQV0KOd0SV9GFz+9Rj/99I+ud3EWZ
/XokUrcrZY+tCyM//4fqeXzfedueLmrCSaXqDE88vLSUvD6NmsUWnK/nd1Hu
vx7pVJZCjYajOHl/BkW7IusiqOmkcrJ/PagX6lP/iEzpRC+fSpNpWu4V2Ymv
L6CexctP2X6wE/3PVNR+W+kAR/L+IsprzfnhFNNO/Lxi6s2OMVVtBeKVUuvu
fXC73t6BXhl16+bJH/MyO9CroNYpDDeK3e1Ar5Ja1isyu+1kB3rVVK7dttF+
tw70aqiYuG22rPnk9TXUvj88mF5M4tdSuraxV2pH2/H9dVTbzNIX+6va8fPq
qVCLM+envW5Hr4FSf7VtiuS1dvQaKM7UmAuih9rRa6TSisvyWS7EE1A5W0x7
3WYTr4kK5VBjn9SJ10zd5CiLrPzehl4z5dtjRBspbkOvhdKUvCye+bwNvVbK
MSOE/iq8Db02alnSvLl5B9rQa6cs8oZ9JdcQr53SexJVstmCeB3UO42BtdXq
xOukehvXSh8bbEWvi+owMWyYX9ZKkT9333hqPYYvW/H13VTbrOszZl5rRb+b
eu5o+cD1CHl/D9Wytn/90w2t+Hk91EBHjx2PIl4vleu9xSuFSbxeao6J8+eD
osTrpfryI72dGlrQ66MMZVYuWZ3egl4f9e6k8vpDd1vQ66fuWdfcfhfYgl4/
xS8+qD9tVwt6/dTfXvsL/nEg3gA1vNAvwXMG8QaoqsLlSfPoxPtGtaykas0H
m9H7Rs2zrJ+2oqIZvW9UYGDoxZA3zeh9p3SZPbKdMc3ofaeYjkMXvE82ozdI
NRvZ6Wh4EG+Qer+x+U6tPfEGKXqSLTvXhHg/qOMh1sFVqsT7QUktu16i8KMJ
vSEquS9naGtVE3pDlMyGDW2V75rQG6I2H9S7cTC2Cb1hKm+di5R5UBN6w1Sd
9VdNuhfxRqgL9QNZ6k7EG6G+vXKVt5pLvBGK67Gj4C894v2kqNHsSV3ixBt/
/qA67nCrAL1RimfGuT4lX4DeKGXA9CsZThKgN0rNOj91cdc1AXpjVFvIoxrx
AAF6YxR9T1qw5Xbi0SDLVNYyfDl5PQ2KkvdWKpkTnwavF6uue6RF3k+Da1kb
7+wWJZ9HA/l5nnEOLY3oiYDYl65lq/Ib0RMB70bGXweSG9ETAa8aR/PXNxrR
E4E1R9fbs081oicCS9+kJj7YRTxRCL7k7eTkRDxRuGBUosm0Ip4ovPbc2iHL
JZ4oWCS+S9KQJ54ofL/ltN5moAE9MRhca5RxqaoBPTH4N/xjvVh6A3picD93
e1hYQgN6YpChp/F2wZUG9MQg+o/Ny1T9iScOKw4u40h4EE8cuDbPDdVWEU8c
aC9XLLD9jXjiYHdJfvVVDvHEwVy73mbSZOJJQEKi5VDE93r0JODZ3L2OS2rr
0ZOAf+bsmKOZVY+eBFwpuxc8KakePQkQLT1sqhFdj54kfPtdW9bmLPEkofOu
xMhFH+JJwvX8UcHoBuJJwteyvvizdsSThD0q9F//Tk94UnBmEd15WId4UnCc
mj1aK008KXib0SbZ1FeHnhRAGOt3ma916ElBdsCTtKWZdehJw+43YnMeJNWh
Jw3z6CfDDGLq0JMGRz+bxLRg4kkDU9rV94gv8aRhgfqJ96u2Em8ScP1WeNus
IN4kWBjww22NJfEmQcbP9j9PGhBvEuw6yL+Yr0q8SfCCfj5otgjxZGBG4cDC
lI5a9GRg1CH9+sbyWvTGnx9/5q39sRY9GbgVUx7zLbEWPRloT4iT7bhZi54s
KG6QDBkLIZ4sBFV2ixkdJp4sNNntWe67nXiyECAlsaTaiXiykDw5umjTAuLJ
gWf0jawRE+LJQVr4cbkkXeLJATehYN8ZOeLJwfa6htoDQzXoycHprjwT/+Ya
9OTBgy5B3S6uQU8e1pUFtNR9qEFPHqzSfeQtE2vQk4d3NcWnH8UQTx544onm
VueJNxkUTTT6G/yINxma3LqiY3cRbzJYHFdWC1hPvMkQ23TY7OBS4k2GlimS
X07OJZ4CvPtwsyTekHgKoMmS1+jSIJ4C3OilfBZLE08BbrONgv/s+oqeAgy4
1pdfyP2KniIsqNwWXBH/deL1DEW4Oj907eqQr+grwuI5xmMiu/D99ooworR8
ZrX9V/w8Rbhy/Eli6zT0/BVh3p0rFnx59CLGvehJh8I6qie8REW4Xmlubpxb
jT2K4GLiyR5KqJ7wBIrwe636WF9oNfYpgbJbXIC2dzX2KYGWeOt2n1XomSlB
A89w56AZevZKsHpanO0DOnruSjCVdvVp0Pcq7FOCvw9O8Ykoq8I+JVhfGTSn
6FUV9inBo/6YD7/dqMI+JdDO356f51+FfUrAmX1K47wbejRlkDAtdD5igx5D
GWas9Fx2yRA9M2UoWRP9tEQOPXtlKC5rWmjTXYl9ysCubUn9+rkS+5Rh/6mT
9beSK7FPGY6uzva4eK0S+5ShzlpC6cGxSuxThqnF8eEdbugJlKF61+EXLnbo
0VRAyWfjzK5p6DFUICHe91O8EnpmKnD0y0yzy98qsE8FLF4Ym92pqMA+FVgS
bXO6KrUC+1SgLiKs0yq2AvtUwFnSxyA9qAL7VGDbdf/+/d7oZauA6eZLysvW
oCdQARGfzPlLrdCjqcL6wKBFXmz0GKrwYWR3xatJ6JmpwuBt93+nd5djnyp0
zLz+4d/icuxThQqf05EnUsqxTxVsQwult90uxz5V+P3JxYI9QeXYpwqPz/94
HrUXvWxVaDrx6HC3C3oCVUit9q7bZo0ejQ6jkacSxAzRY9DB2G3WrTRF9Mzo
sOO+819xg2XYR4c5r8voj2rKsI8OdpoK2uWfyrCPDoFrA9cZPCnDPjrU+8//
9f+QiT46GEzaHjk9AL1sOnhIeHEaPNET0Mf/Ppr+fP07ejQ1OFr3vDXZEj2G
GjSIn0os4qJnpgYp0ja6KgqkTw2Ob/res3ewFPvUYJ7T2EBfbSn2qcGfTt96
LmeVYp8aqC/KjHdOKsU+NXBYUdnzW3Qp9qnB59CYU9RZ9ARq8N5cwtzNBz2a
OsQctkiP3YAeQx1evssSk7FDz0wdqv89kBw8Az17dai66f/ASIf0qUP0Ibjf
JUX61GHlo7/88npLsE8dWJyNfTlVJdinDk4cg6LWTyXYpw7tR160c5+VYJ86
KGT1SR+LLsE+DfjkMmNgIBA9hgZUb5vuF3gAPTMNiHPI8LLcgp79+M+fbI2U
XY6euwY8PHbo0/c56PlrwDWRn59oPNKnAca6bE+eMunTgP5In0PbR4uxTwOG
Z/7zKb21GPs04Er4LI51STH2MUD++Kr1JWnF2MeAT9IpcPZxMfYxICkyKWZ1
FHr2DHi2Y7cDFYieOwOWwBvNBQfR82eA9JqIItet6EUwYIq54/JrK9FLZIDA
4bBL1zzSx4AMK8VPG6aRPgbsm2G+U8AgfZowVyJU46wU6dMc/7TM6EUDX7BP
E7wU2F+06r9gnybMuSoVLFvwBfs0oVKJfVX17Rfs04S8wP05s+K/YJ8meNw6
RNt7Db1ETdh+NUIi/Sx62ZpgXNH2ZIYvegJNOOd5sTzZHT2aFoj2fXN0ckKP
oQWnGlW/yywkfVoQ6lscXjqD9GnBBdEYmTcs0qcFy7r0DV8pkj4t4EWuz8gb
K8I+Lah76vt+uLMI+7Tg6sO+tnnVRdinBUe3RetE5BRhnxZElJ8zlkkpwj5t
OBihPXThIXoMbWg73Opkch09M20QDy/XagxCz14bdtqZcB8fRs9dG0oY85Ze
2IGevzYcG6zfdmYd6dOGWzu0nS4vJX3aYMFkdyT/Rvq04bWmn0z3VNKnDQcM
xc5YaZM+HbAZfgM35UifDkgztBU1fn7GPh0YkIspuN3xGft0YPZblqtN9Wfs
04En0h2RI7mfsU8HZESN1318+xn7dOAfnwS/2MfoJepAQK1lXUQMetk6YLSQ
5hZzET2BDihb7qx+fQI9mi6Yl/JntO9Dj6ELs0YkLI3dSZ8u/B26qPzkGtKn
CwXS09u67EifLlz3C1u2ay7p04VDalNrfk4lfbpwr0QlIEaH9OmCVEKOsrMC
6dMF4+993iwR0qcL1nzLYz97C7FPD/Y2fNBsbyjEPj14/6GX01FciH160LnP
yo+WUYh9evBxqY8U71Uh9ulBaHVEzIZ49Pz1wFyhZHpcNHoRemDskREiehG9
RD3wkRj6e89f6GXrweYbOx06D6An0INVp+bv9/MgfUzYserjN731+Ho5JgSU
Xn5euJz0MsFozPdKuDW+n8cEm7KNR3bMJP1M8NmZYreSjx4wwWL041dbTfJ9
mLCPYcp1kEfPhQlzo9pH3Gnk+zFhi3v7ktC+ggnPhwkbv/v3fRQU4PdlQnfV
nArF8oIJL5gJD0o/CjxyCvD7M6HpZc9QYWrBhBfLhKjXjSNLn6GXOP758sUF
RXHopTLB1tt1065r6GUzYUBRPYgeil4ZE0w+Ds3NCUBPwISIYimnSwfQ62NC
XfOJt56e6NFYEL5hbPOKjejJscA+skNv0Sr0GCyIyDpQb2eDHo8F5t/cLq6b
i54ZC9p72lWOGqEHLChOue4Qr4+ePQuc5FSmt6uh58KCzPqqsDmy6LmzYNtO
2Z2Xx/JxPxZUz7wZNdaXj/uxwDhmlpJvcz7uxwKOfH/CaGU+7seCtMXFm8IK
8nE/Fviu2ik3Kx29RBbsW5wTJXiJXioLTOQMf95+hF42C2RsX6p430avjAWf
5m7OWhKBnoAFLap75GecQ6+PBd1d1Wm8E+jR9OES7X6pgS96cvrg+UFr1uxd
6DH0IW3L1nTHLejx9OHx0Wn7/NagZ6YPvISvhkn26IE+JBa+Khy0Rs9eH5aG
rHVZbIGeiz4oM6vvxk5Dz10fbMr7Hyrqk/30Yf+joBWn1cl++uDKfn9MRp7s
pw+jq/u4kaJkP32IsnbYuLM7D/fTB4mskskpFXm4nz78PJNvzvuYh/vpQ8mV
KQ/jnubhfuPeupgl1jfycD998Oi3a+09i55gvC98zx+v96PXpw8ZRmsvR7mi
R2NDU5CRY/gy9OTYsKPe80jsHPQYbHid9qY5h4sejw0j/BIfOWX0zNjgey7p
h+toLu7HhvkLpzvltubifmxYXrV8y8qSXNyPDZuPxY+0puXifmyIs14ocu1x
Lu7HBhtBrL1rFHr+bPjL/edLq0D0gtnQu7Db0OQgehFsSP4ZdmD2VvRi2eB5
2sFv9Ur0Etmw/WiNXuA89FLHfbkdsz9PRS+bDezm0lhTBnplbHgcZbH4jiR6
AjYoTnUem96fg/uxYfCFd1xWbQ7ux4GQBAWWX14O7scBsd0S9otScnA/Dpie
r5NiPcjB/ThAj06cqRSBnhkH9iaNPaWfRg84sGnGuy1G+9Gz58BvKpu469zQ
c+HAV7Wkimsr0XPnQC0j2b2bQs+HAy+WpMc5G6HnzwERo1D/Im30gjmgd3/a
ZzdZ9CI4cPX02uOSQ9m4Hwc0NuXtet2cjftxwPbnnUMnS7JxPw5IOMcHuqZn
434csFu68/iKpGzcjwM3PM2tV91GT8CBowb34t3D0OvjwGwNwf2QAPRoXOj/
9EProzd6clx4mMtuVHZFj8GFJxyrVq8V6PG4cCrgT8VaCj0zLsy82/+bmzF6
wIV5AQorv+uiZ8+FPFtn86jJ6Llwof6MXK7jaBbux4WQzWI03c4s3I8LodY/
Hg9WZeF+XHA4eye9IScL9+NC6p+z+XUpWbgfFyx7K150x6MXy4W7X7xcFW+g
l8iFTqc4Ketz6KVyweeu0cUAP/SyuVBML28q2o1eGRfCFyd2WWxCT8CFP/SK
zj9wQK+PC1ranU+NAT0aD85JTLf6YIKeHA80crere7DQY/DgR0i1oa4yejwe
TE9dsr5eFD0zHgi4xeEv+jJxPx54f5z5PqohE/fjgTX3ffrFL5m4Hw966A4h
V9IzcT8eyDxLG32QnIn78cD0o6l4/l30/HnQt7volFgkesE8UMhT8LIJQi+C
B60vve5dPYJeLA+s3B5qD+1CL5EHmwz2JXpuQi91/POfWDm3rkAvmwc57q39
vvPRKxvvd73iozYTPcF4fwL37Xsuen08eFzOfH9UHT0aH3QfcrbZTkJPjg9J
SfMj9IYzcD8+/HnHZ4FkRwbux4fQG3UrhqszcD8+FA2Jx4/mZ+B+fGg6unOB
YloG7scHgzsPx4yT0HMZ94Oe5GyIQ8+dDxam7Ij/7pgT+/HheZr58vog9Pz5
cP+CTtlcP/SC+eB9Vpd30wu9CD5YxrezVNzQi+XDi5biF2FO6CXyQcSjpYRl
h14qH2wzn+9ImYte9vjzwr2u7tPRK+NDzDQqVpuJnoAPUbvW8GuU0evjw/Gz
/MLHEujRDGDLLealc4OfyD3k/3sYwP/++8J78YQnvBcL78XCe7HwXiy8Fwvv
xf89hPdi4b1YeC8W3ouF9+KJPuG9WHgvFt6Lhfdi4b1YeC/+7yG8FwvvxcJ7
sfBeLLwXT+wnvBcL78XCe7HwXiy8F0/sJ7wX/1+5F/8PH5I6OQ==
"], {{{},
{GrayLevel[0.5], Opacity[0.2], EdgeForm[None],
GraphicsGroupBox[PolygonBox[CompressedData["
1:eJwl1eN2HQwUBNDc2rZT27bN1G5TK6lt27Zt27btz7bR7qz82GvOC8yc4NDw
kLAoQUFBAT4JROb7QOT9jre84TWveMkLnvOMpzzhMY94yAPuc4+73OE2t7jJ
Da5zjatc4TKXuMgFznOOs5zhNKc4yQmOc4yjHOEwhzjIAfazj73sYTe72MkO
trONrWxhM5vYyAbWs461rGE1q1jJCpazjKUsYTGLWMgC5jOPucxhNrOYyQym
M42pTGEyk5jIBMYzjrGMYTSjGMkIhjOMoQxhMIMYyAD604++hBNGH3rTi570
oDvd6EoXOtOJUDrSgfa0oy1taE0rWtKC5jSjKU1oTAiNaEgD6lOPutShNrWo
SQ2qU42qVKEylahIBcpTjrKUoTSlKEkJilOMohShMIUoSAHyk4+85CE3uchJ
DrKTjaxkITPBZCIjGUhPOtKShtSkIiUpSE4ykpKExCQiIQmITzziEofYxCIm
MYhONKIShUAgspgfxP/8x7/8w9/8xZ/8we/8xq/8ws/8xI/8wPd8x7d8w9d8
xZd8wed8xqcROxCxAbzjLW94zSte8oLnPOMpT3jMIx7ygPvc4y53uM0tbnKD
61zjKle4zCUucoHznOMsZzjNKU5yguMc4yhHOMwhDnKA/exjL3vYzS52soPt
bGMrW9jMJjaygfWsYy1rWM0qVrKC5SxjKUtYzCIWsoD5zGMuc5jNLGYyg+lM
YypTmMwkJjKB8YxjLGMYzShGMoLhDGMoQxjMIAYygP70oy/hhNGH3vSiJz3o
Tje60oXOdCKUjnSgPe1oSxta04qWtKA5zWhKExoTQiMa0oD61KMudahNLWpS
g+pUoypVqEwlKlKB8pSjLGUoTSlKUoLiFKMoRShMIQpSgPzkIy95yE0ucpKD
7GQjK1nITDCZyEgG0pOOtKQhNalISQqSk4ykJCExiUhIAuITj7jEITaxiEkM
ohONqET85I/KPOgU
"]]]}, {}, {}}, {{}, {},
TagBox[
{GrayLevel[0.5], PointSize[0.007333333333333334],
AbsoluteThickness[2], LineBox[CompressedData["
1:eJwl1VN7GAgUBNCktu3Utm0ztdvUTG3btm3btm2u6q5Z7J58+3Dm/oKZGxQS
Gtw9MCAg4KUIu+EITwQiEonIRCEq0YhODGISi9jEIS7xiE8CEpKIxCQhKclI
TgpSkorUpCGItKQjPRnISCYyk4WsZCM7OchJLnKTh7zkIz8FKEghClOEohSj
OCUoSSlKU4aylKM8FahIJSpThapUozo1qEktalOHYOpSj/o0oCGNaEwTmtKM
5rSgJa1oTQhtaEs72tOBjnSiM13oSje6E0oPetKL3vShL/3ozwAGMojBDGEo
wxjOCEYyitGMYSzjGM8EJjKJyUxhKtOYzgxmMovZzGEu85jPAhayiMUsYSnL
WM4KVrKK1axhLetYzwY2sonNbGEr29jODnayi93sYS/72M8BDnKIwxzhKMc4
zglOcorTnOEs5zjPBS5yictc4SrXuM4NbnKL29zhLve4zwMe8ojHPOEpz3jO
C17yDd/yHd/zA694zRve8o73fOAjP/ITP/MLv/Ibv/MHf/IXf/MPn/jMF77y
L2HlDyQc4YlARCIRmShEJRrRiUFMYhGbOMQlHvFJQEISkZgkJCUZyUlBSlKR
mjQEkZZ0pCcDGclEZrKQlWxkJwc5yUVu8pCXfOSnAAUpRGGKUJRiFKcEJSlF
acpQlnKUpwIVqURlqlCValSnBjWpRW3qEExd6lGfBjSkEY1pQlOa0ZwWtKQV
rQmhDW1pR3s60JFOdKYLXelG2HiH0oOe9KI3fehLP/ozgIEMYjBDGMowhjOC
kYxiNGMYyzjGM4GJTGIyU5jKNKYzg5nMYjZzmMs85rOAhSxiMUtYyjKWs4KV
rGI1a1jLOtazgY1sYjNb2Mo2trODnexiN3vYyz72c4CDHOIwRzjKMY5zgpOc
4jRnOMs5znOBi1ziMle4yjWuc4Ob3OI2d7jLPe7zgIc84jFPeMoznvMi8P8f
/B+TN+ZJ
"]]},
Annotation[#, "Charting`Private`Tag#1"]& ]}}],
{"WolframDynamicHighlight", <|
"Label" -> {"XYLabel"}, "Ball" -> {"IndicatedBall"}|>}],
DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {},
Slot["HighlightElements"],
Slot["LayoutOptions"],
Slot["Meta"],
Charting`HighlightActionFunction["DynamicHighlight",
GraphicsComplex[CompressedData["
1:eJzt2flbzHv8//FpT4vWqaZ1ptnK0iIJp7yfWSpEOEXZkxNlK1k6liTH1iJE
FEkHZStLkoOInGhfpL20TvtepFKfjp6vH77f6/MPfK5r5hfXaOY293lcuC7X
k7Vlz6o/RGk02gcxGu2/X//fRzcVGfHrQV0KOd0SV9GFz+9Rj/99I+ud3EWZ
/XokUrcrZY+tCyM//4fqeXzfedueLmrCSaXqDE88vLSUvD6NmsUWnK/nd1Hu
vx7pVJZCjYajOHl/BkW7IusiqOmkcrJ/PagX6lP/iEzpRC+fSpNpWu4V2Ymv
L6CexctP2X6wE/3PVNR+W+kAR/L+IsprzfnhFNNO/Lxi6s2OMVVtBeKVUuvu
fXC73t6BXhl16+bJH/MyO9CroNYpDDeK3e1Ar5Ja1isyu+1kB3rVVK7dttF+
tw70aqiYuG22rPnk9TXUvj88mF5M4tdSuraxV2pH2/H9dVTbzNIX+6va8fPq
qVCLM+envW5Hr4FSf7VtiuS1dvQaKM7UmAuih9rRa6TSisvyWS7EE1A5W0x7
3WYTr4kK5VBjn9SJ10zd5CiLrPzehl4z5dtjRBspbkOvhdKUvCye+bwNvVbK
MSOE/iq8Db02alnSvLl5B9rQa6cs8oZ9JdcQr53SexJVstmCeB3UO42BtdXq
xOukehvXSh8bbEWvi+owMWyYX9ZKkT9333hqPYYvW/H13VTbrOszZl5rRb+b
eu5o+cD1CHl/D9Wytn/90w2t+Hk91EBHjx2PIl4vleu9xSuFSbxeao6J8+eD
osTrpfryI72dGlrQ66MMZVYuWZ3egl4f9e6k8vpDd1vQ66fuWdfcfhfYgl4/
xS8+qD9tVwt6/dTfXvsL/nEg3gA1vNAvwXMG8QaoqsLlSfPoxPtGtaykas0H
m9H7Rs2zrJ+2oqIZvW9UYGDoxZA3zeh9p3SZPbKdMc3ofaeYjkMXvE82ozdI
NRvZ6Wh4EG+Qer+x+U6tPfEGKXqSLTvXhHg/qOMh1sFVqsT7QUktu16i8KMJ
vSEquS9naGtVE3pDlMyGDW2V75rQG6I2H9S7cTC2Cb1hKm+di5R5UBN6w1Sd
9VdNuhfxRqgL9QNZ6k7EG6G+vXKVt5pLvBGK67Gj4C894v2kqNHsSV3ixBt/
/qA67nCrAL1RimfGuT4lX4DeKGXA9CsZThKgN0rNOj91cdc1AXpjVFvIoxrx
AAF6YxR9T1qw5Xbi0SDLVNYyfDl5PQ2KkvdWKpkTnwavF6uue6RF3k+Da1kb
7+wWJZ9HA/l5nnEOLY3oiYDYl65lq/Ib0RMB70bGXweSG9ETAa8aR/PXNxrR
E4E1R9fbs081oicCS9+kJj7YRTxRCL7k7eTkRDxRuGBUosm0Ip4ovPbc2iHL
JZ4oWCS+S9KQJ54ofL/ltN5moAE9MRhca5RxqaoBPTH4N/xjvVh6A3picD93
e1hYQgN6YpChp/F2wZUG9MQg+o/Ny1T9iScOKw4u40h4EE8cuDbPDdVWEU8c
aC9XLLD9jXjiYHdJfvVVDvHEwVy73mbSZOJJQEKi5VDE93r0JODZ3L2OS2rr
0ZOAf+bsmKOZVY+eBFwpuxc8KakePQkQLT1sqhFdj54kfPtdW9bmLPEkofOu
xMhFH+JJwvX8UcHoBuJJwteyvvizdsSThD0q9F//Tk94UnBmEd15WId4UnCc
mj1aK008KXib0SbZ1FeHnhRAGOt3ma916ElBdsCTtKWZdehJw+43YnMeJNWh
Jw3z6CfDDGLq0JMGRz+bxLRg4kkDU9rV94gv8aRhgfqJ96u2Em8ScP1WeNus
IN4kWBjww22NJfEmQcbP9j9PGhBvEuw6yL+Yr0q8SfCCfj5otgjxZGBG4cDC
lI5a9GRg1CH9+sbyWvTGnx9/5q39sRY9GbgVUx7zLbEWPRloT4iT7bhZi54s
KG6QDBkLIZ4sBFV2ixkdJp4sNNntWe67nXiyECAlsaTaiXiykDw5umjTAuLJ
gWf0jawRE+LJQVr4cbkkXeLJATehYN8ZOeLJwfa6htoDQzXoycHprjwT/+Ya
9OTBgy5B3S6uQU8e1pUFtNR9qEFPHqzSfeQtE2vQk4d3NcWnH8UQTx544onm
VueJNxkUTTT6G/yINxma3LqiY3cRbzJYHFdWC1hPvMkQ23TY7OBS4k2GlimS
X07OJZ4CvPtwsyTekHgKoMmS1+jSIJ4C3OilfBZLE08BbrONgv/s+oqeAgy4
1pdfyP2KniIsqNwWXBH/deL1DEW4Oj907eqQr+grwuI5xmMiu/D99ooworR8
ZrX9V/w8Rbhy/Eli6zT0/BVh3p0rFnx59CLGvehJh8I6qie8REW4Xmlubpxb
jT2K4GLiyR5KqJ7wBIrwe636WF9oNfYpgbJbXIC2dzX2KYGWeOt2n1XomSlB
A89w56AZevZKsHpanO0DOnruSjCVdvVp0Pcq7FOCvw9O8Ykoq8I+JVhfGTSn
6FUV9inBo/6YD7/dqMI+JdDO356f51+FfUrAmX1K47wbejRlkDAtdD5igx5D
GWas9Fx2yRA9M2UoWRP9tEQOPXtlKC5rWmjTXYl9ysCubUn9+rkS+5Rh/6mT
9beSK7FPGY6uzva4eK0S+5ShzlpC6cGxSuxThqnF8eEdbugJlKF61+EXLnbo
0VRAyWfjzK5p6DFUICHe91O8EnpmKnD0y0yzy98qsE8FLF4Ym92pqMA+FVgS
bXO6KrUC+1SgLiKs0yq2AvtUwFnSxyA9qAL7VGDbdf/+/d7oZauA6eZLysvW
oCdQARGfzPlLrdCjqcL6wKBFXmz0GKrwYWR3xatJ6JmpwuBt93+nd5djnyp0
zLz+4d/icuxThQqf05EnUsqxTxVsQwult90uxz5V+P3JxYI9QeXYpwqPz/94
HrUXvWxVaDrx6HC3C3oCVUit9q7bZo0ejQ6jkacSxAzRY9DB2G3WrTRF9Mzo
sOO+819xg2XYR4c5r8voj2rKsI8OdpoK2uWfyrCPDoFrA9cZPCnDPjrU+8//
9f+QiT46GEzaHjk9AL1sOnhIeHEaPNET0Mf/Ppr+fP07ejQ1OFr3vDXZEj2G
GjSIn0os4qJnpgYp0ja6KgqkTw2Ob/res3ewFPvUYJ7T2EBfbSn2qcGfTt96
LmeVYp8aqC/KjHdOKsU+NXBYUdnzW3Qp9qnB59CYU9RZ9ARq8N5cwtzNBz2a
OsQctkiP3YAeQx1evssSk7FDz0wdqv89kBw8Az17dai66f/ASIf0qUP0Ibjf
JUX61GHlo7/88npLsE8dWJyNfTlVJdinDk4cg6LWTyXYpw7tR160c5+VYJ86
KGT1SR+LLsE+DfjkMmNgIBA9hgZUb5vuF3gAPTMNiHPI8LLcgp79+M+fbI2U
XY6euwY8PHbo0/c56PlrwDWRn59oPNKnAca6bE+eMunTgP5In0PbR4uxTwOG
Z/7zKb21GPs04Er4LI51STH2MUD++Kr1JWnF2MeAT9IpcPZxMfYxICkyKWZ1
FHr2DHi2Y7cDFYieOwOWwBvNBQfR82eA9JqIItet6EUwYIq54/JrK9FLZIDA
4bBL1zzSx4AMK8VPG6aRPgbsm2G+U8AgfZowVyJU46wU6dMc/7TM6EUDX7BP
E7wU2F+06r9gnybMuSoVLFvwBfs0oVKJfVX17Rfs04S8wP05s+K/YJ8meNw6
RNt7Db1ETdh+NUIi/Sx62ZpgXNH2ZIYvegJNOOd5sTzZHT2aFoj2fXN0ckKP
oQWnGlW/yywkfVoQ6lscXjqD9GnBBdEYmTcs0qcFy7r0DV8pkj4t4EWuz8gb
K8I+Lah76vt+uLMI+7Tg6sO+tnnVRdinBUe3RetE5BRhnxZElJ8zlkkpwj5t
OBihPXThIXoMbWg73Opkch09M20QDy/XagxCz14bdtqZcB8fRs9dG0oY85Ze
2IGevzYcG6zfdmYd6dOGWzu0nS4vJX3aYMFkdyT/Rvq04bWmn0z3VNKnDQcM
xc5YaZM+HbAZfgM35UifDkgztBU1fn7GPh0YkIspuN3xGft0YPZblqtN9Wfs
04En0h2RI7mfsU8HZESN1318+xn7dOAfnwS/2MfoJepAQK1lXUQMetk6YLSQ
5hZzET2BDihb7qx+fQI9mi6Yl/JntO9Dj6ELs0YkLI3dSZ8u/B26qPzkGtKn
CwXS09u67EifLlz3C1u2ay7p04VDalNrfk4lfbpwr0QlIEaH9OmCVEKOsrMC
6dMF4+993iwR0qcL1nzLYz97C7FPD/Y2fNBsbyjEPj14/6GX01FciH160LnP
yo+WUYh9evBxqY8U71Uh9ulBaHVEzIZ49Pz1wFyhZHpcNHoRemDskREiehG9
RD3wkRj6e89f6GXrweYbOx06D6An0INVp+bv9/MgfUzYserjN731+Ho5JgSU
Xn5euJz0MsFozPdKuDW+n8cEm7KNR3bMJP1M8NmZYreSjx4wwWL041dbTfJ9
mLCPYcp1kEfPhQlzo9pH3Gnk+zFhi3v7ktC+ggnPhwkbv/v3fRQU4PdlQnfV
nArF8oIJL5gJD0o/CjxyCvD7M6HpZc9QYWrBhBfLhKjXjSNLn6GXOP758sUF
RXHopTLB1tt1065r6GUzYUBRPYgeil4ZE0w+Ds3NCUBPwISIYimnSwfQ62NC
XfOJt56e6NFYEL5hbPOKjejJscA+skNv0Sr0GCyIyDpQb2eDHo8F5t/cLq6b
i54ZC9p72lWOGqEHLChOue4Qr4+ePQuc5FSmt6uh58KCzPqqsDmy6LmzYNtO
2Z2Xx/JxPxZUz7wZNdaXj/uxwDhmlpJvcz7uxwKOfH/CaGU+7seCtMXFm8IK
8nE/Fviu2ik3Kx29RBbsW5wTJXiJXioLTOQMf95+hF42C2RsX6p430avjAWf
5m7OWhKBnoAFLap75GecQ6+PBd1d1Wm8E+jR9OES7X6pgS96cvrg+UFr1uxd
6DH0IW3L1nTHLejx9OHx0Wn7/NagZ6YPvISvhkn26IE+JBa+Khy0Rs9eH5aG
rHVZbIGeiz4oM6vvxk5Dz10fbMr7Hyrqk/30Yf+joBWn1cl++uDKfn9MRp7s
pw+jq/u4kaJkP32IsnbYuLM7D/fTB4mskskpFXm4nz78PJNvzvuYh/vpQ8mV
KQ/jnubhfuPeupgl1jfycD998Oi3a+09i55gvC98zx+v96PXpw8ZRmsvR7mi
R2NDU5CRY/gy9OTYsKPe80jsHPQYbHid9qY5h4sejw0j/BIfOWX0zNjgey7p
h+toLu7HhvkLpzvltubifmxYXrV8y8qSXNyPDZuPxY+0puXifmyIs14ocu1x
Lu7HBhtBrL1rFHr+bPjL/edLq0D0gtnQu7Db0OQgehFsSP4ZdmD2VvRi2eB5
2sFv9Ur0Etmw/WiNXuA89FLHfbkdsz9PRS+bDezm0lhTBnplbHgcZbH4jiR6
AjYoTnUem96fg/uxYfCFd1xWbQ7ux4GQBAWWX14O7scBsd0S9otScnA/Dpie
r5NiPcjB/ThAj06cqRSBnhkH9iaNPaWfRg84sGnGuy1G+9Gz58BvKpu469zQ
c+HAV7Wkimsr0XPnQC0j2b2bQs+HAy+WpMc5G6HnzwERo1D/Im30gjmgd3/a
ZzdZ9CI4cPX02uOSQ9m4Hwc0NuXtet2cjftxwPbnnUMnS7JxPw5IOMcHuqZn
434csFu68/iKpGzcjwM3PM2tV91GT8CBowb34t3D0OvjwGwNwf2QAPRoXOj/
9EProzd6clx4mMtuVHZFj8GFJxyrVq8V6PG4cCrgT8VaCj0zLsy82/+bmzF6
wIV5AQorv+uiZ8+FPFtn86jJ6Llwof6MXK7jaBbux4WQzWI03c4s3I8LodY/
Hg9WZeF+XHA4eye9IScL9+NC6p+z+XUpWbgfFyx7K150x6MXy4W7X7xcFW+g
l8iFTqc4Ketz6KVyweeu0cUAP/SyuVBML28q2o1eGRfCFyd2WWxCT8CFP/SK
zj9wQK+PC1ranU+NAT0aD85JTLf6YIKeHA80crere7DQY/DgR0i1oa4yejwe
TE9dsr5eFD0zHgi4xeEv+jJxPx54f5z5PqohE/fjgTX3ffrFL5m4Hw966A4h
V9IzcT8eyDxLG32QnIn78cD0o6l4/l30/HnQt7volFgkesE8UMhT8LIJQi+C
B60vve5dPYJeLA+s3B5qD+1CL5EHmwz2JXpuQi91/POfWDm3rkAvmwc57q39
vvPRKxvvd73iozYTPcF4fwL37Xsuen08eFzOfH9UHT0aH3QfcrbZTkJPjg9J
SfMj9IYzcD8+/HnHZ4FkRwbux4fQG3UrhqszcD8+FA2Jx4/mZ+B+fGg6unOB
YloG7scHgzsPx4yT0HMZ94Oe5GyIQ8+dDxam7Ij/7pgT+/HheZr58vog9Pz5
cP+CTtlcP/SC+eB9Vpd30wu9CD5YxrezVNzQi+XDi5biF2FO6CXyQcSjpYRl
h14qH2wzn+9ImYte9vjzwr2u7tPRK+NDzDQqVpuJnoAPUbvW8GuU0evjw/Gz
/MLHEujRDGDLLealc4OfyD3k/3sYwP/++8J78YQnvBcL78XCe7HwXiy8Fwvv
xf89hPdi4b1YeC8W3ouF9+KJPuG9WHgvFt6Lhfdi4b1YeC/+7yG8FwvvxcJ7
sfBeLLwXT+wnvBcL78XCe7HwXiy8F0/sJ7wX/1+5F/8PH5I6OQ==
"], {{{}, {
EdgeForm[],