(* 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[ 36850, 964] NotebookOptionsPosition[ 33605, 902] NotebookOutlinePosition[ 33991, 919] CellTagsIndexPosition[ 33948, 916] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Paternal care, mating, and good genes", "Title", CellChangeTimes->{ 3.813921013748481*^9},ExpressionUUID->"3a99eee4-b9c0-4658-acaa-\ fc8df7c8d417"], Cell[TextData[StyleBox["For Fitzpatrick, CL; Ciresi, CM; and Wade, MJ 2020. \ The evolutionary genetics of paternal care: how good genes and extra-pair \ copulation affect the trade-off between paternal care and mating success. \ Ecology and Evolution.", FontSize->16, FontWeight->"Bold"]], "Text", CellChangeTimes->{ 3.813921040711424*^9},ExpressionUUID->"4f9d397c-94aa-44dd-ac0b-\ 0ed461833475"], Cell[CellGroupData[{ Cell["\<\ This notebook makes a figure that compares the solution for the boundary \ between negative and positive evolution of an allele with pleiotropic effects \ on paternal care and mating success. The point of this figure is to visualize \ the modeled relationship between the indirect genetic effect of paternal care \ and the male mating advantage--this time associated with extra-pair \ copulations at a frequency of 0.3--and the direct genetic effect on that \ relationship. These figures are created using the weak selection approximation equations \ for the change in allele frequency per generation (equation 10b), plugging in \ earlier equations (Eqns 1 and 5) to get deltaq in terms of cmale (indirect \ effect of paternal care), mating advantage (m), and direct genetic effect (s). \ \>", "Subsubsection", CellChangeTimes->{{3.7810192433577223`*^9, 3.781019271631597*^9}, { 3.781019315550785*^9, 3.781019331939583*^9}, {3.781019425230774*^9, 3.781019446516225*^9}, 3.813921105086659*^9, {3.813921276490803*^9, 3.813921282549383*^9}, {3.813921321143654*^9, 3.813921341766429*^9}},ExpressionUUID->"5d452cec-9b58-467b-8d50-\ 1554f0b862af"], Cell[TextData[StyleBox["Eqn 1: Wm = (1 + 2qm)\n\nEqn 5: WNRM = 1 + 2q(s + c\ \[Mars]) + (\[CapitalDelta]qm)(s + 2c\[Mars])\n\nEqn 10b: \ \[CapitalDelta]qtotal ~ (pq){s + [c\[Mars]/2]}/WNRM + (pq)(m/2Wm)\n\n\n\n", FontSize->24]], "Text", CellChangeTimes->{ 3.7810194520974503`*^9, 3.8139211849234543`*^9, {3.8139231280183887`*^9, 3.8139231286174498`*^9}},ExpressionUUID->"a4764c8a-6b28-420f-961f-\ f508b2358f64"], Cell[TextData[StyleBox["First, find \[CapitalDelta]qtotal in terms of q, \ cmale, m, and s. ", FontSize->16]], "Text", CellChangeTimes->{{3.813921362315734*^9, 3.813921452345962*^9}, { 3.8139215122810717`*^9, 3.813921518615109*^9}, {3.8139215870249977`*^9, 3.81392163581843*^9}, 3.813921716382063*^9},ExpressionUUID->"30eb21e3-aaea-4a3f-bdca-\ d87ab1611ebd"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"deltaqtotal", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", "pq", ")"}], RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"(", RowBox[{"cmale", "/", "2"}], ")"}]}], ")"}], "/", "WNRM"}]}], "+", RowBox[{ RowBox[{"(", "pq", ")"}], RowBox[{"(", RowBox[{ RowBox[{"m", "/", "2"}], "Wm"}], ")"}]}]}]}]], "Input", CellChangeTimes->{3.781019522458022*^9}, CellLabel->"In[1]:=",ExpressionUUID->"0356a3e5-f71e-4473-91bb-52f2f7687764"], Cell[BoxData[ RowBox[{ FractionBox[ RowBox[{"m", " ", "pq", " ", "Wm"}], "2"], "+", FractionBox[ RowBox[{"pq", " ", RowBox[{"(", RowBox[{ FractionBox["cmale", "2"], "+", "s"}], ")"}]}], "WNRM"]}]], "Output", CellChangeTimes->{3.782860309258174*^9, 3.803656644275419*^9, 3.804062186662075*^9, 3.8040883927576857`*^9}, CellLabel->"Out[1]=",ExpressionUUID->"bb96d255-f59f-4bb1-8fe3-51a2998e8010"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"deltaqtotal", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], "q"}], ")"}], RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"(", RowBox[{"cmale", "/", "2"}], ")"}]}], ")"}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", RowBox[{"(", RowBox[{"s", "+", "c\[Mars]"}], ")"}]}], "+", RowBox[{ RowBox[{"(", "\[CapitalDelta]qm", ")"}], RowBox[{"(", RowBox[{"s", "+", RowBox[{"2", "c\[Mars]"}]}], ")"}]}]}], ")"}]}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], StyleBox[" ", Background->RGBColor[1, 1, 0]], "q"}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"m", "/", "2"}], RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", " ", "m"}]}], ")"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{3.781019534059299*^9}, CellLabel->"In[2]:=",ExpressionUUID->"37639d64-301e-47ce-8692-54c37a5aef74"], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", "m", " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "m", " ", "q"}]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{ FractionBox["cmale", "2"], "+", "s"}], ")"}]}], RowBox[{"1", "+", RowBox[{"2", " ", "q", " ", RowBox[{"(", RowBox[{"c\[Mars]", "+", "s"}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "c\[Mars]"}], "+", "s"}], ")"}], " ", "\[CapitalDelta]qm"}]}]]}]], "Output", CellChangeTimes->{3.7828603130040627`*^9, 3.803656647538183*^9, 3.804062190746767*^9, 3.804088397628933*^9}, CellLabel->"Out[2]=",ExpressionUUID->"6b001fb7-bca7-4b84-a9f3-d4e690644af1"] }, Open ]], Cell[TextData[StyleBox["\n\n\n\n\nThis time, as opposed to the scenario in \ which extra-pair copulations do not happen, the male mating advantage is \ multiplied by the frequency of extra-pair copulations (e) and the indirect \ effect of paternal care is multiplied by (1-e). ", FontSize->16]], "Text", CellChangeTimes->{ 3.8139217342444897`*^9, {3.8139231360813713`*^9, 3.813923136884953*^9}},ExpressionUUID->"dc5459db-761a-4720-a2b0-\ 4cba0c163738"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"deltaqtotal", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], "q"}], ")"}], RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], RowBox[{"cmale", "/", "2"}]}], ")"}]}], ")"}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", RowBox[{"(", RowBox[{"s", "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], "cmale"}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"e", " ", "m", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], RowBox[{"q", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", " ", "e", " ", "m"}]}], ")"}]}]}], ")"}], RowBox[{"(", RowBox[{"s", "+", RowBox[{"2", RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], "cmale"}]}], ")"}]}]}], ")"}]}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], StyleBox[" ", Background->RGBColor[1, 1, 0]], "q"}], ")"}], RowBox[{"(", RowBox[{"e", " ", RowBox[{"m", "/", "2"}], RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", " ", "e", " ", "m"}]}], ")"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{3.781019546004253*^9, 3.781044282233348*^9}, CellLabel->"In[3]:=",ExpressionUUID->"d8c6c9a6-3c50-4aa6-9e66-093f893cb0c1"], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", "e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}], ")"}]}], "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{ RowBox[{ FractionBox["1", "2"], " ", "cmale", " ", RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}]}], "+", "s"}], ")"}]}], RowBox[{"1", "+", RowBox[{"2", " ", "q", " ", RowBox[{"(", RowBox[{ RowBox[{"cmale", " ", RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}]}], "+", "s"}], ")"}]}], "+", FractionBox[ RowBox[{"e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "cmale", " ", RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}]}], "+", "s"}], ")"}]}], RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}]]}]]}]], "Output", CellChangeTimes->{3.781044422938178*^9, 3.782860315787936*^9, 3.803656649887671*^9, 3.8040621998077593`*^9, 3.804088398620824*^9}, CellLabel->"Out[3]=",ExpressionUUID->"618d6c0d-4024-4840-ad54-11d0b60b89ef"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"FullSimplify", "[", RowBox[{"0", "==", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], "q"}], ")"}], RowBox[{ RowBox[{"(", RowBox[{"s", "+", RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], RowBox[{"cmale", "/", "2"}]}], ")"}]}], ")"}], "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", RowBox[{"(", RowBox[{"s", "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], "cmale"}]}], ")"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"e", " ", "m", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], RowBox[{"q", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", " ", "e", " ", "m"}]}], ")"}]}]}], ")"}], RowBox[{"(", RowBox[{"s", "+", RowBox[{"2", RowBox[{"(", RowBox[{"1", "-", "e"}], ")"}], "cmale"}]}], ")"}]}]}], ")"}]}]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], StyleBox[" ", Background->RGBColor[1, 1, 0]], "q"}], ")"}], RowBox[{"(", RowBox[{"e", " ", RowBox[{"m", "/", "2"}], RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", "q", " ", "e", " ", "m"}]}], ")"}]}], ")"}]}]}]}], "]"}]], "Input", CellChangeTimes->{{3.781019691489245*^9, 3.7810197026152277`*^9}, { 3.78104429057896*^9, 3.7810442945357733`*^9}}, CellLabel->"In[4]:=",ExpressionUUID->"b52236af-bb44-4d95-bdf3-ef23a61e9ec5"], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"e", " ", "m"}], "-", FractionBox[ RowBox[{ RowBox[{"cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}]}], "-", RowBox[{"2", " ", "s"}]}], RowBox[{"1", "-", RowBox[{"2", " ", "cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}]}], "+", RowBox[{"2", " ", "q", " ", "s"}], "+", RowBox[{"e", " ", "m", " ", "q", " ", RowBox[{"(", RowBox[{"2", "+", "s", "+", RowBox[{"3", " ", "q", " ", "s"}]}], ")"}]}]}]]}], ")"}]}], "\[Equal]", "0"}]], "Output", CellChangeTimes->{3.7810197475072823`*^9, 3.781044297902091*^9, 3.78104443947904*^9, 3.7828603209274263`*^9, 3.803656654080372*^9, 3.804062205498034*^9, 3.8040884030195932`*^9}, CellLabel->"Out[4]=",ExpressionUUID->"30877e0c-490d-49af-9101-256156422157"] }, Open ]], Cell["\<\ Then solve for cmale and m in terms of each other when s = - 0.1, 0, and 0.1 \ (when q is set to 0.25 and e is set to 0.3; that is, 30% of males engage in \ extra-pair copulations). \ \>", "Text", CellChangeTimes->{ 3.781020106277916*^9, {3.7810205190194798`*^9, 3.781020530474971*^9}, { 3.8139217636066504`*^9, 3.8139217908799667`*^9}, 3.8139219143962507`*^9, { 3.813923143376781*^9, 3.813923152420474*^9}},ExpressionUUID->"538c1a11-734c-4fd2-8bba-\ efe058c35392"], Cell[BoxData[{ RowBox[{ RowBox[{"s", "=", RowBox[{"-", "0.1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"q", "=", "0.25"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", "0.3"}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.781019895384042*^9, 3.7810199023633823`*^9}, { 3.781020665158094*^9, 3.7810206680478363`*^9}, {3.781044308635157*^9, 3.781044314234817*^9}, {3.782860325222336*^9, 3.78286032549221*^9}, { 3.80406221768554*^9, 3.804062230423661*^9}, 3.8040884106845303`*^9}, CellLabel->"In[1]:=",ExpressionUUID->"f71a69c4-778d-48f7-b1c8-313ae2e529f2"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"e", " ", "m"}], "-", FractionBox[ RowBox[{ RowBox[{"cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}]}], "-", RowBox[{"2", " ", "s"}]}], RowBox[{"1", "-", RowBox[{"2", " ", "cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}]}], "+", RowBox[{"2", " ", "q", " ", "s"}], "+", RowBox[{"e", " ", "m", " ", "q", " ", RowBox[{"(", RowBox[{"2", "+", "s", "+", RowBox[{"3", " ", "q", " ", "s"}]}], ")"}]}]}]]}], ")"}]}], "\[Equal]", "0"}], ",", " ", "cmale"}], "]"}]], "Input", CellChangeTimes->{3.7810206572460814`*^9, 3.7810443342460203`*^9, 3.781044461582417*^9}, CellLabel->"In[4]:=",ExpressionUUID->"6cfa8433-0968-411e-8254-aedf886d623c"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"cmale", "\[Rule]", FractionBox[ RowBox[{"0.014285714285714285`", " ", RowBox[{"(", RowBox[{"3200.`", "\[VeryThinSpace]", "-", RowBox[{"4560.`", " ", "m"}], "-", RowBox[{"657.`", " ", SuperscriptBox["m", "2"]}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]]}], "}"}], "}"}]], "Output", CellChangeTimes->{{3.781020658433914*^9, 3.781020671139035*^9}, { 3.781044342276073*^9, 3.7810443644763536`*^9}, 3.781044464647664*^9, 3.782860335613359*^9, 3.803656716266047*^9, 3.8040622457875566`*^9, 3.804088415361167*^9, 3.808145793998197*^9, 3.808152545669692*^9}, CellLabel->"Out[4]=",ExpressionUUID->"ee81d40c-34f3-4c57-888e-e47cdc02ca1f"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"s", "=", "0"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"q", "=", "0.25"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", "0.3"}], ";"}]}], "Input", CellChangeTimes->{{3.781020677577407*^9, 3.78102068078449*^9}, { 3.808145610994673*^9, 3.8081456113691883`*^9}}, CellLabel->"In[1]:=",ExpressionUUID->"2dfc259d-8fb7-4e4c-941b-28c9fc67ed92"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"e", " ", "m"}], "-", FractionBox[ RowBox[{ RowBox[{"cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}]}], "-", RowBox[{"2", " ", "s"}]}], RowBox[{"1", "-", RowBox[{"2", " ", "cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}]}], "+", RowBox[{"2", " ", "q", " ", "s"}], "+", RowBox[{"e", " ", "m", " ", "q", " ", RowBox[{"(", RowBox[{"2", "+", "s", "+", RowBox[{"3", " ", "q", " ", "s"}]}], ")"}]}]}]]}], ")"}]}], "\[Equal]", "0"}], ",", " ", "cmale"}], "]"}]], "Input", CellChangeTimes->{3.7810207075454693`*^9, 3.781044501093061*^9}, CellLabel->"In[4]:=",ExpressionUUID->"554c05b5-2593-4c40-9b8b-c035a42dd67a"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"cmale", "\[Rule]", RowBox[{"-", FractionBox[ RowBox[{"3.4285714285714284`", " ", "m", " ", RowBox[{"(", RowBox[{"20.`", "\[VeryThinSpace]", "+", RowBox[{"3.`", " ", "m"}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.7810207126259212`*^9, 3.781044502903257*^9, 3.7828603498863983`*^9, 3.8036567253022623`*^9, 3.804062259270372*^9, 3.80408842049975*^9, { 3.8081455868835087`*^9, 3.808145616519434*^9}, 3.808152236079438*^9}, CellLabel->"Out[4]=",ExpressionUUID->"1ba3d858-cc8d-4b1d-8322-5041d21a61e9"] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"s", "=", "0.1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"q", "=", "0.25"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", "0.3"}], ";"}]}], "Input", CellChangeTimes->{{3.781020718780006*^9, 3.7810207213110247`*^9}, { 3.8040622662836018`*^9, 3.804062266669289*^9}, {3.808146558196888*^9, 3.808146561466659*^9}}, CellLabel->"In[1]:=",ExpressionUUID->"33f95984-947c-423d-ad09-dec34fcca496"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ FractionBox["1", "2"], " ", RowBox[{"(", RowBox[{"1", "-", "q"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "e", " ", "m", " ", "q"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"e", " ", "m"}], "-", FractionBox[ RowBox[{ RowBox[{"cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}]}], "-", RowBox[{"2", " ", "s"}]}], RowBox[{"1", "-", RowBox[{"2", " ", "cmale", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "e"}], ")"}], " ", "q", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"e", " ", "m", " ", RowBox[{"(", RowBox[{"1", "+", "q"}], ")"}]}]}], ")"}]}], "+", RowBox[{"2", " ", "q", " ", "s"}], "+", RowBox[{"e", " ", "m", " ", "q", " ", RowBox[{"(", RowBox[{"2", "+", "s", "+", RowBox[{"3", " ", "q", " ", "s"}]}], ")"}]}]}]]}], ")"}]}], "\[Equal]", "0"}], ",", " ", "cmale"}], "]"}]], "Input", CellChangeTimes->{3.7810207285466423`*^9, 3.781044515149741*^9}, CellLabel->"In[4]:=",ExpressionUUID->"dc22fe37-8e9b-4ea7-80bb-c251672ed6a0"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"cmale", "\[Rule]", FractionBox[ RowBox[{"0.014285714285714285`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "3200.`"}], "-", RowBox[{"5040.`", " ", "m"}], "-", RowBox[{"783.`", " ", SuperscriptBox["m", "2"]}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.781020729825838*^9, 3.781044516468658*^9, 3.782860355509798*^9, 3.803656732918436*^9, 3.80406227578855*^9, 3.8040884257430773`*^9, 3.808146571821418*^9, 3.8081526993425694`*^9}, CellLabel->"Out[4]=",ExpressionUUID->"9e495224-87e2-494c-b841-70f5a2d805c6"] }, Open ]], Cell["\<\ Now I have equations in terms of m and cmale when s = -0.1, 0, 0.1. Make a \ plot. Solid line => no viability effects (s = 0). Dashed black line => s = \ -0.1. Gray line => s = 0.1.\ \>", "Text", CellChangeTimes->{{3.781020761286751*^9, 3.7810207932676067`*^9}, { 3.804062285751356*^9, 3.804062289198227*^9}, {3.813922130425007*^9, 3.813922198773364*^9}, {3.81392316409999*^9, 3.813923165255727*^9}},ExpressionUUID->"6b019b18-29b4-431b-bd10-\ 4677866e4910"], Cell[BoxData[ RowBox[{ RowBox[{"PlotViabilityWithEPC", "=", " ", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"0.014285714285714285`", " ", RowBox[{"(", RowBox[{"3200.`", "\[VeryThinSpace]", "-", RowBox[{"4560.`", " ", "m"}], "-", RowBox[{"657.`", " ", SuperscriptBox["m", "2"]}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]], ",", RowBox[{"-", FractionBox[ RowBox[{"3.4285714285714284`", " ", "m", " ", RowBox[{"(", RowBox[{"20.`", "\[VeryThinSpace]", "+", RowBox[{"3.`", " ", "m"}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]]}], ",", FractionBox[ RowBox[{"0.014285714285714285`", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "3200.`"}], "-", RowBox[{"5040.`", " ", "m"}], "-", RowBox[{"783.`", " ", SuperscriptBox["m", "2"]}]}], ")"}]}], RowBox[{"160.`", "\[VeryThinSpace]", "+", RowBox[{"24.`", " ", "m"}], "+", RowBox[{"9.`", " ", SuperscriptBox["m", "2"]}]}]]}], "}"}], ",", RowBox[{"{", RowBox[{"m", ",", RowBox[{"-", "0.5"}], ",", "0.5"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "0.5"}], ",", ".5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.5"}], ",", "0.5"}], "}"}]}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"m", ",", "cmale"}], "}"}]}], ",", RowBox[{"Filling", "\[Rule]", "Top"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Black", ",", RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.04", ",", " ", "0.02"}], "}"}], "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"Black", ",", "Thick"}], "}"}], ",", "Gray"}], "}"}]}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{ 3.7810208382172318`*^9, {3.781044526879842*^9, 3.78104457755939*^9}, { 3.782859918586458*^9, 3.782859920393354*^9}, {3.782860733278945*^9, 3.7828607682809753`*^9}, {3.803656763265592*^9, 3.803656785566168*^9}, { 3.804062317456562*^9, 3.804062355725371*^9}, {3.8040884663741617`*^9, 3.804088494065114*^9}, {3.804245274906399*^9, 3.804245277019252*^9}}, CellLabel->"In[4]:=",ExpressionUUID->"608157d2-c0db-47ad-8059-0e5bf209bbc3"], Cell[CellGroupData[{ Cell[BoxData["PlotViabilityWithEPC"], "Input", CellChangeTimes->{3.804062434282929*^9}, CellLabel->"In[5]:=",ExpressionUUID->"74b5d9fc-7204-40b1-b91b-c807e923ae88"], Cell[BoxData[ GraphicsBox[{GraphicsComplexBox[CompressedData[" 1:eJxl1nk01ekfB3BZi0Kp0CgVUYlfSdrnY4yyZG8h1dgmmkGiBS0Sd2RatUxj 0jITFS2WWzKyfOzu/V4X187l+62ENlspKvKb3+/8nmfO+c1zznOe8/nn8zzn fM77dZ45PsGuO2VlZGS6/9r/OU1g+quxMQ5l/rueQFg1L8i/gMWJb+b5vKjh oL5e5VWRdzvybaas71zDwczC4KnLddrQcXqikjqfhayAUgW19lZ8sSJC4rOA Ba9TS/pboluw8hsTLd7Fdhg+fygOoRn78sTh4ye0Q5Iqq53wtBFn1uXvcAhp gyt2S0Lrfm5AX4h/LpRKITrQO3yBbT0eehvoJHWSwiG5O2GhPbVoVGWskpTd Ct5R1fqNKRL0Nese017QCr/kdyZc+bEGTWJ9No5eaYHH3PyJ9+ZWYed5F0m8 fAvoqVtYbromwgdt20Qj+5uBvzDux2cdAux98apbQ9oEXFnmbHPzcqz47BWX atkEnYt+McVHJTjap2G49mEjTI5MNTcZV4T8887bDms3wiLzvkU33+Vh8jlB hDiuAWKsempywrPx1pVBHvO6HqZtSovfqslHC1k+L9ClHvgfytQE+ndQt43N /LqoDm5KO0+/5F1Db/mpt0v062CHi85DqS8PK6zZgbjztZCjfKiheH8MnMsK vaQhVwvLno2bdcz9OrxVef/Ew18Cev25H1cVpwLXyD4ZbqgB26QWsaM7H1ad Fh3ctqIG1jlMkXuflQ2ue2O4xaur4axnU9Ti7DzImqs1fHhGFVQ4zpzkPKUI 2pKUz4Q8rYSeOdMiay+XwJ/+V1fY54tg/z4j77PTy8Hmg+ZZuesMNLWpu/oJ BKAOSz2O7xKCV8SabiV7EUyXf3xLyVoAJUbLT5QoVIFInccLmVQBBhPyXiTb 1IDMYWHl8xdlMDMvYJ1CnAQMlCrjdjClIPNJNjX0WS0ctZu/uudCCQy55H9X Nq8e2gyXG8eHFoPUzPOrebwGqGNzG3evLII9buqT57Y1wvJf9T5GKhaC9sXn 7vIzm+GW1LjGRpIPadt+v/0uvAUkao+7wv1ywW3tURtBfiuoM0GGTfNyIPO7 oLmrRqQQrjXIXB94BMPK38rlbG6Hzdu/Ly488xBcJavNZyawoFhedntkMx/+ lwcYvLfOa15OOq1JPizGL4hXteegSX/jGw3tTPj//AgDTTumvWUhUedo1KHs TMhO0l2rdIuDm7JTV5QfzPhHP1Xd7FPaZSzY3Y7v+lORDyObk7uUIzhotQiS vBRmQOjA5QtydRxoKtYsSDXK+Md9i8MjvBPTWIh19x7ZasKHWel+uad8OODt tlaLHMiAi3Y7dv2Ux8EU7SHWfGMGOCc0Fxx7yUFicZX2TPl/vqcouXXD9iQW wvfs0zWy5EO6tOKmgxsHdSs+CW0UMsEgrMB1fSYHo43abl/8MuBp6cpUS5YD 63TNnmeaGYABYfWFgxzM2Gbulf4kHZRXqxTlXWVBAYcaNez5MBwyPCPNhQNF JzNHE7VM2D2pr+LW3b/6eUauerYnAxZu+3FWUjMHsfY2X8rmZMB42RTbOX0c 6LUfTd83mA4HGtZ7en3iwDacWfeqks4HE2Y8toqXF1G/Ppj63hwbEVK/PH5Y MnshK6B+TdZ2MfSQVFC/Vq89GdgrKKd+vbEJlF+QWkb98lilxRYlllK/Ft8r NOw/UEL9svXt903aUkz98pzShoq2RdSvKObb5iG9QurXbwajsj/rFVC/RmQu 27SOz6N+2bANz9dhDvVrsEtQ9iYxm/pluiHzyz3bLOqX7s4jusHjHlC/Stf4 jT3lMqhfo1yqmYruferX6ltmO5ZUplC/Tr2Xt2W/SqZ+XTNcOsls3HXql5zE MqzW/xL1K1h+o0Ho/BPUL9Oooht6B0OoXxkBAiO1tD1I/NJq0zcQtp9A4peF RrmOn8YlJH791Bu+K2bPdSR+GWUf/cHdORmJXxccTCYkJaQg8Wv3lTe+19fc R+KX/cCE9Y6CDCR+nVQyONDexUfi1xHf7od9pllI/HLcMHGoZ0s2Er9cTwc+ WJCSg8Qvxjmy1kqSi8SvhzGbxtJ685H4pTV2eVK4TiESvx4WPL3/s3EREr++ yZRb2r6pGIlfpvdcYo0PlCDxK/yGyh9hMaVI/PIIsRNX3CpD4tdNj+eXpj0q R+JXxvhpj/yxAolfPsFXNaNZARK/LC4/m3b/pRBJntXOqjSpaIjQ9bj47lJP DpWOy0m6oxjql+PnDU4HIxm06BhgfkvnUK79jGZpMUP9avDq1PzKjEHiV9rl VZcZAYPLJ/NaK5s45Kkqmd3rYKhfjSmavZ2qDBK/prwK2pScxCDxS2v8UZ2e HgZN5zv4dL3m0El1ccGsQYb6lfP0pyHbESESv1xteEK1cwwSv2xeqG/aLmWQ +FU+e+ztYRkRZjxYYvDuPYdZ4nd/uI0w1K+8zknJn98Kkfglc8J2YdlxBolf Uf52v2jUMkj8qlysYLBsiEHil9O87TGlE0QY/FEv7uMnDhtGwz7vlRFRzx4l OHy/qVeIxDOpqmWpbgyDxLO91+x0g8UMEs/WJtTOP/eWQeJZ9Ibi4tkKIiSe eSQwZrtURWR+6Gbleua1cxYSz7Y8V353wykbiWfeR2J+fXMwB4lng4OGx5qi c5F4lt2bfUFyLB+JZ0c6lJwC7yMSzwwUApadbi1E4lmzw6eAx4NFSDxb9ieM c1AuQeLZysI5uqmTSpF4tkWpfucWwzIknhlZSkbumJQj8cxbxei72UsqkHjG sxbnS20FSDzbW9tbM9FZiMSzu7H6koOBDBLPgs9bDXTEiJB4ZpZU5JgQV4nE M/9l9jp7k8RIPFtY88Eq+m4VEs9iC/w6bwqrkXhmf1VFbn1IDRLPvhz74GM1 WIPEsxn1abw10RIknpWvFtRnK9Yi8axvuEaNF1dLPXu41W3v8ml11DPr38KS dybUUc9U7rZ2BxjUU8/23Z5zIi2znnpWkOz2tf/SBuqZcoTM70+wgXoWfW8j iC0bqWevvlyL7i9upJ5Zu5w+nuvQRD3jG3ZUfapqop6pOQ0nij2aqWe+XUnV U9ubqWdB69StYre3UM/c9f91d3dXC/XshuuRUolfK/Xszu/e1y4OtFLPxGEd skERUuqZj7Zmg3BISj1zej+s6BfVRj2bL0avk2Nt1LMq1dqU6kPt1LPwWHMb U1mWelaoUg+LjrLUs3SrqTGLlTnqmafHwjVazo+QeMYr9bLDUpZ6tqflwOf9 NVlIPBs1VBj++AdLPduSCPMPd7PUMy+dHWX347OQeJZam2t4MZ6lntmsEgpm NbLUs+sT3vk1DbPUs/dTN74U7clC4ln78msHnhxnqWfv9LY4BItZ6tnJtCqV hf0s9SxQErqSN46jni0zvnpS6pWFxLPWca8DYmNY6tlSLcF2TSFLPYN2daO6 1yz17KrOhqfGoyz1LLPfJ7hVgaOeVYq7ozrcs5B4FhGn4rH1GEs9uzK32u11 Ofu3Z+av7t9+yVLPPO9+nyL7maWeGVrz3dPlOOqZT3fMEYXxHLbMlld+bPz3 f/XfK4KOLg== "], {{ {GrayLevel[0], Opacity[0.2], EdgeForm[None], GraphicsGroupBox[ PolygonBox[{{50, 65, 231, 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, 71, 66, 61, 57, 54, 52, 72, 67, 62, 58, 55, 73, 68, 63, 59, 74, 69, 64, 75, 70, 76}}]]}, {}, {GrayLevel[0], Opacity[0.2], EdgeForm[None], GraphicsGroupBox[ PolygonBox[{{126, 56, 60, 77, 147, 141, 136, 132, 129, 127, 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, 148, 142, 137, 133, 130, 128, 149, 143, 138, 134, 131, 150, 144, 139, 135, 151, 145, 140, 152, 146, 153}}]]}, {}, {GrayLevel[0.5], Opacity[0.2], EdgeForm[None], GraphicsGroupBox[ PolygonBox[{{203, 51, 53, 154, 224, 218, 213, 209, 206, 204, 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, 225, 219, 214, 210, 207, 205, 226, 220, 215, 211, 208, 227, 221, 216, 212, 228, 222, 217, 229, 223, 230}}]]}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}, {}}, {{}, \ {}, TagBox[ {GrayLevel[0], AbsoluteThickness[1.6], Opacity[1.], Dashing[{0.04, 0.02}], LineBox[{231, 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, 71, 66, 61, 57, 54, 52, 72, 67, 62, 58, 55, 73, 68, 63, 59, 74, 69, 64, 75, 70, 76, 50}]}, Annotation[#, "Charting`Private`Tag$3176#1"]& ], TagBox[ {GrayLevel[0], Thickness[Large], Opacity[1.], LineBox[{77, 147, 141, 136, 132, 129, 127, 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, 148, 142, 137, 133, 130, 128, 149, 143, 138, 134, 131, 150, 144, 139, 135, 151, 145, 140, 152, 146, 153, 126}]}, Annotation[#, "Charting`Private`Tag$3176#2"]& ], TagBox[ {GrayLevel[0.5], AbsoluteThickness[1.6], Opacity[1.], LineBox[{154, 224, 218, 213, 209, 206, 204, 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, 225, 219, 214, 210, 207, 205, 226, 220, 215, 211, 208, 227, 221, 216, 212, 228, 222, 217, 229, 223, 230, 203}]}, Annotation[#, "Charting`Private`Tag$3176#3"]& ]}}], {}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox[ TagBox["m", HoldForm], TraditionalForm], FormBox["cmale", TraditionalForm]}, AxesOrigin->{0, 0}, DisplayFunction->Identity, Frame->{{False, False}, {False, False}}, FrameLabel->{{None, None}, {None, None}}, FrameTicks->{{Automatic, Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, Charting`ScaledFrameTicks[{Identity, Identity}]}}, GridLines->{None, None}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None, "CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Identity[#]& )[ Part[#, 1]], (Identity[#]& )[ Part[#, 2]]}& )}, "AxesInFront" -> True}, PlotRange->{{-0.5, 0.5}, {-0.5, 0.5}}, PlotRangeClipping->True, PlotRangePadding->{{0, 0}, {0, 0}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.8040624356874533`*^9, 3.8040885035317383`*^9, 3.804245285858976*^9}, CellLabel->"Out[5]=",ExpressionUUID->"9b3de18a-a332-41ed-bab4-6201e9a0895b"] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1501, 972}, WindowMargins->{{Automatic, 355}, {Automatic, 22}}, Magnification:>1.5 Inherited, FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \ 2018)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 156, 3, 146, "Title",ExpressionUUID->"3a99eee4-b9c0-4658-acaa-fc8df7c8d417"], Cell[739, 27, 402, 8, 90, "Text",ExpressionUUID->"4f9d397c-94aa-44dd-ac0b-0ed461833475"], Cell[CellGroupData[{ Cell[1166, 39, 1173, 20, 414, "Subsubsection",ExpressionUUID->"5d452cec-9b58-467b-8d50-1554f0b862af"], Cell[2342, 61, 420, 7, 473, "Text",ExpressionUUID->"a4764c8a-6b28-420f-961f-f508b2358f64"], Cell[2765, 70, 372, 7, 54, "Text",ExpressionUUID->"30eb21e3-aaea-4a3f-bdca-d87ab1611ebd"], Cell[CellGroupData[{ Cell[3162, 81, 513, 16, 46, "Input",ExpressionUUID->"0356a3e5-f71e-4473-91bb-52f2f7687764"], Cell[3678, 99, 427, 11, 79, "Output",ExpressionUUID->"bb96d255-f59f-4bb1-8fe3-51a2998e8010"] }, Open ]], Cell[CellGroupData[{ Cell[4142, 115, 1125, 37, 46, "Input",ExpressionUUID->"37639d64-301e-47ce-8692-54c37a5aef74"], Cell[5270, 154, 900, 27, 83, "Output",ExpressionUUID->"6b001fb7-bca7-4b84-a9f3-d4e690644af1"] }, Open ]], Cell[6185, 184, 460, 8, 267, "Text",ExpressionUUID->"dc5459db-761a-4720-a2b0-4cba0c163738"], Cell[CellGroupData[{ Cell[6670, 196, 1650, 54, 109, "Input",ExpressionUUID->"d8c6c9a6-3c50-4aa6-9e66-093f893cb0c1"], Cell[8323, 252, 1352, 39, 95, "Output",ExpressionUUID->"618d6c0d-4024-4840-ad54-11d0b60b89ef"] }, Open ]], Cell[CellGroupData[{ Cell[9712, 296, 1784, 56, 109, "Input",ExpressionUUID->"b52236af-bb44-4d95-bdf3-ef23a61e9ec5"], Cell[11499, 354, 1368, 38, 73, "Output",ExpressionUUID->"30877e0c-490d-49af-9101-256156422157"] }, Open ]], Cell[12882, 395, 494, 15, 260, "Text",ExpressionUUID->"538c1a11-734c-4fd2-8bba-efe058c35392"], Cell[13379, 412, 660, 14, 171, "Input",ExpressionUUID->"f71a69c4-778d-48f7-b1c8-313ae2e529f2"], Cell[CellGroupData[{ Cell[14064, 430, 1406, 39, 117, "Input",ExpressionUUID->"6cfa8433-0968-411e-8254-aedf886d623c"], Cell[15473, 471, 872, 19, 79, "Output",ExpressionUUID->"ee81d40c-34f3-4c57-888e-e47cdc02ca1f"] }, Open ]], Cell[16360, 493, 484, 13, 171, "Input",ExpressionUUID->"2dfc259d-8fb7-4e4c-941b-28c9fc67ed92"], Cell[CellGroupData[{ Cell[16869, 510, 1379, 38, 117, "Input",ExpressionUUID->"554c05b5-2593-4c40-9b8b-c035a42dd67a"], Cell[18251, 550, 786, 18, 71, "Output",ExpressionUUID->"1ba3d858-cc8d-4b1d-8322-5041d21a61e9"] }, Open ]], Cell[19052, 571, 458, 10, 109, "Input",ExpressionUUID->"33f95984-947c-423d-ad09-dec34fcca496"], Cell[CellGroupData[{ Cell[19535, 585, 1379, 38, 117, "Input",ExpressionUUID->"dc22fe37-8e9b-4ea7-80bb-c251672ed6a0"], Cell[20917, 625, 799, 19, 79, "Output",ExpressionUUID->"9e495224-87e2-494c-b841-70f5a2d805c6"] }, Open ]], Cell[21731, 647, 483, 16, 329, "Text",ExpressionUUID->"6b019b18-29b4-431b-bd10-4677866e4910"], Cell[22217, 665, 2879, 74, 225, "Input",ExpressionUUID->"608157d2-c0db-47ad-8059-0e5bf209bbc3"], Cell[CellGroupData[{ Cell[25121, 743, 166, 2, 46, "Input",ExpressionUUID->"74b5d9fc-7204-40b1-b91b-c807e923ae88"], Cell[25290, 747, 8275, 150, 372, "Output",ExpressionUUID->"9b3de18a-a332-41ed-bab4-6201e9a0895b"] }, Open ]] }, Open ]] }, Open ]] } ] *)