-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathllio_stat_analysis_defense_split2.nb
8522 lines (8409 loc) · 435 KB
/
llio_stat_analysis_defense_split2.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 9.0' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 157, 7]
NotebookDataLength[ 445026, 8513]
NotebookOptionsPosition[ 441581, 8409]
NotebookOutlinePosition[ 441924, 8424]
CellTagsIndexPosition[ 441881, 8421]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell["DPS equation", "Section",
CellChangeTimes->{{3.6008619935060844`*^9, 3.60086199885638*^9}}],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"mit", "[", "d_", "]"}], " ", ":=", " ",
RowBox[{"1", "/",
RowBox[{"(",
RowBox[{"100", "/",
RowBox[{"(",
RowBox[{"100", "+",
RowBox[{"(",
RowBox[{"d", "/", "20"}], ")"}]}], ")"}]}], ")"}]}]}], "\n",
RowBox[{"ehp", " ", "=", " ",
RowBox[{"FullSimplify", "[",
RowBox[{
RowBox[{"(",
RowBox[{"hp", "/",
RowBox[{"(",
RowBox[{"66", "/", "25"}], ")"}]}], ")"}], " ",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"mit", "[", "ar", "]"}], "*", "scale"}], ")"}], " ", "+", " ",
RowBox[{"(",
RowBox[{
RowBox[{"mit", "[", "mr", "]"}], "*",
RowBox[{"(",
RowBox[{"1", "-", "scale"}], ")"}]}], ")"}]}], ")"}]}],
"]"}]}], "\n",
RowBox[{
RowBox[{"s", " ", "=", " ", "0.3"}], ";"}], "\n",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"mit", "[", "ar", "]"}], "*", "scale"}], ")"}], " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"hp", " ", "\[Rule]", " ",
RowBox[{"2000", "*",
RowBox[{"(",
RowBox[{"66", "/", "25"}], ")"}]}]}], ",", " ",
RowBox[{"ar", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"mr", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"scale", " ", "\[Rule]", " ", "s"}]}], "}"}]}], "\n",
RowBox[{
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"mit", "[", "mr", "]"}], "*",
RowBox[{"(",
RowBox[{"1", "-", "scale"}], ")"}]}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}]}], "}"}], "}"}], ",", " ", "0"}],
"]"}], " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"hp", " ", "\[Rule]", " ",
RowBox[{"2000", "*",
RowBox[{"(",
RowBox[{"66", "/", "25"}], ")"}]}]}], ",", " ",
RowBox[{"ar", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"mr", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"scale", " ", "\[Rule]", " ", "s"}]}], "}"}]}], "\n",
RowBox[{"ehp", " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"hp", " ", "\[Rule]", " ",
RowBox[{"2000", "*",
RowBox[{"(",
RowBox[{"66", "/", "25"}], ")"}]}]}], ",", " ",
RowBox[{"ar", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"mr", " ", "\[Rule]", " ", "2000"}], ",", " ",
RowBox[{"scale", " ", "\[Rule]", " ", "0.5"}]}], "}"}]}], "\n",
RowBox[{"Clear", "[",
RowBox[{"hp", ",", "ar", ",", "mr", ",", "scale"}], "]"}], "\n",
RowBox[{
RowBox[{"g", " ", "=", " ",
RowBox[{"Grad", "[",
RowBox[{"ehp", ",", " ",
RowBox[{"{",
RowBox[{"hp", ",", "ar", ",", "mr"}], "}"}]}], "]"}]}], "\n"}], "\n",
RowBox[{"Manipulate", "[",
RowBox[{
RowBox[{"VectorPlot3D", "[",
RowBox[{
RowBox[{"g", " ", "/.", " ",
RowBox[{"scale", " ", "\[Rule]", " ", "ss"}]}], ",", " ",
RowBox[{"{",
RowBox[{"hp", ",", " ", "0", ",", " ", "10000"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"ar", ",", " ", "0", ",", " ", "10000"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"mr", ",", " ", "0", ",", " ", "10000"}], "}"}], ",", " ",
RowBox[{"VectorPoints", " ", "->", " ", "Coarse"}], ",", " ",
RowBox[{"AxesLabel", "\[Rule]",
RowBox[{"{",
RowBox[{"hp", ",", "ar", ",", "mr"}], "}"}]}]}], "]"}], ",", " ",
RowBox[{"{",
RowBox[{"ss", ",", " ", "0", ",", " ", "1"}], "}"}]}], "]"}]}], "Code",
CellChangeTimes->{{3.6008582880445685`*^9, 3.600858317031201*^9}, {
3.6008583606546545`*^9, 3.600858405198162*^9}, {3.6008585196806087`*^9,
3.6008585279720755`*^9}, {3.6008587146705756`*^9, 3.600858726157227*^9}, {
3.6008587775481143`*^9, 3.6008588250217857`*^9}, {3.6008588760766635`*^9,
3.600858914301814*^9}, {3.600858949211779*^9, 3.6008590917868013`*^9}, {
3.6008591925964737`*^9, 3.6008592011679573`*^9}, {3.6008592866647677`*^9,
3.6008592931361265`*^9}, 3.60085955540089*^9, {3.600859660742817*^9,
3.6008596651970673`*^9}, {3.600859705605336*^9, 3.600859754452085*^9},
3.600859825116067*^9, {3.6008599745644765`*^9, 3.600860019812018*^9}, {
3.6008601696604548`*^9, 3.6008601942908354`*^9}, {3.6008611483235254`*^9,
3.600861184066537*^9}, 3.6008619190028877`*^9, 3.6008619928730445`*^9,
3.600862700412858*^9, {3.601051006153372*^9, 3.601051016539962*^9}, {
3.601051140944957*^9, 3.6010511647382956`*^9}, {3.6010514263150144`*^9,
3.6010514684563923`*^9}, {3.6010516368658676`*^9,
3.6010516405210733`*^9}, {3.601053253155814*^9, 3.6010532801833344`*^9}, {
3.601055114159534*^9, 3.60105515308673*^9}, {3.601055227406912*^9,
3.601055269782296*^9}, {3.601055904006985*^9, 3.6010559065831304`*^9}, {
3.6010561467876415`*^9, 3.601056192587224*^9}, {3.601315280565883*^9,
3.601315350984846*^9}, {3.6013156399451065`*^9, 3.6013156413771877`*^9}, {
3.601315845048648*^9, 3.6013158778644943`*^9}, {3.601315918193764*^9,
3.6013159432731805`*^9}, {3.6014124835197706`*^9,
3.6014124880110254`*^9}, {3.601412519451806*^9, 3.6014125256191564`*^9}, {
3.6014126851221924`*^9, 3.601412695338772*^9}, {3.6014129210825605`*^9,
3.601412923970725*^9}, {3.601413099091646*^9, 3.601413126185181*^9}, {
3.6014134244180765`*^9, 3.6014135609068108`*^9}, {3.601413869515295*^9,
3.601413949970848*^9}, {3.601413980691593*^9, 3.6014139858508797`*^9}, {
3.601415063779956*^9, 3.601415092515579*^9}, {3.6014151448505487`*^9,
3.6014151977625465`*^9}, 3.601415272113759*^9, {3.601418900912349*^9,
3.6014189036635003`*^9}, {3.601418948561048*^9, 3.60141895070317*^9}, {
3.6014295386580267`*^9, 3.6014296848003116`*^9}, {3.6014298144596586`*^9,
3.601430130370556*^9}, {3.6014301805153923`*^9, 3.6014302679473505`*^9}, {
3.601430330091871*^9, 3.6014304138266096`*^9}, {3.601430459084174*^9,
3.601430555771657*^9}, {3.6014307026579795`*^9, 3.6014307079702797`*^9}, {
3.6014307430552673`*^9, 3.6014307491736145`*^9}, {3.601430974742394*^9,
3.601430990076263*^9}, {3.601503279048199*^9, 3.6015032892157717`*^9}, {
3.601504285114812*^9, 3.6015043720087013`*^9}, {3.601504445884859*^9,
3.6015044688201485`*^9}, {3.6015045031740828`*^9,
3.6015045888389025`*^9}, {3.6015046252539525`*^9, 3.601504676164817*^9}, {
3.6015047216533766`*^9, 3.60150474944394*^9}, {3.601504789927219*^9,
3.601504794381469*^9}, {3.6015048360938168`*^9, 3.601504836326829*^9}}],
Cell[BoxData[
FractionBox[
RowBox[{"hp", " ",
RowBox[{"(",
RowBox[{"2000", "+", "mr", "+",
RowBox[{"ar", " ", "scale"}], "-",
RowBox[{"mr", " ", "scale"}]}], ")"}]}], "5280"]], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.6015048374178915`*^9}],
Cell[BoxData["0.6`"], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.6015048374208913`*^9}],
Cell[BoxData["1.4`"], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.601504837422891*^9}],
Cell[BoxData["4000.`"], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.6015048374248905`*^9}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
FractionBox[
RowBox[{"2000", "+", "mr", "+",
RowBox[{"ar", " ", "scale"}], "-",
RowBox[{"mr", " ", "scale"}]}], "5280"], ",",
FractionBox[
RowBox[{"hp", " ", "scale"}], "5280"], ",",
FractionBox[
RowBox[{"hp", " ",
RowBox[{"(",
RowBox[{"1", "-", "scale"}], ")"}]}], "5280"]}], "}"}]], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.601504837425891*^9}],
Cell[BoxData[
TagBox[
StyleBox[
DynamicModuleBox[{$CellContext`ss$$ = 0., Typeset`show$$ = True,
Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu",
Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ =
"\"untitled\"", Typeset`specs$$ = {{
Hold[$CellContext`ss$$], 0, 1}}, Typeset`size$$ = {297., {166., 170.}},
Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ =
True, $CellContext`ss$16786$$ = 0},
DynamicBox[Manipulate`ManipulateBoxes[
1, StandardForm, "Variables" :> {$CellContext`ss$$ = 0},
"ControllerVariables" :> {
Hold[$CellContext`ss$$, $CellContext`ss$16786$$, 0]},
"OtherVariables" :> {
Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$,
Typeset`animator$$, Typeset`animvar$$, Typeset`name$$,
Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
Typeset`skipInitDone$$}, "Body" :> VectorPlot3D[
ReplaceAll[$CellContext`g, $CellContext`scale -> $CellContext`ss$$], \
{$CellContext`hp, 0, 10000}, {$CellContext`ar, 0, 10000}, {$CellContext`mr, 0,
10000}, VectorPoints -> Coarse,
AxesLabel -> {$CellContext`hp, $CellContext`ar, $CellContext`mr}],
"Specifications" :> {{$CellContext`ss$$, 0, 1}}, "Options" :> {},
"DefaultOptions" :> {}],
ImageSizeCache->{348., {212., 217.}},
SingleEvaluation->True],
Deinitialization:>None,
DynamicModuleValues:>{},
SynchronousInitialization->True,
UnsavedVariables:>{Typeset`initDone$$},
UntrackedVariables:>{Typeset`size$$}], "Manipulate",
Deployed->True,
StripOnInput->False],
Manipulate`InterpretManipulate[1]]], "Output",
CellChangeTimes->{
3.6014301322306566`*^9, {3.6014301837695813`*^9, 3.6014302683313675`*^9}, {
3.6014303418475323`*^9, 3.6014304142376337`*^9}, {3.6014304531908455`*^9,
3.6014305753657618`*^9}, 3.6014307083652973`*^9, 3.601430749810651*^9,
3.6014308805640593`*^9, 3.601430991221322*^9, 3.6014450156978865`*^9,
3.601503292487956*^9, {3.601504314901493*^9, 3.6015043724987354`*^9}, {
3.6015044512291594`*^9, 3.6015044691871705`*^9}, {3.601504566610652*^9,
3.6015045890889177`*^9}, {3.6015046267410355`*^9, 3.6015046765788403`*^9},
3.601504749809961*^9, 3.6015047947154884`*^9, 3.6015048376219034`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{"constraints", " ", "=", " ",
RowBox[{"{",
RowBox[{
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"hp", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"ar", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"mr", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{
RowBox[{"hp", " ", "+", " ", "ar", " ", "+", " ", "mr"}], " ", "==",
" ", "3000"}]}], "}"}]}], ";"}], "\n",
RowBox[{
RowBox[{"{",
RowBox[{"knownBestEHP", ",", "knownBestDistro"}], "}"}], " ", "=", " ",
RowBox[{"Assuming", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}]}], "}"}], ",", " ",
RowBox[{"FullSimplify", "[",
RowBox[{"Maximize", "[",
RowBox[{
RowBox[{"{",
RowBox[{"ehp", ",", "constraints"}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"hp", ",", "ar", ",", "mr"}], "}"}]}], "]"}], "]"}]}],
"]"}]}]}], "Code",
CellChangeTimes->{{3.6008582880445685`*^9, 3.600858317031201*^9}, {
3.6008583606546545`*^9, 3.600858405198162*^9}, {3.6008585196806087`*^9,
3.6008585279720755`*^9}, {3.6008587146705756`*^9, 3.600858726157227*^9}, {
3.6008587775481143`*^9, 3.6008588250217857`*^9}, {3.6008588760766635`*^9,
3.600858914301814*^9}, {3.600858949211779*^9, 3.6008590917868013`*^9}, {
3.6008591925964737`*^9, 3.6008592011679573`*^9}, {3.6008592866647677`*^9,
3.6008592931361265`*^9}, 3.60085955540089*^9, {3.600859660742817*^9,
3.6008596651970673`*^9}, {3.600859705605336*^9, 3.600859754452085*^9},
3.600859825116067*^9, {3.6008599745644765`*^9, 3.600860019812018*^9}, {
3.6008601696604548`*^9, 3.6008601942908354`*^9}, {3.6008611483235254`*^9,
3.600861184066537*^9}, {3.6008619190028877`*^9, 3.6008619311325693`*^9}, {
3.6008627409331427`*^9, 3.6008627412881575`*^9}, {3.6008668001335545`*^9,
3.600866805786877*^9}, {3.6013153554981003`*^9, 3.6013153764172773`*^9}, {
3.6013156572420797`*^9, 3.601315662449378*^9}, {3.601412956090539*^9,
3.6014129631379433`*^9}, {3.601413034578991*^9, 3.601413120714865*^9}, {
3.601414034873663*^9, 3.6014140715857377`*^9}, {3.6014152918018694`*^9,
3.6014153070177364`*^9}, {3.601415345681926*^9, 3.601415354273408*^9}, {
3.6014153956027555`*^9, 3.601415405586321*^9}, {3.601418929874985*^9,
3.6014189768326454`*^9}, {3.601424353445258*^9, 3.6014243535332627`*^9}, {
3.6014310066102*^9, 3.6014310143616333`*^9}, {3.6014911412131352`*^9,
3.601491168068657*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
FractionBox["12500", "11"],
RowBox[{
FractionBox["1", "3"], "\[LessEqual]", "scale", "\[LessEqual]",
FractionBox["2", "3"]}]},
{
RowBox[{"-",
FractionBox[
RowBox[{"3125", " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"-", "5"}], "+",
RowBox[{"3", " ", "scale"}]}], ")"}], "2"]}],
RowBox[{"66", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+", "scale"}], ")"}]}]]}],
RowBox[{"scale", "<",
FractionBox["1", "3"]}]},
{
FractionBox[
RowBox[{"3125", " ",
SuperscriptBox[
RowBox[{"(",
RowBox[{"2", "+",
RowBox[{"3", " ", "scale"}]}], ")"}], "2"]}],
RowBox[{"66", " ", "scale"}]],
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], ",",
RowBox[{"{",
RowBox[{
RowBox[{"hp", "\[Rule]",
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{"3000",
RowBox[{
FractionBox["1", "3"], "\[LessEqual]", "scale", "\[LessEqual]",
FractionBox["2", "3"]}]},
{
FractionBox[
RowBox[{"500", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "5"}], "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}],
RowBox[{
RowBox[{"-", "1"}], "+", "scale"}]],
RowBox[{
RowBox[{"3", " ", "scale"}], "<", "1"}]},
{
FractionBox[
RowBox[{"500", " ",
RowBox[{"(",
RowBox[{"2", "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}], "scale"],
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]}], ",",
RowBox[{"ar", "\[Rule]",
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
FractionBox[
RowBox[{"500", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "2"}], "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}], "scale"],
RowBox[{
RowBox[{"3", " ", "scale"}], ">", "2"}]},
{"0",
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]}], ",",
RowBox[{"mr", "\[Rule]",
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
FractionBox[
RowBox[{"500", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}],
RowBox[{
RowBox[{"-", "1"}], "+", "scale"}]],
RowBox[{
RowBox[{"3", " ", "scale"}], "<", "1"}]},
{"0",
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.6014131119113674`*^9, 3.601413121470914*^9},
3.601413434782665*^9, 3.601413593882674*^9, 3.60141404197406*^9,
3.6014140720777655`*^9, {3.601415285384511*^9, 3.6014153137261105`*^9},
3.601415361918847*^9, 3.601415415835902*^9, {3.6014189246876907`*^9,
3.6014189802138367`*^9}, 3.601424356980459*^9, 3.6014301651005187`*^9,
3.6014305795579996`*^9, {3.6014307777042255`*^9, 3.601430782713515*^9}, {
3.6014308792189827`*^9, 3.6014308851693163`*^9}, {3.6014309969136505`*^9,
3.6014310154406943`*^9}, 3.6014414131297827`*^9, 3.601441448362773*^9,
3.601445020687169*^9, {3.6014911455713773`*^9, 3.6014911695077333`*^9},
3.601505027404582*^9}]
}, Open ]]
}, Open ]],
Cell[CellGroupData[{
Cell["Path of steepest ascent ", "Section",
CellChangeTimes->{{3.600861976437124*^9, 3.600861982000432*^9},
3.600862015768332*^9}],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"g", " ", "=", " ",
RowBox[{"Grad", "[",
RowBox[{
RowBox[{"ehp", " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"hp", " ", "\[Rule]", " ", "x"}], ",", " ",
RowBox[{"ar", " ", "\[Rule]", " ", "y"}], ",", " ",
RowBox[{"mr", " ", "\[Rule]", " ", "z"}]}], "}"}]}], ",", " ",
RowBox[{"{",
RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], "]"}]}]], "Code",
CellChangeTimes->{{3.6008582880445685`*^9, 3.600858317031201*^9}, {
3.6008583606546545`*^9, 3.600858405198162*^9}, {3.6008585196806087`*^9,
3.6008585279720755`*^9}, {3.6008587146705756`*^9, 3.600858726157227*^9}, {
3.6008587775481143`*^9, 3.6008588250217857`*^9}, {3.6008588760766635`*^9,
3.600858914301814*^9}, {3.600858949211779*^9, 3.6008590917868013`*^9}, {
3.6008591925964737`*^9, 3.6008592480765963`*^9}, 3.6008592981444135`*^9, {
3.600859361409969*^9, 3.600859390363598*^9}, {3.600859468987023*^9,
3.6008595268242817`*^9}, {3.600859843032069*^9, 3.6008598535276604`*^9}, {
3.600860042906317*^9, 3.6008600814234843`*^9}, {3.6008603562669497`*^9,
3.6008605080604973`*^9}, 3.6008605564452143`*^9, {3.600860664480299*^9,
3.6008607054555993`*^9}, {3.6008607416336403`*^9,
3.6008607849750795`*^9}, {3.6008608287525377`*^9, 3.600860853552938*^9}, {
3.600860894119221*^9, 3.600860904605811*^9}, {3.6008610451407194`*^9,
3.600861046743804*^9}, {3.6008614013037605`*^9, 3.6008614307954197`*^9}, {
3.600861479250142*^9, 3.6008615097028556`*^9}, {3.6008616865738096`*^9,
3.6008617041587987`*^9}, 3.6008619403410883`*^9, 3.6008619756150737`*^9,
3.6008620189805126`*^9, {3.60105620620499*^9, 3.6010562172286105`*^9}, {
3.601315386441842*^9, 3.60131540435185*^9}, {3.60131597090473*^9,
3.601315970992735*^9}, {3.6014155462242837`*^9, 3.601415550200514*^9}}],
Cell[BoxData[
RowBox[{"{",
RowBox[{
FractionBox[
RowBox[{"2000", "+",
RowBox[{"scale", " ", "y"}], "+", "z", "-",
RowBox[{"scale", " ", "z"}]}], "5280"], ",",
FractionBox[
RowBox[{"scale", " ", "x"}], "5280"], ",",
FractionBox[
RowBox[{
RowBox[{"(",
RowBox[{"1", "-", "scale"}], ")"}], " ", "x"}], "5280"]}],
"}"}]], "Output",
CellChangeTimes->{3.6008619420441895`*^9, 3.6008627438162994`*^9,
3.6010510461936255`*^9, 3.601051188374627*^9, 3.6010514397577715`*^9,
3.601051476479843*^9, 3.601051652776764*^9, 3.601053291789987*^9,
3.6010552815339575`*^9, 3.6010559133165092`*^9, 3.601056217672634*^9,
3.601216735674362*^9, 3.6013154802581205`*^9, 3.6013156738950167`*^9,
3.601315971349761*^9, 3.601413137697833*^9, 3.6014155086471543`*^9,
3.6014155504415226`*^9, 3.6014190153588285`*^9, 3.6014306134449196`*^9,
3.6014307920890408`*^9, 3.6014308921307135`*^9, 3.6014310248602343`*^9,
3.6014414633936305`*^9, 3.6014450236053348`*^9, 3.6015050296727085`*^9}]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"hp2arUnbounded1", ",", "hp2mrUnbounded1"}], "}"}], "}"}], " ",
"=", " ",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"y", "[", "x", "]"}], ",",
RowBox[{"z", "[", "x", "]"}]}], "}"}], " ", "/.", " ",
RowBox[{"DSolve", "[",
RowBox[{
RowBox[{"{", "\n", "\t",
RowBox[{
RowBox[{
RowBox[{
RowBox[{"y", "'"}], "[", "x", "]"}], " ", "\[Equal]", " ",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"g", "[",
RowBox[{"[", "1", "]"}], "]"}], " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"y", " ", "\[Rule]", " ",
RowBox[{"y", "[", "x", "]"}]}], ",", " ",
RowBox[{"z", " ", "\[Rule]", " ",
RowBox[{"z", "[", "x", "]"}]}]}], "}"}]}], ")"}], " ", "/",
" ",
RowBox[{"(",
RowBox[{"g", "[",
RowBox[{"[", "2", "]"}], "]"}], ")"}]}]}], ",", " ", "\n", "\t",
RowBox[{
RowBox[{
RowBox[{"z", "'"}], "[", "x", "]"}], " ", "\[Equal]", " ",
RowBox[{
RowBox[{"(",
RowBox[{
RowBox[{"g", "[",
RowBox[{"[", "1", "]"}], "]"}], " ", "/.", " ",
RowBox[{"{",
RowBox[{
RowBox[{"y", " ", "\[Rule]", " ",
RowBox[{"y", "[", "x", "]"}]}], ",", " ",
RowBox[{"z", " ", "\[Rule]", " ",
RowBox[{"z", "[", "x", "]"}]}]}], "}"}]}], ")"}], " ", "/",
" ",
RowBox[{"(",
RowBox[{"g", "[",
RowBox[{"[", "3", "]"}], "]"}], ")"}]}]}], ",", " ", "\n", "\t",
RowBox[{
RowBox[{"y", "[",
RowBox[{"hp", " ", "/.", " ", "knownBestDistro"}], "]"}], " ",
"\[Equal]", " ",
RowBox[{"(",
RowBox[{"ar", " ", "/.", " ", "knownBestDistro"}], ")"}]}], ",",
" ", "\n", "\t",
RowBox[{
RowBox[{"z", "[",
RowBox[{"hp", " ", "/.", " ", "knownBestDistro"}], "]"}], " ",
"\[Equal]", " ",
RowBox[{"(",
RowBox[{"mr", " ", "/.", " ", "knownBestDistro"}], ")"}]}]}],
"}"}], ",", " ", "\n", "\t",
RowBox[{"{",
RowBox[{
RowBox[{"y", "[", "x", "]"}], ",",
RowBox[{"z", "[", "x", "]"}]}], "}"}], ",", " ", "x"}], "]"}]}]}],
";"}], "\n"}], "\n",
RowBox[{"hp2arUnbounded", " ", "=", " ",
RowBox[{"Assuming", "[",
RowBox[{
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"FullSimplify", "[", "hp2arUnbounded1", "]"}]}], "]"}]}], "\n",
RowBox[{
RowBox[{"hp2mrUnbounded", " ", "=", " ",
RowBox[{"Assuming", "[",
RowBox[{
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"FullSimplify", "[", "hp2mrUnbounded1", "]"}]}], "]"}]}],
"\n"}], "\n",
RowBox[{"hp2arCondition", " ", "=", " ",
RowBox[{"Assuming", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", " ", ">=", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}]}], "}"}], ",", " ",
RowBox[{"FullSimplify", "@",
RowBox[{"Reduce", "[",
RowBox[{"{",
RowBox[{
RowBox[{"hp2arUnbounded", " ", ">=", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"x", " ", ">=", " ", "0"}]}], "}"}], "]"}]}]}], "]"}]}], "\n",
RowBox[{
RowBox[{"hp2mrCondition", " ", "=", " ",
RowBox[{"Assuming", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"x", " ", ">=", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}]}], "}"}], ",", " ",
RowBox[{"FullSimplify", "@",
RowBox[{"Reduce", "[",
RowBox[{"{",
RowBox[{
RowBox[{"hp2mrUnbounded", " ", ">=", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"x", " ", ">=", " ", "0"}]}], "}"}], "]"}]}]}], "]"}]}],
"\n"}], "\n",
RowBox[{
RowBox[{"hp2ar", "[", "hp_", "]"}], " ", ":=", " ",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"hp2arUnbounded", ",", " ", "hp2arCondition"}], "}"}], "}"}],
" ", "/.", " ",
RowBox[{"x", " ", "\[Rule]", " ", "hp"}]}], ",", " ", "0"}],
"]"}]}], "\n",
RowBox[{
RowBox[{
RowBox[{"hp2mr", "[", "hp_", "]"}], " ", ":=", " ",
RowBox[{"Piecewise", "[",
RowBox[{
RowBox[{
RowBox[{"{",
RowBox[{"{",
RowBox[{"hp2mrUnbounded", ",", " ", "hp2mrCondition"}], "}"}], "}"}],
" ", "/.", " ",
RowBox[{"x", " ", "\[Rule]", " ", "hp"}]}], ",", " ", "0"}], "]"}]}],
"\n"}], "\n",
RowBox[{
RowBox[{"(*",
RowBox[{"FindInstance", "[",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"(", "hp2mrUnbounded", ")"}], " ", ">=", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "\[GreaterEqual]", " ", "0"}], ",", " ",
RowBox[{"scale", " ", "<", " ", "1"}], ",", " ",
RowBox[{"x", " ", ">", " ", "0"}]}], "}"}], ",", " ",
RowBox[{"{",
RowBox[{"x", ",", "scale"}], "}"}]}], "]"}], "*)"}]}]}], "Code",
CellChangeTimes->{{3.6008582880445685`*^9, 3.600858317031201*^9}, {
3.6008583606546545`*^9, 3.600858405198162*^9}, {3.6008585196806087`*^9,
3.6008585279720755`*^9}, {3.6008587146705756`*^9, 3.600858726157227*^9}, {
3.6008587775481143`*^9, 3.6008588250217857`*^9}, {3.6008588760766635`*^9,
3.600858914301814*^9}, {3.600858949211779*^9, 3.6008590917868013`*^9}, {
3.6008591925964737`*^9, 3.6008592480765963`*^9}, 3.6008592981444135`*^9, {
3.600859361409969*^9, 3.600859390363598*^9}, {3.600859468987023*^9,
3.6008595268242817`*^9}, {3.600859843032069*^9, 3.6008598535276604`*^9}, {
3.600860042906317*^9, 3.6008600814234843`*^9}, {3.6008603562669497`*^9,
3.6008605080604973`*^9}, 3.6008605564452143`*^9, {3.600860664480299*^9,
3.6008607054555993`*^9}, {3.6008607416336403`*^9,
3.6008607849750795`*^9}, {3.6008608287525377`*^9, 3.600860853552938*^9}, {
3.600860894119221*^9, 3.600860904605811*^9}, {3.6008610451407194`*^9,
3.600861046743804*^9}, {3.6008614013037605`*^9, 3.6008614307954197`*^9}, {
3.600861479250142*^9, 3.6008615097028556`*^9}, {3.6008616865738096`*^9,
3.6008617041587987`*^9}, 3.6008619403410883`*^9, 3.600862041899803*^9, {
3.601051211349924*^9, 3.601051228659898*^9}, {3.6010513775892735`*^9,
3.601051379611392*^9}, {3.601051489108554*^9, 3.601051493181777*^9}, {
3.601056284804413*^9, 3.601056376189551*^9}, {3.6013154115232525`*^9,
3.601315475223838*^9}, {3.6013157473481503`*^9, 3.601315751378377*^9}, {
3.601415561251135*^9, 3.601415611779002*^9}, {3.601415642149723*^9,
3.6014156822949924`*^9}, {3.601418138326144*^9, 3.6014181430784135`*^9}, {
3.6014183337762127`*^9, 3.6014184067263517`*^9}, {3.6014184712950087`*^9,
3.6014184716990323`*^9}, {3.601418554654726*^9, 3.601418590927787*^9}, {
3.601419044124462*^9, 3.6014191016857185`*^9}, {3.601419585183117*^9,
3.601419604659215*^9}, {3.6014196529089537`*^9, 3.6014196603893776`*^9}, {
3.601419902735102*^9, 3.601419907973405*^9}, {3.601419955279084*^9,
3.6014199993435755`*^9}, {3.6014315629997168`*^9, 3.601431616161729*^9}, {
3.6014316528558073`*^9, 3.60143174919127*^9}, {3.6014317824811516`*^9,
3.6014317936087866`*^9}, {3.6014318731052856`*^9, 3.601431888673174*^9}, {
3.6014320291441317`*^9, 3.6014320969129705`*^9}, {3.601432147169818*^9,
3.6014321591764984`*^9}, {3.6014322985603943`*^9,
3.6014323331363487`*^9}, {3.601433938209285*^9, 3.6014339446406493`*^9}, {
3.6014356655991497`*^9, 3.601435669200354*^9}, {3.6014359942077675`*^9,
3.6014359951268244`*^9}, {3.6014360468207483`*^9,
3.6014360684749804`*^9}, {3.6014361796602745`*^9, 3.601436251531351*^9}, {
3.6014363144389105`*^9, 3.601436369006007*^9}, {3.6014364069811583`*^9,
3.6014364071901674`*^9}, {3.6014364953971677`*^9, 3.601436534164364*^9}, {
3.601436608350567*^9, 3.601436608991603*^9}, {3.6014367516366854`*^9,
3.6014367532727776`*^9}, {3.6014368188454924`*^9, 3.601436827067958*^9}, {
3.6014368909465733`*^9, 3.601437019229845*^9}, {3.6014371261439023`*^9,
3.6014371280870075`*^9}, {3.601437269198003*^9, 3.6014372959295216`*^9}, {
3.6014373496755667`*^9, 3.6014374239297743`*^9}, {3.601437457308665*^9,
3.6014377747836466`*^9}, {3.6014378221273336`*^9, 3.601437848263815*^9}, {
3.601437901863846*^9, 3.6014379137195234`*^9}, {3.6014380082018757`*^9,
3.601438022412681*^9}, {3.6014902593701687`*^9, 3.601490267215618*^9}, {
3.601491196616274*^9, 3.6014912003034773`*^9}}],
Cell[BoxData[
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
FractionBox[
RowBox[{
RowBox[{"-", "9000000"}], "+",
SuperscriptBox["x", "2"]}],
RowBox[{"9000", " ", "scale"}]],
RowBox[{
FractionBox["1", "3"], "\[LessEqual]", "scale", "\[LessEqual]",
FractionBox["2", "3"]}]},
{
RowBox[{
RowBox[{"250", " ",
RowBox[{"(",
RowBox[{"3", "-",
FractionBox["5", "scale"]}], ")"}]}], "-",
FractionBox[
RowBox[{
SuperscriptBox[
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+", "scale"}], ")"}], "2"], " ",
SuperscriptBox["x", "2"]}],
RowBox[{"1000", " ", "scale", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "5"}], "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}]]}],
RowBox[{
RowBox[{"3", " ", "scale"}], "<", "1"}]},
{
RowBox[{
FractionBox[
RowBox[{"750", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "2"}], "+", "scale"}], ")"}]}], "scale"], "+",
FractionBox[
RowBox[{"scale", " ",
SuperscriptBox["x", "2"]}],
RowBox[{"1000", " ",
RowBox[{"(",
RowBox[{"2", "+",
RowBox[{"3", " ", "scale"}]}], ")"}]}]]}],
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.601431575027398*^9, 3.601431616995781*^9, {3.601431672510921*^9,
3.6014317593578415`*^9}, {3.601431790021579*^9, 3.6014317944768314`*^9}, {
3.6014318792616353`*^9, 3.601431889180197*^9}, {3.6014320796559925`*^9,
3.6014320972329845`*^9}, {3.601432149846965*^9, 3.6014321597255297`*^9}, {
3.6014322992714305`*^9, 3.601432333819393*^9}, {3.601433941220455*^9,
3.6014339452256823`*^9}, 3.6014356698243895`*^9, 3.601435995587846*^9, {
3.6014360492858934`*^9, 3.601436068924*^9}, {3.6014361808073387`*^9,
3.6014362090839415`*^9}, {3.601436240051696*^9, 3.601436251965372*^9}, {
3.6014363203092427`*^9, 3.601436369525031*^9}, 3.6014364076151953`*^9, {
3.6014364958731894`*^9, 3.601436534680394*^9}, 3.6014366095836368`*^9,
3.6014367538258038`*^9, {3.601436821271624*^9, 3.6014368277859936`*^9}, {
3.6014368915646076`*^9, 3.601436959132436*^9}, {3.6014369899821835`*^9,
3.601437019799872*^9}, 3.601437128611039*^9, 3.601437317597745*^9, {
3.6014373524687195`*^9, 3.6014373792752385`*^9}, 3.6014374248088183`*^9, {
3.601437462897977*^9, 3.601437466393175*^9}, {3.601437504499334*^9,
3.6014375452796435`*^9}, 3.601437590529207*^9, 3.6014376219189854`*^9, {
3.601437671862816*^9, 3.6014377752896748`*^9}, {3.6014378230243797`*^9,
3.6014378491238575`*^9}, {3.601437902592887*^9, 3.601437915064595*^9}, {
3.6014380103849945`*^9, 3.6014380235557404`*^9}, 3.6014414669318256`*^9,
3.6014450279385743`*^9, 3.6014902679106536`*^9, 3.601491202217586*^9,
3.6015050324508667`*^9}],
Cell[BoxData[
TagBox[GridBox[{
{"\[Piecewise]", GridBox[{
{
RowBox[{"-",
FractionBox[
RowBox[{
RowBox[{"-", "9000000"}], "+",
SuperscriptBox["x", "2"]}],
RowBox[{"9000", " ",
RowBox[{"(",
RowBox[{
RowBox[{"-", "1"}], "+", "scale"}], ")"}]}]]}],
RowBox[{
FractionBox["1", "3"], "\[LessEqual]", "scale", "\[LessEqual]",
FractionBox["2", "3"]}]},
{
FractionBox[
RowBox[{
RowBox[{"250000", " ",