-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathBoltzmann.nb
More file actions
5534 lines (5431 loc) · 217 KB
/
Copy pathBoltzmann.nb
File metadata and controls
5534 lines (5431 loc) · 217 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 11.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 216049, 5526]
NotebookOptionsPosition[ 208147, 5385]
NotebookOutlinePosition[ 208517, 5401]
CellTagsIndexPosition[ 208474, 5398]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["Boltzmann Notebook", "Title",
CellChangeTimes->{{3.780741896999395*^9,
3.7807419186545143`*^9}},ExpressionUUID->"59410aec-699f-45e3-8f8f-\
884cce0d2e8a"],
Cell["\<\
Integrates even polynomials of speed inside Boltzmann collision integrals\
\>", "Subtitle",
CellChangeTimes->{{3.780741963853879*^9, 3.7807420295299797`*^9}, {
3.8498422719166317`*^9,
3.8498422837167673`*^9}},ExpressionUUID->"3937187c-3fda-4dc7-9eb1-\
b6b9a3e3e6d3"],
Cell[CellGroupData[{
Cell["INPUT", "Section",
CellChangeTimes->{{3.780742049264736*^9, 3.780742052960974*^9}, {
3.780744008851202*^9,
3.780744009731817*^9}},ExpressionUUID->"786f6e70-94eb-4cd6-b59e-\
26d985318e8b"],
Cell[TextData[StyleBox["Integral to be evaluated:",
FontColor->RGBColor[0.5019607843137255, 0., 0.]]], "Text",
CellChangeTimes->{{3.780742064415801*^9,
3.780742074028367*^9}},ExpressionUUID->"55514076-15bf-45fd-be23-\
ff51b05628c1"],
Cell[BoxData[
RowBox[{
RowBox[{"Integral", "=", "dres"}], ";"}]], "Input",
CellChangeTimes->{{3.7807420812646456`*^9, 3.780742094215147*^9},
3.780742339474197*^9, 3.7809093011050177`*^9, 3.78090949882063*^9,
3.780909750172062*^9, 3.7809098074784164`*^9, 3.780910464611392*^9,
3.7809105222985*^9, 3.7854843796968613`*^9, 3.7855728094126153`*^9,
3.78557286642575*^9, 3.7855729681540303`*^9, 3.7855732375175424`*^9,
3.78557333554687*^9, {3.793244414472411*^9, 3.793244415161874*^9},
3.793246702775085*^9, 3.7932470247568884`*^9, {3.7932561957165847`*^9,
3.7932562011302204`*^9}, 3.793256244261176*^9, 3.801627327863653*^9, {
3.801627404638936*^9, 3.8016274184463954`*^9}, {3.801627454168588*^9,
3.8016274876565723`*^9}, {3.8016278276175027`*^9,
3.8016278393643255`*^9}, {3.801628560598565*^9, 3.8016285634399953`*^9}, {
3.801628861607648*^9, 3.8016288646725454`*^9}, {3.8016288978900433`*^9,
3.801628938792036*^9}, {3.801632142491333*^9, 3.8016321497317686`*^9}, {
3.8016322902718525`*^9, 3.801632291045686*^9}, {3.8016330342869625`*^9,
3.8016330583762493`*^9}, {3.8016380765480995`*^9,
3.8016380784274683`*^9}, {3.848716604433584*^9, 3.848716605150857*^9}, {
3.848721753383397*^9, 3.848721773702955*^9}, {3.848721825179243*^9,
3.8487218304456806`*^9}, {3.848722150087372*^9, 3.848722168180901*^9}, {
3.8487222617746043`*^9, 3.8487222642567205`*^9}, {3.8487223572872543`*^9,
3.8487223649826655`*^9}, {3.84872244579508*^9, 3.8487224474344416`*^9}, {
3.848722482822032*^9, 3.84872248315419*^9}, {3.8487232891556053`*^9,
3.848723289542025*^9}, {3.8488991557351737`*^9, 3.8488991558513136`*^9}, {
3.849069068504829*^9, 3.8490690689276056`*^9}, 3.849069773482706*^9,
3.8490701696613913`*^9, {3.8490705514684134`*^9, 3.8490705517854137`*^9}, {
3.8498418422087917`*^9, 3.8498418437182207`*^9}},
CellLabel->
"In[153]:=",ExpressionUUID->"a7ac6e7f-cd5d-4c7e-85a1-977a75aad117"],
Cell[TextData[{
StyleBox["Put",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" ",
StyleBox["True",
FontWeight->"Bold"],
" ",
StyleBox["to show only diagonal terms",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" (i = j)"
}], "Text",
CellChangeTimes->{{3.780742134102508*^9, 3.780742153021083*^9}, {
3.785568333310904*^9, 3.7855683688107967`*^9}, {3.8016284610596313`*^9,
3.801628474348863*^9}, {3.801636812401534*^9,
3.8016368131353197`*^9}},ExpressionUUID->"bafb30d6-f64e-4dca-8874-\
62fb386d9f3f"],
Cell[BoxData[
RowBox[{
RowBox[{"diagonal", " ", "=", " ", "False"}], ";"}]], "Input",
CellLabel->
"In[154]:=",ExpressionUUID->"47eaf164-ed00-48ab-ae3e-8c4fd9f0dfb4"],
Cell[TextData[{
StyleBox["Put",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" ",
StyleBox["True",
FontWeight->"Bold"],
StyleBox[" to show equal temperature case",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" (i",
Cell[BoxData[
SubscriptBox["\[Theta]", "i"]],
CellChangeTimes->{{3.7770013311846304`*^9, 3.7770014031530986`*^9}, {
3.7770014479532614`*^9, 3.7770015758575025`*^9}, {3.7770016421087933`*^9,
3.7770017446245146`*^9}, {3.777001805143985*^9, 3.7770018234143457`*^9}, {
3.7770018576373453`*^9, 3.777001868701682*^9}, {3.777001960960461*^9,
3.777002046364682*^9}, {3.777002085201646*^9, 3.777002234055195*^9}, {
3.777002276610862*^9, 3.777002277305982*^9}, {3.7770023369077373`*^9,
3.77700234647639*^9}, {3.777004730111148*^9, 3.7770047367489705`*^9}, {
3.777004775713829*^9, 3.777004794393256*^9}, {3.7770048782432184`*^9,
3.7770048867193995`*^9}, {3.777005135356024*^9, 3.777005138503594*^9}, {
3.7770062715015907`*^9, 3.7770062725989447`*^9}, {3.777006388216528*^9,
3.777006388494771*^9}, {3.7770065328895864`*^9, 3.7770065333972597`*^9}, {
3.77700781314746*^9, 3.7770079815532436`*^9}, 3.777008100965679*^9, {
3.777008249647401*^9, 3.7770083143522425`*^9}, {3.777008496651662*^9,
3.7770085508750143`*^9}, {3.777008733992568*^9, 3.7770087984972696`*^9}, {
3.777008944898369*^9, 3.7770089526989555`*^9}, {3.777008994553643*^9,
3.7770090420542374`*^9}, {3.7770100865381565`*^9,
3.7770101090553946`*^9}, {3.7770102036221323`*^9, 3.777010224623185*^9}, {
3.7770103190307875`*^9, 3.7770103367719975`*^9}, {3.7770103799261312`*^9,
3.7770103895797234`*^9}, {3.7770104279021797`*^9,
3.7770104668006506`*^9}, {3.777010654236351*^9, 3.777010701877531*^9}, {
3.7770107342543592`*^9, 3.777010744769043*^9}, {3.777010776904726*^9,
3.7770107779093328`*^9}, {3.7770108212382236`*^9, 3.777010822369748*^9}, {
3.7770108546333265`*^9, 3.7770108653657227`*^9}, {3.7770116249534693`*^9,
3.7770116287053666`*^9}, {3.777011780344282*^9, 3.777011784597994*^9}, {
3.777011840045618*^9, 3.777011850205907*^9}, {3.777011886060293*^9,
3.777011899008877*^9}, {3.7770119986424265`*^9, 3.7770119988123016`*^9}, {
3.777012461300623*^9, 3.7770125885732007`*^9}, 3.7770126235147543`*^9, {
3.7770127558962736`*^9, 3.7770128074602027`*^9}, {3.777013115598998*^9,
3.777013129585382*^9}, {3.77701338321282*^9, 3.777013426523162*^9}, {
3.777013501281876*^9, 3.7770135014224634`*^9}, {3.7770135420890017`*^9,
3.777013551428228*^9}, {3.7770135883321295`*^9, 3.777013599872896*^9}, {
3.777013898564535*^9, 3.7770139001422243`*^9}, {3.7770148684195857`*^9,
3.777014899341736*^9}, {3.7770149301506357`*^9, 3.7770149321735177`*^9}, {
3.777014974464608*^9, 3.777014989319892*^9}, {3.777015078648814*^9,
3.777015105074156*^9}, {3.777015138841237*^9, 3.7770152518462467`*^9}, {
3.777015428301572*^9, 3.7770155043138466`*^9}, {3.7770157862401342`*^9,
3.777015856140811*^9}, {3.7770158895212708`*^9, 3.777015966786131*^9}, {
3.777016078005509*^9, 3.777016081899946*^9}, {3.7770174394982657`*^9,
3.7770174801726165`*^9}, {3.7770175363251467`*^9,
3.7770175889580936`*^9}, {3.7770176851513367`*^9,
3.7770177307735434`*^9}, {3.7770177809950657`*^9, 3.777017790951309*^9}, {
3.7770178493678446`*^9, 3.7770178796944556`*^9}, {3.777018150556406*^9,
3.77701815689841*^9}, {3.7770184764393396`*^9, 3.777018497824078*^9}, {
3.777018783772696*^9, 3.7770189402198563`*^9}, {3.777019077844077*^9,
3.777019093371052*^9}, {3.7770191626527576`*^9, 3.777019206255988*^9},
3.7770193795090685`*^9, 3.7770199771218805`*^9, 3.77702004620735*^9,
3.7770203335764894`*^9, {3.777020373131014*^9, 3.7770204167059307`*^9}, {
3.777020452412899*^9, 3.777020458414177*^9}, 3.777020580158021*^9,
3.7770208604724264`*^9, {3.7770214239901314`*^9, 3.7770215587137823`*^9},
3.77702235543779*^9, {3.7770224889440002`*^9, 3.7770225509347467`*^9}, {
3.777022594285115*^9, 3.777022616208768*^9}, {3.7770226959211955`*^9,
3.777022724936646*^9}, 3.7770227641543627`*^9, {3.777022826108985*^9,
3.7770230346753197`*^9}, {3.7770230729548206`*^9,
3.7770231418733463`*^9}, {3.7770231793220196`*^9, 3.777023317276226*^9}, {
3.7770238133730392`*^9, 3.7770238175017853`*^9}, {3.7770238502622967`*^9,
3.777023852922696*^9}, {3.7770248992730246`*^9, 3.777024954407259*^9},
3.7770249998332295`*^9, {3.777025123671203*^9, 3.777025184243395*^9}, {
3.7770254800369577`*^9, 3.777025496144828*^9}, {3.7770255781193275`*^9,
3.77702559318647*^9}, 3.7770256844670024`*^9, 3.777025790958192*^9,
3.777025823845621*^9, {3.7770260285500045`*^9, 3.777026047321204*^9},
3.7770260901824503`*^9, 3.7770263240340824`*^9, {3.7770288035891037`*^9,
3.777028806417425*^9}, {3.777028836542922*^9, 3.777028838958272*^9}, {
3.7770306662364397`*^9, 3.777030671737853*^9}, {3.7770310688904696`*^9,
3.7770310744367933`*^9}, {3.7770334816802864`*^9, 3.777033494148533*^9}, {
3.7770335742358685`*^9, 3.7770335836170654`*^9}, 3.7770336203881903`*^9, {
3.7770336652801285`*^9, 3.7770336719325943`*^9}, {3.7770337440514107`*^9,
3.7770337622553825`*^9}, {3.777034374130868*^9, 3.777034385032163*^9}, {
3.777101223083253*^9, 3.7771012581692953`*^9}, 3.777185891378599*^9,
3.7771859723458595`*^9, 3.7771860561144*^9, 3.7771861201130857`*^9, {
3.777186305096449*^9, 3.7771863105111217`*^9}, {3.777186985656376*^9,
3.7771869890912657`*^9}, {3.7771871579795165`*^9,
3.7771871618791203`*^9}, {3.7771876152725973`*^9, 3.7771876271682916`*^9},
3.7771892445617137`*^9, {3.7771899235496063`*^9,
3.7771900314103465`*^9}, {3.7771909323189116`*^9,
3.7771909510833845`*^9}, {3.7771938372914205`*^9, 3.777193915792196*^9}, {
3.777193954969129*^9, 3.777193974695012*^9}, 3.777194178867585*^9, {
3.7771945208125277`*^9, 3.7771945343602047`*^9}, {3.777194620203471*^9,
3.7771946586358624`*^9}, {3.777194823489591*^9, 3.7771948537420063`*^9}, {
3.7771949178206434`*^9, 3.7771949429272194`*^9}, {3.777194994396601*^9,
3.77719500347902*^9}, {3.7771950441500373`*^9, 3.7771950478577595`*^9}, {
3.7771950984896555`*^9, 3.7771950996749277`*^9}, 3.7771951436958847`*^9, {
3.7771952412769413`*^9, 3.777195266360532*^9}, 3.777195673281556*^9, {
3.777196023254566*^9, 3.7771960257550306`*^9}, 3.7771962108413315`*^9,
3.7771964348207383`*^9, {3.777196571470793*^9, 3.777196603700994*^9}, {
3.77719671640376*^9, 3.777196744655404*^9}, {3.777197041016348*^9,
3.7771970857498674`*^9}, {3.7771979803363595`*^9,
3.7771979809319925`*^9}, {3.7771985013183146`*^9,
3.7771985038037896`*^9}, {3.77719854132374*^9, 3.7771985416865177`*^9},
3.778066242768161*^9, 3.7784834449847374`*^9, {3.7784837002916636`*^9,
3.77848370204801*^9}, {3.780736502102478*^9, 3.7807365167849236`*^9},
3.780736938371971*^9, {3.780737007460211*^9, 3.7807370188105755`*^9}, {
3.780737373989303*^9, 3.78073738205152*^9}, {3.78073744060362*^9,
3.7807374451267796`*^9}, {3.780737768503789*^9, 3.780737788173134*^9},
3.7807378333061094`*^9, {3.7807408103030705`*^9, 3.780740813742545*^9}, {
3.7807417728658915`*^9, 3.7807417796989555`*^9}, {3.7807421629929233`*^9,
3.7807421689971156`*^9}, {3.7807423441213446`*^9,
3.7807423727735453`*^9}, {3.780742403029793*^9, 3.7807424804653244`*^9}, {
3.7807426625065765`*^9, 3.7807428612785363`*^9}, {3.780742911085021*^9,
3.7807429341225705`*^9}, 3.78090925493522*^9, {3.780909432760846*^9,
3.7809094520279646`*^9}, {3.7809097142728496`*^9, 3.780909744665778*^9}, {
3.7809099036772923`*^9, 3.780909904270878*^9}, {3.780910146375532*^9,
3.7809101467318783`*^9}, {3.785573324541329*^9, 3.7855733300553274`*^9}, {
3.7855736052321033`*^9, 3.7855736253168097`*^9}, 3.7856501558932457`*^9, {
3.801636489784712*^9, 3.801636610571928*^9}, {3.8016366543281145`*^9,
3.801636696693453*^9}},ExpressionUUID->
"695dbd1e-72d5-4d0d-88bb-d5c8f8f9e445"],
" = j",
Cell[BoxData[
SubscriptBox["\[Theta]", "j"]],
CellChangeTimes->{{3.7770013311846304`*^9, 3.7770014031530986`*^9}, {
3.7770014479532614`*^9, 3.7770015758575025`*^9}, {3.7770016421087933`*^9,
3.7770017446245146`*^9}, {3.777001805143985*^9, 3.7770018234143457`*^9}, {
3.7770018576373453`*^9, 3.777001868701682*^9}, {3.777001960960461*^9,
3.777002046364682*^9}, {3.777002085201646*^9, 3.777002234055195*^9}, {
3.777002276610862*^9, 3.777002277305982*^9}, {3.7770023369077373`*^9,
3.77700234647639*^9}, {3.777004730111148*^9, 3.7770047367489705`*^9}, {
3.777004775713829*^9, 3.777004794393256*^9}, {3.7770048782432184`*^9,
3.7770048867193995`*^9}, {3.777005135356024*^9, 3.777005138503594*^9}, {
3.7770062715015907`*^9, 3.7770062725989447`*^9}, {3.777006388216528*^9,
3.777006388494771*^9}, {3.7770065328895864`*^9, 3.7770065333972597`*^9}, {
3.77700781314746*^9, 3.7770079815532436`*^9}, 3.777008100965679*^9, {
3.777008249647401*^9, 3.7770083143522425`*^9}, {3.777008496651662*^9,
3.7770085508750143`*^9}, {3.777008733992568*^9, 3.7770087984972696`*^9}, {
3.777008944898369*^9, 3.7770089526989555`*^9}, {3.777008994553643*^9,
3.7770090420542374`*^9}, {3.7770100865381565`*^9,
3.7770101090553946`*^9}, {3.7770102036221323`*^9, 3.777010224623185*^9}, {
3.7770103190307875`*^9, 3.7770103367719975`*^9}, {3.7770103799261312`*^9,
3.7770103895797234`*^9}, {3.7770104279021797`*^9,
3.7770104668006506`*^9}, {3.777010654236351*^9, 3.777010701877531*^9}, {
3.7770107342543592`*^9, 3.777010744769043*^9}, {3.777010776904726*^9,
3.7770107779093328`*^9}, {3.7770108212382236`*^9, 3.777010822369748*^9}, {
3.7770108546333265`*^9, 3.7770108653657227`*^9}, {3.7770116249534693`*^9,
3.7770116287053666`*^9}, {3.777011780344282*^9, 3.777011784597994*^9}, {
3.777011840045618*^9, 3.777011850205907*^9}, {3.777011886060293*^9,
3.777011899008877*^9}, {3.7770119986424265`*^9, 3.7770119988123016`*^9}, {
3.777012461300623*^9, 3.7770125885732007`*^9}, 3.7770126235147543`*^9, {
3.7770127558962736`*^9, 3.7770128074602027`*^9}, {3.777013115598998*^9,
3.777013129585382*^9}, {3.77701338321282*^9, 3.777013426523162*^9}, {
3.777013501281876*^9, 3.7770135014224634`*^9}, {3.7770135420890017`*^9,
3.777013551428228*^9}, {3.7770135883321295`*^9, 3.777013599872896*^9}, {
3.777013898564535*^9, 3.7770139001422243`*^9}, {3.7770148684195857`*^9,
3.777014899341736*^9}, {3.7770149301506357`*^9, 3.7770149321735177`*^9}, {
3.777014974464608*^9, 3.777014989319892*^9}, {3.777015078648814*^9,
3.777015105074156*^9}, {3.777015138841237*^9, 3.7770152518462467`*^9}, {
3.777015428301572*^9, 3.7770155043138466`*^9}, {3.7770157862401342`*^9,
3.777015856140811*^9}, {3.7770158895212708`*^9, 3.777015966786131*^9}, {
3.777016078005509*^9, 3.777016081899946*^9}, {3.7770174394982657`*^9,
3.7770174801726165`*^9}, {3.7770175363251467`*^9,
3.7770175889580936`*^9}, {3.7770176851513367`*^9,
3.7770177307735434`*^9}, {3.7770177809950657`*^9, 3.777017790951309*^9}, {
3.7770178493678446`*^9, 3.7770178796944556`*^9}, {3.777018150556406*^9,
3.77701815689841*^9}, {3.7770184764393396`*^9, 3.777018497824078*^9}, {
3.777018783772696*^9, 3.7770189402198563`*^9}, {3.777019077844077*^9,
3.777019093371052*^9}, {3.7770191626527576`*^9, 3.777019206255988*^9},
3.7770193795090685`*^9, 3.7770199771218805`*^9, 3.77702004620735*^9,
3.7770203335764894`*^9, {3.777020373131014*^9, 3.7770204167059307`*^9}, {
3.777020452412899*^9, 3.777020458414177*^9}, 3.777020580158021*^9,
3.7770208604724264`*^9, {3.7770214239901314`*^9, 3.7770215587137823`*^9},
3.77702235543779*^9, {3.7770224889440002`*^9, 3.7770225509347467`*^9}, {
3.777022594285115*^9, 3.777022616208768*^9}, {3.7770226959211955`*^9,
3.777022724936646*^9}, 3.7770227641543627`*^9, {3.777022826108985*^9,
3.7770230346753197`*^9}, {3.7770230729548206`*^9,
3.7770231418733463`*^9}, {3.7770231793220196`*^9, 3.777023317276226*^9}, {
3.7770238133730392`*^9, 3.7770238175017853`*^9}, {3.7770238502622967`*^9,
3.777023852922696*^9}, {3.7770248992730246`*^9, 3.777024954407259*^9},
3.7770249998332295`*^9, {3.777025123671203*^9, 3.777025184243395*^9}, {
3.7770254800369577`*^9, 3.777025496144828*^9}, {3.7770255781193275`*^9,
3.77702559318647*^9}, 3.7770256844670024`*^9, 3.777025790958192*^9,
3.777025823845621*^9, {3.7770260285500045`*^9, 3.777026047321204*^9},
3.7770260901824503`*^9, 3.7770263240340824`*^9, {3.7770288035891037`*^9,
3.777028806417425*^9}, {3.777028836542922*^9, 3.777028838958272*^9}, {
3.7770306662364397`*^9, 3.777030671737853*^9}, {3.7770310688904696`*^9,
3.7770310744367933`*^9}, {3.7770334816802864`*^9, 3.777033494148533*^9}, {
3.7770335742358685`*^9, 3.7770335836170654`*^9}, 3.7770336203881903`*^9, {
3.7770336652801285`*^9, 3.7770336719325943`*^9}, {3.7770337440514107`*^9,
3.7770337622553825`*^9}, {3.777034374130868*^9, 3.777034385032163*^9}, {
3.777101223083253*^9, 3.7771012581692953`*^9}, 3.777185891378599*^9,
3.7771859723458595`*^9, 3.7771860561144*^9, 3.7771861201130857`*^9, {
3.777186305096449*^9, 3.7771863105111217`*^9}, {3.777186985656376*^9,
3.7771869890912657`*^9}, {3.7771871579795165`*^9,
3.7771871618791203`*^9}, {3.7771876152725973`*^9, 3.7771876271682916`*^9},
3.7771892445617137`*^9, {3.7771899235496063`*^9,
3.7771900314103465`*^9}, {3.7771909323189116`*^9,
3.7771909510833845`*^9}, {3.7771938372914205`*^9, 3.777193915792196*^9}, {
3.777193954969129*^9, 3.777193974695012*^9}, 3.777194178867585*^9, {
3.7771945208125277`*^9, 3.7771945343602047`*^9}, {3.777194620203471*^9,
3.7771946586358624`*^9}, {3.777194823489591*^9, 3.7771948537420063`*^9}, {
3.7771949178206434`*^9, 3.7771949429272194`*^9}, {3.777194994396601*^9,
3.77719500347902*^9}, {3.7771950441500373`*^9, 3.7771950478577595`*^9}, {
3.7771950984896555`*^9, 3.7771950996749277`*^9}, 3.7771951436958847`*^9, {
3.7771952412769413`*^9, 3.777195266360532*^9}, 3.777195673281556*^9, {
3.777196023254566*^9, 3.7771960257550306`*^9}, 3.7771962108413315`*^9,
3.7771964348207383`*^9, {3.777196571470793*^9, 3.777196603700994*^9}, {
3.77719671640376*^9, 3.777196744655404*^9}, {3.777197041016348*^9,
3.7771970857498674`*^9}, {3.7771979803363595`*^9,
3.7771979809319925`*^9}, {3.7771985013183146`*^9,
3.7771985038037896`*^9}, {3.77719854132374*^9, 3.7771985416865177`*^9},
3.778066242768161*^9, 3.7784834449847374`*^9, {3.7784837002916636`*^9,
3.77848370204801*^9}, {3.780736502102478*^9, 3.7807365167849236`*^9},
3.780736938371971*^9, {3.780737007460211*^9, 3.7807370188105755`*^9}, {
3.780737373989303*^9, 3.78073738205152*^9}, {3.78073744060362*^9,
3.7807374451267796`*^9}, {3.780737768503789*^9, 3.780737788173134*^9},
3.7807378333061094`*^9, {3.7807408103030705`*^9, 3.780740813742545*^9}, {
3.7807417728658915`*^9, 3.7807417796989555`*^9}, {3.7807421629929233`*^9,
3.7807421689971156`*^9}, {3.7807423441213446`*^9,
3.7807423727735453`*^9}, {3.780742403029793*^9, 3.7807424804653244`*^9}, {
3.7807426625065765`*^9, 3.7807428612785363`*^9}, {3.780742911085021*^9,
3.7807429341225705`*^9}, 3.78090925493522*^9, {3.780909432760846*^9,
3.7809094520279646`*^9}, {3.7809097142728496`*^9, 3.780909744665778*^9}, {
3.7809099036772923`*^9, 3.780909904270878*^9}, {3.780910146375532*^9,
3.7809101467318783`*^9}, {3.785573324541329*^9, 3.7855733300553274`*^9}, {
3.7855736052321033`*^9, 3.7855736253168097`*^9}, 3.7856501558932457`*^9, {
3.801636489784712*^9, 3.801636610571928*^9}, {3.8016366543281145`*^9,
3.801636696693453*^9}},ExpressionUUID->
"a19a0059-353a-4eee-b6cb-6009617622b8"],
" = T)"
}], "Text",
CellChangeTimes->{{3.780742134102508*^9, 3.780742153021083*^9}, {
3.785568333310904*^9, 3.7855683688107967`*^9}, {3.8016284610596313`*^9,
3.801628474348863*^9}, {3.8016367912162495`*^9, 3.8016368215922403`*^9}, {
3.8016378909168043`*^9,
3.8016378921215553`*^9}},ExpressionUUID->"4197169a-a930-48cc-80c6-\
583a54e8b7e5"],
Cell[BoxData[
RowBox[{
RowBox[{"equaltemp", " ", "=", " ", "False"}], ";"}]], "Input",
CellChangeTimes->{{3.801636827215705*^9, 3.801636835902751*^9}},
CellLabel->
"In[155]:=",ExpressionUUID->"d634fc4b-464e-49c8-9858-1dc2f28a2578"],
Cell[TextData[{
StyleBox["Put ",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
StyleBox["True",
FontWeight->"Bold"],
" ",
StyleBox["to show equal flow speed case",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" (",
Cell[BoxData[
SubscriptBox[
OverscriptBox["u", "\[RightVector]"], "i"]],
CellChangeTimes->{{3.7770013311846304`*^9, 3.7770014031530986`*^9}, {
3.7770014479532614`*^9, 3.7770015758575025`*^9}, {3.7770016421087933`*^9,
3.7770017446245146`*^9}, {3.777001805143985*^9, 3.7770018234143457`*^9}, {
3.7770018576373453`*^9, 3.777001868701682*^9}, {3.777001960960461*^9,
3.777002046364682*^9}, {3.777002085201646*^9, 3.777002234055195*^9}, {
3.777002276610862*^9, 3.777002277305982*^9}, {3.7770023369077373`*^9,
3.77700234647639*^9}, {3.777004730111148*^9, 3.7770047367489705`*^9}, {
3.777004775713829*^9, 3.777004794393256*^9}, {3.7770048782432184`*^9,
3.7770048867193995`*^9}, {3.777005135356024*^9, 3.777005138503594*^9}, {
3.7770062715015907`*^9, 3.7770062725989447`*^9}, {3.777006388216528*^9,
3.777006388494771*^9}, {3.7770065328895864`*^9, 3.7770065333972597`*^9}, {
3.77700781314746*^9, 3.7770079815532436`*^9}, 3.777008100965679*^9, {
3.777008249647401*^9, 3.7770083143522425`*^9}, {3.777008496651662*^9,
3.7770085508750143`*^9}, {3.777008733992568*^9, 3.7770087984972696`*^9}, {
3.777008944898369*^9, 3.7770089526989555`*^9}, {3.777008994553643*^9,
3.7770090420542374`*^9}, {3.7770100865381565`*^9,
3.7770101090553946`*^9}, {3.7770102036221323`*^9, 3.777010224623185*^9}, {
3.7770103190307875`*^9, 3.7770103367719975`*^9}, {3.7770103799261312`*^9,
3.7770103895797234`*^9}, {3.7770104279021797`*^9,
3.7770104668006506`*^9}, {3.777010654236351*^9, 3.777010701877531*^9}, {
3.7770107342543592`*^9, 3.777010744769043*^9}, {3.777010776904726*^9,
3.7770107779093328`*^9}, {3.7770108212382236`*^9, 3.777010822369748*^9}, {
3.7770108546333265`*^9, 3.7770108653657227`*^9}, {3.7770116249534693`*^9,
3.7770116287053666`*^9}, {3.777011780344282*^9, 3.777011784597994*^9}, {
3.777011840045618*^9, 3.777011850205907*^9}, {3.777011886060293*^9,
3.777011899008877*^9}, {3.7770119986424265`*^9, 3.7770119988123016`*^9}, {
3.777012461300623*^9, 3.7770125885732007`*^9}, 3.7770126235147543`*^9, {
3.7770127558962736`*^9, 3.7770128074602027`*^9}, {3.777013115598998*^9,
3.777013129585382*^9}, {3.77701338321282*^9, 3.777013426523162*^9}, {
3.777013501281876*^9, 3.7770135014224634`*^9}, {3.7770135420890017`*^9,
3.777013551428228*^9}, {3.7770135883321295`*^9, 3.777013599872896*^9}, {
3.777013898564535*^9, 3.7770139001422243`*^9}, {3.7770148684195857`*^9,
3.777014899341736*^9}, {3.7770149301506357`*^9, 3.7770149321735177`*^9}, {
3.777014974464608*^9, 3.777014989319892*^9}, {3.777015078648814*^9,
3.777015105074156*^9}, {3.777015138841237*^9, 3.7770152518462467`*^9}, {
3.777015428301572*^9, 3.7770155043138466`*^9}, {3.7770157862401342`*^9,
3.777015856140811*^9}, {3.7770158895212708`*^9, 3.777015966786131*^9}, {
3.777016078005509*^9, 3.777016081899946*^9}, {3.7770174394982657`*^9,
3.7770174801726165`*^9}, {3.7770175363251467`*^9,
3.7770175889580936`*^9}, {3.7770176851513367`*^9,
3.7770177307735434`*^9}, {3.7770177809950657`*^9, 3.777017790951309*^9}, {
3.7770178493678446`*^9, 3.7770178796944556`*^9}, {3.777018150556406*^9,
3.77701815689841*^9}, {3.7770184764393396`*^9, 3.777018497824078*^9}, {
3.777018783772696*^9, 3.7770189402198563`*^9}, {3.777019077844077*^9,
3.777019093371052*^9}, {3.7770191626527576`*^9, 3.777019206255988*^9},
3.7770193795090685`*^9, 3.7770199771218805`*^9, 3.77702004620735*^9,
3.7770203335764894`*^9, {3.777020373131014*^9, 3.7770204167059307`*^9}, {
3.777020452412899*^9, 3.777020458414177*^9}, 3.777020580158021*^9,
3.7770208604724264`*^9, {3.7770214239901314`*^9, 3.7770215587137823`*^9},
3.77702235543779*^9, {3.7770224889440002`*^9, 3.7770225509347467`*^9}, {
3.777022594285115*^9, 3.777022616208768*^9}, {3.7770226959211955`*^9,
3.777022724936646*^9}, 3.7770227641543627`*^9, {3.777022826108985*^9,
3.7770230346753197`*^9}, {3.7770230729548206`*^9,
3.7770231418733463`*^9}, {3.7770231793220196`*^9, 3.777023317276226*^9}, {
3.7770238133730392`*^9, 3.7770238175017853`*^9}, {3.7770238502622967`*^9,
3.777023852922696*^9}, {3.7770248992730246`*^9, 3.777024954407259*^9},
3.7770249998332295`*^9, {3.777025123671203*^9, 3.777025184243395*^9}, {
3.7770254800369577`*^9, 3.777025496144828*^9}, {3.7770255781193275`*^9,
3.77702559318647*^9}, 3.7770256844670024`*^9, 3.777025790958192*^9,
3.777025823845621*^9, {3.7770260285500045`*^9, 3.777026047321204*^9},
3.7770260901824503`*^9, 3.7770263240340824`*^9, {3.7770288035891037`*^9,
3.777028806417425*^9}, {3.777028836542922*^9, 3.777028838958272*^9}, {
3.7770306662364397`*^9, 3.777030671737853*^9}, {3.7770310688904696`*^9,
3.7770310744367933`*^9}, {3.7770334816802864`*^9, 3.777033494148533*^9}, {
3.7770335742358685`*^9, 3.7770335836170654`*^9}, 3.7770336203881903`*^9, {
3.7770336652801285`*^9, 3.7770336719325943`*^9}, {3.7770337440514107`*^9,
3.7770337622553825`*^9}, {3.777034374130868*^9, 3.777034385032163*^9}, {
3.777101223083253*^9, 3.7771012581692953`*^9}, 3.777185891378599*^9,
3.7771859723458595`*^9, 3.7771860561144*^9, 3.7771861201130857`*^9, {
3.777186305096449*^9, 3.7771863105111217`*^9}, {3.777186985656376*^9,
3.7771869890912657`*^9}, {3.7771871579795165`*^9,
3.7771871618791203`*^9}, {3.7771876152725973`*^9, 3.7771876271682916`*^9},
3.7771892445617137`*^9, {3.7771899235496063`*^9,
3.7771900314103465`*^9}, {3.7771909323189116`*^9,
3.7771909510833845`*^9}, {3.7771938372914205`*^9, 3.777193915792196*^9}, {
3.777193954969129*^9, 3.777193974695012*^9}, 3.777194178867585*^9, {
3.7771945208125277`*^9, 3.7771945343602047`*^9}, {3.777194620203471*^9,
3.7771946586358624`*^9}, {3.777194823489591*^9, 3.7771948537420063`*^9}, {
3.7771949178206434`*^9, 3.7771949429272194`*^9}, {3.777194994396601*^9,
3.77719500347902*^9}, {3.7771950441500373`*^9, 3.7771950478577595`*^9}, {
3.7771950984896555`*^9, 3.7771950996749277`*^9}, 3.7771951436958847`*^9, {
3.7771952412769413`*^9, 3.777195266360532*^9}, 3.777195673281556*^9, {
3.777196023254566*^9, 3.7771960257550306`*^9}, 3.7771962108413315`*^9,
3.7771964348207383`*^9, {3.777196571470793*^9, 3.777196603700994*^9}, {
3.77719671640376*^9, 3.777196744655404*^9}, {3.777197041016348*^9,
3.7771970857498674`*^9}, {3.7771979803363595`*^9,
3.7771979809319925`*^9}, {3.7771985013183146`*^9,
3.7771985038037896`*^9}, {3.77719854132374*^9, 3.7771985416865177`*^9},
3.778066242768161*^9, 3.7784834449847374`*^9, {3.7784837002916636`*^9,
3.77848370204801*^9}, {3.780736502102478*^9, 3.7807365167849236`*^9},
3.780736938371971*^9, {3.780737007460211*^9, 3.7807370188105755`*^9}, {
3.780737373989303*^9, 3.78073738205152*^9}, {3.78073744060362*^9,
3.7807374451267796`*^9}, {3.780737768503789*^9, 3.780737788173134*^9},
3.7807378333061094`*^9, {3.7807408103030705`*^9, 3.780740813742545*^9}, {
3.7807417728658915`*^9, 3.7807417796989555`*^9}, {3.7807421629929233`*^9,
3.7807421689971156`*^9}, {3.7807423441213446`*^9,
3.7807423727735453`*^9}, {3.780742403029793*^9, 3.7807424804653244`*^9}, {
3.7807426625065765`*^9, 3.7807428612785363`*^9}, {3.780742911085021*^9,
3.7807429341225705`*^9}, 3.78090925493522*^9, {3.780909432760846*^9,
3.7809094520279646`*^9}, {3.7809097142728496`*^9, 3.780909744665778*^9}, {
3.7809099036772923`*^9, 3.780909904270878*^9}, {3.780910146375532*^9,
3.7809101467318783`*^9}, {3.785573324541329*^9, 3.7855733300553274`*^9}, {
3.7855736052321033`*^9, 3.7855736253168097`*^9}, 3.7856501558932457`*^9, {
3.801636489784712*^9, 3.801636610571928*^9}, {3.8016366543281145`*^9,
3.801636696693453*^9}},ExpressionUUID->
"023b91d8-59d7-4a80-b07c-cfb61b3afd61"],
" = ",
Cell[BoxData[
SubscriptBox[
OverscriptBox["u", "\[RightVector]"], "j"]],
CellChangeTimes->{{3.7770013311846304`*^9, 3.7770014031530986`*^9}, {
3.7770014479532614`*^9, 3.7770015758575025`*^9}, {3.7770016421087933`*^9,
3.7770017446245146`*^9}, {3.777001805143985*^9, 3.7770018234143457`*^9}, {
3.7770018576373453`*^9, 3.777001868701682*^9}, {3.777001960960461*^9,
3.777002046364682*^9}, {3.777002085201646*^9, 3.777002234055195*^9}, {
3.777002276610862*^9, 3.777002277305982*^9}, {3.7770023369077373`*^9,
3.77700234647639*^9}, {3.777004730111148*^9, 3.7770047367489705`*^9}, {
3.777004775713829*^9, 3.777004794393256*^9}, {3.7770048782432184`*^9,
3.7770048867193995`*^9}, {3.777005135356024*^9, 3.777005138503594*^9}, {
3.7770062715015907`*^9, 3.7770062725989447`*^9}, {3.777006388216528*^9,
3.777006388494771*^9}, {3.7770065328895864`*^9, 3.7770065333972597`*^9}, {
3.77700781314746*^9, 3.7770079815532436`*^9}, 3.777008100965679*^9, {
3.777008249647401*^9, 3.7770083143522425`*^9}, {3.777008496651662*^9,
3.7770085508750143`*^9}, {3.777008733992568*^9, 3.7770087984972696`*^9}, {
3.777008944898369*^9, 3.7770089526989555`*^9}, {3.777008994553643*^9,
3.7770090420542374`*^9}, {3.7770100865381565`*^9,
3.7770101090553946`*^9}, {3.7770102036221323`*^9, 3.777010224623185*^9}, {
3.7770103190307875`*^9, 3.7770103367719975`*^9}, {3.7770103799261312`*^9,
3.7770103895797234`*^9}, {3.7770104279021797`*^9,
3.7770104668006506`*^9}, {3.777010654236351*^9, 3.777010701877531*^9}, {
3.7770107342543592`*^9, 3.777010744769043*^9}, {3.777010776904726*^9,
3.7770107779093328`*^9}, {3.7770108212382236`*^9, 3.777010822369748*^9}, {
3.7770108546333265`*^9, 3.7770108653657227`*^9}, {3.7770116249534693`*^9,
3.7770116287053666`*^9}, {3.777011780344282*^9, 3.777011784597994*^9}, {
3.777011840045618*^9, 3.777011850205907*^9}, {3.777011886060293*^9,
3.777011899008877*^9}, {3.7770119986424265`*^9, 3.7770119988123016`*^9}, {
3.777012461300623*^9, 3.7770125885732007`*^9}, 3.7770126235147543`*^9, {
3.7770127558962736`*^9, 3.7770128074602027`*^9}, {3.777013115598998*^9,
3.777013129585382*^9}, {3.77701338321282*^9, 3.777013426523162*^9}, {
3.777013501281876*^9, 3.7770135014224634`*^9}, {3.7770135420890017`*^9,
3.777013551428228*^9}, {3.7770135883321295`*^9, 3.777013599872896*^9}, {
3.777013898564535*^9, 3.7770139001422243`*^9}, {3.7770148684195857`*^9,
3.777014899341736*^9}, {3.7770149301506357`*^9, 3.7770149321735177`*^9}, {
3.777014974464608*^9, 3.777014989319892*^9}, {3.777015078648814*^9,
3.777015105074156*^9}, {3.777015138841237*^9, 3.7770152518462467`*^9}, {
3.777015428301572*^9, 3.7770155043138466`*^9}, {3.7770157862401342`*^9,
3.777015856140811*^9}, {3.7770158895212708`*^9, 3.777015966786131*^9}, {
3.777016078005509*^9, 3.777016081899946*^9}, {3.7770174394982657`*^9,
3.7770174801726165`*^9}, {3.7770175363251467`*^9,
3.7770175889580936`*^9}, {3.7770176851513367`*^9,
3.7770177307735434`*^9}, {3.7770177809950657`*^9, 3.777017790951309*^9}, {
3.7770178493678446`*^9, 3.7770178796944556`*^9}, {3.777018150556406*^9,
3.77701815689841*^9}, {3.7770184764393396`*^9, 3.777018497824078*^9}, {
3.777018783772696*^9, 3.7770189402198563`*^9}, {3.777019077844077*^9,
3.777019093371052*^9}, {3.7770191626527576`*^9, 3.777019206255988*^9},
3.7770193795090685`*^9, 3.7770199771218805`*^9, 3.77702004620735*^9,
3.7770203335764894`*^9, {3.777020373131014*^9, 3.7770204167059307`*^9}, {
3.777020452412899*^9, 3.777020458414177*^9}, 3.777020580158021*^9,
3.7770208604724264`*^9, {3.7770214239901314`*^9, 3.7770215587137823`*^9},
3.77702235543779*^9, {3.7770224889440002`*^9, 3.7770225509347467`*^9}, {
3.777022594285115*^9, 3.777022616208768*^9}, {3.7770226959211955`*^9,
3.777022724936646*^9}, 3.7770227641543627`*^9, {3.777022826108985*^9,
3.7770230346753197`*^9}, {3.7770230729548206`*^9,
3.7770231418733463`*^9}, {3.7770231793220196`*^9, 3.777023317276226*^9}, {
3.7770238133730392`*^9, 3.7770238175017853`*^9}, {3.7770238502622967`*^9,
3.777023852922696*^9}, {3.7770248992730246`*^9, 3.777024954407259*^9},
3.7770249998332295`*^9, {3.777025123671203*^9, 3.777025184243395*^9}, {
3.7770254800369577`*^9, 3.777025496144828*^9}, {3.7770255781193275`*^9,
3.77702559318647*^9}, 3.7770256844670024`*^9, 3.777025790958192*^9,
3.777025823845621*^9, {3.7770260285500045`*^9, 3.777026047321204*^9},
3.7770260901824503`*^9, 3.7770263240340824`*^9, {3.7770288035891037`*^9,
3.777028806417425*^9}, {3.777028836542922*^9, 3.777028838958272*^9}, {
3.7770306662364397`*^9, 3.777030671737853*^9}, {3.7770310688904696`*^9,
3.7770310744367933`*^9}, {3.7770334816802864`*^9, 3.777033494148533*^9}, {
3.7770335742358685`*^9, 3.7770335836170654`*^9}, 3.7770336203881903`*^9, {
3.7770336652801285`*^9, 3.7770336719325943`*^9}, {3.7770337440514107`*^9,
3.7770337622553825`*^9}, {3.777034374130868*^9, 3.777034385032163*^9}, {
3.777101223083253*^9, 3.7771012581692953`*^9}, 3.777185891378599*^9,
3.7771859723458595`*^9, 3.7771860561144*^9, 3.7771861201130857`*^9, {
3.777186305096449*^9, 3.7771863105111217`*^9}, {3.777186985656376*^9,
3.7771869890912657`*^9}, {3.7771871579795165`*^9,
3.7771871618791203`*^9}, {3.7771876152725973`*^9, 3.7771876271682916`*^9},
3.7771892445617137`*^9, {3.7771899235496063`*^9,
3.7771900314103465`*^9}, {3.7771909323189116`*^9,
3.7771909510833845`*^9}, {3.7771938372914205`*^9, 3.777193915792196*^9}, {
3.777193954969129*^9, 3.777193974695012*^9}, 3.777194178867585*^9, {
3.7771945208125277`*^9, 3.7771945343602047`*^9}, {3.777194620203471*^9,
3.7771946586358624`*^9}, {3.777194823489591*^9, 3.7771948537420063`*^9}, {
3.7771949178206434`*^9, 3.7771949429272194`*^9}, {3.777194994396601*^9,
3.77719500347902*^9}, {3.7771950441500373`*^9, 3.7771950478577595`*^9}, {
3.7771950984896555`*^9, 3.7771950996749277`*^9}, 3.7771951436958847`*^9, {
3.7771952412769413`*^9, 3.777195266360532*^9}, 3.777195673281556*^9, {
3.777196023254566*^9, 3.7771960257550306`*^9}, 3.7771962108413315`*^9,
3.7771964348207383`*^9, {3.777196571470793*^9, 3.777196603700994*^9}, {
3.77719671640376*^9, 3.777196744655404*^9}, {3.777197041016348*^9,
3.7771970857498674`*^9}, {3.7771979803363595`*^9,
3.7771979809319925`*^9}, {3.7771985013183146`*^9,
3.7771985038037896`*^9}, {3.77719854132374*^9, 3.7771985416865177`*^9},
3.778066242768161*^9, 3.7784834449847374`*^9, {3.7784837002916636`*^9,
3.77848370204801*^9}, {3.780736502102478*^9, 3.7807365167849236`*^9},
3.780736938371971*^9, {3.780737007460211*^9, 3.7807370188105755`*^9}, {
3.780737373989303*^9, 3.78073738205152*^9}, {3.78073744060362*^9,
3.7807374451267796`*^9}, {3.780737768503789*^9, 3.780737788173134*^9},
3.7807378333061094`*^9, {3.7807408103030705`*^9, 3.780740813742545*^9}, {
3.7807417728658915`*^9, 3.7807417796989555`*^9}, {3.7807421629929233`*^9,
3.7807421689971156`*^9}, {3.7807423441213446`*^9,
3.7807423727735453`*^9}, {3.780742403029793*^9, 3.7807424804653244`*^9}, {
3.7807426625065765`*^9, 3.7807428612785363`*^9}, {3.780742911085021*^9,
3.7807429341225705`*^9}, 3.78090925493522*^9, {3.780909432760846*^9,
3.7809094520279646`*^9}, {3.7809097142728496`*^9, 3.780909744665778*^9}, {
3.7809099036772923`*^9, 3.780909904270878*^9}, {3.780910146375532*^9,
3.7809101467318783`*^9}, {3.785573324541329*^9, 3.7855733300553274`*^9}, {
3.7855736052321033`*^9, 3.7855736253168097`*^9}, 3.7856501558932457`*^9, {
3.801636489784712*^9, 3.801636610571928*^9}, {3.8016366543281145`*^9,
3.801636696693453*^9}},ExpressionUUID->
"14a5c3e4-7761-464f-b417-41884354734f"],
")"
}], "Text",
CellChangeTimes->{{3.780742134102508*^9, 3.780742153021083*^9}, {
3.785568333310904*^9, 3.7855683688107967`*^9}, {3.8016284610596313`*^9,
3.801628474348863*^9}, {3.8016367912162495`*^9, 3.8016368500141296`*^9}, {
3.801637663926503*^9,
3.8016377206169195`*^9}},ExpressionUUID->"87b809eb-ceeb-48eb-a8c5-\
8d3e3b21d8fd"],
Cell[BoxData[
RowBox[{
RowBox[{"equalflow", " ", "=", " ", "False"}], ";"}]], "Input",
CellChangeTimes->{{3.801636827215705*^9, 3.801636835902751*^9}, {
3.8016377381280966`*^9, 3.80163773861451*^9}, {3.801638168410508*^9,
3.8016381709190054`*^9}, {3.8016383444024067`*^9, 3.8016383452908263`*^9}},
CellLabel->
"In[156]:=",ExpressionUUID->"c1d7c452-cb42-411b-95b7-387752606bf9"],
Cell[TextData[{
StyleBox["Integral names to evaluate (see ",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
StyleBox["Readme.md",
FontColor->RGBColor[0., 0., 0.]],
StyleBox[" for descriptions):",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
"\n",
StyleBox["Just put the name in the r.h.s. of the",
FontColor->RGBColor[0.5019607843137255, 0., 0.]],
" ",
StyleBox["Integral = ...;",
FontWeight->"Bold"],
" ",
StyleBox["definition above.",
FontColor->RGBColor[0.5019607843137255, 0., 0.]]
}], "Text",
CellChangeTimes->{{3.780742194274277*^9, 3.780742239456069*^9}, {
3.7807423116479416`*^9, 3.780742331124953*^9}, {3.7807424974772024`*^9,
3.7807425679381075`*^9}, {3.7807426178802958`*^9, 3.7807426449684772`*^9},
3.7854849303319407`*^9, {3.8016376367579803`*^9, 3.8016376377637725`*^9}, {
3.8498423185326757`*^9, 3.849842327662681*^9}, {3.849842455415328*^9,
3.8498424561521273`*^9}},ExpressionUUID->"8272f064-5333-4daf-b534-\
815de67bddac"],
Cell[BoxData[{
RowBox[{
RowBox[{"c", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"csonin", "=",
RowBox[{"1", "+",
RowBox[{
SubscriptBox["\[Alpha]", "i"],
RowBox[{"(",
RowBox[{
RowBox[{"1", "/", "2"}], "-",
RowBox[{
RowBox[{
RowBox[{
SubscriptBox["v", "i"], ".",
SubscriptBox["v", "i"]}], "/", "3"}], "/",
SubscriptBox["\[Theta]", "i"]}], "+",
RowBox[{
RowBox[{"1", "/", "30"}], "*",
RowBox[{"(",
RowBox[{
SubscriptBox["v", "i"], ".",
SubscriptBox["v", "i"]}], ")"}], "*",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
SubscriptBox["v", "i"], ".",
SubscriptBox["v", "i"]}], ")"}], "/",
SubscriptBox["\[Theta]", "i"]}], "/",
SubscriptBox["\[Theta]", "i"]}]}]}], ")"}]}], "+",
RowBox[{"0", "*",
SubscriptBox["\[Alpha]", "j"],
RowBox[{"(",
RowBox[{
RowBox[{"1", "/", "2"}], "-",
RowBox[{
RowBox[{
RowBox[{
SubscriptBox["v", "j"], ".",
SubscriptBox["v", "j"]}], "/", "3"}], "/",
SubscriptBox["\[Theta]", "j"]}], "+",
RowBox[{
RowBox[{"1", "/", "30"}], "*",
RowBox[{"(",
RowBox[{
SubscriptBox["v", "j"], ".",
SubscriptBox["v", "j"]}], ")"}], "*",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
SubscriptBox["v", "j"], ".",
SubscriptBox["v", "j"]}], ")"}], "/",
SubscriptBox["\[Theta]", "j"]}], "/",
SubscriptBox["\[Theta]", "j"]}]}]}], ")"}]}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"p", "=",
RowBox[{
RowBox[{"TraditionalForm", "[",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"i", "*",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"]}], "+",
RowBox[{"j", "*",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}]}], ")"}], "/",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}], "]"}], "+",
RowBox[{
RowBox[{
RowBox[{"TraditionalForm", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "/",
RowBox[{"TraditionalForm", "[",
SuperscriptBox[
RowBox[{"Abs", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "2"],
"]"}]}], "*",
RowBox[{
RowBox[{"(",
RowBox[{
SubscriptBox["v", "k"], "-", "ud"}], ")"}], ".",
RowBox[{"(",
RowBox[{"wd", "-", "w"}], ")"}]}]}]}]}], ";"}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"r", "=",
RowBox[{
RowBox[{"TraditionalForm", "[",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "]"}], "+",
RowBox[{
RowBox[{
RowBox[{"TraditionalForm", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "/",
RowBox[{"TraditionalForm", "[",
SuperscriptBox[
RowBox[{"Abs", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "2"],
"]"}]}], "*",
RowBox[{
SubscriptBox["v", "i"], ".",
RowBox[{"(",
RowBox[{"wd", "-", "w"}], ")"}]}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"rres", "=",
RowBox[{
RowBox[{
RowBox[{"TraditionalForm", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "/",
RowBox[{"TraditionalForm", "[",
SuperscriptBox[
RowBox[{"Abs", "[",
RowBox[{
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "i"], "-",
SubscriptBox[
RowBox[{"OverVector", "[", "u", "]"}], "j"]}], "]"}], "2"], "]"}]}],
"*",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"-",
SuperscriptBox[
SubscriptBox["v", "i"], "\[Dagger]"]}], "+",
SubscriptBox["v", "i"]}], ")"}], ".",
RowBox[{"(",
RowBox[{"wd", "-", "w"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"b", "=",
RowBox[{
RowBox[{
SubscriptBox["v", "k"], ".",
SubscriptBox["v", "k"]}], "/", "3"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"bsonin", "=",
RowBox[{"b", "*", "csonin"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dagg", "=",
RowBox[{
RowBox[{
SubscriptBox["v", "i"], ".",
SubscriptBox["v", "i"]}], "/", "3"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dres", "=",
RowBox[{
RowBox[{
RowBox[{"-",
RowBox[{
SuperscriptBox[
SubscriptBox["v", "i"], "\[Dagger]"], ".",
SuperscriptBox[
SubscriptBox["v", "i"], "\[Dagger]"]}]}], "/", "3"}], "+",
RowBox[{
RowBox[{
SubscriptBox["v", "i"], ".",
SubscriptBox["v", "i"]}], "/", "3"}]}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"daggsonin", "=",
RowBox[{"dagg", "*", "csonin"}]}], ";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"dressonin", "=",
RowBox[{"dres", "*", "csonin"}]}], ";"}], "\n",
RowBox[{
RowBox[{"bp", "=",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "*",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "/", "15"}], "/",
SubscriptBox["\[Theta]", "i"]}], "/",
SubscriptBox["\[Theta]", "i"]}]}], "-",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "*",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "/", "15"}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "/",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}], "/",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"i", "*",
SubscriptBox["\[Theta]", "i"]}], "-",
RowBox[{"j", "*",
SubscriptBox["\[Theta]", "j"]}]}], ")"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{"i", "*",
SubscriptBox["\[Theta]", "i"]}], "-",
RowBox[{"j", "*",
SubscriptBox["\[Theta]", "j"]}]}], ")"}]}], "+",
RowBox[{"2", "*",
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "/", "9"}], "*",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"i", "*", "i"}], "-",
RowBox[{"j", "*", "j", "*",
RowBox[{
SubscriptBox["\[Theta]", "j"], "/",
SubscriptBox["\[Theta]", "i"]}]}], "+",
RowBox[{"i", "*", "j"}], "-",
RowBox[{"i", "*", "j", "*",
RowBox[{
SubscriptBox["\[Theta]", "j"], "/",
SubscriptBox["\[Theta]", "i"]}]}]}], ")"}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "/",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}], "/",
RowBox[{"(",
RowBox[{"i", "+", "j"}], ")"}]}]}]}]}], ";"}], "\n",
RowBox[{
RowBox[{"daggpi", "=",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "/", "15"}]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "/", "15"}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}]}], "+",
RowBox[{"2", "*",
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "/", "9"}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}]}]}], ")"}], "/",
SubscriptBox["\[Theta]", "i"]}], "/",
SubscriptBox["\[Theta]", "i"]}]}], ";"}], "\n",
RowBox[{
RowBox[{"daggpj", "=",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"u", ".", "u"}], ")"}], "/", "15"}]}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "*",
RowBox[{
RowBox[{"(",
RowBox[{"w", ".", "w"}], ")"}], "/", "15"}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "j"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "i"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}], "*",
RowBox[{
SubscriptBox["\[Theta]", "j"], "/",
RowBox[{"(",
RowBox[{
SubscriptBox["\[Theta]", "i"], "+",
SubscriptBox["\[Theta]", "j"]}], ")"}]}]}], "-",
RowBox[{"2", "*",
RowBox[{"(",