-
Notifications
You must be signed in to change notification settings - Fork 0
/
VM15D.nb
1652 lines (1620 loc) · 77.7 KB
/
VM15D.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 11.2' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 79416, 1644]
NotebookOptionsPosition[ 77388, 1601]
NotebookOutlinePosition[ 77781, 1617]
CellTagsIndexPosition[ 77738, 1614]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["\<\
Hamiltonian fluid reduction of the 1.5D Vlasov-Maxwell equations\
\>", "Title",
CellChangeTimes->{{3.753546115869246*^9, 3.7535463066624527`*^9}, {
3.7535463938212347`*^9, 3.75354641630978*^9}, {3.753546468332477*^9,
3.753546516482757*^9}, {3.753546568015329*^9, 3.753546591777028*^9}, {
3.75354662412742*^9, 3.753546630879004*^9}, 3.7535467420299587`*^9, {
3.753546819042474*^9, 3.7535469607739477`*^9}, {3.753547004470717*^9,
3.7535473086410027`*^9}, {3.753547339057282*^9, 3.753547431758519*^9}, {
3.753547679649687*^9, 3.753547841298791*^9}, {3.753547873892469*^9,
3.753547997760531*^9}, {3.753548064804661*^9, 3.7535480965215063`*^9}, {
3.753548222471622*^9, 3.75354822469851*^9}, 3.753548255573922*^9, {
3.753548297843718*^9, 3.753548327473954*^9}, {3.75354840657583*^9,
3.7535486788413763`*^9}, {3.7535508958785954`*^9, 3.75355103797435*^9}, {
3.7535510680542307`*^9, 3.753551199675138*^9}, {3.753552123616346*^9,
3.75355213157292*^9}, {3.7535522117292957`*^9, 3.753552263591567*^9}, {
3.7535524666206207`*^9, 3.75355246972407*^9}, {3.75355358085892*^9,
3.753553585199996*^9}, {3.753553706330567*^9, 3.75355371716817*^9}, {
3.7535537739296637`*^9, 3.7535537767674723`*^9}, {3.753554133422422*^9,
3.753554136646284*^9}, {3.753554371302622*^9, 3.753554373821686*^9}, {
3.753554746331951*^9, 3.753554766218598*^9}, {3.753559924350161*^9,
3.753559926412731*^9}, {3.753560119208755*^9, 3.753560121022291*^9}, {
3.7535601606797237`*^9, 3.7535602822091703`*^9}, {3.7535620858282833`*^9,
3.753562087899042*^9}, {3.753562170990518*^9, 3.753562172083952*^9}, {
3.753736134375054*^9, 3.753736155861287*^9}, {3.753736193644183*^9,
3.753736195394608*^9}, {3.753737247821797*^9, 3.75373725097248*^9}, {
3.753737336035337*^9, 3.753737337993511*^9}, {3.753737377308753*^9,
3.753737380210651*^9}, {3.753737762336728*^9, 3.753737765414063*^9},
3.753800351347273*^9, {3.753800406903344*^9, 3.75380041809588*^9}, {
3.753800522264441*^9, 3.753800723785987*^9}, {3.753800792890429*^9,
3.753800821791657*^9}, {3.753800892025881*^9, 3.75380090191398*^9},
3.753805444868209*^9, {3.753805709254031*^9, 3.753805739333407*^9},
3.75380582195673*^9, {3.760958929067404*^9, 3.760959007138435*^9}, {
3.76095904856314*^9, 3.760959132234791*^9}, {3.760959171155839*^9,
3.760959181266307*^9}, {3.760959235713327*^9, 3.7609592459042377`*^9}, {
3.760959353279406*^9, 3.760959427864913*^9}, {3.760959476739066*^9,
3.760959478648314*^9}, {3.76095961790515*^9, 3.7609596200067987`*^9}, {
3.760960897949101*^9, 3.760960904482552*^9}, {3.760961459704638*^9,
3.7609615735241957`*^9}, {3.761021997278446*^9, 3.761022042876995*^9}, {
3.761024008372117*^9, 3.7610240228501053`*^9}, {3.761024086882766*^9,
3.7610241870963373`*^9}, {3.761024227817196*^9, 3.761024281303578*^9}, {
3.761024318393301*^9, 3.761024320247088*^9}, {3.7610243522867517`*^9,
3.761024354662118*^9}, {3.761024462376025*^9, 3.761024464565569*^9}, {
3.761026840595697*^9, 3.761026843488748*^9}, {3.761030403686119*^9,
3.761030457435652*^9}, {3.7613045448736258`*^9, 3.7613045553290997`*^9}, {
3.7670211648020906`*^9, 3.76702117672124*^9}, {3.771647938741407*^9,
3.771647938919251*^9}, {3.773732905956273*^9, 3.7737329066989183`*^9}, {
3.773733036708457*^9, 3.773733106032875*^9}, {3.773733190511716*^9,
3.773733319054948*^9}, {3.773733352909101*^9, 3.7737333557086563`*^9}, {
3.773733565211183*^9, 3.773733622308074*^9}, {3.7737336523942947`*^9,
3.773733920744495*^9}, {3.773733968559547*^9, 3.773734164409007*^9},
3.773734705027666*^9, {3.773736156430729*^9, 3.773736163646825*^9}, {
3.7737362798767548`*^9, 3.773736280289803*^9}, {3.773736314894017*^9,
3.7737363198613577`*^9}, {3.77373639125883*^9, 3.7737363969178247`*^9}, {
3.773736482196463*^9, 3.773736489200789*^9}, {3.773737514911749*^9,
3.773737520249723*^9}, {3.7738110292098703`*^9, 3.77381121179672*^9}, {
3.7738123478864202`*^9, 3.773812360278975*^9}, {3.7739752469919024`*^9,
3.773975267384288*^9}, {3.773976219831649*^9, 3.773976238171723*^9}, {
3.773978245258353*^9, 3.7739782456642113`*^9}, {3.7739800121428623`*^9,
3.773980016744731*^9}, {3.773983885141548*^9, 3.773983908873831*^9}, {
3.7739862090961943`*^9, 3.773986213855732*^9}, {3.773986536605913*^9,
3.773986542163089*^9}, {3.7740652177172813`*^9, 3.774065457264765*^9}, {
3.774072776815504*^9, 3.77407281609195*^9}, {3.77414838068399*^9,
3.774148389264089*^9}, 3.774510053828025*^9, {3.774693162498599*^9,
3.774693175696795*^9}, {3.7747635317512827`*^9, 3.774763771492153*^9}, {
3.774763834665825*^9, 3.7747638854819117`*^9}, {3.7747639525927687`*^9,
3.774763958011119*^9}, {3.774764058573048*^9, 3.774764063877684*^9}, {
3.774765514332431*^9, 3.7747655261794453`*^9}, {3.774765581557225*^9,
3.7747655988267527`*^9}, 3.774771508525646*^9, {3.77676401428086*^9,
3.776764143911002*^9}, {3.776764180839233*^9, 3.776764234553648*^9},
3.776764264783798*^9, {3.776764318517859*^9, 3.77676433689266*^9}, {
3.776764418171701*^9, 3.7767645982133427`*^9}, {3.776764643851592*^9,
3.7767647617530203`*^9}, {3.776765596483079*^9, 3.776765599913987*^9}, {
3.7767656696352386`*^9, 3.77676567821799*^9}, {3.776765730153578*^9,
3.776765752896085*^9}, {3.776765873425028*^9, 3.776765888575511*^9},
3.776766625043023*^9, {3.7767672074698763`*^9, 3.77676724256505*^9}, {
3.777026066937914*^9, 3.7770262182070293`*^9}, {3.8285989552487097`*^9,
3.828598961888832*^9}, {3.828598992127674*^9, 3.8285990151913652`*^9}, {
3.828599050387712*^9, 3.828599168233653*^9}, {3.8285992937054157`*^9,
3.828599310159523*^9}},ExpressionUUID->"df979f07-c71f-4966-b27b-\
19072cd565d9"],
Cell[TextData[{
"based on [C. Chandre, B.A. Shadwick, ",
StyleBox["Hamiltonian fluid reduction of the 1.5D Vlasov-Maxwell equations,",
FontSlant->"Italic"],
" ",
ButtonBox["Physics of Plasmas",
BaseStyle->"Hyperlink",
ButtonData->{
URL["https://doi.org/10.1063/5.0056155"], None},
ButtonNote->"https://doi.org/10.1063/5.0056155"],
" 28, 092114 (2021)]\n"
}], "Text",
CellChangeTimes->{{3.753546115869246*^9, 3.7535463066624527`*^9}, {
3.7535463938212347`*^9, 3.75354641630978*^9}, {3.753546468332477*^9,
3.753546516482757*^9}, {3.753546568015329*^9, 3.753546591777028*^9}, {
3.75354662412742*^9, 3.753546630879004*^9}, 3.7535467420299587`*^9, {
3.753546819042474*^9, 3.7535469607739477`*^9}, {3.753547004470717*^9,
3.7535473086410027`*^9}, {3.753547339057282*^9, 3.753547431758519*^9}, {
3.753547679649687*^9, 3.753547841298791*^9}, {3.753547873892469*^9,
3.753547997760531*^9}, {3.753548064804661*^9, 3.7535480965215063`*^9}, {
3.753548222471622*^9, 3.75354822469851*^9}, 3.753548255573922*^9, {
3.753548297843718*^9, 3.753548327473954*^9}, {3.75354840657583*^9,
3.7535486788413763`*^9}, {3.7535508958785954`*^9, 3.75355103797435*^9}, {
3.7535510680542307`*^9, 3.753551199675138*^9}, {3.753552123616346*^9,
3.75355213157292*^9}, {3.7535522117292957`*^9, 3.753552263591567*^9}, {
3.7535524666206207`*^9, 3.75355246972407*^9}, {3.75355358085892*^9,
3.753553585199996*^9}, {3.753553706330567*^9, 3.75355371716817*^9}, {
3.7535537739296637`*^9, 3.7535537767674723`*^9}, {3.753554133422422*^9,
3.753554136646284*^9}, {3.753554371302622*^9, 3.753554373821686*^9}, {
3.753554746331951*^9, 3.753554766218598*^9}, {3.753559924350161*^9,
3.753559926412731*^9}, {3.753560119208755*^9, 3.753560121022291*^9}, {
3.7535601606797237`*^9, 3.7535602822091703`*^9}, {3.7535620858282833`*^9,
3.753562087899042*^9}, {3.753562170990518*^9, 3.753562172083952*^9}, {
3.753736134375054*^9, 3.753736155861287*^9}, {3.753736193644183*^9,
3.753736195394608*^9}, {3.753737247821797*^9, 3.75373725097248*^9}, {
3.753737336035337*^9, 3.753737337993511*^9}, {3.753737377308753*^9,
3.753737380210651*^9}, {3.753737762336728*^9, 3.753737765414063*^9},
3.753800351347273*^9, {3.753800406903344*^9, 3.75380041809588*^9}, {
3.753800522264441*^9, 3.753800723785987*^9}, {3.753800792890429*^9,
3.753800821791657*^9}, {3.753800892025881*^9, 3.75380090191398*^9},
3.753805444868209*^9, {3.753805709254031*^9, 3.753805739333407*^9},
3.75380582195673*^9, {3.760958929067404*^9, 3.760959007138435*^9}, {
3.76095904856314*^9, 3.760959132234791*^9}, {3.760959171155839*^9,
3.760959181266307*^9}, {3.760959235713327*^9, 3.7609592459042377`*^9}, {
3.760959353279406*^9, 3.760959427864913*^9}, {3.760959476739066*^9,
3.760959478648314*^9}, {3.76095961790515*^9, 3.7609596200067987`*^9}, {
3.760960897949101*^9, 3.760960904482552*^9}, {3.760961459704638*^9,
3.7609615735241957`*^9}, {3.761021997278446*^9, 3.761022042876995*^9}, {
3.761024008372117*^9, 3.7610240228501053`*^9}, {3.761024086882766*^9,
3.7610241870963373`*^9}, {3.761024227817196*^9, 3.761024281303578*^9}, {
3.761024318393301*^9, 3.761024320247088*^9}, {3.7610243522867517`*^9,
3.761024354662118*^9}, {3.761024462376025*^9, 3.761024464565569*^9}, {
3.761026840595697*^9, 3.761026843488748*^9}, {3.761030403686119*^9,
3.761030457435652*^9}, {3.7613045448736258`*^9, 3.7613045553290997`*^9}, {
3.7670211648020906`*^9, 3.76702117672124*^9}, {3.771647938741407*^9,
3.771647938919251*^9}, {3.773732905956273*^9, 3.7737329066989183`*^9}, {
3.773733036708457*^9, 3.773733106032875*^9}, {3.773733190511716*^9,
3.773733319054948*^9}, {3.773733352909101*^9, 3.7737333557086563`*^9}, {
3.773733565211183*^9, 3.773733622308074*^9}, {3.7737336523942947`*^9,
3.773733920744495*^9}, {3.773733968559547*^9, 3.773734164409007*^9},
3.773734705027666*^9, {3.773736156430729*^9, 3.773736163646825*^9}, {
3.7737362798767548`*^9, 3.773736280289803*^9}, {3.773736314894017*^9,
3.7737363198613577`*^9}, {3.77373639125883*^9, 3.7737363969178247`*^9}, {
3.773736482196463*^9, 3.773736489200789*^9}, {3.773737514911749*^9,
3.773737520249723*^9}, {3.7738110292098703`*^9, 3.77381121179672*^9}, {
3.7738123478864202`*^9, 3.773812360278975*^9}, {3.7739752469919024`*^9,
3.773975267384288*^9}, {3.773976219831649*^9, 3.773976238171723*^9}, {
3.773978245258353*^9, 3.7739782456642113`*^9}, {3.7739800121428623`*^9,
3.773980016744731*^9}, {3.773983885141548*^9, 3.773983908873831*^9}, {
3.7739862090961943`*^9, 3.773986213855732*^9}, {3.773986536605913*^9,
3.773986542163089*^9}, {3.7740652177172813`*^9, 3.774065457264765*^9}, {
3.774072776815504*^9, 3.77407281609195*^9}, {3.77414838068399*^9,
3.774148389264089*^9}, 3.774510053828025*^9, {3.774693162498599*^9,
3.774693175696795*^9}, {3.7747635317512827`*^9, 3.774763771492153*^9}, {
3.774763834665825*^9, 3.7747638854819117`*^9}, {3.7747639525927687`*^9,
3.774763958011119*^9}, {3.774764058573048*^9, 3.774764063877684*^9}, {
3.774765514332431*^9, 3.7747655261794453`*^9}, {3.774765581557225*^9,
3.7747655988267527`*^9}, 3.774771508525646*^9, {3.77676401428086*^9,
3.776764143911002*^9}, {3.776764180839233*^9, 3.776764234553648*^9},
3.776764264783798*^9, {3.776764318517859*^9, 3.77676433689266*^9}, {
3.776764418171701*^9, 3.7767645982133427`*^9}, {3.776764643851592*^9,
3.7767647617530203`*^9}, {3.776765596483079*^9, 3.776765599913987*^9}, {
3.7767656696352386`*^9, 3.77676567821799*^9}, {3.776765730153578*^9,
3.776765752896085*^9}, {3.776765873425028*^9, 3.776765888575511*^9},
3.776766625043023*^9, {3.7767672074698763`*^9, 3.77676724256505*^9}, {
3.777026066937914*^9, 3.7770262182070293`*^9}, {3.8285989552487097`*^9,
3.828598961888832*^9}, {3.828598992127674*^9, 3.8285990151913652`*^9}, {
3.828599050387712*^9, 3.828599168233653*^9}, {3.8285992937054157`*^9,
3.8285993541755877`*^9}, {3.829192500565709*^9, 3.829192547569335*^9}, {
3.829192588673379*^9, 3.829192609319477*^9}, {3.8413099414194393`*^9,
3.8413099670489597`*^9}, {3.841310007525659*^9,
3.841310007525931*^9}},ExpressionUUID->"45049cb9-be61-4b15-a57e-\
05ecd587aef1"],
Cell[CellGroupData[{
Cell["\<\
Definition of the fluid variables, the bracket and the closure functions\
\>", "Subtitle",
CellChangeTimes->{{3.753546115869246*^9, 3.7535463066624527`*^9}, {
3.7535463938212347`*^9, 3.75354641630978*^9}, {3.753546468332477*^9,
3.753546516482757*^9}, {3.753546568015329*^9, 3.753546591777028*^9}, {
3.75354662412742*^9, 3.753546630879004*^9}, 3.7535467420299587`*^9, {
3.753546819042474*^9, 3.7535469607739477`*^9}, {3.753547004470717*^9,
3.7535473086410027`*^9}, {3.753547339057282*^9, 3.753547431758519*^9}, {
3.753547679649687*^9, 3.753547841298791*^9}, {3.753547873892469*^9,
3.753547997760531*^9}, {3.753548064804661*^9, 3.7535480965215063`*^9}, {
3.753548222471622*^9, 3.75354822469851*^9}, 3.753548255573922*^9, {
3.753548297843718*^9, 3.753548327473954*^9}, {3.75354840657583*^9,
3.7535486788413763`*^9}, {3.7535508958785954`*^9, 3.75355103797435*^9}, {
3.7535510680542307`*^9, 3.753551199675138*^9}, {3.753552123616346*^9,
3.75355213157292*^9}, {3.7535522117292957`*^9, 3.753552263591567*^9}, {
3.7535524666206207`*^9, 3.75355246972407*^9}, {3.75355358085892*^9,
3.753553585199996*^9}, {3.753553706330567*^9, 3.75355371716817*^9}, {
3.7535537739296637`*^9, 3.7535537767674723`*^9}, {3.753554133422422*^9,
3.753554136646284*^9}, {3.753554371302622*^9, 3.753554373821686*^9}, {
3.753554746331951*^9, 3.753554766218598*^9}, {3.753559924350161*^9,
3.753559926412731*^9}, {3.753560119208755*^9, 3.753560121022291*^9}, {
3.7535601606797237`*^9, 3.7535602822091703`*^9}, {3.7535620858282833`*^9,
3.753562087899042*^9}, {3.753562170990518*^9, 3.753562172083952*^9}, {
3.753736134375054*^9, 3.753736155861287*^9}, {3.753736193644183*^9,
3.753736195394608*^9}, {3.753737247821797*^9, 3.75373725097248*^9}, {
3.753737336035337*^9, 3.753737337993511*^9}, {3.753737377308753*^9,
3.753737380210651*^9}, {3.753737762336728*^9, 3.753737765414063*^9},
3.753800351347273*^9, {3.753800406903344*^9, 3.75380041809588*^9}, {
3.753800522264441*^9, 3.753800723785987*^9}, {3.753800792890429*^9,
3.753800821791657*^9}, {3.753800892025881*^9, 3.75380090191398*^9},
3.753805444868209*^9, {3.753805709254031*^9, 3.753805739333407*^9},
3.75380582195673*^9, {3.760958929067404*^9, 3.760959007138435*^9}, {
3.76095904856314*^9, 3.760959132234791*^9}, {3.760959171155839*^9,
3.760959181266307*^9}, {3.760959235713327*^9, 3.7609592459042377`*^9}, {
3.760959353279406*^9, 3.760959427864913*^9}, {3.760959476739066*^9,
3.760959478648314*^9}, {3.76095961790515*^9, 3.7609596200067987`*^9}, {
3.760960897949101*^9, 3.760960904482552*^9}, {3.760961459704638*^9,
3.7609615735241957`*^9}, {3.761021997278446*^9, 3.761022042876995*^9}, {
3.761024008372117*^9, 3.7610240228501053`*^9}, {3.761024086882766*^9,
3.7610241870963373`*^9}, {3.761024227817196*^9, 3.761024281303578*^9}, {
3.761024318393301*^9, 3.761024320247088*^9}, {3.7610243522867517`*^9,
3.761024354662118*^9}, {3.761024462376025*^9, 3.761024464565569*^9}, {
3.761026840595697*^9, 3.761026843488748*^9}, {3.761030403686119*^9,
3.761030457435652*^9}, {3.7613045448736258`*^9, 3.7613045553290997`*^9}, {
3.7670211648020906`*^9, 3.76702117672124*^9}, {3.771647938741407*^9,
3.771647938919251*^9}, {3.773732905956273*^9, 3.7737329066989183`*^9}, {
3.773733036708457*^9, 3.773733106032875*^9}, {3.773733190511716*^9,
3.773733319054948*^9}, {3.773733352909101*^9, 3.7737333557086563`*^9}, {
3.773733565211183*^9, 3.773733622308074*^9}, {3.7737336523942947`*^9,
3.773733920744495*^9}, {3.773733968559547*^9, 3.773734164409007*^9},
3.773734705027666*^9, {3.773736156430729*^9, 3.773736163646825*^9}, {
3.7737362798767548`*^9, 3.773736280289803*^9}, {3.773736314894017*^9,
3.7737363198613577`*^9}, {3.77373639125883*^9, 3.7737363969178247`*^9}, {
3.773736482196463*^9, 3.773736489200789*^9}, {3.773737514911749*^9,
3.773737520249723*^9}, {3.7738110292098703`*^9, 3.77381121179672*^9}, {
3.7738123478864202`*^9, 3.773812360278975*^9}, {3.7739752469919024`*^9,
3.773975267384288*^9}, {3.773976219831649*^9, 3.773976238171723*^9}, {
3.773978245258353*^9, 3.7739782456642113`*^9}, {3.7739800121428623`*^9,
3.773980016744731*^9}, {3.773983885141548*^9, 3.773983908873831*^9}, {
3.7739862090961943`*^9, 3.773986213855732*^9}, {3.773986536605913*^9,
3.773986542163089*^9}, {3.7740652177172813`*^9, 3.774065457264765*^9}, {
3.774072776815504*^9, 3.77407281609195*^9}, {3.77414838068399*^9,
3.774148389264089*^9}, 3.774510053828025*^9, {3.774693162498599*^9,
3.774693175696795*^9}, {3.7747635317512827`*^9, 3.774763771492153*^9}, {
3.774763834665825*^9, 3.7747638854819117`*^9}, {3.7747639525927687`*^9,
3.774763958011119*^9}, {3.774764058573048*^9, 3.774764063877684*^9}, {
3.774765514332431*^9, 3.7747655261794453`*^9}, {3.774765581557225*^9,
3.7747655988267527`*^9}, 3.774771508525646*^9, {3.77676401428086*^9,
3.776764143911002*^9}, {3.776764180839233*^9, 3.776764234553648*^9},
3.776764264783798*^9, {3.776764318517859*^9, 3.77676433689266*^9}, {
3.776764418171701*^9, 3.7767645982133427`*^9}, {3.776764643851592*^9,
3.7767647617530203`*^9}, {3.776765596483079*^9, 3.776765599913987*^9}, {
3.7767656696352386`*^9, 3.77676567821799*^9}, {3.776765730153578*^9,
3.776765752896085*^9}, {3.776765873425028*^9, 3.776765888575511*^9},
3.776766625043023*^9, {3.7767672074698763`*^9, 3.77676724256505*^9}, {
3.777026066937914*^9, 3.7770262182070293`*^9}, {3.8285989552487097`*^9,
3.828598961888832*^9}, {3.828598992127674*^9, 3.8285990151913652`*^9}, {
3.828599050387712*^9, 3.828599168233653*^9}, {3.8285992937054157`*^9,
3.8285993541755877`*^9}, {3.829192500565709*^9, 3.829192547569335*^9}, {
3.8291925920395937`*^9,
3.829192593223858*^9}},ExpressionUUID->"976dbdcd-3506-4f9a-b7ed-\
7278a2731541"],
Cell[BoxData[{
RowBox[{
RowBox[{"Clear", "[", "\"\<Global\[OpenCurlyQuote]*\>\"", "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{"<<", "VariationalMethods`"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"SetOptions", "[",
RowBox[{"Simplify", ",",
RowBox[{"TimeConstraint", "\[Rule]", "5000"}]}], "]"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"Var", "=",
RowBox[{"{",
RowBox[{
RowBox[{"\[Rho]", "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["P", "x"], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["P", "p"], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"0", ",", "2"}]], "[", "z", "]"}]}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"VarR", "=",
RowBox[{"{",
RowBox[{
RowBox[{
SubscriptBox["P", "x"], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}], ",",
RowBox[{
SubscriptBox["S",
RowBox[{"0", ",", "2"}]], "[", "z", "]"}]}], "}"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"3", ",", "0"}]], "_"], "[",
RowBox[{"p_", ",", "s_", ",", "t_", ",", "y_"}], "]"}], ":=",
RowBox[{
RowBox[{"s",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}], "-",
RowBox[{
RowBox[{"s", "^", "2"}], "/",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"2", ",", "1"}]], "_"], "[",
RowBox[{"p_", ",", "s_", ",", "t_", ",", "y_"}], "]"}], ":=",
RowBox[{
RowBox[{"t",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}], "-",
RowBox[{"s", " ",
RowBox[{"t", "/",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"1", ",", "2"}]], "_"], "[",
RowBox[{"p_", ",", "s_", ",", "t_", ",", "y_"}], "]"}], ":=",
RowBox[{
RowBox[{"y",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}], "-",
RowBox[{
RowBox[{"t", "^", "2"}], "/",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"0", ",", "3"}]], "_"], "[",
RowBox[{"p_", ",", "s_", ",", "t_", ",", "y_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"t", "/", "s"}],
RowBox[{"(",
RowBox[{
RowBox[{"3", "y"}], "-",
RowBox[{"2",
RowBox[{
RowBox[{"t", "^", "2"}], "/", "s"}]}]}], ")"}],
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}], "-",
RowBox[{
RowBox[{"t", "^", "3"}], "/",
RowBox[{"(",
RowBox[{"s",
RowBox[{"(",
RowBox[{"\[Alpha]", "-", "p"}], ")"}]}], ")"}]}], "+",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"y", "-",
RowBox[{
RowBox[{"t", "^", "2"}], "/", "s"}]}], ")"}], "^",
RowBox[{"(",
RowBox[{"4", "/", "3"}], ")"}]}], "/",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"\[CapitalPsi]", "'"}], "[",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"y", "-",
RowBox[{
RowBox[{"t", "^", "2"}], "/", "s"}]}], ")"}], "^",
RowBox[{"(",
RowBox[{"1", "/", "3"}], ")"}]}],
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"p", "-", "\[Alpha]"}], ")"}], "^", "2"}], "+", "s"}],
")"}], "/", "s"}]}], "]"}], ")"}], "^", "2"}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"Bracket", "[",
RowBox[{"F_", ",", "G_"}], "]"}], ":=",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "1", "]"}], "]"}], ",", "z"}], "]"}],
RowBox[{"G", "[",
RowBox[{"[", "3", "]"}], "]"}]}], "-",
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "3", "]"}], "]"}],
RowBox[{"D", "[",
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "1", "]"}], "]"}], ",", "z"}], "]"}]}], "-",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["P", "x"], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "2", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"G", "[",
RowBox[{"[", "3", "]"}], "]"}]}], "-",
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "3", "]"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "2", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{"2",
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}]}], "+",
RowBox[{"4",
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}]}], "+",
RowBox[{"2",
RowBox[{
SubscriptBox["S",
RowBox[{"0", ",", "2"}]], "[", "z", "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}]}]}], ")"}]}], "+",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "2", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{"2",
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}],
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{"2",
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}],
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "2", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}], "+",
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "3", "]"}], "]"}],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"0", ",", "2"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"2", ",", "0"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"1", ",", "1"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "+",
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
SubscriptBox["S",
RowBox[{"0", ",", "2"}]], "[", "z", "]"}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}]}], ")"}],
RowBox[{"G", "[",
RowBox[{"[", "3", "]"}], "]"}]}], "+",
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"3", ",", "0"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}], "+",
RowBox[{"2",
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"2", ",", "1"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "4", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}], "+", "\[IndentingNewLine]",
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"1", ",", "2"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}], "+",
RowBox[{"2",
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"0", ",", "3"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}], "+",
RowBox[{
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"2", ",", "1"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}], "+",
RowBox[{"2",
RowBox[{
OverscriptBox[
SubscriptBox["S",
RowBox[{"1", ",", "2"}]], "_"], "@@", "VarR"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}],
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}]}], "-",
RowBox[{
RowBox[{
RowBox[{"F", "[",
RowBox[{"[", "6", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}],
RowBox[{"D", "[",
RowBox[{
RowBox[{
RowBox[{"G", "[",
RowBox[{"[", "5", "]"}], "]"}], "/",
RowBox[{"\[Rho]", "[", "z", "]"}]}], ",", "z"}], "]"}]}]}],
")"}]}]}]}], ";"}]}], "Input",
CellChangeTimes->{{3.753546115869246*^9, 3.7535463066624527`*^9}, {
3.7535463938212347`*^9, 3.75354641630978*^9}, {3.753546468332477*^9,
3.753546516482757*^9}, {3.753546568015329*^9, 3.753546591777028*^9}, {
3.75354662412742*^9, 3.753546630879004*^9}, 3.7535467420299587`*^9, {
3.753546819042474*^9, 3.7535469607739477`*^9}, {3.753547004470717*^9,
3.7535473086410027`*^9}, {3.753547339057282*^9, 3.753547431758519*^9}, {
3.753547679649687*^9, 3.753547841298791*^9}, {3.753547873892469*^9,
3.753547997760531*^9}, {3.753548064804661*^9, 3.7535480965215063`*^9}, {
3.753548222471622*^9, 3.75354822469851*^9}, 3.753548255573922*^9, {
3.753548297843718*^9, 3.753548327473954*^9}, {3.75354840657583*^9,
3.7535486788413763`*^9}, {3.7535508958785954`*^9, 3.75355103797435*^9}, {
3.7535510680542307`*^9, 3.753551199675138*^9}, {3.753552123616346*^9,
3.75355213157292*^9}, {3.7535522117292957`*^9, 3.753552263591567*^9}, {
3.7535524666206207`*^9, 3.75355246972407*^9}, {3.75355358085892*^9,
3.753553585199996*^9}, {3.753553706330567*^9, 3.75355371716817*^9}, {
3.7535537739296637`*^9, 3.7535537767674723`*^9}, {3.753554133422422*^9,
3.753554136646284*^9}, {3.753554371302622*^9, 3.753554373821686*^9}, {
3.753554746331951*^9, 3.753554766218598*^9}, {3.753559924350161*^9,
3.753559926412731*^9}, {3.753560119208755*^9, 3.753560121022291*^9}, {
3.7535601606797237`*^9, 3.7535602822091703`*^9}, {3.7535620858282833`*^9,
3.753562087899042*^9}, {3.753562170990518*^9, 3.753562172083952*^9}, {
3.753736134375054*^9, 3.753736155861287*^9}, {3.753736193644183*^9,
3.753736195394608*^9}, {3.753737247821797*^9, 3.75373725097248*^9}, {
3.753737336035337*^9, 3.753737337993511*^9}, {3.753737377308753*^9,
3.753737380210651*^9}, {3.753737762336728*^9, 3.753737765414063*^9},
3.753800351347273*^9, {3.753800406903344*^9, 3.75380041809588*^9}, {
3.753800522264441*^9, 3.753800723785987*^9}, {3.753800792890429*^9,
3.753800821791657*^9}, {3.753800892025881*^9, 3.75380090191398*^9},
3.753805444868209*^9, {3.753805709254031*^9, 3.753805739333407*^9},
3.75380582195673*^9, {3.760958929067404*^9, 3.760959007138435*^9}, {
3.76095904856314*^9, 3.760959132234791*^9}, {3.760959171155839*^9,
3.760959181266307*^9}, {3.760959235713327*^9, 3.7609592459042377`*^9}, {
3.760959353279406*^9, 3.760959427864913*^9}, {3.760959476739066*^9,
3.760959478648314*^9}, {3.76095961790515*^9, 3.7609596200067987`*^9}, {
3.760960897949101*^9, 3.760960904482552*^9}, {3.760961459704638*^9,
3.7609615735241957`*^9}, {3.761021997278446*^9, 3.761022042876995*^9}, {
3.761024008372117*^9, 3.7610240228501053`*^9}, {3.761024086882766*^9,
3.7610241870963373`*^9}, {3.761024227817196*^9, 3.761024281303578*^9}, {
3.761024318393301*^9, 3.761024320247088*^9}, {3.7610243522867517`*^9,
3.761024354662118*^9}, {3.761024462376025*^9, 3.761024464565569*^9}, {
3.761026840595697*^9, 3.761026843488748*^9}, {3.761030403686119*^9,
3.761030457435652*^9}, {3.7613045448736258`*^9, 3.7613045553290997`*^9}, {
3.7670211648020906`*^9, 3.76702117672124*^9}, {3.771647938741407*^9,
3.771647938919251*^9}, {3.773732905956273*^9, 3.7737329066989183`*^9}, {
3.773733036708457*^9, 3.773733106032875*^9}, {3.773733190511716*^9,
3.773733319054948*^9}, {3.773733352909101*^9, 3.7737333557086563`*^9}, {
3.773733565211183*^9, 3.773733622308074*^9}, {3.7737336523942947`*^9,
3.773733920744495*^9}, {3.773733968559547*^9, 3.773734164409007*^9},
3.773734705027666*^9, {3.773736156430729*^9, 3.773736163646825*^9}, {
3.7737362798767548`*^9, 3.773736280289803*^9}, {3.773736314894017*^9,
3.7737363198613577`*^9}, {3.77373639125883*^9, 3.7737363969178247`*^9}, {
3.773736482196463*^9, 3.773736489200789*^9}, {3.773737514911749*^9,
3.773737520249723*^9}, {3.7738110292098703`*^9, 3.77381121179672*^9}, {
3.7738123478864202`*^9, 3.773812360278975*^9}, {3.7739752469919024`*^9,
3.773975267384288*^9}, {3.773976219831649*^9, 3.773976238171723*^9}, {
3.773978245258353*^9, 3.7739782456642113`*^9}, {3.7739800121428623`*^9,
3.773980016744731*^9}, {3.773983885141548*^9, 3.773983908873831*^9}, {
3.7739862090961943`*^9, 3.773986213855732*^9}, {3.773986536605913*^9,
3.773986542163089*^9}, {3.7740652177172813`*^9, 3.774065457264765*^9}, {
3.774072776815504*^9, 3.77407281609195*^9}, {3.77414838068399*^9,
3.774148389264089*^9}, 3.774510053828025*^9, {3.774693162498599*^9,
3.774693175696795*^9}, {3.7747635317512827`*^9, 3.774763771492153*^9}, {
3.774763834665825*^9, 3.7747638854819117`*^9}, {3.7747639525927687`*^9,
3.774763958011119*^9}, {3.774764058573048*^9, 3.774764063877684*^9}, {
3.774765514332431*^9, 3.7747655261794453`*^9}, {3.774765581557225*^9,
3.7747655988267527`*^9}, 3.774771508525646*^9, {3.77676401428086*^9,
3.776764143911002*^9}, {3.776764180839233*^9, 3.776764234553648*^9},
3.776764264783798*^9, {3.776764318517859*^9, 3.77676433689266*^9}, {
3.776764418171701*^9, 3.7767645982133427`*^9}, {3.776764643851592*^9,
3.7767647617530203`*^9}, {3.776765596483079*^9, 3.776765599913987*^9}, {
3.7767656696352386`*^9, 3.77676567821799*^9}, {3.776765730153578*^9,
3.776765752896085*^9}, {3.776765873425028*^9, 3.776765888575511*^9},
3.776766625043023*^9, {3.7767672074698763`*^9, 3.77676724256505*^9}, {
3.777026066937914*^9, 3.7770262182070293`*^9}, {3.8285989552487097`*^9,
3.828598961888832*^9}, {3.828598992127674*^9, 3.8285990151913652`*^9}, {
3.828599050387712*^9, 3.828599063573513*^9}, {3.828599190283894*^9,
3.828599224031007*^9}, 3.828599258876816*^9, {3.8286001074891376`*^9,
3.828600113856056*^9}, {3.828600144729143*^9, 3.828600407662435*^9}, {
3.8286011336986637`*^9, 3.828601194284034*^9}, {3.8286022358143578`*^9,
3.828602253002366*^9}, {3.828602716499908*^9, 3.8286027460968857`*^9}, {
3.828603076153716*^9, 3.8286030879613867`*^9}, {3.8286040771041803`*^9,
3.828604252882765*^9}, 3.828604654597836*^9, {3.8286048303539457`*^9,
3.828604831407155*^9}, {3.828606473669828*^9, 3.82860648767556*^9}, {
3.8286065179970827`*^9, 3.828606776888225*^9}, {3.828608269979114*^9,
3.828608279794272*^9}, {3.828608364033588*^9, 3.828608397625924*^9}, {
3.828695614073001*^9, 3.828695634739829*^9}, {3.828695770966875*^9,
3.8286957894462013`*^9}, {3.828695952480435*^9, 3.828696024600217*^9}, {
3.828696057119451*^9, 3.82869611507728*^9}, {3.828696653278467*^9,
3.828696655869323*^9}, {3.8290429587098103`*^9, 3.829042978188632*^9}, {
3.829185687545992*^9, 3.8291856883129787`*^9}, 3.82918572506422*^9, {
3.829185826398035*^9, 3.829185845468186*^9}, {3.829191429555284*^9,
3.8291914488461847`*^9}, {3.829191497736092*^9, 3.8291915004482107`*^9}, {
3.829197688041815*^9, 3.829197688367358*^9}, {3.829199446142827*^9,
3.829199447534585*^9}},
CellLabel->"In[51]:=",ExpressionUUID->"4a861d29-6513-41f3-bf2b-f7a3db9eacbd"]
}, Open ]],
Cell[CellGroupData[{
Cell["\<\
Jacobi identity (the calculation takes about an hour)\
\>", "Subtitle",
CellChangeTimes->{{3.753546115869246*^9, 3.7535463066624527`*^9}, {
3.7535463938212347`*^9, 3.75354641630978*^9}, {3.753546468332477*^9,
3.753546516482757*^9}, {3.753546568015329*^9, 3.753546591777028*^9}, {
3.75354662412742*^9, 3.753546630879004*^9}, 3.7535467420299587`*^9, {
3.753546819042474*^9, 3.7535469607739477`*^9}, {3.753547004470717*^9,
3.7535473086410027`*^9}, {3.753547339057282*^9, 3.753547431758519*^9}, {
3.753547679649687*^9, 3.753547841298791*^9}, {3.753547873892469*^9,
3.753547997760531*^9}, {3.753548064804661*^9, 3.7535480965215063`*^9}, {
3.753548222471622*^9, 3.75354822469851*^9}, 3.753548255573922*^9, {
3.753548297843718*^9, 3.753548327473954*^9}, {3.75354840657583*^9,
3.7535486788413763`*^9}, {3.7535508958785954`*^9, 3.75355103797435*^9}, {
3.7535510680542307`*^9, 3.753551199675138*^9}, {3.753552123616346*^9,
3.75355213157292*^9}, {3.7535522117292957`*^9, 3.753552263591567*^9}, {
3.7535524666206207`*^9, 3.75355246972407*^9}, {3.75355358085892*^9,
3.753553585199996*^9}, {3.753553706330567*^9, 3.75355371716817*^9}, {
3.7535537739296637`*^9, 3.7535537767674723`*^9}, {3.753554133422422*^9,
3.753554136646284*^9}, {3.753554371302622*^9, 3.753554373821686*^9}, {
3.753554746331951*^9, 3.753554766218598*^9}, {3.753559924350161*^9,
3.753559926412731*^9}, {3.753560119208755*^9, 3.753560121022291*^9}, {
3.7535601606797237`*^9, 3.7535602822091703`*^9}, {3.7535620858282833`*^9,
3.753562087899042*^9}, {3.753562170990518*^9, 3.753562172083952*^9}, {
3.753736134375054*^9, 3.753736155861287*^9}, {3.753736193644183*^9,
3.753736195394608*^9}, {3.753737247821797*^9, 3.75373725097248*^9}, {
3.753737336035337*^9, 3.753737337993511*^9}, {3.753737377308753*^9,
3.753737380210651*^9}, {3.753737762336728*^9, 3.753737765414063*^9},
3.753800351347273*^9, {3.753800406903344*^9, 3.75380041809588*^9}, {
3.753800522264441*^9, 3.753800723785987*^9}, {3.753800792890429*^9,
3.753800821791657*^9}, {3.753800892025881*^9, 3.75380090191398*^9},
3.753805444868209*^9, {3.753805709254031*^9, 3.753805739333407*^9},
3.75380582195673*^9, {3.760958929067404*^9, 3.760959007138435*^9}, {
3.76095904856314*^9, 3.760959132234791*^9}, {3.760959171155839*^9,
3.760959181266307*^9}, {3.760959235713327*^9, 3.7609592459042377`*^9}, {
3.760959353279406*^9, 3.760959427864913*^9}, {3.760959476739066*^9,
3.760959478648314*^9}, {3.76095961790515*^9, 3.7609596200067987`*^9}, {
3.760960897949101*^9, 3.760960904482552*^9}, {3.760961459704638*^9,
3.7609615735241957`*^9}, {3.761021997278446*^9, 3.761022042876995*^9}, {
3.761024008372117*^9, 3.7610240228501053`*^9}, {3.761024086882766*^9,
3.7610241870963373`*^9}, {3.761024227817196*^9, 3.761024281303578*^9}, {
3.761024318393301*^9, 3.761024320247088*^9}, {3.7610243522867517`*^9,
3.761024354662118*^9}, {3.761024462376025*^9, 3.761024464565569*^9}, {
3.761026840595697*^9, 3.761026843488748*^9}, {3.761030403686119*^9,
3.761030457435652*^9}, {3.7613045448736258`*^9, 3.7613045553290997`*^9}, {
3.7670211648020906`*^9, 3.76702117672124*^9}, {3.771647938741407*^9,
3.771647938919251*^9}, {3.773732905956273*^9, 3.7737329066989183`*^9}, {
3.773733036708457*^9, 3.773733106032875*^9}, {3.773733190511716*^9,
3.773733319054948*^9}, {3.773733352909101*^9, 3.7737333557086563`*^9}, {
3.773733565211183*^9, 3.773733622308074*^9}, {3.7737336523942947`*^9,
3.773733920744495*^9}, {3.773733968559547*^9, 3.773734164409007*^9},
3.773734705027666*^9, {3.773736156430729*^9, 3.773736163646825*^9}, {
3.7737362798767548`*^9, 3.773736280289803*^9}, {3.773736314894017*^9,
3.7737363198613577`*^9}, {3.77373639125883*^9, 3.7737363969178247`*^9}, {
3.773736482196463*^9, 3.773736489200789*^9}, {3.773737514911749*^9,
3.773737520249723*^9}, {3.7738110292098703`*^9, 3.77381121179672*^9}, {
3.7738123478864202`*^9, 3.773812360278975*^9}, {3.7739752469919024`*^9,