-
Notifications
You must be signed in to change notification settings - Fork 0
/
FAST3D.nb
3442 lines (3391 loc) · 161 KB
/
FAST3D.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.3' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 160837, 3434]
NotebookOptionsPosition[ 156620, 3369]
NotebookOutlinePosition[ 156962, 3384]
CellTagsIndexPosition[ 156919, 3381]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"hullcurve", " ", "=", " ",
RowBox[{"2.", "*",
RowBox[{
RowBox[{"Abs", "[",
RowBox[{"x", "/", "2."}], "]"}], "^", "3"}]}]}], "\[IndentingNewLine]",
RowBox[{"hull", " ", "=", " ",
RowBox[{"(",
RowBox[{
RowBox[{"x", "\[GreaterEqual]",
RowBox[{"-", "2."}]}], "&&",
RowBox[{"x", "\[LessEqual]", "2."}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0."}], "&&",
RowBox[{"y", "\[LessEqual]", "10."}], "&&",
RowBox[{"z", "\[GreaterEqual]", "hullcurve"}], "&&",
RowBox[{"z", "\[LessEqual]", "2."}]}], ")"}]}], "\[IndentingNewLine]",
RowBox[{"hullreg", " ", "=", " ",
RowBox[{"ImplicitRegion", "[",
RowBox[{"hull", ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{"hulldensity", " ", "=", " ", "300."}], "\[IndentingNewLine]",
RowBox[{"waterdensity", " ", "=", " ", "1000."}], "\[IndentingNewLine]",
RowBox[{"ballastdensity", " ", "=", " ",
RowBox[{"5000000.", "/",
RowBox[{"(",
RowBox[{"81.", "*", "\[Pi]"}], ")"}]}]}], "\[IndentingNewLine]",
RowBox[{"ballast", " ", "=", " ",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"x", "^", "2"}], "+",
RowBox[{
RowBox[{"(",
RowBox[{"z", "-", "0.09"}], ")"}], "^", "2"}]}], "\[LessEqual]",
RowBox[{"0.09", "^", "2"}]}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0."}], "&&",
RowBox[{"y", "\[LessEqual]", "10."}]}], ",", " ",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{"boatdens", " ", "=", " ",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"ballastdensity", "+", "hulldensity"}], ",", " ",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
"ballast"}]}], "}"}], "}"}], ",", " ", "hulldensity"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{"Show", "[",
RowBox[{
RowBox[{"RegionPlot3D", "[",
RowBox[{"hullreg", ",",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Yellow", ",",
RowBox[{"Opacity", "[", "0.5", "]"}]}], "]"}]}], ",",
RowBox[{"Mesh", "\[Rule]", "None"}]}], "]"}], ",", " ",
RowBox[{"RegionPlot3D", "[",
RowBox[{"ballast", ",", " ",
RowBox[{"PlotStyle", "\[Rule]",
RowBox[{"Directive", "[",
RowBox[{"Blue", ",",
RowBox[{"Opacity", "[", "0.5", "]"}]}], "]"}]}], ",",
RowBox[{"Mesh", "\[Rule]", "None"}]}], "]"}]}], "]"}]}], "Input",
CellChangeTimes->{{3.7590709082786827`*^9, 3.759070976680941*^9}, {
3.7590710164343243`*^9, 3.759071086100382*^9}},
CellLabel->"In[1]:=",ExpressionUUID->"2f6ff26d-771a-4b03-a8f4-d7ffe4270b88"],
Cell[BoxData[
RowBox[{"0.25`", " ",
SuperscriptBox[
RowBox[{"Abs", "[", "x", "]"}], "3"]}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529499645`*^9},
CellLabel->"Out[1]=",ExpressionUUID->"30f5de73-8294-4afc-9cc2-20e243ffe95c"],
Cell[BoxData[
RowBox[{
RowBox[{"x", "\[GreaterEqual]",
RowBox[{"-", "2.`"}]}], "&&",
RowBox[{"x", "\[LessEqual]", "2.`"}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.25`", " ",
SuperscriptBox[
RowBox[{"Abs", "[", "x", "]"}], "3"]}]}], "&&",
RowBox[{"z", "\[LessEqual]", "2.`"}]}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529499645`*^9},
CellLabel->"Out[2]=",ExpressionUUID->"fad5957c-29a8-4e1f-8450-ca9eb2b8f180"],
Cell[BoxData[
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{"x", "\[GreaterEqual]",
RowBox[{"-", "2.`"}]}], "&&",
RowBox[{"x", "\[LessEqual]", "2.`"}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "10.`"}], "&&",
RowBox[{"z", "\[GreaterEqual]",
RowBox[{"0.25`", " ",
SuperscriptBox[
RowBox[{"Abs", "[", "x", "]"}], "3"]}]}], "&&",
RowBox[{"z", "\[LessEqual]", "2.`"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529499645`*^9},
CellLabel->"Out[3]=",ExpressionUUID->"51032e18-4cf9-4290-b3b1-24843425d81f"],
Cell[BoxData["300.`"], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529499645`*^9},
CellLabel->"Out[4]=",ExpressionUUID->"e3c3d0c1-604f-4de0-806a-d9a244d33388"],
Cell[BoxData["1000.`"], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529656177`*^9},
CellLabel->"Out[5]=",ExpressionUUID->"61edbee9-fe39-4f6e-9df9-a18db0c9882f"],
Cell[BoxData["19648.758406406832`"], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529656177`*^9},
CellLabel->"Out[6]=",ExpressionUUID->"737c93fb-ee98-47e3-9644-7ca4ee0ca538"],
Cell[BoxData[
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"-", "0.09`"}], "+", "z"}], ")"}], "2"]}], "\[LessEqual]",
"0.0081`"}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "10.`"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.7591925529656177`*^9},
CellLabel->"Out[7]=",ExpressionUUID->"f24be689-1806-4151-95d5-0dc7cf81ac31"],
Cell[BoxData[
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{"19948.758406406832`",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{
RowBox[{
SuperscriptBox["x", "2"], "+",
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"-", "0.09`"}], "+", "z"}], ")"}], "2"]}],
"\[LessEqual]", "0.0081`"}], "&&",
RowBox[{"y", "\[GreaterEqual]", "0.`"}], "&&",
RowBox[{"y", "\[LessEqual]", "10.`"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}]},
{"300.`",
TagBox["True",
"PiecewiseDefault",
AutoDelete->True]}
},
AllowedDimensions->{2, Automatic},
Editable->True,
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.84]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}},
Selectable->True]}
},
GridBoxAlignment->{
"Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}},
"RowsIndexed" -> {}},
GridBoxItemSize->{
"Columns" -> {{Automatic}}, "ColumnsIndexed" -> {}, "Rows" -> {{1.}},
"RowsIndexed" -> {}},
GridBoxSpacings->{"Columns" -> {
Offset[0.27999999999999997`], {
Offset[0.35]},
Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
Offset[0.2], {
Offset[0.4]},
Offset[0.2]}, "RowsIndexed" -> {}}],
"Piecewise",
DeleteWithContents->True,
Editable->False,
SelectWithContents->True,
Selectable->False]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.759192553012484*^9},
CellLabel->"Out[8]=",ExpressionUUID->"be25423b-2559-4e3f-9ec6-3de59b1a085f"],
Cell[BoxData[
Graphics3DBox[{GraphicsComplex3DBox[CompressedData["
1:eJyFWrtuJFUQHYM0EcnAqCWkltgeGn8DWXeICfcTkHajXZaMT9jMAotoA6d8
AGiQCGBl8QWdgBgSWmiRmGRkTGu9XmM8ffvUdZ2qEp2M5syZ03WrTt17+1F9
8uT+gzdms9kPB7PZm7efv60fb9aPr36cTcfHJ/XRyeFF890efnot+NEery+b
8fPDG8H3tM36umGd6Wi/3P/twb+CbxJPdMbzbIYGOhPteaST4jhoWQc46wCH
DuJkHPopnrnoQwc460T4B+MJ3mqRT+gD5/iBc/zAI31/XAtzXuA8LuCRjl/H
pdEBzjoRnvJWmDiB+/kvQh3gY15Osw/rcSAX4s9PR+F3zXmB47zM9/NfBvUt
jQ5w1kmf74nOk/0fn72W+NdTf+g8T317c9PoelVGBzh0tB+qwA9VC309rlXg
21Xg25XoMO7Xq265XsD9/NQyLuYD1/NY3WIe036uzbiAY1zTPKXnvft/NmP8
Dy9p3tslnXsvad7bNof7H45fCT6Gc3Ru9CMceQOePl80jPP8pvlzg/N8onUW
LcbLOPPT16XhA/f5teFzXwPnftTx5/7S+Ud9t1SvukW9dF3qFvl/f/9x+DfV
cdNMv+t6/f6zj9/qMA6fM474gXN9NT4EeK4vcK4vcK7jzf545x+Kv6c4l0aH
6wWc66XjKVvOG9cROOZJP8+14Af7r9/mfllP+dN1hB82jR4v6t6TDuaTF9M8
sxX9cXyHPxkcdWdcz29binMwOPIAnP0wfv1+J35I5+0azR8MH/MD83ne0Ppz
w/f3IVvjKz0u3g9sxT+cT3/d34qvmO+v11vxFfP9dXkrfgNf560K8lYZPtY7
P28rw9c+7MgPtYn///jwuc/vm8/2vx7/Rb49a75Kv4vfPtr/8MfXgoMf4fA/
4/A/4/A/41ivGUe9OE4fH0z83BfA4X/mc19onbnh8zwMHP3CfO4jfd5FEM/C
8NPXpeEDBx91R3/5dS+CcRVGB33HOsBZh/sUOmpduKMD3/p5K42OWi/u6ABn
vu7rs0bHszI44oniRH8xjnWEce3zTvz8aMzT54LD+MD1vIG+6wwf+x/G07Ez
uO7fLujfLujfLuhfjufK4Nhv+3HC513Q7530NefN7/dO+pr5fr930te+Pvd7
J33NfL/fO+lrPx7u90762vfJMoi/COIvAv3a8Ll/9bhKg3Pf6XFVgqPQqG/y
4duomxwRDp8zno5dwL8yOHzo69g4fXww8bMPgcNvzGcfap254bMPgcNvzGcf
6vMugngWhp+OpeED9+MvgvgLw9frS+azD/W4ymC8peHr9SLz2bc6D1WQH+tn
vY5wPCvDT0cd5DPv63z+ucGh/8W0PvD6Ahz/AD5dHwmOdY356DvG07HD/Rbi
Xxk++s7XsfH7+GDGxX0HnOd/xn39ueFz3wHn+Z9xP55FEM/C8NPXpeED9+Mv
gviLQL82fO47Pa7S4NxHelxVoLMK4/H9k+/XMa51evEP9pPAYbQI1z7vjc8Z
1/Ew38bj44OJ0/dzL75lvu/nXvLj67Of+8C3feDPPvBhH/itD3zVB77qA1/1
ga9646vx70+vJQ/pf5cN58evY/bVo3Rd0Oh9y1kDHAUH/svkI/BxXQ8c/Ag/
mfzNeDp2Bkf8Pn/Wcjz+uIbGj38w8bM/oaP3OVkHOOuwb3U8c6MD3I9nbnT0
vijrsM91PAuDp69Lg/N9LcaZz/e1GPfjL0388CfrcB9Bh6+Xwef+0nWsgjpW
QZyrIM8rw0f80bjAn65/id81Gs/9gv2R6tN8P1nzb/Fx//ON+FjuJ4/z5b2X
Zl34la4n8NwK8QLHczE+L/qUcfwfuOqjO/GrvriDK5/f4tALniOIz1mH/a/j
rA0/4eZ5gfE/8OB5gfEtcPanjjP7cLr/T/uHTaPx/HxnrfMiz/sYhz72yXo+
6ck/+f6P1sn5Ac75Ac75AY78cJycN8Sp+70XPnDW4TxDR88DPY2rNjpcF319
ke+b4Xfgqo4PL816N/7t+BWtv+eCgwmc+5T1dV+w/mD09bp2Tvy54Jg39Hsg
20aPKz+XhI72VY6f3zPR81J+7qnzsDQ67EMdf2H47E+tXxs+v3fBOHT0uErD
Zz8D5/coGGc+fMh89qc+78rwMV7gYx1Ps6/GNJxcpPXo2WvxJ/oBOPjA12pf
mu+j6uto896C6GMdAQ4d4Bgf7/cQP3zuxz9I/ODD/8wHrvr/Nh7giEefd250
gPN54X/mA+fz8ntWOv5FEP8iiH9hdNBHrAOcdVBHrgv3o85PEeSnEH2tU/N7
OC33L/TRj6wPnPW5f3VdyqAupdHhvoYO+pd1gLMO97uubxXUtzI67v7hNM8D
/rhWxp/Iv++HOqiXHPj+HPwIR33/AyIm0OA=
"], {
{RGBColor[1, 1, 0], Opacity[0.5], EdgeForm[None],
StyleBox[GraphicsGroup3DBox[Polygon3DBox[CompressedData["
1:eJxdmQm4V3Max8+hhCwTuYSWq1KNJu2RKLRJbosWKlqEpLp1W1EpKSmVbgtN
WiSRJSk01TALM2bMvhmDMcZYmsxYBmMwmPfz/D7nufeZ5/F9zu/3/t71+/7+
V+c9pWPKB046JMuyWYdmWfyXNQzkgXMCyDsEmgVqukbWPHBYVqXPeRv3h6lb
S71CxvoIdfHzzcDhniFvpU5Ln9gfqc5R2tTxrNCto21zdY82n2Ot4SjP65jT
Yfo7Rv3axqntvrk+jgsc7/oY7b+hv6NcdxA19Uncusbm/ATPeZ7kukGgXuBk
9RsFTtGuoetTA6Xqn1gtF3Ir0b6VeicpR1ZfvyXatfKc+K0DTQMHlZdq3yRw
hrEam0MbdZqod4L1IzvNOA3Ua2YfmytvLU4zXiv9N1H/dOM0t64SYzQxB3Ra
qNdB228F+gS6By7Qlhgt1cWmrfrtzOcMY2N7caBvYJtn7bThrGM1O57tRQe5
ammsM0UrZW20aauPg9rkxuK383SebFp7/kyeZGebU2f9Ie+UJRzUX2GDfkf1
OtkffHfJ0r07N1BDf+cY+2x9d7H/p5rbeeqeKxf4/3ugW5buf9fAWXINz+dn
6e53U9bZOLXUxeZYY1J3T3vURxl3mvvcI0v38UJxvL7rqnescTtr2ytL97iH
8Q4ovyhL96p3lu4Hto29A32zqrvFupF6DfTXVPuDnvU1t232sJ1or6yPcXuq
20I9ZHMCc9U9qP4W5fcFRgaGWs9Z7ln3tyf9PB8VuMSzftV6OsAeDZb7QYGB
ygbYm5H6nhWYHpgaGCKv2F0m10Ose4S5jrRH1HW5HF9mLsMDZeqMNsdLs3Rv
iD9MLi/Xx0hzuTJwhT5GKx/j/gr9jrAfw5QPVbewudLnGDnerg75jDf+uMB1
gesDV6lTpqxM+XXVzlhfLbdj5fZaORznk/019uZqezFWm+Hyi+8JcjFeGXxP
CUx0T8xyuZ8kv2Pkokw98quQ78nuJ2o72V7hY5ocV+h7ijr8FrbJzyxt8T1D
Xqcpn5mlO7E2cFdgtfyjd1CdG9S7yf2N9mqWdS+z1qX2Zpbn9Ghe4Gbls7Ut
9pzPUTZb2Y3as5/r+WjrwNft9mFxYL69u8pe3Go/FtmfhdY1X50F6s/3bKw2
tyhbK5Yoo55xxqxpTGLflqV7QBx+k6PkFts7tGc/1POVWfo7uEKdysCdrle7
X6tsuXoT5HVVlu7ISnW4S+sCd2fpHqwJfNuzde7LtVtr3VNcLze/9Vm6L/i4
177yN+mewKbABtcz3W9SD/nGLN2d9cabru7mLN2ZjeqtU47f+5RtcD1Hf1uU
4f8GZQ3lFr7pP3eF3jeO/2nUCjQK3J+lu3BAcM/5f9OD9u0B19vt+VDXDwUe
s287Ao/a/0fcL1TWWd2HtWseMf8Zz+/ar2cCO7N0D/C3y3wfD+zO0j3ZZR7c
hScCT2bpHu12Tx+eMnZnYz1lj74T2Gv/97jnbuwL7M/S/djrvtK8apoTOTyd
pfuC7veydBeQ/dY8fh34fpbuzQ8Dz2bpLvwg8Jy9fdb9Gn1Q813abdAO+5/I
EXxvzVJvDiin/ucDP8rSvXhOO+wPyVNtS80JHz+V9+fl6Ff24wP7s8M1nP/c
en8mhx9U04XLFzx7TNkL+v6p9dOHXxoD3n/jnr79QjxunOWR63cCy/KkRz9+
J+jBHwIv2pPfu6cnfwy8ZB9edA+Pfwq8bE9ecg+vrwRele+X3cPXnwOv2Y9X
lf1F+ev26jXrpdYPA/8y9w/V/1vgr4E37ccbyt/y7O3Ajz1H/o5nP9SOeP+w
N++63l6t/+/J7fs+35Pz983pDfv/unE2m8dHck2+NfN0XiOen9gPzj+2P/8O
fGr/PnFPD/4T+CxL9/tT9/Tg88AX9uQz9/Tgv4Ev7ckX7unBV4Gv7cmX2u3R
LzX/ylz4Byzcc4cPzVNP8jxxdIhn1EBPOP+r8r8Y4xV9/E0554fliSM4eEs5
5/zNe9Pzd5S/ZbyXzZffHX8n52Xpt3VAvcPz9HfxiHgelacaarvebg+PyVOv
jo7nkXm6B7U8/6fy99X7QHs4OC6exyr7ZjyPz1Ov6sSzfaB1oF2e9OhT3TyB
3pTkqX6eX9mXEs9OylOvToxnvTz1Bxl7entCnmJwX76Rp95h+7ln3AtyOSVP
3JycJy6KmOzpbT1jl+gDX/xWqOnUPPF7inbY9/LZ1fXpgaaBJoFm8tXU/z/R
q4Z54g4OW7intw0Cp9mbI/RRWz9N9UN/Gnt+pD1EvyLQIU91ts1TT1rYg7b2
4MxAyzzVc0aglb6/FajvfT3VM2qld208b+0efwP0VUd/reWinbp1lVFPDX03
sd+cd1S3sOnovRhi7IHukXcOnJWnPncKnGOfznZ/onX3ydP/i5t5Vk87Ypwf
OM/+dHTdzVrPVdbFWB21Oc+zU4x5nn66itJAeaC7dV6gXY/AhYGe9rW7dhdq
V+K6kPVQfoE8kVdv70Fj1430h26van56KSv1HlTIA/fvorzqfhz093pJ7MsC
/fK0Z93fnq7J0znysfHsK5/4a6ENei3tE/1vZRziXez6UvvHXRkgsOuvTh/1
2vjEx+WBoe7b6r+D+/b2pYPyTsoHGav9//kb7NmZ6uNvmHfhssBw78cw953U
u8gcBnsXLjLuCPs/XFkX+b/SXl0RGGn/kFXY85Hqd7W+0dY/Slxs/gPswVXy
dpX50Ie+2vXUL3dhZuDGwA3q4GusNY/S52DXg6xplLKxrq8PXG2/Cxm+Rpvr
KO3GqHdd4Fo5uyYwTj6QjZejcfI5RJ0ixlifE/QzUc6wK1fO+WT5m+i+q+eT
3F9vPmXmTY5TPB+pfYXyy7Wf4r6rPqZUk1+t7Rj9sZ+q/5kCjmcFpgem5ek3
QK9m2Kdp6hU6ZfqY6tl09WbIyyhjXq/NDfosM4+b9Ffmepp6/dS90fVMz2cr
4zknMNf9HP3drN7N6tzkGet56s9xPd8cH1B/jP2a7b7c9S3qz3ONjwX2HB+3
5unfxPx7eIX53hVYlKf7skCd2wILA4vzdHc4n2Cs69S5Xb3Z7m/xfIL5IFuq
3h15ujtLlC/zbLJnFca7xfNNgXsDW8yRXFcG1gdWBTYH7gxUKudv2Er3q9yz
XuN+tWue/cXawHd9znB9p3zcLSqVrdbPem1Xu96iziZzvNN8sF0X+LaoVEYt
91Q7Y71BVCobYk79PaeHG9XZaKxN6m9Qb72czJW3jcZaZz3wPt5+Usd96jNr
ZObG7A/brdV4Z32/qFS2TdtV5knsB8RDgYcD29W/X1l/7fD/oH7vU8YZ3xz4
9sC3FOzn6w8/jwQeFZXKVrku1f8O4/Jc4Pp+Y23X9rHATu12VrN9XDuei1yj
y53dpY+F5rFT2W5Rqgw/T+gD+ZOBp0SpMu74njzdc3wsVY93U36P+/SzR1+c
LfP8KZ97Ramyh/S3RyzRllxvs469+t6fp+9jfNPiG9Y2eawQD9mP4vsf3/72
Gwsft+pjnzoN1SPeYmvnuxrfEY4WzP/59sP8jG+DNXzWdM0Z3wtKBN8C6gnm
+fUF3wJaG5O5f6MsfSfhewnfdviuwr74dsO3lGcE/7biewrfTJirni26KKvh
3S++DXazBnC+NTBfZ5bbQzCzLb5d9DJX0NtcAbN4vkXwLYEZV19j9reOwu+J
xr80q/rewnOQ60HmBAabExhSzb6HuTB3XKE9eQ0zJ2alzMeYqzKv5h2TmXU/
Zf3FNeZWzM0HivHmd6kxivx4TnI9yfxAufmByXIHKuRus79/fvvkN0K+mB0z
z2SeyCx5jntmicyXi9kze2bM87M0M1yg32nWTE2LrAcw6yu+CdxubWCptYFl
1ob9DHPi30G8Y/DvJ2pdaZ3Ut8bamM0yS2O+NtU9M1tmqBs8Q+9u6wfrzXW6
elOtmVnqTM82mgf13+s5cy/mJ0/KzRbP4WWrfaXWndYLNw9mVXNUZmjFrP1R
78S8LM39thpvc5Z+M+g9ou4G4x3QL3O2xWK3nML3Y8bH7w5l1HOPtTJPfcC8
brcOesB8lLnccnn9gXzB3ffllv4wo1vh2bPyiA0zQGayzFP36Y+56ff0hz2z
u9X28GnP8bVfGT6YET2h3+fsAbX+0jrJ/SdyypzyAzm6Q7sl8vdjOWQG+Qt9
4OtHnsPTz+wTds/rm/3PtYOXPdZN314wJpz9NquaLf/G3JlbMSN5RQ5+l1XN
ol+0TvCS9WP/e8/R/be+mH3+QR/04TW5hstXlTFXfMXe4OvlrGrO/Ed9cPYn
Zdi/Lqdw/K41v6MvZjPs/yG/cPOhXMDl2/IJf28qK+aVL+j3Dc/pFXOkX1vP
x8akDx95jt9/KVttvpXafCq/+GBW86k1fZ5VzYq/lEfwtfXD5Reeo/cfeYWD
/+oDDpgjvSGXzHNek2PmMK/ri/nVq3L8lT7I6TPjYM+c6k3rYM7zkTwd6f/z
4JhZFJwyI2TeBF9/V4f5FLzzm+b/8cyumHm9J8fMEd+WP2Zc78ofc6pP7A9z
KuLzO2AehYx7xBzpYzlg1oSMv5v85srlifnPl3LJvOhr62bmAy/UykwJXf4G
87e4m3rMfNBlrnaa/36BD+YwtayLmniXhFfmPzX1W18Z+4bawUVjfcEls4nj
5YR5FXMWuGHOAo/U3VJd5n7MV9Cr7R7e4eMMdWvpH97hgxlbMWtjLlLM+Vrp
6zj3ta2Z9+4LzIl5CfeS+V177eGJ+Uk9wfzkZOvrKS/M4Tqpe6v87FeP+Qmc
dpSzLvahm3zVN34D0V3f2Jyrbok51ZXH3vaFvIaZE7Ml5irF7Ioncw1mUQPl
63T1mpl3b7lDr589oP4O1oTf4dYARpgTtV5mfPwOsWcH/A2UWtNI6+Hfv7xz
FHOw0ebH/KeY/VDHtcYrlcdybXh3nGHMa9RlTsCMoZ89nGy84t14tvnyXnaF
fscpYz9RO+yLd1xkd+iLmuYbDywwP7BIX8wb5uqjmEcwn+A9e6Fn7JllMBdZ
bm77rJEerjLHJeaE38XKuqvTU7015snfTfrD75h3WN5r1xnvHvPhPZP32I2e
8X58l+A9ekVe9d6LDe81vM89aN0PW/NCZcX7Je9+W8UW49CbtXKGzQ45oqY9
csr+cWvj/Y/3p9vcP6Fu8Q5ZvFM+al67PCtsHjEv3uF4l+N9r3iXe1Lszqve
BXcZE973qvs/z9KoNw==
"]]],
Lighting->Automatic]}, {}, {}, {}, {}}],
GraphicsComplex3DBox[CompressedData["
1:eJx1WkuIHFUU7fipmI+LwqbAUDF0WYoICoKIIM5rd7pQjAgiiIgYVyYRNRtX
URClERMEQVpEsDGSWUiwzYzE0cRMTFyoqASMSVyIooPfXQcFwZ6pPre45947
DIScOZx73/289+q913t09307Lup0Oj+s63Qunv574+6XH3n76PhoZ/Zz54/P
dt974EAaH7z55/PbDgl+wxpvOe245Ir/3uwuHF3ceW5h57lf0okD3+5/+q2P
04Y9967fs/K+8NNHe8/cv+V4gj74H171wu/bn/mUdCaio/kT4V+47qF9x79c
IH/GonPXq/X0d12fdYCzXehMvZ7+HhL9xm7WZ33gvp+Z0b9mzfDmPusDZ33g
0Nfjyo0OcOhs+m705O0bF2VcW0/fsfX0hkWjw/oNu2v0gfv63T702e6F7y89
cWDjosnXr3OrBXHE2GV/mngWxh/g7A9w35+6D39Yn+3uWoUXrjR2gbNd4Bxn
4Kzf5LcM6qE0+sBZHzj0T1127b+3ffMB1fMx4Z9fdWfnNmMXONsFznaBs12t
f4zqrSd2oQ/cr89eOK6mfuYpDlUQzyqIZxXEswrs1n3Y1fNk3cc8qeu5Dvqo
Fn9YJ8I3Nn4bHb+PajMvAce8tL3xS+bhxw8v3fT8wTcER5wZxzwMHPM85hPM
86zzcP+J6y9/bZnWha7wgaMfGUcfMY76Zxx1yDjqinG9jozEf9SDH586xJFH
xpFHxpFH7U8dxLPuI54rTT1KHje/dPc7vRf3C45xMQ7+ank/+M8Jwieig3Uf
ONZ95rMO1lnWAc46wH1/MqOj12X2JzM6WE9ZR6+/y4bvjysPxpUH48qNDvLo
56srfN53NXUypPEWho9+YT5w5qOPwGfcj09pdNB3rKPXl2WKW8/wgTMf/ej7
WRm+7q+hiT/6kXH0I8//WNdYn+MDHPH56ovpz9K7lMfBbH4eST9ueeqv8abP
nxMcfODQAb+8+vWVvx/7jHQmyedPDB99x3bRR4yjLyI/p+KHl/5k/dzwgTMf
9e/rdw0f9e/7XwgfcdB90cYfOOsAZx3dLwPys+5z/HVfDATX9T8g/crguj6t
XfYTOPLe1OdQ6me6uEw93SU48sg46odx1I/PzwyO+onsroVp30nSzwN+N7Bb
GBx59P0pDY58+fwq8KcWnNdx7Csa/kDwtX/+OGl0dBwGEh/NH0icGUecGUc8
59YmvnsMznz4w3zEH/xZHSauw4Z/9hPU2ys/3Tr9vWUvcAwJOP7P/N9m+278
Hetmh36a+fxUYh3g2p9J8u1OjJ/Aff2J0UdfsD7jzPf9zIwO48z3/cyMPuqE
9Rlnvj/ePBhvHow3NzqoK9ZhnPl+3IogbkUQt8LoYN5gHcaZ78e5DOJcBnEu
jQ7mJdZhnPl+vnpBvnpBvnrBuKpgXFUwrop15iI/G7U6yHttdGBf8wcJ8wDN
Y9IfjIO/9keZD8/O+fMJ+BPSB5/nh0FCv7O+Pw9AP+v7+tzXg4Q+Zb7fv/CH
1wX4w32KuHUN3+9H+F8E/nPfDRL6iPl+f2G8ZTBe7qNBQr1x3lHnjKNu/frh
dRxx4PocSr3p/UO7vmv9odQJ7TdEh/ZLKdh3pWDfJXn39c3+KgX7qxTsr1Kw
v0rB/kriz+NCPH0/2/0S9iHB95TELcLp+0L6mnHkxednBg++mwT39c33keTF
t1sYnL9ftN3S4MiL708v8Kf93vm6+e6Q+OM7BTjiDBz7Oj6/wvkZ+HRuKTjO
wVgH391sF99BvJ9kPvLu+9mej4GPemA+n4Np/czw+bwLfNQJ8/lcS/uTB/7k
ho+64nyhrnw/C6PD39fgAwefccRfj6s0fNSn738dxKcyOuBHOqiT2TkP1bOc
O1E9D6U+g3N14QOfneeTzhFjl861hM+4rlv2pz23Zz7r6HpudYCzDp/rav3M
6AD3/cmMjq7/ocF9/3PD53NaHc+u4QP3x1UYvr6XbPl8Hqv9L018gHN8+HwV
fOB+Xsz5atLnS2y3kvsFHYfa8IH7dV6HOojP+uachOp8JLiu85GsC+Dj3o3X
C+DcX+CxXdyzsF3Gdb+wPxPxB3zdFy0fOPuP+mcd4OyPrv8RxaE2OuCzju6L
kcH9+HQNH7jPrw2f7/vgP+of6z5w1DnjqHPGUc/Aed+O/TavC1iPOJ7Q4XHB
/9l9Aa0X84LreV7urwWn+3HB6T2D4PSeQezSewPh07228HGvwX4yrteXebI7
SfCfcbar15dWBzjHATjd1wvujysz+sDpXp7WHXknIDjb5fc/sKvXqdZu8P6H
+nds+H7c8iBu5l0Q9fWY4pkH8cyNXb0+ct2a90I0P4ypbrvGLnA/j0WQR/Mu
iNZfzmMR5LEwdvV63doN3gWl4F1QCt4FJX4XxHa57/Q+gevKvCMSnP0J3hEJ
7tdzafzR+5DWn+B9UQreFyV+X6TrxLzXSsG7I9rncF/0gr7oBX1h3iMJHtn1
81UF+TLvlAT389W+U9L9Zd7/CO73b2381Drc1+b9Eq3j3Ne1eU/F6w6/L9Xr
cvteJbo31/eY7f1RdA7sn1vOftrzOrmfAp/fDWJcHfqhc2k5H/Zxc68k570+
39zvyPmtr2PuWeQ81tc39yBzwb2GnKMyHtw7zKEvfB1zLyDnpf69T4wjX/8D
bjlL1Q==
"], {
{RGBColor[0, 0, 1], Opacity[0.5], EdgeForm[None],
StyleBox[GraphicsGroup3DBox[Polygon3DBox[CompressedData["
1:eJxNmgu4llMWx79Xp3KpDkqU6khR0v0yaqKSoiR0UkQpRlLUUbqPM04qdDSR
CqUbpUgXp0IulRJmmMidaRhmXOY+wzCDY8ysn/f3zvmex3r2Xmv912Wvvc/X
u9fW9MqS4gmH5HK5x6vlcvFfLgmC7xh0aNAR8gXyzOs5b5BLbWoENQxqFlQU
1EK7juL3JSkOzCn6ej5J59XUN1TOvIHzDF8jL85hyguDagUdGXS4uJrmVisv
90Pla4qp7nhYnr/aebFqK6un71p5NuSM/ijxdZQfbqxC18waCtSdIr62c2x6
Oj9GXf2gukHHSk2CGqs7xlj1xJDT0dp0VHacdvg8Xpsi5Q3ENhLfSGxjfTYQ
1zyoZdAZ2jZ0flLQiUEnBzXV7kT9NNVnoTl1NIcTzONo862rr+PV4aOZfqAW
ef5PyvPTXHwTc2om31x/x+gz27uW6grcb/atQ1DboGnGaB10qtg2+mml7lRz
KRLXRr9txbQzt9Zi2ovpoJ92+ujsHrfVvov4XyWpfJo6cJ20a+e8o36LtAP7
I2N00aa9cnwNCuqq3WnSWmXdnJ9mXu217+b4Y3Wd8rDIuqsnX/jTnZPTucYk
dg/zYO0DgnoFnamsh3kXiK0j30M8vwuJdgXyfZwf7v4ViZ2mvI66Avf2ef/O
sDsrl56dVtaiQBmxi4P6B/VzfnZQX8dzxJBHb9dwlvqeebjMzzl5dv2U97cu
PbUt1q5YeX/n/fLs+ufFPc/YxdYWm7uDKoNudF/Yk/ODLnQ/qMkFyuGHWMMX
khSzOWiq/rvqY5Axplkz8BcHXSRujDRY2WbjTdMv+BeTFNPLnEcEjdIf8mFB
Q/V7WdClzi/Nwww1Bv4uUY98uPaX6reP4xXGGGOeZzq/Mmik+MvFjlE23BzH
iMf+J/pkfpU89bguqMT1jNUG2Qzlo5WN1eZq9dc4wo/Lk413jmxC0PX6uFZ+
tLEmKBvrONpaL7L2vbUFNyWX7v8Fyjgzk5xPystxtHUh94nmX2ze2TonKmO8
wflUfV+ov4FBk6UF+p0hfqKYycqRTdfPdNdQrO7nQeXaLlLPmeZszwz6qSMY
zkJp0E3ip+sb/WAxA+Qr5X+mD+a3BM3JpeerPFd1xmbn0nM7xLFUWan8bGOW
aoePm4MWK0c/XN0wx1lBZc7LzOMy8Yuly5WV62+WewM/Ms/XUPGcySvUlRtz
sTJ88W8Av6H8Js5VN9u9oi5rtbvNWswTc6sjZ/Nq9+QO9zXTzdX2du2RjdOG
+Xj5xfq4yvl893W+/Dz9luftPWseLd2hnPO22ZzvDFqoHeNduaozt1D9DHNY
rPxOfY0Wv9D5AvkZYheJu1Mdv6uclyXya81jijnBV8SZ2Bu0J0lx7O09Qfe6
h0uD1uWtmz24zxrOdb5OHPu+TOJcLZfYnxUS9V6pzXxruMr8Viuf516tcyQO
fwPPBD0p9n7xzNcEPaBsgfMpSTrf7PorrccCddjw97LeNc12fp/E/CHrUGlN
yPt2812q/F51lcZZ5nyp9g/rr9JabHBdyB/R50ZrsynoQfWrla8St0n5KuUr
Xe9GsZu1/dr1b1RP3C1Bj1qnjerXiV+knJpslH9Qf1u12x+6x639U9ZmW9Av
ktT/Y8bYrmyb61xhPGQvJ6ldhXj8IXtCf4xPOyfeFrG71e2QsNnuup7U11Oe
jZ3WHD+7HHc7fz7oILmIezZoj7Jd4iqVP+PIGp4Tx7jPvdirv+dc3wvu5y51
+5S9KH9Q/CZHcngt6E1r/ZqYjEf/sva/DHrJOL+SR/dKHvagspedvyL+VW32
Oz+g7qAxN+rvBeevi3tF/YPaIH8DWZLWeqf4t9T9HpvQvRr025h/IP+msd4J
ejvo3aD3XOc7+mf8dd78oDz+f5Nnw/w3ntF3jQ3ufW2Rc4bBvJ+kOnx9mkv/
Jv7gGt4yl18G5qMYf2eunLkK+W3qX1XGWeVcvpSkcjCfBH0sbdcO/IEkrccn
xkb/YS71/4GyHcoP+HdFTT+Tp8Z8C2XfSru1If96oTsmqEbQ36ljjG8l6fr+
JObP7EnI3kjSkb3ao353Vockle2zRuj+Yp3+5vygI/xftXnTWH/RpmaS1rzE
PDivL2mL3T/0+Q/5t5NUv18M8i+CPs+l9uj/Ke4dab9nED9fiv8ql+reNSfW
QWz4A+7zQbFfOv9SH/8S9570hufh37mU/3XQwcQz53q/0Qdn+1v3c0oeFgxn
7hv1W8VwDr5T9p9cet44Q5/Ify3mW32vVw7mP56hj8R9Yb2+86x8r5yzcGjQ
f917LpqJ5+OQJOXZv8R8Wet72v8xl9qBYw+qJSmxBwVJWpPqSVqzr82NeNXd
d/Yf3Fd5uOryrOmwmB8ubVXGt1ulaz0q5Ed7no9wzUWOtWKsrbyWsu/1+Z15
oK+TpDLWUeh6C+WRf2vNjkzSeDW0OdJ5klTZkSO5Ngw63rGBNQFziPnW1bZA
eTVrUqgMfT0x/L3Wd14o7jhlxyqvqay6MQrNi33dbO0PFQ/f0DWTG98eZbn0
O4T8Lw/ZSH2cEtRC7ElBzZK07o2tJ3RCktacPWmkvsiaH6FtU2Od6NhUfZF7
0MT5CeKbiSlS1z6ouTJ0nDt+51437lZjn2TeLbSto/8i9wu+lXrW1lKbDH+U
OHycKuYU563EHC3uZHPJYrYWV2S+U4PaKGutvK724Nsqb+fYVuqgjHlH+SLP
AWeiU1CXoM5Bg4IuCLpQfKH6to4dHTkbJY6dxBV5Huo7r28M8v6Rsq5Jeo5P
k9i/HkHd1CGj9t2TtP7MOTvdtcdP9h0I/seeD7CNpdODzgjqqb8e2mJ3rJhe
4pqIZc7ZY36Ca+rsetur76Uvzs2ZQX0d+wSdFXS22H5BQ4POVX+ymF7i2pvD
2fLoODfnKOuX5+dc5+g4L/2d9zd+H/m+7hn5FosdIA77890HZJyhgfo+T31v
13amuZ+nvrE59tLXAGMN0AdzzksH4w+U7+u8nfz5+sx8n2/cwdoVmx/yS4Mm
BF0vHt1FYhmH5M07y48IGh50mTW4SDt0nO2LHUucg+0mnlyvDbpGOWeLc3WJ
9heb1zBzm+oIf7n44frBx1VBP9H3lUGjHa9y7O18VNAVSXp2wfZU31d9hrta
H2Plr9FulPuD/bigMc7HqBvlPPu74DydrWycMnrp9N7puxN3vLVo5B5cJz/e
+EPdl2HOx4m5TvxQ8y/RpkR+kD6v0wYfE4Mm5eHG6mOi+muMPdb9LtHPzCT9
TaD2/JtfoQ6baUHTlbNPu9QNdj8ni7tB/XRthimHnxH0U+PcoHySMWeqm6rN
dGUjtOO9gTeJlv5ugZ0dVBp0Y9CtQbcFlQXNsSY3Bc3SnnNzs/zPxM1ybKS/
m/U5Rzm4W8wPeXdjzBHL2ZmrfI54aro+6H51nLF5+igPul2f8+Tnuwby/7lx
b1M+Uvw8/S7QhvEOdYx36gObJeY61zjktlDft4udr6xMOb7vyvPLfJFYxsVB
d4tfKP4ecejuDVpqrLvFLBKDbpm4pfLLg1aIv0+bpcpWKitTvsy1LTHucv12
N3dqtUr7lcofDFotlelvrHuyXNlE932J61ys7WrxK8Q/oF/GNdrif51xS82x
TBl+12qzxnmpWPxzRm7SH+fkoaANYh8W80hejI36XKccfqvjI/rGvsI1bsrz
y3yzuu5i16jLYm3Vz8NitwQ9Kp/NK4xBzEnmsz3oyaCngh4X+5g+we9UXqEP
cttmTlu0z/BP5K2J+Q51+H7avLZp84w65Ddos0P/T4rj94C/nVvFV5gLOH63
duv72STthUIl4p7T9159Pqvfp8XtFVOhb96ZeYvnXY41TzEG7+u8H/N+vlvZ
VOtSIQ45v5GTXfdOcdlbOPa8I/P+1zSXvm/zPs27N2/bvCM2kY5VNtl14hue
N2rewxuKJ1/els9Q1sIYzcTWNzb/drUydmPjNxLX2Jjc77OeE++w7bTh3yr+
TeTfQd4su6ZXrv+/0XaSeKPN9APE0POjZ0IfkLfkDv7u806LfRexvA2cp5y3
wux9F/50+bXGQ84bEO9MvMXQCx6ij0n64R2C9zbeNXgDGJVL3yegEerxe6G+
wF2sv+y9jTe4q+XxscAYa9WfJWagOt6Z6OfzfjXePHj36u1eTTMmeU/RD7Gn
ymdvP7xh3mRO5EM/nt40PfqZ1ot3z8HyzEvFYl9uPvS7y3LpPZO3AXr5y/U3
K5feO8v0u1QZbw702nn7okdLj4s+Lb3/FfqgJ7zdvaV3Sj/6Pn1Xmusi1zvd
XFaZT/aWdqP8yjzfG/RBPHp3W9xb1rtenl4cfW3u1fQO3zcmMr4vOWf0yugf
0wt+ynUgo+9FL+jZXNpbQ79L2V7nrIc+9C/0v0lfD5nDBu2JmfWX4bMeMT1m
epL0hugL0aOlv/WCMR4xxkbtN6jbJ8/f30vasVZ80Xukb0r/hx4o/dT3lP3Q
a01SGf3f/foAR3/sdWtFnbLeKiP8R7m0J0rPCtwBfaxXTi/1A9fxofl8aJ05
T0vcS3D0ULebF301+lWcHfqiWW+V3hd14ffgU/eE3iQ9rj9r/7F69u017b4W
T0+NvgPYP8r/znXsyVX1TbFD/5n7nMk4U69aU+pN/xEbeof0y/5l7elh0tOk
lvQq306q+pufa0ePjrzBIfu7WPpn/7SW9Cp/sM+l9SDOi9q8aY2mGBss8b9y
38Ai4/+F4be4tfsMnvpyT6cnUdd6/NeakEdiTaGc+ppJVR+COxTfzyOs0/fW
iryruR7u//StahizunmSX0FeLQ6xFvSO6G3Rg8v6Tcf7byT9BPoR9Gy449MX
oHfTUr5QPvOBHWesqVh80SdqpQ19nTbyYLuKz/qGrAcZ91bup1mvp4014H5L
36K266Ofwb9tZyRVvQ36HfQ3uOtzv+TOmfUy4LmrD5JvLo486X+01/5kbZqb
dx9l5DHA/LmDc6/Pek8D1XcwRjtx5+iD2FeaD3aDzSO7317hGrnDcR8f4nyw
+Z1mzbqo66ysp/WCv8i8fvh/c9z/np6bbsbhO4/v3uGeI+6Y3PO4/w8VN9z6
Z3d7zhz3cr7DuEfxncp3MN/GfINO0DffZNm321jXw3cbdzTuddz/xhmHdfNt
RJ7TJO6hfDfxbT7LPLkncFe5xzxmmDc00/z5nuY7M/veLHWdo42NHXcL7hjc
ucrF3WKe81zXSPV8r3Jvm639PLHIsrsfsuzuV25+3LXu0tdc45P74jzdHfrI
7nXw/Dv2jb/vYJeZK3/b3/nbm91/Vlh37jmrlXG34V7DXYFvfO4I94ulftxD
HlDGnWON/CZ59Ev1QezN+sqwG4yxwdqu1Xd2j9ymzSb3m296vtW5E3A3mKSM
+0J2X9meVN0ltmj/RFJ1rwLLdz33hZ3q2FPO1B7jbFH/qLKnjcO3+46k6g4z
2TzQP6ff/wHhHjsQ
"]]],
Lighting->Automatic]}, {}, {}, {}, {}}]},
DisplayFunction->Identity,
FaceGridsStyle->Automatic,
Method->{"DefaultBoundaryStyle" -> Directive[
GrayLevel[0.3]]},
PlotRange->{All, All, All},
PlotRangePadding->{Automatic, Automatic, Automatic},
Ticks->{Automatic, Automatic, Automatic},
ViewPoint->{1.3, -2.4, 2.},
ViewVertical->{0., 0., 1.}]], "Output",
CellChangeTimes->{{3.759070955384362*^9, 3.759070977652513*^9}, {
3.759071021266911*^9, 3.75907105987184*^9}, 3.7590710948405046`*^9,
3.759071884142231*^9, 3.7590727720663548`*^9, 3.759192553280632*^9},
CellLabel->"Out[9]=",ExpressionUUID->"080d0719-e0d3-4b32-ac66-96ffe52139ab"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"Mass", "[",
RowBox[{"boatregion_", ",", " ", "boatdensity_"}], "]"}], ":=",
RowBox[{"Integrate", "[",
RowBox[{"boatdensity", ",",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", " ", "y", ",", " ", "z"}], "}"}], "\[Element]",
"boatregion"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"boatmass", " ", "=", " ",
RowBox[{"Mass", "[",
RowBox[{"hullreg", ",", " ", "boatdens"}],
"]"}]}], "\[IndentingNewLine]"}], "Input",
CellChangeTimes->{{3.75907084715088*^9, 3.7590709012747087`*^9},
3.759071138603026*^9, {3.759071259499156*^9, 3.759071266928605*^9}, {
3.759071300947072*^9, 3.759071307409132*^9}, 3.759071376767318*^9},
CellLabel->"In[10]:=",ExpressionUUID->"148f209f-1c4c-4a67-bd1b-e73e94ce68a7"],
Cell[BoxData["22999.99934311148`"], "Output",
CellChangeTimes->{3.75907134693904*^9, 3.7590713798961782`*^9,
3.759071886737342*^9, 3.759072775328446*^9, 3.7591925571537666`*^9},
CellLabel->"Out[11]=",ExpressionUUID->"7785f581-47bb-474f-aebe-b937a02b0b5a"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"CenterOfMass", "[",
RowBox[{"boatregion_", ",", " ", "boatdensity_"}], "]"}], ":=", " ",
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{"boatdensity", " ", "*", " ",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], ",", " ",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
"boatregion"}]}], "]"}], "/", "boatmass"}]}], "\[IndentingNewLine]",
RowBox[{"boatcom", " ", "=", " ",
RowBox[{"CenterOfMass", "[",
RowBox[{"hullreg", ",", " ", "boatdens"}], "]"}]}]}], "Input",
CellChangeTimes->{{3.759071318573318*^9, 3.759071330976748*^9},
3.759071392873054*^9},
CellLabel->"In[12]:=",ExpressionUUID->"bbe2e87c-4594-4617-88eb-689f8882b028"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-", "3.02406055700285`*^-17"}], ",", "5.0000000000000036`", ",",
"0.9139751432480252`"}], "}"}]], "Output",
CellChangeTimes->{
3.759071228139221*^9, 3.759071278282048*^9, {3.759071338256373*^9,
3.759071395475768*^9}, 3.759071889439376*^9, 3.759072778598921*^9,
3.759192558719369*^9},
CellLabel->"Out[13]=",ExpressionUUID->"2b49b75c-d297-4833-8c6a-7446296b6b6c"]
}, Open ]],
Cell[BoxData[{
RowBox[{
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]_", ",", " ", "\[Phi]_"}], "]"}], ":=", " ",
RowBox[{"Quiet", "[",
RowBox[{"ImplicitRegion", "[",
RowBox[{
RowBox[{
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"y", ",", " ",
RowBox[{"\[Phi]", "\[Equal]", "90"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"x", ",", " ",
RowBox[{"\[Theta]", "==", "90"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "z"}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}]}], "}"}], ",", "z"}],
"]"}], "\[LessEqual]",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"draft", ",", " ",
RowBox[{
RowBox[{"\[Theta]", "\[Equal]", "90"}], "||",
RowBox[{"\[Phi]", "\[Equal]", "90"}]}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-", "1"}], "*",
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "x"}], "+",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "y"}], "+",
"draft"}], ")"}]}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}]}], "}"}], ",",
RowBox[{
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "x"}], "+",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "*", "y"}], "+",
"draft"}]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",",
RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"y", ",",
RowBox[{"-", "1"}], ",", "11"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"z", ",",
RowBox[{"-", "1"}], ",", "3"}], "}"}]}], "}"}]}], "]"}],
"]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CalculateDraft3D", "[",
RowBox[{
"boatregion_", ",", " ", "boatdensity_", ",", " ", "\[Theta]_", ",", " ",
"\[Phi]_"}], "]"}], ":=",
RowBox[{"Quiet", "[",
RowBox[{
RowBox[{"FindRoot", "[",
RowBox[{
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{"1.", ",", " ",
RowBox[{
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}], "\[Element]",
RowBox[{"RegionIntersection", "[",
RowBox[{"boatregion", ",", " ",
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]", ",", " ", "\[Phi]"}], "]"}]}], "]"}]}]}],
"]"}], "\[Equal]",
RowBox[{"boatmass", "/", "1000."}]}], ",", " ",
RowBox[{"{",
RowBox[{"draft", ",", " ",
RowBox[{"1", "-",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "/", "3"}], "-",
RowBox[{
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}], "/", "3"}]}]}], "}"}]}],
"]"}], "[",
RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"DisplacedWater3D", "[",
RowBox[{
"boatregion_", ",", " ", "boatdensity_", ",", " ", "\[Theta]_", ",", " ",
"\[Phi]_"}], "]"}], ":=", " ",
RowBox[{"Quiet", "[",
RowBox[{
RowBox[{"RegionIntersection", "[",
RowBox[{"boatregion", ",", " ",
RowBox[{"Waterline3D", "[",
RowBox[{"\[Theta]", ",", " ", "\[Phi]"}], "]"}]}], "]"}], "/.",
RowBox[{"CalculateDraft3D", "[",
RowBox[{
"boatregion", ",", " ", "boatdensity", ",", " ", "\[Theta]", ",", " ",
"\[Phi]"}], "]"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"CenterOfBuoyancy3D", "[",
RowBox[{
"boatregion_", ",", " ", "boatdensity_", ",", " ", "\[Theta]_", ",", " ",
"\[Phi]_"}], "]"}], ":=",
RowBox[{"Quiet", "[",
RowBox[{"RegionCentroid", "[",
RowBox[{"DisplacedWater3D", "[",
RowBox[{
"boatregion", ",", " ", "boatdensity", ",", " ", "\[Theta]", ",", " ",
"\[Phi]"}], "]"}], "]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"RightingMoment3D", "[",
RowBox[{
"boatregion_", ",", " ", "boatdensity_", ",", " ", "\[Theta]_", ",", " ",
"\[Phi]_"}], "]"}], ":=",
RowBox[{"Quiet", "[",
RowBox[{"Cross", "[",
RowBox[{
RowBox[{
RowBox[{"CenterOfBuoyancy3D", "[",
RowBox[{
"boatregion", ",", " ", "boatdensity", ",", " ", "\[Theta]", ",", " ",
"\[Phi]"}], "]"}], "-", "boatcom"}], ",", " ",
RowBox[{"boatmass", "*", "9.8", "*",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", " ",
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", "1"}], "}"}], "/",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", " ",
RowBox[{
RowBox[{"\[Theta]", ">", "90"}], "||",
RowBox[{"\[Phi]", ">", "90"}]}]}], "}"}], "}"}], ",", "1"}],
"]"}]}], "/",
RowBox[{"Norm", "[",
RowBox[{"{",
RowBox[{
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Theta]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", " ",
RowBox[{"-",
RowBox[{"Tan", "[",
RowBox[{"\[Phi]", "*",
RowBox[{"\[Pi]", "/", "180"}]}], "]"}]}], ",", "1"}], "}"}],
"]"}]}]}]}], "]"}], "]"}]}]}], "Input",
CellChangeTimes->{{3.75907150036087*^9, 3.75907151420553*^9}, {
3.759071689998314*^9, 3.759071692052773*^9}, {3.759071729200437*^9,
3.7590717906234617`*^9}, {3.7590718315538177`*^9, 3.75907184764012*^9}},
CellLabel->"In[14]:=",ExpressionUUID->"0440646f-917b-4090-ad63-896091dbb2a8"],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"RightingMoment3D", "[",
RowBox[{"hullreg", ",", " ", "boatdens", ",", "89", ",", " ", "0"}], "]"}],
"//", "AbsoluteTiming"}]], "Input",
CellChangeTimes->{{3.759071905531032*^9, 3.759071907593807*^9}},
CellLabel->"In[19]:=",ExpressionUUID->"ea65a17a-7f9c-4f1d-9057-0e171d1fee1f"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0.5075956513202365`", ",",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1.0481674990928554`*^-11"}], ",",
RowBox[{"-", "64060.320286829585`"}], ",",
RowBox[{"-", "6.004947580563875`*^-10"}]}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{
3.7590714312247763`*^9, 3.759071802545567*^9, 3.7590718608981047`*^9, {
3.7590718944794893`*^9, 3.7590719083857803`*^9}, 3.7591925593308125`*^9},
CellLabel->"Out[19]=",ExpressionUUID->"c35b9f1e-5c72-49ee-9c29-19c6fadc1625"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"RightingMoment3D", "[",
RowBox[{"hullreg", ",", " ", "boatdens", ",", "89", ",", "10"}], "]"}], "//",
"AbsoluteTiming"}]], "Input",
CellChangeTimes->{{3.759071449664168*^9, 3.7590714927464*^9}, {
3.759071545919116*^9, 3.759071584836463*^9}, {3.759071619381941*^9,
3.7590716496829844`*^9}, {3.7590719297996807`*^9, 3.7590719330211697`*^9}},
CellLabel->"In[20]:=",ExpressionUUID->"78c14fa7-2771-443b-9bc2-26805db7c946"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"3.0628179504916737`", ",",
RowBox[{"{",
RowBox[{"272.503137033268`", ",",
RowBox[{"-", "64059.56161180969`"}], ",", "4316.265180379211`"}],
"}"}]}], "}"}]], "Output",
CellChangeTimes->{{3.759071465261001*^9, 3.759071472382268*^9}, {
3.7590715406538877`*^9, 3.7590715878343143`*^9}, {3.759071635507448*^9,
3.7590716506601763`*^9}, {3.759071918704681*^9, 3.75907193461478*^9},
3.759192562400545*^9},
CellLabel->"Out[20]=",ExpressionUUID->"2151b666-6472-4bc9-848a-61f460a5bdfb"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"f", "[",
RowBox[{"angle1_", ",", " ", "angle2_"}], "]"}], ":=",
RowBox[{"N", "[",
RowBox[{
RowBox[{"Norm", "[",
RowBox[{"RightingMoment3D", "[",
RowBox[{
"hullreg", ",", " ", "boatdens", ",", "angle1", ",", " ", "angle2"}],
"]"}], "]"}], "*",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"-", "1"}], ",", " ",
RowBox[{
RowBox[{
RowBox[{"RightingMoment3D", "[",
RowBox[{
"hullreg", ",", " ", "boatdens", ",", "angle1", ",", " ",
"angle2"}], "]"}], "[",
RowBox[{"[", "2", "]"}], "]"}], ">", "0."}]}], "}"}], "}"}], ",",
"1"}], "]"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable1", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "0.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "0."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable2", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "2.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "2."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable3", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "4.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "4."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable4", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "6.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "6."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable5", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "8.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "8."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable6", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "10.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "10."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable7", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "12.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "12."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable8", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "14.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "14."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable9", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "16.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "16."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable10", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "18.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "18."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",
RowBox[{"angleTable11", "=",
RowBox[{"Quiet", "[",
RowBox[{"Parallelize", "[",
RowBox[{"Table", "[",
RowBox[{
RowBox[{"{",
RowBox[{"angle1", ",", " ", "20.", ",",
RowBox[{"f", "[",
RowBox[{"angle1", ",", "20."}], "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{"angle1", ",", "1.", ",", "176.", ",", "5."}], "}"}]}], "]"}],
"]"}], "]"}]}], "\[IndentingNewLine]",