(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 11.2' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 145139, 3167] NotebookOptionsPosition[ 141244, 3097] NotebookOutlinePosition[ 141582, 3112] CellTagsIndexPosition[ 141539, 3109] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Supplementary notebook Heterogeneity in behaviour and movement can influence the stability of \ predator-prey periodic travelling waves - Renato Andrade*\.b9 and Christina \ A. Cobbold\.b2\ \>", "Title", CellChangeTimes->{{3.8705177461778717`*^9, 3.87051776170647*^9}, 3.870517926333305*^9, {3.870770141813418*^9, 3.870770161460511*^9}, { 3.871045329182517*^9, 3.871045358670982*^9}, {3.871045424887229*^9, 3.871045427087481*^9}, 3.871045464459895*^9, {3.871045503855434*^9, 3.871045515591259*^9}, {3.871045610468348*^9, 3.8710456160299253`*^9}, 3.8710458157541428`*^9, {3.871046448883108*^9, 3.8710464596551456`*^9}},ExpressionUUID->"553ed30a-6730-473d-8682-\ f67e849b952c"], Cell[CellGroupData[{ Cell["\<\ Bulletin of Mathematical Biology \.b9 School of Mathematics and Statistics, University of Glasgow, Glasgow, \ G12 8QQ, United Kingdom \.b2 Boyd Orr Centre for Population and Ecosystem Health, University of \ Glasgow, Glasgow, G12 8QW, United Kingdom *Corresponding author: Renato.AntunesCostadeAndrade@glasgow.ac.uk\ \>", "Section", CellChangeTimes->{{3.871045369500484*^9, 3.8710453854416924`*^9}, { 3.871045431981436*^9, 3.871045459614421*^9}, {3.871045520425824*^9, 3.8710455413276176`*^9}, {3.871046210132594*^9, 3.871046212992589*^9}, 3.871046284036085*^9, {3.871046471287299*^9, 3.8710464714033213`*^9}},ExpressionUUID->"a747abc8-3b73-4236-b76d-\ 9f7e5a389046"], Cell[BoxData[ RowBox[{ RowBox[{"(*", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ "Computation", " ", "and", " ", "evaluation", " ", "of", " ", "\[Omega]1", " ", "according", " ", "to", " ", "the", " ", "procedure", " ", "described", " ", "by", " ", "Sherrat", " ", "et", " ", RowBox[{"al", ".", " ", RowBox[{"(", "2003", ")"}]}], " ", "The", " ", "Effects", " ", "of", " ", "the", " ", "Size", " ", "and", " ", "Shape", " ", "of", " ", "Landscape", " ", "Features", " ", "on", " ", "the", " ", "Formation", " ", "of", " ", "Traveling", " ", "Waves", " ", "in", " ", "Cyclic", " ", RowBox[{ RowBox[{"Populations", "."}], " ", "\[IndentingNewLine]", "\[IndentingNewLine]", "\\", "partial_T"}], " ", "h"}], "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\\", RowBox[{"partial", "^", "2"}]}], "_X", " ", "h"}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "h"}], ")"}], "h"}], "-", " ", RowBox[{ RowBox[{"\\", RowBox[{"Big", "[", RowBox[{ RowBox[{ RowBox[{"\\", "frac"}], RowBox[{"{", RowBox[{"C", " ", "h"}], "}"}], RowBox[{"{", RowBox[{"1", "+", RowBox[{"C", " ", "h"}]}], "}"}]}], "+", " ", RowBox[{ RowBox[{"\\", "frac"}], RowBox[{"{", RowBox[{ RowBox[{"\\", "beta"}], " ", RowBox[{"C", "\\", "eta"}], " ", "h"}], "}"}], RowBox[{ RowBox[{"{", RowBox[{"1", "+", RowBox[{ RowBox[{"C", "\\", "eta"}], " ", "h"}]}], "}"}], "\\", "Big"}]}]}], "]"}]}], RowBox[{"p", "\n", " ", "\\", "partial_T"}], " ", "p"}]}], "=", " ", RowBox[{ RowBox[{ RowBox[{"\\", "delta", "\\", RowBox[{"partial", "^", "2"}]}], "_X", " ", "p"}], "-", " ", RowBox[{ RowBox[{"\\", "frac"}], RowBox[{"{", "p", "}"}], RowBox[{"{", RowBox[{"E_", "1", "B_", "1"}], "}"}]}], "+", " ", RowBox[{ RowBox[{"\\", RowBox[{"Big", "[", RowBox[{ RowBox[{ RowBox[{"\\", "frac"}], RowBox[{"{", RowBox[{"C", " ", "h"}], "}"}], RowBox[{"{", RowBox[{"B_", "1", RowBox[{"(", RowBox[{"1", "+", RowBox[{"C", " ", "h"}]}], ")"}]}], "}"}]}], "+", " ", RowBox[{ RowBox[{"\\", "frac"}], RowBox[{"{", RowBox[{ RowBox[{"\\", "beta"}], " ", RowBox[{"C", "\\", "eta"}], " ", "h"}], "}"}], RowBox[{ RowBox[{"{", RowBox[{"B_", "1", RowBox[{"(", RowBox[{"1", "+", RowBox[{ RowBox[{"C", "\\", "eta"}], " ", "h"}]}], ")"}]}], "}"}], "\\", "Big"}]}]}], "]"}]}], "p"}]}]}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", "*)"}], "\[IndentingNewLine]", " "}]], "Input", CellChangeTimes->{{3.81417891404978*^9, 3.8141789143097677`*^9}, 3.816256573055167*^9, {3.8162599972433023`*^9, 3.816259997708579*^9}, { 3.816263781663917*^9, 3.8162637819055634`*^9}, {3.816269141377301*^9, 3.816269141997353*^9}, {3.8453812271665688`*^9, 3.845381264945526*^9}, { 3.8677486445031347`*^9, 3.8677486488376713`*^9}, {3.867748793861174*^9, 3.867748814926955*^9}, {3.8677494337739353`*^9, 3.867749563904695*^9}, { 3.867749669543584*^9, 3.867749705930317*^9}, {3.867749739983571*^9, 3.867749755681912*^9}, {3.870170912759418*^9, 3.870170932642469*^9}, 3.8701710166117563`*^9, {3.870338097458365*^9, 3.870338116262569*^9}, { 3.870339865705101*^9, 3.870339865816143*^9}, {3.870517760828545*^9, 3.8705177888400517`*^9}},ExpressionUUID->"64d0e723-f54c-4c8f-907e-\ 2dd284af16b7"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Section 0: Define the reaction functions and compute the steady states\ \>", "Section", CellChangeTimes->{{3.870517804342021*^9, 3.870517832958167*^9}, { 3.870769989639295*^9, 3.8707699919567547`*^9}, {3.870770199363674*^9, 3.870770215121654*^9}, 3.8710461785341454`*^9},ExpressionUUID->"8a8cbf01-963c-4b63-a4d5-\ d3c0812c0320"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{"ClearAll", "[", RowBox[{ "\[Beta]", ",", "c", ",", "\[Eta]", ",", "B1", ",", "B2", ",", "E1", ",", "F", ",", "G", ",", "p", ",", "h", ",", "w", ",", " ", "fullsol", ",", " ", "hs", ",", " ", "ps"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"$Assumptions", "=", RowBox[{ RowBox[{"\[Eta]", " ", ">", "0"}], " ", "&&", " ", RowBox[{"c", " ", ">", "0"}], " ", "&&", RowBox[{"B1", " ", ">", " ", "0"}], " ", "&&", RowBox[{"B2", " ", ">", " ", "0"}], " ", "&&", " ", RowBox[{"A", " ", ">", " ", "0"}], " ", "&&", RowBox[{"\[Beta]", ">", "0"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"F", "[", RowBox[{"h_", ",", "p_"}], "]"}], ":=", " ", RowBox[{ RowBox[{"c", " ", "p", " ", RowBox[{"h", "/", RowBox[{"(", RowBox[{"B1", RowBox[{"(", RowBox[{"1", "+", " ", RowBox[{"c", " ", "h"}]}], ")"}]}], ")"}]}]}], "+", RowBox[{"c", " ", "\[Beta]", " ", "\[Eta]", " ", "p", " ", RowBox[{"h", "/", RowBox[{"(", RowBox[{"B1", RowBox[{"(", RowBox[{"1", "+", RowBox[{"c", " ", "\[Eta]", " ", "h"}]}], ")"}]}], ")"}]}]}], " ", "-", " ", RowBox[{"p", "/", RowBox[{"(", RowBox[{"E1", " ", "B1"}], ")"}]}]}]}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"G", "[", RowBox[{"h_", ",", "p_"}], "]"}], ":=", " ", RowBox[{ RowBox[{ RowBox[{"-", " ", "c"}], " ", "p", " ", RowBox[{"h", "/", RowBox[{"(", RowBox[{"1", "+", " ", RowBox[{"c", " ", "h"}]}], ")"}]}]}], "-", RowBox[{"\[Beta]", " ", "c", " ", "\[Eta]", " ", "p", " ", RowBox[{"h", "/", RowBox[{"(", RowBox[{"1", "+", " ", RowBox[{"c", " ", "\[Eta]", " ", "h"}]}], ")"}]}]}], " ", "+", " ", RowBox[{"h", RowBox[{"(", RowBox[{"1", " ", "-", " ", "h"}], ")"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"w", "[", "h_", "]"}], ":=", " ", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "h"}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"c", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"c", " ", "h"}]}], ")"}]}], " ", "+", " ", RowBox[{"\[Beta]", " ", "c", " ", RowBox[{"\[Eta]", "/", RowBox[{"(", RowBox[{"1", " ", "+", " ", RowBox[{"c", " ", "\[Eta]", " ", "h"}]}], ")"}]}]}]}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"fullsol", " ", "=", " ", RowBox[{"FullSimplify", "[", RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], " ", "\[Equal]", " ", "0"}], " ", "&&", " ", RowBox[{ RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}], " ", "\[Equal]", " ", "0"}]}], ",", " ", RowBox[{"{", RowBox[{"h", ",", "p"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"coexistence", " ", "steady", " ", "state"}], "*)"}], "\[IndentingNewLine]", RowBox[{"hs", " ", "=", " ", RowBox[{"fullsol", "[", RowBox[{"[", RowBox[{"3", ",", "1", ",", "2"}], "]"}], "]"}]}], "\[IndentingNewLine]", RowBox[{"ps", " ", "=", " ", RowBox[{"fullsol", "[", RowBox[{"[", RowBox[{"3", ",", "2", ",", "2"}], "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.81417891404978*^9, 3.8141789143097677`*^9}, 3.816256573055167*^9, {3.8162599972433023`*^9, 3.816259997708579*^9}, { 3.816263781663917*^9, 3.8162637819055634`*^9}, {3.816269141377301*^9, 3.816269141997353*^9}, {3.8453812271665688`*^9, 3.845381264945526*^9}, { 3.8677486445031347`*^9, 3.8677486488376713`*^9}, {3.867748793861174*^9, 3.867748814926955*^9}, {3.8677494337739353`*^9, 3.867749563904695*^9}, { 3.867749669543584*^9, 3.867749705930317*^9}, {3.867749739983571*^9, 3.867749755681912*^9}, {3.870170912759418*^9, 3.870170932642469*^9}, 3.8701710166117563`*^9, {3.870338097458365*^9, 3.870338116262569*^9}, { 3.870339865705101*^9, 3.870339865816143*^9}, {3.870517760828545*^9, 3.870517817597672*^9}},ExpressionUUID->"44499988-2afd-4ab3-bec9-\ 71c5368b14bd"], Cell[BoxData[ FractionBox["2", RowBox[{"c", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "E1", "-", "\[Eta]", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}]]], "Output", CellChangeTimes->{ 3.816260007673848*^9, 3.816263786465261*^9, 3.816269157561713*^9, 3.816620923614305*^9, 3.819990395149886*^9, 3.819990809910742*^9, 3.8199946792954206`*^9, 3.819995280068986*^9, 3.819995683395197*^9, 3.8199961827248287`*^9, 3.819996729713945*^9, 3.8199998318620567`*^9, 3.820004741691141*^9, 3.820007594656899*^9, 3.820007777358328*^9, 3.820007939913474*^9, 3.820008092843512*^9, 3.820008219914859*^9, 3.820011056558919*^9, 3.820014223345696*^9, 3.8200167981356363`*^9, 3.820055027956944*^9, 3.820057775097437*^9, 3.820060691073681*^9, 3.820065181693656*^9, 3.820068311837648*^9, 3.8200785327777042`*^9, 3.820082788383857*^9, 3.8200829424604797`*^9, 3.8200867924535427`*^9, 3.82008682426169*^9, 3.820090470148801*^9, 3.8200943528089046`*^9, { 3.820231892267478*^9, 3.820231906539508*^9}, 3.8202331290342293`*^9, 3.8202334365161*^9, 3.820237702435618*^9, 3.82024313463249*^9, 3.820247691595439*^9, 3.8206676355097313`*^9, 3.820667766657905*^9, 3.820683546949781*^9, 3.820683937586443*^9, 3.820684270768063*^9, 3.8206843389322357`*^9, 3.820684441623336*^9, 3.820687759119823*^9, 3.820689815695136*^9, 3.820693050107286*^9, 3.820696790485519*^9, 3.820699123347419*^9, 3.820702324395194*^9, 3.8207046551083097`*^9, 3.8207454765532627`*^9, 3.8207456125370827`*^9, 3.8207522933832073`*^9, 3.820753767563983*^9, 3.820924145052561*^9, 3.820924840493137*^9, 3.8209279594886217`*^9, 3.820934785373683*^9, 3.82093681473761*^9, 3.820937042309434*^9, 3.8209403541566677`*^9, 3.821205686841371*^9, 3.821206263316573*^9, 3.821207191607079*^9, 3.821207340621151*^9, 3.821207476857024*^9, 3.821211821023712*^9, 3.821211997828783*^9, 3.821274104968186*^9, 3.821276185891083*^9, 3.821295982223572*^9, { 3.821296225054063*^9, 3.8212962379156322`*^9}, 3.821296413362687*^9, 3.821303895664366*^9, 3.821349858640416*^9, 3.8213517997944307`*^9, 3.821362282576243*^9, 3.821363411658312*^9, 3.8214406287364197`*^9, 3.822395405097415*^9, 3.822395636715575*^9, 3.822405577659626*^9, 3.8225697940579453`*^9, 3.8225699735076113`*^9, 3.822570774146659*^9, 3.822578204074615*^9, 3.822581457769534*^9, 3.822581571348942*^9, 3.822581760629026*^9, 3.8225818897678823`*^9, 3.822582092106943*^9, 3.8225821912707043`*^9, 3.82258307771313*^9, 3.82258451501829*^9, 3.822585778995179*^9, 3.822587794670021*^9, 3.822589254385124*^9, { 3.822589560928432*^9, 3.822589589107071*^9}, 3.822589651892469*^9, 3.822589716730098*^9, 3.822589781441945*^9, 3.822590065251359*^9, 3.822591962132984*^9, 3.82319899928161*^9, 3.823200382408862*^9, 3.823292329220593*^9, 3.823292414137678*^9, 3.823296016216352*^9, 3.823297391816389*^9, 3.823555258607255*^9, 3.823557082278533*^9, 3.823558486357008*^9, 3.8236002485928593`*^9, 3.823600318167491*^9, 3.824298143255495*^9, 3.82429820695124*^9, 3.841324870614168*^9, 3.8413250635928917`*^9, 3.841326150889401*^9, 3.84155995108565*^9, 3.842334300719336*^9, 3.842334333515789*^9, 3.842336800012926*^9, 3.842339435044715*^9, {3.8423394973027353`*^9, 3.842339516068049*^9}, 3.84303282058718*^9, 3.8430342615753317`*^9, 3.843035936495295*^9, 3.8430434933237*^9, 3.8436355840781603`*^9, 3.843635640180191*^9, 3.843636466680826*^9, 3.8436381482155857`*^9, 3.844171863357952*^9, 3.844176215910716*^9, 3.844177025157576*^9, 3.84417908325461*^9, 3.8442374499387836`*^9, 3.844239977229495*^9, 3.844258182375104*^9, 3.8444959709745293`*^9, 3.844860375306011*^9, 3.845381281407596*^9, 3.845389195422049*^9, 3.8466851742182903`*^9, 3.846685327975951*^9, 3.8487571795391903`*^9, 3.851178147241667*^9, 3.8679199803860407`*^9, 3.8679234561146994`*^9, 3.870171034333346*^9, 3.870172164471468*^9, 3.870176237111171*^9, 3.870179396705922*^9, 3.8702459047620773`*^9, 3.870338740452734*^9, 3.870519803830413*^9, 3.8705199969808083`*^9, 3.870522178053238*^9, 3.870522755453874*^9, 3.8705243400722427`*^9, 3.8707750261843853`*^9, 3.8710462135374126`*^9},ExpressionUUID->"4c270777-89bc-4627-805d-\ 42019592b15e"], Cell[BoxData[ FractionBox[ RowBox[{"2", " ", "E1"}], RowBox[{"2", "+", RowBox[{"c", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "E1", "-", "\[Eta]", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", RowBox[{"c", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], "+", FractionBox[ RowBox[{"c", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"c", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "E1", "-", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "c"}], ")"}], " ", "\[Eta]"}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"c", "+", "\[Beta]", "+", RowBox[{"c", " ", "\[Beta]"}]}], ")"}], " ", "\[Eta]"}]}], ")"}]}]}]], "+", FractionBox[ RowBox[{"2", "+", RowBox[{"c", " ", RowBox[{"(", RowBox[{"1", "+", "\[Eta]", "-", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}]}], ")"}]}]}], RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"c", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "E1", "-", RowBox[{ RowBox[{"(", RowBox[{"1", "+", "c"}], ")"}], " ", "\[Eta]"}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"c", "+", "\[Beta]", "+", RowBox[{"c", " ", "\[Beta]"}]}], ")"}], " ", "\[Eta]"}]}], ")"}]}]}]]}]]], "Output", CellChangeTimes->{ 3.816260007673848*^9, 3.816263786465261*^9, 3.816269157561713*^9, 3.816620923614305*^9, 3.819990395149886*^9, 3.819990809910742*^9, 3.8199946792954206`*^9, 3.819995280068986*^9, 3.819995683395197*^9, 3.8199961827248287`*^9, 3.819996729713945*^9, 3.8199998318620567`*^9, 3.820004741691141*^9, 3.820007594656899*^9, 3.820007777358328*^9, 3.820007939913474*^9, 3.820008092843512*^9, 3.820008219914859*^9, 3.820011056558919*^9, 3.820014223345696*^9, 3.8200167981356363`*^9, 3.820055027956944*^9, 3.820057775097437*^9, 3.820060691073681*^9, 3.820065181693656*^9, 3.820068311837648*^9, 3.8200785327777042`*^9, 3.820082788383857*^9, 3.8200829424604797`*^9, 3.8200867924535427`*^9, 3.82008682426169*^9, 3.820090470148801*^9, 3.8200943528089046`*^9, { 3.820231892267478*^9, 3.820231906539508*^9}, 3.8202331290342293`*^9, 3.8202334365161*^9, 3.820237702435618*^9, 3.82024313463249*^9, 3.820247691595439*^9, 3.8206676355097313`*^9, 3.820667766657905*^9, 3.820683546949781*^9, 3.820683937586443*^9, 3.820684270768063*^9, 3.8206843389322357`*^9, 3.820684441623336*^9, 3.820687759119823*^9, 3.820689815695136*^9, 3.820693050107286*^9, 3.820696790485519*^9, 3.820699123347419*^9, 3.820702324395194*^9, 3.8207046551083097`*^9, 3.8207454765532627`*^9, 3.8207456125370827`*^9, 3.8207522933832073`*^9, 3.820753767563983*^9, 3.820924145052561*^9, 3.820924840493137*^9, 3.8209279594886217`*^9, 3.820934785373683*^9, 3.82093681473761*^9, 3.820937042309434*^9, 3.8209403541566677`*^9, 3.821205686841371*^9, 3.821206263316573*^9, 3.821207191607079*^9, 3.821207340621151*^9, 3.821207476857024*^9, 3.821211821023712*^9, 3.821211997828783*^9, 3.821274104968186*^9, 3.821276185891083*^9, 3.821295982223572*^9, { 3.821296225054063*^9, 3.8212962379156322`*^9}, 3.821296413362687*^9, 3.821303895664366*^9, 3.821349858640416*^9, 3.8213517997944307`*^9, 3.821362282576243*^9, 3.821363411658312*^9, 3.8214406287364197`*^9, 3.822395405097415*^9, 3.822395636715575*^9, 3.822405577659626*^9, 3.8225697940579453`*^9, 3.8225699735076113`*^9, 3.822570774146659*^9, 3.822578204074615*^9, 3.822581457769534*^9, 3.822581571348942*^9, 3.822581760629026*^9, 3.8225818897678823`*^9, 3.822582092106943*^9, 3.8225821912707043`*^9, 3.82258307771313*^9, 3.82258451501829*^9, 3.822585778995179*^9, 3.822587794670021*^9, 3.822589254385124*^9, { 3.822589560928432*^9, 3.822589589107071*^9}, 3.822589651892469*^9, 3.822589716730098*^9, 3.822589781441945*^9, 3.822590065251359*^9, 3.822591962132984*^9, 3.82319899928161*^9, 3.823200382408862*^9, 3.823292329220593*^9, 3.823292414137678*^9, 3.823296016216352*^9, 3.823297391816389*^9, 3.823555258607255*^9, 3.823557082278533*^9, 3.823558486357008*^9, 3.8236002485928593`*^9, 3.823600318167491*^9, 3.824298143255495*^9, 3.82429820695124*^9, 3.841324870614168*^9, 3.8413250635928917`*^9, 3.841326150889401*^9, 3.84155995108565*^9, 3.842334300719336*^9, 3.842334333515789*^9, 3.842336800012926*^9, 3.842339435044715*^9, {3.8423394973027353`*^9, 3.842339516068049*^9}, 3.84303282058718*^9, 3.8430342615753317`*^9, 3.843035936495295*^9, 3.8430434933237*^9, 3.8436355840781603`*^9, 3.843635640180191*^9, 3.843636466680826*^9, 3.8436381482155857`*^9, 3.844171863357952*^9, 3.844176215910716*^9, 3.844177025157576*^9, 3.84417908325461*^9, 3.8442374499387836`*^9, 3.844239977229495*^9, 3.844258182375104*^9, 3.8444959709745293`*^9, 3.844860375306011*^9, 3.845381281407596*^9, 3.845389195422049*^9, 3.8466851742182903`*^9, 3.846685327975951*^9, 3.8487571795391903`*^9, 3.851178147241667*^9, 3.8679199803860407`*^9, 3.8679234561146994`*^9, 3.870171034333346*^9, 3.870172164471468*^9, 3.870176237111171*^9, 3.870179396705922*^9, 3.8702459047620773`*^9, 3.870338740452734*^9, 3.870519803830413*^9, 3.8705199969808083`*^9, 3.870522178053238*^9, 3.870522755453874*^9, 3.8705243400722427`*^9, 3.8707750261843853`*^9, 3.871046213539187*^9},ExpressionUUID->"f3a2ff0b-d3d3-44b9-a6ad-\ c625cd4ac009"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Section 1: Compute the jacobian and C_crit", "Section", CellChangeTimes->{{3.7988347922346067`*^9, 3.798834797097238*^9}, { 3.867749949175293*^9, 3.867749949414136*^9}, {3.870170923268626*^9, 3.8701709234331903`*^9}, {3.870517841927836*^9, 3.870517861015704*^9}, { 3.870518104782816*^9, 3.870518104929731*^9}, 3.8705182872660513`*^9, { 3.870769998305819*^9, 3.870769999905858*^9}, {3.870770202865884*^9, 3.870770219372768*^9}},ExpressionUUID->"9a1e7251-ec19-4baa-86eb-\ cecbf78ffcc2"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ClearAll", "[", RowBox[{ "\[Beta]", ",", "J", ",", "c", ",", "\[Eta]", ",", "B1", ",", "B2", ",", "E1", ",", "\[Mu]", ",", "trace", ",", "det", ",", "r", ",", "i", ",", "\[CapitalLambda]0", ",", "\[CapitalOmega]0", ",", "\[CapitalOmega]1", ",", "p", ",", "h", ",", " ", "ccrit1", ",", "ccrit2", ",", "\[Xi]1", ",", "\[Xi]2"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"J", " ", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "p"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "h"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "p"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "h"}], "]"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"trace", " ", "=", " ", RowBox[{ RowBox[{"Tr", "[", "J", "]"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"h", " ", "\[Rule]", " ", "hs"}], ",", RowBox[{"p", "\[Rule]", " ", "ps"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"det", " ", "=", " ", RowBox[{ RowBox[{"Det", "[", "J", "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"h", " ", "\[Rule]", " ", "hs"}], ",", RowBox[{"p", "\[Rule]", " ", "ps"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ccrit1", " ", "=", " ", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{"trace", "\[Equal]", "0"}], ",", "c"}], "]"}], "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"equivalent", " ", "way", " ", "to", " ", "solve", " ", RowBox[{"Tr", "[", "J", "]"}]}], " ", "=", " ", "0"}], ",", " ", RowBox[{ RowBox[{"but", " ", "more", " ", "efficient"}], " ", "-", " ", RowBox[{ "See", " ", "stability", " ", "analysis", " ", "in", " ", "Cobbold", " ", "et", " ", "al", " ", RowBox[{"(", "2022", ")"}], " ", "for", " ", "details"}]}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{"ccrit2", " ", "=", RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{"Solve", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"w", "'"}], "[", "hs", "]"}], "\[Equal]", "0"}], ",", "c"}], "]"}], "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Simplify", "[", RowBox[{"ccrit2", "/.", " ", RowBox[{"\[Eta]", " ", "\[Rule]", " ", "1"}]}], "]"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"FortranForm", "[", RowBox[{"ccrit2", " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"\[Eta]", "\[Rule]", " ", "eta"}], ",", " ", RowBox[{"\[Beta]", " ", "\[Rule]", " ", "beta"}]}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.7987342234154997`*^9, 3.798734231613965*^9}, { 3.798734261816148*^9, 3.798734377619363*^9}, {3.798734428142951*^9, 3.798734688256218*^9}, {3.798734785033711*^9, 3.798734815958508*^9}, { 3.798734898870325*^9, 3.798734932680476*^9}, {3.798735055896419*^9, 3.7987350807379427`*^9}, {3.798735379767914*^9, 3.798735409429327*^9}, { 3.7987356794800463`*^9, 3.798735724674944*^9}, {3.798735774777719*^9, 3.7987359603264503`*^9}, {3.798736095842292*^9, 3.7987361374053802`*^9}, { 3.798736849448769*^9, 3.7987368609479237`*^9}, {3.798736916737837*^9, 3.7987369189394913`*^9}, {3.798737883546331*^9, 3.798737943045732*^9}, { 3.798737983267951*^9, 3.7987380417562313`*^9}, {3.798738107227046*^9, 3.798738134723774*^9}, {3.798738313205209*^9, 3.7987383348834677`*^9}, { 3.79873936455468*^9, 3.798739392178103*^9}, {3.79873953774231*^9, 3.798739597158895*^9}, {3.798739693609086*^9, 3.7987397673084183`*^9}, { 3.798740011522107*^9, 3.798740011718122*^9}, {3.798740192116337*^9, 3.798740199511115*^9}, {3.798740276535761*^9, 3.798740285584331*^9}, { 3.798740498459552*^9, 3.7987405015225487`*^9}, {3.798740662441784*^9, 3.798740681889448*^9}, {3.798740801073596*^9, 3.7987408011624126`*^9}, { 3.798741549973358*^9, 3.7987415791581783`*^9}, {3.798741711383276*^9, 3.798741764741913*^9}, {3.798741964503049*^9, 3.7987419937323847`*^9}, { 3.798742026373466*^9, 3.7987420288416*^9}, {3.798742103384511*^9, 3.798742103669217*^9}, {3.798742373271051*^9, 3.7987423862248583`*^9}, { 3.798743334370755*^9, 3.7987433490221767`*^9}, {3.798743385934226*^9, 3.798743569443079*^9}, {3.798743601993536*^9, 3.798743635056342*^9}, { 3.798743679253008*^9, 3.7987436859660683`*^9}, {3.798743733982664*^9, 3.7987437363037663`*^9}, {3.798743781503202*^9, 3.79874383437232*^9}, { 3.7987485628756723`*^9, 3.798748563498851*^9}, {3.798748595418427*^9, 3.798748677190031*^9}, 3.798748715901243*^9, 3.798748949627234*^9, { 3.798819247413789*^9, 3.798819282015173*^9}, {3.798819330066749*^9, 3.798819330223613*^9}, {3.798819420978258*^9, 3.798819440256166*^9}, { 3.7988194953892593`*^9, 3.798819621349909*^9}, {3.7988196581787567`*^9, 3.798819747880352*^9}, {3.7988206212259073`*^9, 3.798820629960053*^9}, { 3.798821922445283*^9, 3.798821989651827*^9}, {3.7988226092308493`*^9, 3.798822640219384*^9}, 3.798834812659114*^9, {3.799018009180895*^9, 3.799018097714417*^9}, {3.799018133337532*^9, 3.79901813365396*^9}, { 3.799018164712735*^9, 3.799018184643711*^9}, 3.799018404078745*^9, 3.799095147055685*^9, {3.7990959816296043`*^9, 3.799095999046344*^9}, 3.79935294656887*^9, {3.799353005230137*^9, 3.7993531653261538`*^9}, { 3.799353226283279*^9, 3.79935332453985*^9}, {3.799353368971451*^9, 3.799353387529275*^9}, 3.7993534514079*^9, {3.799427218647582*^9, 3.7994272233122063`*^9}, {3.7995152202878923`*^9, 3.7995152824914618`*^9}, {3.799515317934978*^9, 3.799515426110711*^9}, { 3.79952160096795*^9, 3.7995216275478153`*^9}, {3.799526075170764*^9, 3.799526088662923*^9}, {3.799526447468687*^9, 3.799526460354454*^9}, { 3.799579501807143*^9, 3.799579574687928*^9}, {3.799609254763646*^9, 3.799609274873742*^9}, {3.801504308248908*^9, 3.8015043216274137`*^9}, { 3.801516281122512*^9, 3.80151632308362*^9}, {3.8045249250037603`*^9, 3.804524977548204*^9}, {3.8045250494701443`*^9, 3.804525071616247*^9}, 3.804526170254887*^9, {3.8045263425450907`*^9, 3.804526358049371*^9}, { 3.804608753519327*^9, 3.804608798348752*^9}, {3.8046153486403637`*^9, 3.804615350145711*^9}, {3.8046154901122923`*^9, 3.804615492992744*^9}, { 3.8046155818387203`*^9, 3.804615609453364*^9}, {3.804615642131398*^9, 3.804615651262536*^9}, {3.81160608964854*^9, 3.811606098914525*^9}, { 3.8116061658769617`*^9, 3.811606195832057*^9}, {3.811606250173057*^9, 3.811606270060581*^9}, {3.811606364734374*^9, 3.811606392145708*^9}, { 3.8116637841275473`*^9, 3.811663791219623*^9}, {3.811663949841942*^9, 3.81166395117531*^9}, {3.811664015681662*^9, 3.811664056260593*^9}, { 3.8116846500229883`*^9, 3.811684692212419*^9}, {3.8116986099636908`*^9, 3.811698657357214*^9}, {3.811698776057794*^9, 3.8116987804587193`*^9}, 3.8122104788933983`*^9, {3.812296732964787*^9, 3.812296788104147*^9}, { 3.8127307975248013`*^9, 3.8127308236870327`*^9}, {3.8127308701054373`*^9, 3.812730880804455*^9}, {3.812730941703039*^9, 3.812730978434101*^9}, { 3.812732674242885*^9, 3.812732722551504*^9}, {3.812735047826624*^9, 3.812735058190894*^9}, {3.812901367908091*^9, 3.812901371261791*^9}, { 3.8129016591200867`*^9, 3.812901659460806*^9}, {3.812902248599148*^9, 3.812902249696342*^9}, 3.812902344596095*^9, {3.812902375145059*^9, 3.812902376068527*^9}, {3.812902517211801*^9, 3.812902523526093*^9}, { 3.812902799940284*^9, 3.812902802116323*^9}, {3.812903014570705*^9, 3.812903014866517*^9}, 3.8129034487084723`*^9, {3.813158531031769*^9, 3.813158538294826*^9}, {3.81315861866444*^9, 3.813158634454507*^9}, { 3.8132471380545397`*^9, 3.8132471642299013`*^9}, {3.8132474987587843`*^9, 3.813247500232027*^9}, {3.813248515616477*^9, 3.813248541322818*^9}, { 3.813335951812253*^9, 3.813336019632717*^9}, {3.8133361354400663`*^9, 3.813336184654544*^9}, {3.813336260786762*^9, 3.813336266680437*^9}, { 3.813502776963819*^9, 3.8135028212315073`*^9}, {3.813504385752734*^9, 3.81350438959519*^9}, {3.813504529732304*^9, 3.813504532731739*^9}, { 3.8135047245804873`*^9, 3.813504724924851*^9}, {3.813505312532651*^9, 3.813505325823781*^9}, {3.813507499935174*^9, 3.813507636946794*^9}, { 3.8135076764114637`*^9, 3.813507740562434*^9}, {3.813507777159842*^9, 3.813507806639048*^9}, {3.813507843145596*^9, 3.813507845680951*^9}, { 3.81350789629276*^9, 3.81350790674308*^9}, {3.813508001862009*^9, 3.813508096396082*^9}, {3.813508133864497*^9, 3.8135081578117447`*^9}, { 3.813508191802782*^9, 3.813508222410645*^9}, {3.813509678957717*^9, 3.813509680305846*^9}, {3.813678604240458*^9, 3.81367860456435*^9}, { 3.813678739891315*^9, 3.8136787757554483`*^9}, {3.813679976254689*^9, 3.813680028917473*^9}, {3.8136800929169607`*^9, 3.813680150909163*^9}, { 3.813680202229609*^9, 3.8136802079190407`*^9}, {3.813680330760276*^9, 3.813680367461582*^9}, {3.81368044653712*^9, 3.813680449435545*^9}, { 3.813680513403907*^9, 3.813680519310547*^9}, {3.813680564452283*^9, 3.813680591864854*^9}, {3.8137639723874702`*^9, 3.81376397387121*^9}, { 3.814111295960573*^9, 3.814111302157703*^9}, {3.815475781986075*^9, 3.815475822207179*^9}, {3.820924139986253*^9, 3.820924140139441*^9}, 3.844189684238967*^9, {3.8441897145010033`*^9, 3.8441897146569443`*^9}, { 3.844237500708506*^9, 3.844237509170483*^9}, {3.845390073479212*^9, 3.84539007376473*^9}, {3.8453901118380623`*^9, 3.8453902072030687`*^9}, { 3.845390735798182*^9, 3.845390790813398*^9}, {3.846685300901208*^9, 3.846685361781104*^9}, {3.851178142001629*^9, 3.851178159674836*^9}, { 3.851178448257082*^9, 3.851178450426318*^9}, {3.851178617957581*^9, 3.851178640559816*^9}, 3.851178690787167*^9, 3.867748672635831*^9, { 3.8677487611953773`*^9, 3.867748761505496*^9}, {3.8677497642159967`*^9, 3.8677497855044622`*^9}, {3.870170965372242*^9, 3.8701710109652147`*^9}, { 3.8701710980175753`*^9, 3.8701711013918447`*^9}, {3.8701720859322863`*^9, 3.870172089535202*^9}, 3.870172908113249*^9, {3.870179547193363*^9, 3.870179550851925*^9}, {3.870245724219845*^9, 3.870245814105981*^9}, { 3.870247035731124*^9, 3.870247067598001*^9}, {3.870338425664592*^9, 3.870338431120647*^9}, 3.870517849377664*^9, {3.870519760601693*^9, 3.870519780140018*^9}, {3.870519819406116*^9, 3.870519822195179*^9}, 3.870519898460595*^9},ExpressionUUID->"c218540f-34cc-4e4f-bfe4-\ 3d27be67ba86"], Cell[BoxData[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["E1", "3"], " ", SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}], "4"]}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]", "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]", "-", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}]}], ")"}]}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Eta]", " ", RowBox[{"(", RowBox[{"\[Beta]", "+", "\[Eta]", "-", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "2"]}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Beta]", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}]}], ")"}]}]}], ")"}]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Eta]", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"\[Beta]", " ", RowBox[{"(", RowBox[{"2", "-", "\[Eta]", "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "3"], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "2"]}]}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]", "+", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "4"]}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "-", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "3"], " ", RowBox[{"(", RowBox[{"1", "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}]}], ")"}]}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"3", "-", RowBox[{"2", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], "+", RowBox[{"\[Eta]", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "4"}], "+", "\[Eta]", "+", RowBox[{"3", " ", SuperscriptBox["\[Beta]", "2"], " ", SuperscriptBox["\[Eta]", "4"]}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "-", RowBox[{"\[Beta]", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], "-", RowBox[{"2", " ", SuperscriptBox["\[Beta]", "2"], " ", SuperscriptBox["\[Eta]", "3"], " ", RowBox[{"(", RowBox[{"2", "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}], "+", RowBox[{"\[Beta]", " ", SuperscriptBox["\[Eta]", "2"], " ", RowBox[{"(", RowBox[{"\[Beta]", "-", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Beta]", " ", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}]}], ")"}]}]}], ")"}]}]}], ")"}]}], "+", RowBox[{ SuperscriptBox["E1", "2"], " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "3"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Eta]", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"\[Beta]", " ", RowBox[{"(", RowBox[{ SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Eta]", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "2"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]], "+", RowBox[{"\[Beta]", " ", "\[Eta]", " ", RowBox[{"(", RowBox[{"2", "-", RowBox[{"3", " ", "\[Eta]"}], "+", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], "2"], "-", RowBox[{"2", " ", "E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", "\[Eta]"}], ")"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", RowBox[{"\[Beta]", " ", "\[Eta]"}]}], ")"}]}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"E1", "+", RowBox[{"E1", " ", "\[Beta]", " ", "\[Eta]"}]}], ")"}], "2"]}]]}], ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]}], ")"}]}]], "Output", CellChangeTimes->{ 3.846685347542964*^9, 3.848757185270575*^9, 3.851178456780199*^9, 3.8511786471076717`*^9, 3.8511786950986557`*^9, 3.867919984572547*^9, 3.8679234602739735`*^9, 3.870171078548855*^9, 3.8701720672355223`*^9, 3.870172207642228*^9, 3.870176326700933*^9, 3.870179660652628*^9, 3.870245928388256*^9, 3.870338754179783*^9, {3.870519816761382*^9, 3.8705198289742928`*^9}, 3.870519903413982*^9, 3.870520010879369*^9, 3.8705221930197897`*^9, 3.870522280284779*^9, 3.870522771336607*^9, 3.87052436379018*^9, 3.870775040534955*^9, 3.8710462289111*^9},ExpressionUUID->"e230c138-949c-4822-a4a0-95eb5f08a269"], Cell[BoxData[ FractionBox[ RowBox[{"1", "+", "E1", "+", RowBox[{"E1", " ", "\[Beta]"}]}], RowBox[{ RowBox[{"-", "1"}], "+", "E1", "+", RowBox[{"E1", " ", "\[Beta]"}]}]]], "Output", CellChangeTimes->{ 3.846685347542964*^9, 3.848757185270575*^9, 3.851178456780199*^9, 3.8511786471076717`*^9, 3.8511786950986557`*^9, 3.867919984572547*^9, 3.8679234602739735`*^9, 3.870171078548855*^9, 3.8701720672355223`*^9, 3.870172207642228*^9, 3.870176326700933*^9, 3.870179660652628*^9, 3.870245928388256*^9, 3.870338754179783*^9, {3.870519816761382*^9, 3.8705198289742928`*^9}, 3.870519903413982*^9, 3.870520010879369*^9, 3.8705221930197897`*^9, 3.870522280284779*^9, 3.870522771336607*^9, 3.87052436379018*^9, 3.870775040534955*^9, 3.871046228917008*^9},ExpressionUUID->"194cbf59-492e-415c-9c27-\ 26abcbcad0f7"] }, Open ]], Cell[BoxData[{ RowBox[{"ClearAll", "[", RowBox[{ "\[Beta]", ",", "J", ",", "c", ",", "\[Eta]", ",", "B1", ",", "B2", ",", "E1", ",", "r", ",", "i", ",", "\[CapitalLambda]0", ",", "\[CapitalOmega]0", ",", "\[CapitalOmega]1", ",", "p", ",", "h", ",", " ", "ccrit"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ccrit", " ", "=", " ", "ccrit2"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Relevant", " ", "parameters", " ", "for", " ", "the", " ", "normal", " ", "form"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", " ", "=", " ", RowBox[{"trace", "/", "2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"i", " ", "=", RowBox[{"Sqrt", "[", RowBox[{"det", "-", RowBox[{"trace", "*", RowBox[{"trace", "/", "4"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[CapitalLambda]0", " ", "=", " ", RowBox[{ RowBox[{"D", "[", RowBox[{"r", ",", "c"}], "]"}], " ", "/.", " ", RowBox[{"c", " ", "\[Rule]", " ", "ccrit"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[CapitalOmega]0", " ", "=", " ", RowBox[{"i", " ", "/.", " ", RowBox[{"c", " ", "\[Rule]", " ", "ccrit"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[CapitalOmega]1", " ", "=", " ", RowBox[{ RowBox[{"D", "[", RowBox[{"i", ",", "c"}], "]"}], " ", "/.", " ", RowBox[{"c", " ", "\[Rule]", " ", "ccrit"}]}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.813680758001863*^9, 3.8136807586853523`*^9}, { 3.813680808101639*^9, 3.813680820183529*^9}, 3.81376353706448*^9, { 3.8137639843127937`*^9, 3.813763985370322*^9}, {3.814193648924788*^9, 3.814193657343454*^9}, {3.8154758256004953`*^9, 3.8154758468677588`*^9}, { 3.8209282464851637`*^9, 3.820928260105897*^9}, {3.845381283680087*^9, 3.845381385632813*^9}, {3.845389189623674*^9, 3.8453892238188334`*^9}, { 3.846685367158064*^9, 3.8466854225054703`*^9}, {3.8677497980082808`*^9, 3.867749810296505*^9}, {3.867749864839696*^9, 3.8677498745217953`*^9}, { 3.8701720964640713`*^9, 3.870172096623859*^9}, {3.870176198623065*^9, 3.870176198703475*^9}, {3.8702458986477203`*^9, 3.8702458988814297`*^9}, { 3.8703384573814917`*^9, 3.8703384590143147`*^9}},ExpressionUUID->"287e0f4c-0953-4b8b-810f-\ 69f8de8594ad"] }, Open ]], Cell[CellGroupData[{ Cell["Set the fixed parameters", "Subchapter", CellChangeTimes->{{3.81411125487113*^9, 3.814111264266161*^9}, 3.8206842028251*^9, {3.870517905850609*^9, 3.870517909616067*^9}, { 3.870518114102515*^9, 3.870518114881226*^9}},ExpressionUUID->"6ff83796-23ad-4212-8639-\ 61206155d192"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ClearAll", "[", RowBox[{"\[Beta]", ",", "E1", ",", "B1", ",", "\[Eta]"}], " ", "]"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Fixed", " ", "paramaters"}], "*)"}]}], "\[IndentingNewLine]", "\"\<\[Beta] = \>\"", \ "\[IndentingNewLine]", RowBox[{"\[Beta]", " ", "=", "1.0"}], "\[IndentingNewLine]", "\"\<\[Eta] = \>\"", "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Eta]", "=", "7.0"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "Tests", " ", "and", " ", "constraints", " ", "for", " ", "existence", " ", "of", " ", "the", " ", "Hopf"}], "-", "bifurcation", " ", "-", " ", RowBox[{"See", " ", "Cobbold", " ", "et", " ", RowBox[{"al", ".", " ", RowBox[{"(", "2022", ")"}]}]}]}], "*)"}]}], "\[IndentingNewLine]", "\"\\"", "\ \[IndentingNewLine]", RowBox[{ RowBox[{"1", "/", RowBox[{"(", RowBox[{"1", " ", "+", " ", "\[Beta]"}], ")"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Assuming", "[", RowBox[{ RowBox[{"E1", ">", "0"}], ",", RowBox[{"Simplify", "[", "ccrit", "]"}]}], "]"}], "\[IndentingNewLine]", "\"\\"", "\ \[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[Beta]", " ", "\[Eta]"}], " ", "+", " ", "1"}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"\[Beta]", " ", RowBox[{"\[Eta]", "^", "2"}]}], " ", "+", "1"}], ")"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.813847853068179*^9, 3.81384788416678*^9}, 3.813848089010256*^9, {3.813848929843068*^9, 3.8138489375473127`*^9}, { 3.813849116691184*^9, 3.813849229963903*^9}, {3.813849274557877*^9, 3.813849404011653*^9}, {3.813849493841926*^9, 3.8138495732796164`*^9}, { 3.81384961757297*^9, 3.813849665649061*^9}, {3.8138496969413424`*^9, 3.813849844826907*^9}, {3.813849893647382*^9, 3.813849938407579*^9}, { 3.813849971011075*^9, 3.813850020658668*^9}, {3.813850094185698*^9, 3.8138501068948717`*^9}, {3.813850193559969*^9, 3.813850193696293*^9}, { 3.813850228651712*^9, 3.8138502316161213`*^9}, {3.81385035557405*^9, 3.81385037554002*^9}, {3.8138543319895573`*^9, 3.813854382740534*^9}, { 3.813854447122476*^9, 3.813854632428308*^9}, {3.813854664099441*^9, 3.813854836585938*^9}, {3.8138552857225533`*^9, 3.813855292585596*^9}, { 3.81385535459562*^9, 3.81385539408436*^9}, {3.813855444241312*^9, 3.813855453732688*^9}, {3.813855590036703*^9, 3.813855602050605*^9}, { 3.8138557566693077`*^9, 3.813855766776437*^9}, {3.813855815078727*^9, 3.813855883970489*^9}, {3.8138627457136393`*^9, 3.813862785736258*^9}, { 3.8138670675825033`*^9, 3.8138670779395523`*^9}, {3.813867140738015*^9, 3.8138671690339117`*^9}, {3.813917147611693*^9, 3.813917148679953*^9}, { 3.8139172374457073`*^9, 3.813917246639744*^9}, {3.813917368606563*^9, 3.813917371032949*^9}, {3.813917567966646*^9, 3.813917656021612*^9}, { 3.813917860400943*^9, 3.813917906387423*^9}, {3.814016062176503*^9, 3.814016072368937*^9}, {3.8140167824211817`*^9, 3.814016910787592*^9}, { 3.8140170323411083`*^9, 3.814017050196142*^9}, {3.8140171577558203`*^9, 3.81401721457301*^9}, {3.814017283937005*^9, 3.814017295991026*^9}, { 3.81401732735608*^9, 3.814017328365224*^9}, {3.814018541955846*^9, 3.8140185447230873`*^9}, {3.81401858350347*^9, 3.81401858464624*^9}, { 3.814018794041625*^9, 3.814018841429224*^9}, {3.814019088684404*^9, 3.8140191804739532`*^9}, {3.814019313804188*^9, 3.814019314224742*^9}, { 3.81401948942514*^9, 3.8140194971850033`*^9}, {3.8140196711865664`*^9, 3.814019692289475*^9}, {3.814019833930777*^9, 3.8140198370236692`*^9}, { 3.8140251460453587`*^9, 3.814025226361556*^9}, {3.814025277880512*^9, 3.8140253065867863`*^9}, {3.81402547669005*^9, 3.81402548415982*^9}, { 3.8140255672318697`*^9, 3.814025613192268*^9}, {3.814025699721527*^9, 3.814025733384205*^9}, {3.814025779880883*^9, 3.814025833746921*^9}, { 3.8140260851726313`*^9, 3.8140261466396112`*^9}, {3.8140261925051193`*^9, 3.814026256449338*^9}, {3.8140263004909153`*^9, 3.814026302738495*^9}, { 3.8140276953855753`*^9, 3.814027747090849*^9}, {3.814089238937214*^9, 3.814089245845149*^9}, {3.814089332706204*^9, 3.814089338258224*^9}, { 3.814089451443385*^9, 3.8140894819478827`*^9}, {3.8140993755286617`*^9, 3.8140993819464817`*^9}, {3.814111233827095*^9, 3.814111234585766*^9}, { 3.8141937352989483`*^9, 3.8141937354208307`*^9}, {3.814261444878549*^9, 3.814261445448695*^9}, 3.8142638822664003`*^9, 3.8144385934982243`*^9, { 3.814438690490912*^9, 3.814438739303619*^9}, {3.8144389062214327`*^9, 3.814438906332591*^9}, {3.814441073479191*^9, 3.814441073692215*^9}, { 3.814446045628806*^9, 3.814446127392247*^9}, {3.814446235928728*^9, 3.8144462362331047`*^9}, 3.814463036566448*^9, {3.8144631932764997`*^9, 3.814463198849606*^9}, 3.8144633417275543`*^9, {3.814525244370717*^9, 3.814525244456441*^9}, {3.81452891029209*^9, 3.814528910393132*^9}, { 3.814726881246932*^9, 3.814726926456852*^9}, {3.81473119049043*^9, 3.81473119075737*^9}, {3.814955642172789*^9, 3.8149556521052628`*^9}, { 3.814955731111578*^9, 3.814955732748191*^9}, {3.814955770165793*^9, 3.814955770703888*^9}, {3.814958012277232*^9, 3.81495801802678*^9}, { 3.814963351670281*^9, 3.814963352308837*^9}, {3.8149636132573433`*^9, 3.814963618729459*^9}, {3.814971929590786*^9, 3.814971931029499*^9}, { 3.8149772264248667`*^9, 3.814977230036426*^9}, {3.814980964727211*^9, 3.814980965599844*^9}, {3.814986347919033*^9, 3.8149863619952307`*^9}, { 3.814993743733018*^9, 3.814993746475277*^9}, 3.815040446523501*^9, { 3.8150427185404177`*^9, 3.815042739444206*^9}, 3.815068339736587*^9, { 3.8150781203016872`*^9, 3.81507812124226*^9}, 3.815078295869111*^9, { 3.815216117454565*^9, 3.815216206708527*^9}, {3.815223181591799*^9, 3.815223197012802*^9}, {3.8152259801526737`*^9, 3.8152259845960417`*^9}, { 3.815232672911747*^9, 3.815232824725021*^9}, {3.815235791612417*^9, 3.815235802240562*^9}, {3.81523807518472*^9, 3.815238078522222*^9}, { 3.815240629978879*^9, 3.815240630041697*^9}, {3.815243010902005*^9, 3.8152430163596487`*^9}, {3.8152443806909027`*^9, 3.815244396457423*^9}, { 3.81524538362494*^9, 3.81524542639746*^9}, {3.815247871792701*^9, 3.8152478742026176`*^9}, 3.815253721754489*^9, {3.815257555839304*^9, 3.815257558989148*^9}, {3.815296706578437*^9, 3.8152967072615347`*^9}, { 3.815303477759705*^9, 3.815303478901589*^9}, {3.815309569997998*^9, 3.815309573127572*^9}, 3.815317877141842*^9, {3.815323899307622*^9, 3.8153239140408993`*^9}, {3.815323955295681*^9, 3.815323959418511*^9}, { 3.815326328892902*^9, 3.8153263314450903`*^9}, {3.815329646147807*^9, 3.815329649503935*^9}, 3.815332008549094*^9, {3.815336238939065*^9, 3.815336242929007*^9}, {3.815340372906146*^9, 3.815340381187139*^9}, { 3.8153429784528217`*^9, 3.8153429829636583`*^9}, {3.8154684337393637`*^9, 3.815468437636644*^9}, {3.815472468824298*^9, 3.815472470981311*^9}, { 3.815472994658319*^9, 3.81547299579732*^9}, 3.8154736252480288`*^9, { 3.815473939148929*^9, 3.815473945568803*^9}, {3.815474093893311*^9, 3.815474095084651*^9}, {3.815475089234436*^9, 3.815475108372904*^9}, { 3.815475180975233*^9, 3.815475203266609*^9}, {3.8154753006150227`*^9, 3.815475301035687*^9}, {3.815475851153833*^9, 3.815475854303007*^9}, { 3.815478941441476*^9, 3.8154789425356627`*^9}, {3.8154813185778933`*^9, 3.815481327694879*^9}, {3.815483414481146*^9, 3.815483418661022*^9}, { 3.815486096663032*^9, 3.815486098006947*^9}, {3.81549091013527*^9, 3.8154909112058973`*^9}, {3.8154931089640083`*^9, 3.815493110016161*^9}, { 3.815495196745026*^9, 3.815495199957196*^9}, {3.81549745604943*^9, 3.815497456277731*^9}, {3.815502230000288*^9, 3.8155022309750233`*^9}, { 3.815504266276462*^9, 3.8155042682055798`*^9}, {3.815506734665807*^9, 3.8155067390055447`*^9}, 3.8155114500774403`*^9, {3.8155142328610764`*^9, 3.815514236947043*^9}, {3.815566735931329*^9, 3.815566742279497*^9}, { 3.815569228376494*^9, 3.815569229573752*^9}, {3.815571798097188*^9, 3.81557181465698*^9}, {3.815574425029833*^9, 3.815574426840082*^9}, 3.815578656169038*^9, {3.8155807929320097`*^9, 3.815580797085114*^9}, { 3.8155837573953857`*^9, 3.815583759350883*^9}, {3.815585790818914*^9, 3.815585797102417*^9}, 3.815588337241083*^9, {3.815594860985737*^9, 3.815594865280905*^9}, {3.8155970599631243`*^9, 3.815597060677361*^9}, { 3.815600596215267*^9, 3.8156005980830193`*^9}, {3.8162518848030376`*^9, 3.816251888347568*^9}, {3.816252183744873*^9, 3.816252186351049*^9}, { 3.816256841593628*^9, 3.816256843854102*^9}, {3.8162570288467484`*^9, 3.8162570351428347`*^9}, {3.816263811049534*^9, 3.816263918027423*^9}, { 3.816264051079234*^9, 3.816264053085997*^9}, {3.816264100992408*^9, 3.8162641230377283`*^9}, {3.816621705469274*^9, 3.8166217067531843`*^9}, { 3.819990347265308*^9, 3.8199903753188553`*^9}, {3.819990483726313*^9, 3.8199904907960377`*^9}, {3.819999790686673*^9, 3.819999804099369*^9}, { 3.820004704611774*^9, 3.820004708579886*^9}, {3.820007553301672*^9, 3.8200075555168858`*^9}, {3.820007723792789*^9, 3.820007727207649*^9}, { 3.8200079182222757`*^9, 3.82000792015942*^9}, {3.8200080739081573`*^9, 3.820008074217718*^9}, {3.8200110130127068`*^9, 3.820011023689701*^9}, { 3.8200141839651117`*^9, 3.820014186834499*^9}, {3.8200167541796103`*^9, 3.820016764008779*^9}, {3.820054989971693*^9, 3.8200549919657183`*^9}, { 3.820057744767446*^9, 3.820057745917789*^9}, {3.8200606607351637`*^9, 3.820060662676178*^9}, {3.820065147961714*^9, 3.820065150827683*^9}, { 3.820068280579616*^9, 3.820068282139679*^9}, {3.82007502272544*^9, 3.8200750373067617`*^9}, {3.8200784763842087`*^9, 3.820078479850442*^9}, { 3.820082753480589*^9, 3.82008275957128*^9}, 3.8200829122399054`*^9, { 3.820086758459858*^9, 3.820086761259838*^9}, {3.820090442770599*^9, 3.82009044313873*^9}, {3.820094323336689*^9, 3.8200943261514997`*^9}, { 3.820231363709322*^9, 3.820231364911353*^9}, {3.820231448980515*^9, 3.8202316586881647`*^9}, {3.8202317665463343`*^9, 3.82023179479694*^9}, { 3.8202318676842003`*^9, 3.820231869181456*^9}, {3.820233088260557*^9, 3.820233088343128*^9}, {3.8202333774814653`*^9, 3.820233390511446*^9}, { 3.82023767993183*^9, 3.820237682425991*^9}, {3.820243109387423*^9, 3.820243112084889*^9}, {3.8202476637875023`*^9, 3.820247665549646*^9}, { 3.820667535853363*^9, 3.8206675727657547`*^9}, {3.820684183665489*^9, 3.820684205656611*^9}, {3.820684235767518*^9, 3.820684246489502*^9}, { 3.820684317683414*^9, 3.820684319064251*^9}, {3.8206844046414337`*^9, 3.8206844051977*^9}, {3.820687675556394*^9, 3.82068767789845*^9}, { 3.8206897462957783`*^9, 3.8206897499109173`*^9}, {3.8206930158189373`*^9, 3.820693017331942*^9}, {3.820696719878152*^9, 3.820696723881123*^9}, { 3.820699104252255*^9, 3.820699106178296*^9}, {3.820702286715938*^9, 3.820702294027163*^9}, {3.820704623117988*^9, 3.820704625496578*^9}, { 3.820752256916325*^9, 3.820752268297263*^9}, {3.820935073104765*^9, 3.8209350735634823`*^9}, {3.821211778288844*^9, 3.8212117827185*^9}, { 3.821296204646468*^9, 3.8212962076337767`*^9}, {3.821349644394158*^9, 3.8213496467786283`*^9}, {3.8213517714599857`*^9, 3.821351773461198*^9}, { 3.8213622462608747`*^9, 3.82136225015648*^9}, {3.82136337627873*^9, 3.8213633783876143`*^9}, {3.822395373671977*^9, 3.822395374080752*^9}, { 3.8224055384543467`*^9, 3.822405538943552*^9}, {3.82256976116735*^9, 3.822569764095625*^9}, 3.822569941284091*^9, 3.822578176175564*^9, { 3.8225814109161673`*^9, 3.822581423920046*^9}, {3.8225815413929853`*^9, 3.822581541836155*^9}, {3.822581734173427*^9, 3.822581734495254*^9}, { 3.8225818623627577`*^9, 3.822581862513247*^9}, {3.822582058756192*^9, 3.822582059063018*^9}, {3.8225844873896627`*^9, 3.8225844876510353`*^9}, { 3.822585748894499*^9, 3.8225857497139053`*^9}, {3.822587769069633*^9, 3.822587770064171*^9}, {3.8225892219754543`*^9, 3.8225892227103*^9}, { 3.822589541463584*^9, 3.822589544052847*^9}, {3.822589634114785*^9, 3.822589634635572*^9}, {3.822589698088089*^9, 3.822589698396234*^9}, { 3.8225897556410513`*^9, 3.822589764728793*^9}, {3.822590017771882*^9, 3.822590018294026*^9}, {3.8225919313799467`*^9, 3.822591931623008*^9}, { 3.823198948936804*^9, 3.823198949458517*^9}, {3.823200343746475*^9, 3.823200343890626*^9}, {3.823288448622332*^9, 3.823288448965056*^9}, { 3.823288482448771*^9, 3.823288601222763*^9}, {3.823288666388864*^9, 3.8232889910948877`*^9}, {3.823291400932473*^9, 3.823291401249125*^9}, { 3.82329599250668*^9, 3.823295994343355*^9}, {3.823297148165262*^9, 3.8232972717869377`*^9}, {3.8232973083471746`*^9, 3.823297369180789*^9}, { 3.8235570605732317`*^9, 3.823557062651964*^9}, {3.823558465416039*^9, 3.823558466961543*^9}, {3.8236002463490343`*^9, 3.823600286499044*^9}, 3.824298117566319*^9, 3.824298186422042*^9, {3.824298394352096*^9, 3.824298394748021*^9}, {3.838386131568082*^9, 3.838386184514793*^9}, { 3.841324611766222*^9, 3.841324847738583*^9}, {3.841324895964373*^9, 3.8413248963384447`*^9}, {3.841325018315042*^9, 3.84132502249907*^9}, { 3.8413261274271803`*^9, 3.841326127588489*^9}, {3.8423342326069593`*^9, 3.842334232690463*^9}, {3.842336719830854*^9, 3.84233671991934*^9}, { 3.842339418400468*^9, 3.842339418478218*^9}, {3.8430359142557707`*^9, 3.843035914310397*^9}, {3.8430434732807913`*^9, 3.843043473710072*^9}, { 3.8436355220718927`*^9, 3.843635544439931*^9}, {3.843636448787025*^9, 3.843636448845724*^9}, {3.84363812735776*^9, 3.843638127841715*^9}, { 3.844171863721094*^9, 3.84417186451587*^9}, {3.844258183413237*^9, 3.844258184100439*^9}, {3.844258242097528*^9, 3.844258243538804*^9}, 3.84486037069027*^9, {3.867749889058716*^9, 3.867749936179331*^9}, { 3.8679199261026363`*^9, 3.8679199274306736`*^9}, {3.8679231011136265`*^9, 3.8679231016291246`*^9}, {3.870171114504943*^9, 3.8701711207071257`*^9}, { 3.870338144641303*^9, 3.87033815756642*^9}, 3.8705181482874193`*^9, { 3.870519560121092*^9, 3.870519576315165*^9}, {3.870522238568857*^9, 3.87052226256144*^9}, {3.87052272699298*^9, 3.870522730585904*^9}, { 3.8705243103357563`*^9, 3.870524316114233*^9}, {3.871045958467011*^9, 3.871045963910185*^9}},ExpressionUUID->"4cb1ff4a-a751-482f-b501-\ ee2009ee5be2"], Cell[BoxData["\<\"\[Beta] = \"\>"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.8710462289988937`*^9},ExpressionUUID->"76f778b7-ad19-4460-9570-\ 4122290b7c2f"], Cell[BoxData["1.`"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.871046228999864*^9},ExpressionUUID->"168eb232-a29e-437c-a1a2-\ 915461a092a1"], Cell[BoxData["\<\"\[Eta] = \"\>"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.871046229000695*^9},ExpressionUUID->"8178d135-110f-48d2-8d24-\ bb2b83c38d68"], Cell[BoxData["7.`"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.87104622900154*^9},ExpressionUUID->"eeacdccc-f38a-49f5-bda9-\ 4a4ccea9eb1b"], Cell[BoxData["\<\"E1 must be bigger than\"\>"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.871046229002355*^9},ExpressionUUID->"f19230cd-6fff-41d6-9914-\ 8ac5d1b7e7a9"], Cell[BoxData["0.5`"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.871046229003243*^9},ExpressionUUID->"133a27c2-1a89-447b-b1a3-\ 55ee45f38e5e"], Cell[BoxData[ FractionBox[ RowBox[{"0.5625000000000001`", "\[VeryThinSpace]", "+", RowBox[{"1.28`", " ", SuperscriptBox["E1", "3"]}], "-", RowBox[{"0.09`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "0.72`"}], "+", RowBox[{"0.035`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}]}], ")"}]}], "+", RowBox[{ SuperscriptBox["E1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "0.44`"}], "+", RowBox[{"0.16`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}]}], ")"}]}]}], RowBox[{ RowBox[{"-", "3.87`"}], "+", RowBox[{"8.`", " ", SuperscriptBox["E1", "3"]}], "+", RowBox[{"0.64125`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}], "+", RowBox[{"E1", " ", RowBox[{"(", RowBox[{"12.870000000000003`", "\[VeryThinSpace]", "-", RowBox[{"1.50125`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}]}], ")"}]}], "+", RowBox[{ SuperscriptBox["E1", "2"], " ", RowBox[{"(", RowBox[{ RowBox[{"-", "16.509999999999998`"}], "+", RowBox[{"1.`", " ", SqrtBox[ RowBox[{"36.`", "\[VeryThinSpace]", "-", RowBox[{"72.`", " ", "E1"}], "+", RowBox[{"64.`", " ", SuperscriptBox["E1", "2"]}]}]]}]}], ")"}]}]}]]], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.8710462291512947`*^9},ExpressionUUID->"e36be28a-ae22-4597-b095-\ 8cce1c19edc1"], Cell[BoxData["\<\"must be bigger than\"\>"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.871046229153404*^9},ExpressionUUID->"e77b46f3-a040-42aa-86d0-\ 4a825e3c2005"], Cell[BoxData["0.16`"], "Output", CellChangeTimes->{{3.844258251966549*^9, 3.8442582535639257`*^9}, 3.8448603810195293`*^9, 3.867919984884941*^9, 3.8679234606020203`*^9, 3.8701722590168257`*^9, 3.870176327705471*^9, 3.870338754268497*^9, 3.870520010973961*^9, 3.8705221931446533`*^9, 3.870522771448851*^9, 3.870524363878355*^9, 3.87077504062361*^9, 3.8710462291543694`*^9},ExpressionUUID->"a01c009f-05c5-48ea-a6b9-\ 28e15682a134"] }, Open ]], Cell[CellGroupData[{ Cell["Section 2: Normal form transformation", "Section", CellChangeTimes->{{3.79883476781522*^9, 3.798834770810676*^9}, 3.798834834560101*^9, {3.8705181266021976`*^9, 3.870518146645891*^9}, { 3.870770009689802*^9, 3.870770012472755*^9}, {3.870770223959663*^9, 3.870770224310103*^9}},ExpressionUUID->"eca8fcbe-e1d7-4292-acc1-\ 708515933e00"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ClearAll", "[", RowBox[{ "c", ",", "CurlyF", ",", "CurlyG", ",", "p", ",", "h", ",", "S", ",", "\[Mu]", ",", "J"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"J", " ", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "p"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "h"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "p"}], "]"}], ",", RowBox[{"D", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}], ",", "h"}], "]"}]}], "}"}]}], "}"}], " ", "/.", RowBox[{"{", RowBox[{ RowBox[{"h", " ", "\[Rule]", " ", "hs"}], ",", RowBox[{"p", "\[Rule]", " ", "ps"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Construct", " ", "the", " ", "transformation", " ", "matrix"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"S", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"i", "/", RowBox[{"J", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"J", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "-", "r"}], ")"}]}], "/", RowBox[{"J", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"c", "=", RowBox[{"ccrit", " ", "+", " ", "\[Mu]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"$Assumptions", " ", "=", " ", RowBox[{"\[Mu]", ">", "0"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"To", " ", "obtain", " ", "the", " ", "non"}], "-", RowBox[{ "linear", " ", "part", " ", "of", " ", "the", " ", "transformed", " ", "equations"}]}], ",", " ", RowBox[{ RowBox[{ "subtract", " ", "the", " ", "linear", " ", "part", " ", "with", " ", "\[Mu]"}], " ", "=", " ", "0"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"p", ",", "h"}], "}"}], " ", "=", " ", RowBox[{ RowBox[{"S", ".", RowBox[{"{", RowBox[{"P", ",", "H"}], "}"}]}], "+", RowBox[{"{", RowBox[{"ps", ",", "hs"}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"normal", " ", "form", " ", "non"}], "-", RowBox[{"linear", " ", "part"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"CurlyF", ",", "CurlyG"}], "}"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Inverse", "[", "S", "]"}], ".", RowBox[{"{", RowBox[{ RowBox[{"F", "[", RowBox[{"h", ",", "p"}], "]"}], ",", RowBox[{"G", "[", RowBox[{"h", ",", "p"}], "]"}]}], "}"}]}], " ", "-", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "\[CapitalOmega]0"}], " ", "H"}], ",", RowBox[{"\[CapitalOmega]0", " ", "P"}]}], "}"}]}], " ", "/.", " ", RowBox[{"\[Mu]", " ", "\[Rule]", " ", "0"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"testing", " ", "-", " ", RowBox[{ "all", " ", "first", " ", "derivatives", " ", "should", " ", "be", " ", "0", " ", "at", " ", RowBox[{"(", RowBox[{"H", ",", "P"}], ")"}]}]}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"0", ",", "0"}], ")"}], ".", " ", "Use"}], " ", "random", " ", "values", " ", "for", " ", "E1"}]}], ",", RowBox[{ "B1", " ", "because", " ", "it", " ", "takes", " ", "large", " ", "amounts", " ", "of", " ", "time", " ", "to", " ", "calculate", " ", "in", " ", "the", " ", "general", " ", "model"}]}], "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Chop", "[", RowBox[{ RowBox[{"{", RowBox[{"CurlyF", ",", "CurlyG"}], "}"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}], ",", " ", RowBox[{"E1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}], ",", " ", RowBox[{"B1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}]}], "}"}]}], "]"}], "\n", RowBox[{"Chop", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", "H"}], "]"}], ",", RowBox[{"D", "[", RowBox[{"CurlyF", ",", "P"}], "]"}]}], "}"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}], ",", " ", RowBox[{"E1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}], ",", " ", RowBox[{"B1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}]}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Chop", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", "H"}], "]"}], ",", RowBox[{"D", "[", RowBox[{"CurlyG", ",", "P"}], "]"}]}], "}"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}], ",", " ", RowBox[{"E1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}], ",", " ", RowBox[{"B1", " ", "\[Rule]", " ", RowBox[{"RandomReal", "[", "]"}]}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.7988348487424717`*^9, 3.798835058317482*^9}, { 3.7988351146992407`*^9, 3.798835127315126*^9}, {3.798835167865234*^9, 3.798835171332934*^9}, {3.798835910492838*^9, 3.798835910735024*^9}, { 3.799089357389271*^9, 3.7990893613624372`*^9}, {3.799089494167211*^9, 3.799089494679369*^9}, {3.799089532238475*^9, 3.799089540145364*^9}, { 3.7990895782853394`*^9, 3.799089641146096*^9}, {3.799095082815982*^9, 3.7990951288206987`*^9}, {3.799095167891527*^9, 3.799095185425871*^9}, { 3.799095245484201*^9, 3.7990953238312902`*^9}, {3.799095360044105*^9, 3.7990954235622*^9}, {3.799095714923703*^9, 3.7990958135048103`*^9}, { 3.799095911875805*^9, 3.799095966110196*^9}, {3.799096079122109*^9, 3.799096094571516*^9}, {3.7990961922298822`*^9, 3.79909626022867*^9}, { 3.799096304144013*^9, 3.79909635795863*^9}, {3.7990964045000687`*^9, 3.799096433764577*^9}, {3.799096522034423*^9, 3.799096525982604*^9}, { 3.7991756819372187`*^9, 3.7991757372789183`*^9}, {3.799175800104061*^9, 3.799175803921924*^9}, {3.799175858712253*^9, 3.799175886754517*^9}, { 3.7993387694242268`*^9, 3.799338783626523*^9}, {3.799340963783258*^9, 3.79934098016286*^9}, {3.799341017545733*^9, 3.799341059865057*^9}, { 3.7993411116901627`*^9, 3.7993411195365953`*^9}, {3.799341155821512*^9, 3.7993411559383287`*^9}, {3.7993412017996397`*^9, 3.799341201917316*^9}, { 3.7993413101345863`*^9, 3.799341368057558*^9}, {3.799341398316296*^9, 3.799341421220695*^9}, {3.799341856231319*^9, 3.799341863809473*^9}, { 3.799341956006513*^9, 3.7993419901563168`*^9}, {3.799342043087694*^9, 3.7993420546829643`*^9}, {3.799345378463757*^9, 3.799345411140834*^9}, 3.7993454475953283`*^9, {3.7993455592593412`*^9, 3.799345570136397*^9}, { 3.799345907815031*^9, 3.799345910576844*^9}, {3.799345965669057*^9, 3.799345973306756*^9}, {3.7993461028902893`*^9, 3.799346114918655*^9}, { 3.7993461613473663`*^9, 3.79934616280412*^9}, {3.7993463326362553`*^9, 3.799346366078745*^9}, 3.7993468408627872`*^9, {3.799347365776374*^9, 3.799347373470913*^9}, {3.799347411066841*^9, 3.79934741915621*^9}, { 3.799352443777791*^9, 3.799352517746533*^9}, {3.799352553905456*^9, 3.799352561603738*^9}, {3.799353998602828*^9, 3.7993540846574793`*^9}, { 3.799354147093877*^9, 3.7993541472216587`*^9}, {3.799355354053388*^9, 3.799355358129335*^9}, 3.799426903423232*^9, 3.7994270193947487`*^9, { 3.799427071335475*^9, 3.799427132602598*^9}, {3.799427170339905*^9, 3.799427175107684*^9}, 3.799427450412916*^9, {3.7994277844449368`*^9, 3.799427784977961*^9}, {3.799428099604631*^9, 3.799428202654277*^9}, { 3.799428257737465*^9, 3.799428285275159*^9}, {3.799428371150613*^9, 3.799428371944509*^9}, {3.799432972609164*^9, 3.799432973428636*^9}, { 3.799433195476212*^9, 3.799433241881061*^9}, {3.799433704817164*^9, 3.799433753994651*^9}, {3.799433824445882*^9, 3.7994338290509644`*^9}, { 3.799433887657426*^9, 3.7994339094446917`*^9}, 3.799433958651613*^9, 3.799434005948946*^9, {3.799434097575324*^9, 3.799434097867982*^9}, { 3.7994341501401777`*^9, 3.7994341637718897`*^9}, {3.7994343171328297`*^9, 3.7994343218248262`*^9}, {3.799434545355991*^9, 3.799434611326601*^9}, { 3.7994346547788267`*^9, 3.799434655351083*^9}, {3.7994346873006372`*^9, 3.799434726789886*^9}, {3.79943654072611*^9, 3.799436561166198*^9}, { 3.799436679535039*^9, 3.799436704532559*^9}, {3.799436755435421*^9, 3.79943677885997*^9}, {3.79943686608679*^9, 3.799436870942277*^9}, { 3.799436931420692*^9, 3.799436947060602*^9}, {3.799437021258567*^9, 3.799437032834518*^9}, {3.799437171602499*^9, 3.7994371769389973`*^9}, { 3.7994372163917007`*^9, 3.799437267688345*^9}, {3.799437470491337*^9, 3.799437480242732*^9}, {3.799437515084627*^9, 3.799437526251411*^9}, { 3.799437598391882*^9, 3.799437603193061*^9}, {3.799437638294631*^9, 3.799437655167232*^9}, {3.7994376919733257`*^9, 3.799437692169332*^9}, { 3.799437728381845*^9, 3.79943778304117*^9}, {3.799437896444227*^9, 3.799437905190357*^9}, {3.799438728159216*^9, 3.799438747646502*^9}, { 3.7994387779116077`*^9, 3.7994387780447407`*^9}, {3.7995154611705647`*^9, 3.799515495222445*^9}, {3.7995181731959877`*^9, 3.799518372124672*^9}, { 3.799518433904194*^9, 3.7995184409794807`*^9}, {3.799520981571788*^9, 3.799521076431209*^9}, {3.799521648914073*^9, 3.799521666672291*^9}, { 3.799521730356344*^9, 3.7995217643979*^9}, 3.79952188246239*^9, 3.799521935965744*^9, {3.799522057493698*^9, 3.7995220805286217`*^9}, { 3.799522165023143*^9, 3.799522179180606*^9}, {3.799522229041395*^9, 3.799522245064209*^9}, {3.799522363330195*^9, 3.79952239415147*^9}, { 3.799522428674479*^9, 3.799522472521854*^9}, {3.7995225253933496`*^9, 3.799522534849473*^9}, {3.799522592330235*^9, 3.799522593111589*^9}, { 3.7995226605833397`*^9, 3.799522663460053*^9}, {3.799522749455667*^9, 3.799522753274911*^9}, 3.799522813097053*^9, {3.799522853370618*^9, 3.799522862425549*^9}, {3.799522895973949*^9, 3.79952294079729*^9}, { 3.799522997403302*^9, 3.799523013138096*^9}, {3.7995231220637817`*^9, 3.7995231488245068`*^9}, {3.799523213460495*^9, 3.7995232859664583`*^9}, { 3.799523318977768*^9, 3.79952344534625*^9}, {3.7995236048521214`*^9, 3.799523612834495*^9}, {3.799523839401*^9, 3.799523840281034*^9}, { 3.799523904363489*^9, 3.7995239093145027`*^9}, {3.79952394076747*^9, 3.799523941328424*^9}, {3.7995245015780373`*^9, 3.799524530445471*^9}, { 3.799524674776546*^9, 3.7995246791407957`*^9}, {3.7995248623218193`*^9, 3.799524863277532*^9}, {3.799525646766768*^9, 3.7995256562962313`*^9}, { 3.799525780800481*^9, 3.7995258423179283`*^9}, {3.799525880241434*^9, 3.79952592755272*^9}, {3.799526200777349*^9, 3.7995262066368103`*^9}, { 3.799526260758526*^9, 3.799526294862406*^9}, {3.7995264700024023`*^9, 3.7995264804760857`*^9}, {3.799526674062392*^9, 3.7995267449218197`*^9}, { 3.799526785724243*^9, 3.799526838065874*^9}, 3.799526936296583*^9, { 3.799526990905416*^9, 3.799527012173915*^9}, {3.7995270564255247`*^9, 3.7995270709659367`*^9}, {3.799529891874312*^9, 3.7995299342250423`*^9}, { 3.799529981296877*^9, 3.799529982927733*^9}, {3.799530035298167*^9, 3.799530064391821*^9}, {3.799530112790209*^9, 3.7995301319878483`*^9}, { 3.799530188464777*^9, 3.799530194949585*^9}, {3.799532845019743*^9, 3.799532874935936*^9}, {3.799532906684589*^9, 3.799532920320079*^9}, { 3.799533109214205*^9, 3.7995331172680273`*^9}, {3.7995332921024942`*^9, 3.799533348859406*^9}, {3.7995335181146584`*^9, 3.799533542122404*^9}, { 3.7995337631396847`*^9, 3.799533786452848*^9}, {3.799534017199758*^9, 3.799534021084608*^9}, {3.799534160201983*^9, 3.7995341655081253`*^9}, 3.799578782186886*^9, {3.799578832291832*^9, 3.799578878062407*^9}, 3.799578922368538*^9, {3.799578971353388*^9, 3.799578976368635*^9}, { 3.799579135781185*^9, 3.799579163018886*^9}, {3.799579674573395*^9, 3.7995796749720993`*^9}, {3.799579717886251*^9, 3.799579740674549*^9}, { 3.799579865621243*^9, 3.799579894685594*^9}, {3.799609106577841*^9, 3.799609116385035*^9}, {3.799609184274444*^9, 3.7996091930780487`*^9}, { 3.79960930458004*^9, 3.7996093163673773`*^9}, {3.799609351508379*^9, 3.799609380579658*^9}, {3.799609418941764*^9, 3.799609421248654*^9}, { 3.7996098616348867`*^9, 3.799609919025742*^9}, {3.799610141119808*^9, 3.799610205713431*^9}, {3.799610274307171*^9, 3.799610294682085*^9}, { 3.799611656862503*^9, 3.799611682303583*^9}, {3.7996117302654743`*^9, 3.799611854072398*^9}, {3.799611981821328*^9, 3.799611996072226*^9}, { 3.799612298988312*^9, 3.79961237842865*^9}, {3.799613217032671*^9, 3.799613260734578*^9}, {3.799873561365934*^9, 3.799873583684581*^9}, { 3.8007346329921703`*^9, 3.800734645957387*^9}, {3.80073475642328*^9, 3.800734756798697*^9}, 3.800893217392641*^9, {3.8010864372876587`*^9, 3.8010864534559317`*^9}, {3.802025201791235*^9, 3.802025250757512*^9}, 3.802025282739533*^9, {3.802028667967813*^9, 3.80202867069632*^9}, { 3.811664279943062*^9, 3.811664368521779*^9}, {3.811684784066297*^9, 3.811684790849947*^9}, {3.811684845407648*^9, 3.811684851508335*^9}, { 3.811684890619267*^9, 3.8116849294083347`*^9}, {3.811685012307527*^9, 3.811685017078315*^9}, {3.811685063030916*^9, 3.811685073464663*^9}, { 3.8116851323240557`*^9, 3.811685196387958*^9}, {3.8116852514166117`*^9, 3.811685252489002*^9}, {3.81168554627309*^9, 3.811685547381469*^9}, { 3.8116855886608953`*^9, 3.811685591500963*^9}, {3.811685631061051*^9, 3.811685692388916*^9}, {3.811698673121139*^9, 3.8116986751184998`*^9}, { 3.811698734820423*^9, 3.81169875492728*^9}, {3.81169879167634*^9, 3.811698803750761*^9}, 3.811700196822193*^9, {3.811700266633748*^9, 3.811700270311401*^9}, {3.812210536021146*^9, 3.812210536728833*^9}, { 3.812224973222921*^9, 3.8122249930275373`*^9}, {3.812296809817294*^9, 3.812296827424548*^9}, {3.812732744680821*^9, 3.812732746289981*^9}, { 3.812732797745035*^9, 3.812732799685396*^9}, {3.8127328368572407`*^9, 3.8127328376096983`*^9}, {3.8129013998091183`*^9, 3.812901403365157*^9}, { 3.812902178293212*^9, 3.812902183083437*^9}, {3.8132485475776663`*^9, 3.813248552755896*^9}, {3.813680712806982*^9, 3.813680722245483*^9}, 3.8136807553143587`*^9, 3.813680802323978*^9, {3.813680897491807*^9, 3.813680903961872*^9}, {3.813763545202949*^9, 3.813763648571319*^9}, { 3.813830183689724*^9, 3.813830195947248*^9}, {3.814111456450548*^9, 3.814111476165648*^9}, {3.8141115133237877`*^9, 3.814111522628132*^9}, { 3.8141115647774982`*^9, 3.814111573468211*^9}, {3.8141116203724127`*^9, 3.8141116300214453`*^9}, {3.814111767063436*^9, 3.8141118027032213`*^9}, { 3.814111862344685*^9, 3.814111868257221*^9}, 3.814111905843955*^9, { 3.8141122425241137`*^9, 3.8141123712809057`*^9}, {3.8141125694999866`*^9, 3.814112580300042*^9}, {3.814179260650024*^9, 3.8141792611823063`*^9}, { 3.814193790738351*^9, 3.814193814095142*^9}, {3.814193856718955*^9, 3.814193859755105*^9}, {3.814446527779943*^9, 3.814446554983561*^9}, { 3.814955790237856*^9, 3.814955796624888*^9}, {3.8154750306588573`*^9, 3.815475034502404*^9}, {3.815475116457958*^9, 3.815475127116906*^9}, { 3.815475310317337*^9, 3.815475311263208*^9}, {3.8154757127805357`*^9, 3.815475717378522*^9}, {3.815475764313862*^9, 3.815475765406542*^9}, { 3.8199904131395617`*^9, 3.819990451203085*^9}, {3.819990811157269*^9, 3.819990813756776*^9}, {3.819994610091831*^9, 3.8199946315758266`*^9}, { 3.8209241807031193`*^9, 3.820924201500821*^9}, {3.867749992193771*^9, 3.8677500794320717`*^9}, {3.870171135547574*^9, 3.870171139955401*^9}, { 3.870246830914426*^9, 3.87024684233675*^9}, {3.870246968231255*^9, 3.870246993860992*^9}, 3.870518130720254*^9, {3.870519613933207*^9, 3.870519649800346*^9}, {3.8705205124931307`*^9, 3.870520570257195*^9}, { 3.870522013195073*^9, 3.870522017077323*^9}, {3.870522049987616*^9, 3.8705220504311333`*^9}, {3.870522183932538*^9, 3.870522225762966*^9}, { 3.870522311856592*^9, 3.870522338604567*^9}, {3.8705223689548473`*^9, 3.8705223899346523`*^9}, {3.8705227027947683`*^9, 3.870522722028512*^9}, { 3.8705242872421*^9, 3.8705243037731123`*^9}, {3.871045985882107*^9, 3.871046034510119*^9}, {3.8710460705304728`*^9, 3.871046107057857*^9}},ExpressionUUID->"0f5a2bd7-6372-4d3a-a9e5-\ 51ea3d601cab"], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.8703387612520943`*^9, 3.870520018138721*^9, 3.870521995863604*^9, 3.870522034194953*^9, 3.8705222105463123`*^9, 3.8705222907811213`*^9, 3.870522322257926*^9, 3.870522373322958*^9, 3.870522412416808*^9, 3.870522694945424*^9, 3.870524276524688*^9, 3.870524370938775*^9, 3.8707750477161913`*^9, 3.871046236188785*^9},ExpressionUUID->"949a2f6a-e00e-4b68-9303-\ 24e8d157f9a4"], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.8703387612520943`*^9, 3.870520018138721*^9, 3.870521995863604*^9, 3.870522034194953*^9, 3.8705222105463123`*^9, 3.8705222907811213`*^9, 3.870522322257926*^9, 3.870522373322958*^9, 3.870522412416808*^9, 3.870522694945424*^9, 3.870524276524688*^9, 3.870524370938775*^9, 3.8707750477161913`*^9, 3.871046236688994*^9},ExpressionUUID->"ed50d2d5-ea87-41e9-a58b-\ 5587d38d352f"], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{3.8703387612520943`*^9, 3.870520018138721*^9, 3.870521995863604*^9, 3.870522034194953*^9, 3.8705222105463123`*^9, 3.8705222907811213`*^9, 3.870522322257926*^9, 3.870522373322958*^9, 3.870522412416808*^9, 3.870522694945424*^9, 3.870524276524688*^9, 3.870524370938775*^9, 3.8707750477161913`*^9, 3.8710462368697653`*^9},ExpressionUUID->"5b8e8da4-2463-40ca-9e58-\ 466494018b78"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Section 4: Evaluate \[Omega]_1 and generate the values for the E_1xB_1 plot \ (Figure 5)\ \>", "Section", CellChangeTimes->{{3.814117174619775*^9, 3.814117179560883*^9}, { 3.814118123554666*^9, 3.814118131276762*^9}, {3.815474985680491*^9, 3.815474991934511*^9}, {3.8206678287695513`*^9, 3.820667833264654*^9}, { 3.86775013008704*^9, 3.867750136951311*^9}, {3.870339528194578*^9, 3.8703395282718678`*^9}, {3.870518174005309*^9, 3.8705181954597197`*^9}, { 3.870522131788577*^9, 3.870522135311233*^9}, {3.870770016234364*^9, 3.870770017570657*^9}, {3.8707700628976297`*^9, 3.87077007117458*^9}, { 3.870770227417651*^9, 3.8707702276031857`*^9}, {3.871045872205563*^9, 3.87104589012372*^9}},ExpressionUUID->"21293870-4c70-42fd-883f-\ aec614437a87"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ClearAll", "[", RowBox[{ "E1", ",", "B1", ",", "list\[Omega]1", ",", "stabilitytest", ",", "Nx", ",", "Ny", ",", "Maxy", ",", "Maxx", ",", "Miny", ",", "Minx", ",", "p", ",", "h", ",", "\[Omega]1", ",", "HS", ",", "PS", ",", "Listccrit", ",", "listomega0"}], "]"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"grid", " ", "size"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Nx", " ", "=", " ", "20"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Ny", " ", "=", "20"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"parameter", " ", "range"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Maxy", " ", "=", " ", "10"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Miny", " ", "=", "1.00"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Maxx", " ", "=", " ", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Minx", " ", "=", " ", "0.55"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Mu]", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{"\"\\"", "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"AbsoluteTiming", "[", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"testing", " ", "if", " ", "hs"}], ",", " ", RowBox[{ "ps", " ", "and", " ", "ccrit", " ", "are", " ", "within", " ", "the", " ", "correct", " ", "range"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"HS", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"hs", " ", "<", " ", "0"}], " ", "||", " ", RowBox[{"hs", " ", ">", " ", "1"}]}], ",", RowBox[{"-", "1"}], ",", " ", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"PS", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"ps", " ", "<", " ", "0"}], " ", ",", RowBox[{"-", "1"}], ",", " ", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Listccrit", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"ccrit", " ", "<", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[Beta]", " ", "\[Eta]"}], " ", "+", " ", "1"}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"\[Beta]", " ", RowBox[{"\[Eta]", "^", "2"}]}], " ", "+", "1"}], ")"}]}]}], " ", ",", RowBox[{"-", "1"}], ",", " ", "1"}], "]"}], ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"ArrayPlot", "[", RowBox[{"HS", ",", RowBox[{"AxesLabel", "\[Rule]", "Automatic"}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Miny", ",", "Maxy"}], "}"}], ",", RowBox[{"{", RowBox[{"Minx", ",", "Maxx"}], "}"}]}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"B1", ",", "E1"}], "}"}]}], ",", RowBox[{"ColorRules", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"1", " ", "\[Rule]", " ", "White"}], ",", RowBox[{ RowBox[{"-", "1"}], "\[Rule]", " ", "Black"}]}], "}"}]}], ",", RowBox[{"DataReversed", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"ArrayPlot", "[", RowBox[{"PS", ",", RowBox[{"AxesLabel", "\[Rule]", "Automatic"}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Miny", ",", "Maxy"}], "}"}], ",", RowBox[{"{", RowBox[{"Minx", ",", "Maxx"}], "}"}]}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"B1", ",", "E1"}], "}"}]}], ",", RowBox[{"ColorRules", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"1", " ", "\[Rule]", " ", "White"}], ",", RowBox[{ RowBox[{"-", "1"}], "\[Rule]", " ", "Black"}]}], "}"}]}], ",", RowBox[{"DataReversed", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"ArrayPlot", "[", RowBox[{"Listccrit", ",", RowBox[{"AxesLabel", "\[Rule]", "Automatic"}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Miny", ",", "Maxy"}], "}"}], ",", RowBox[{"{", RowBox[{"Minx", ",", "Maxx"}], "}"}]}], "}"}]}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"B1", ",", "E1"}], "}"}]}], ",", RowBox[{"ColorRules", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"1", " ", "\[Rule]", " ", "White"}], ",", RowBox[{ RowBox[{"-", "1"}], "\[Rule]", " ", "Black"}]}], "}"}]}], ",", RowBox[{"DataReversed", "\[Rule]", "True"}]}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Computing", " ", "\[CapitalLambda]1", " ", "and", " ", "\[CapitalOmega]2", " ", "calculating", " ", "each", " ", "derivative", " ", "separately", " ", "for", " ", "efficiency"}], "*)"}], "\[IndentingNewLine]", RowBox[{"ClearAll", "[", RowBox[{ "dfppp", ",", "dfphh", ",", "dgpph", ",", "dfpp", ",", "dgpp", ",", "dfhh", ",", "dghh", ",", "dfph", ",", "dgph", ",", "dgppp", ",", "dfpph", ",", "dfhhh", ",", "list\[CapitalOmega]0", ",", "Lambda1", ",", "Omega2", ",", "dummy"}], "]"}], "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", " ", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", RowBox[{"{", RowBox[{"P", ",", "3"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfppp", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", " ", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", "P", ",", "H", ",", "H"}], "]"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfphh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", "P", ",", "P", ",", "H"}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dgpph", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", RowBox[{"{", RowBox[{"H", ",", "3"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dghhh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", RowBox[{"{", RowBox[{"P", ",", "2"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfpp", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", RowBox[{"{", RowBox[{"P", ",", "2"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dgpp", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", RowBox[{"{", RowBox[{"H", ",", "2"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfhh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", RowBox[{"{", RowBox[{"H", ",", "2"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dghh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", "P", ",", "H"}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfph", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", "P", ",", "H"}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dgph", " ", "=", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", "=", " ", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", RowBox[{"{", RowBox[{"P", ",", "3"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dgppp", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyG", ",", "P", ",", "H", ",", "H"}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dgphh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", "P", ",", "P", ",", "H"}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfpph", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", RowBox[{ RowBox[{"D", "[", RowBox[{"CurlyF", ",", RowBox[{"{", RowBox[{"H", ",", "3"}], "}"}]}], "]"}], "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"dfhhh", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"dummy", " ", "=", " ", RowBox[{"\[CapitalOmega]0", "/.", " ", RowBox[{"{", RowBox[{ RowBox[{"H", "\[Rule]", " ", "0"}], ",", RowBox[{"P", " ", "\[Rule]", " ", "0"}]}], "}"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"list\[CapitalOmega]0", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"dummy", ",", RowBox[{"{", RowBox[{"B1", ",", "Miny", ",", "Maxy", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxy", "-", "Miny"}], ")"}], "/", RowBox[{"(", RowBox[{"Ny", "-", "1"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"E1", ",", "Minx", ",", "Maxx", ",", RowBox[{ RowBox[{"(", RowBox[{"Maxx", "-", "Minx"}], ")"}], "/", RowBox[{"(", RowBox[{"Nx", "-", "1"}], ")"}]}]}], "}"}]}], "]"}]}], ";", " ", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Lambda1", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ "dfppp", "+", "dfphh", "+", "dgpph", " ", "+", " ", "dghhh"}], ")"}]}], "/", "16"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"dfpp", " ", "dgpp"}], " ", "-", " ", RowBox[{"dfhh", " ", "dghh"}], " ", "-", " ", RowBox[{"dfph", RowBox[{"(", RowBox[{"dfpp", " ", "+", " ", "dfhh"}], ")"}]}], "+", " ", RowBox[{"dgph", RowBox[{"(", RowBox[{"dghh", " ", "+", " ", "dgpp"}], ")"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"16", " ", "list\[CapitalOmega]0"}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Omega2", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ "dgppp", " ", "+", " ", "dgphh", " ", "-", " ", "dfpph", " ", "-", " ", "dfhhh"}], ")"}], "/", "16"}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"dfpp", RowBox[{"(", RowBox[{"dgph", " ", "-", " ", "dfhh"}], ")"}]}], " ", "+", " ", RowBox[{"dghh", RowBox[{"(", RowBox[{"dfph", " ", "-", " ", "dgpp"}], ")"}]}], " ", "-", " ", RowBox[{"3", "dfhh", RowBox[{"(", RowBox[{"dgph", "+", "dfhh"}], ")"}]}], "-", RowBox[{"3", "dgpp", RowBox[{"(", RowBox[{"dfph", "+", "dgpp"}], ")"}]}], " ", "-", " ", RowBox[{"2", RowBox[{ RowBox[{"(", RowBox[{"dfpp", "+", "dfhh", " ", "-", " ", "dgph"}], ")"}], "^", "2"}]}], "-", RowBox[{"2", RowBox[{ RowBox[{"(", RowBox[{"dgpp", "+", "dghh", " ", "-", " ", "dfph"}], ")"}], "^", "2"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{"48", " ", "list\[CapitalOmega]0"}], ")"}]}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"\[Omega]1", " ", "=", " ", RowBox[{ RowBox[{"-", "\[CapitalOmega]2"}], "/", "\[CapitalLambda]1"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"\[Omega]1", " ", "=", " ", RowBox[{ RowBox[{"-", " ", "Omega2"}], "/", "Lambda1"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"stability", " ", RowBox[{"criterion", ":", " ", RowBox[{"|", "\[Omega]1", "|", " ", RowBox[{"<", " ", RowBox[{"1.110468", " ", "-", " ", "0"}]}]}], ":", " ", "stable"}]}], ",", " ", RowBox[{"1", ":", " ", "unstable"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"stabilitytest", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Im", "[", RowBox[{"\[Omega]1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "]"}], " ", "\[Equal]", "0"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Omega]1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "^", "2"}], "<", "1.233"}], ",", "0", ",", "1"}], "]"}], ",", " ", RowBox[{"-", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "Ny"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "Nx"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", "plot", "*)"}], "\[IndentingNewLine]", RowBox[{"ArrayPlot", "[", RowBox[{"stabilitytest", ",", RowBox[{"AxesLabel", "\[Rule]", "Automatic"}], ",", RowBox[{"ColorRules", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"1", " ", "\[Rule]", " ", "Blue"}], ",", " ", RowBox[{"0", " ", "\[Rule]", " ", "Red"}], ",", RowBox[{ RowBox[{"-", "1"}], "\[Rule]", " ", "Black"}]}], "}"}]}], ",", RowBox[{"FrameTicks", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Miny", ",", "Maxy"}], "}"}], ",", RowBox[{"{", RowBox[{"Minx", ",", "Maxx"}], "}"}]}], "}"}]}], ",", RowBox[{"DataReversed", "\[Rule]", " ", "True"}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"B1", ",", "E1"}], "}"}]}]}], "]"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", "\n", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{ "save", " ", "data", " ", "to", " ", "generate", " ", "the", " ", "E1xB1", " ", RowBox[{"plot", ":", " ", "E1B1Plots_Minx"}]}], "=", "*"}], ",", RowBox[{"Maxx", "=", "*"}], ",", RowBox[{"Miny", "=", "*"}], ",", RowBox[{"Maxy", "=", "*"}], ",", RowBox[{"beta", "=", "*"}], ",", RowBox[{"eta", "=", RowBox[{"*", RowBox[{".", "txt"}]}]}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Export", "[", RowBox[{ RowBox[{"\"\\"", "<>", RowBox[{"ToString", "[", "Minx", "]"}], "<>", "\"\<,Maxx=\>\"", "<>", RowBox[{"ToString", "[", "Maxx", "]"}], "<>", "\"\<,Miny=\>\"", "<>", RowBox[{"ToString", "[", "Miny", "]"}], "<>", "\"\<,Maxy=\>\"", "<>", RowBox[{"ToString", "[", "Maxy", "]"}], "<>", "\"\<,beta=\>\"", "<>", RowBox[{"ToString", "[", "\[Beta]", "]"}], "<>", "\"\<,eta=\>\"", "<>", RowBox[{"ToString", "[", "\[Eta]", "]"}], "<>", "\"\<.txt\>\""}], " ", ",", "stabilitytest"}], "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.8141172604613867`*^9, 3.8141172880342493`*^9}, { 3.814117344804783*^9, 3.8141173457022657`*^9}, {3.814117383929709*^9, 3.81411745305208*^9}, 3.814178882578498*^9, {3.81419387545437*^9, 3.814193906338717*^9}, {3.814193990564268*^9, 3.814194029833171*^9}, { 3.814194154058144*^9, 3.814194160717815*^9}, {3.814194227721949*^9, 3.814194227944449*^9}, {3.814194574472636*^9, 3.8141945917447433`*^9}, { 3.8142608332790728`*^9, 3.81426087708571*^9}, {3.8142609208643637`*^9, 3.81426092923744*^9}, {3.814260969260674*^9, 3.814260984641901*^9}, { 3.814261505617228*^9, 3.814261537533218*^9}, {3.81426192230507*^9, 3.814261936678113*^9}, {3.814287623076413*^9, 3.814287644202793*^9}, { 3.814287695730712*^9, 3.814287750385285*^9}, {3.814438245766509*^9, 3.814438257003475*^9}, {3.814438381958399*^9, 3.814438387886157*^9}, { 3.814438631753172*^9, 3.814438672396487*^9}, {3.814439830151916*^9, 3.814439831321644*^9}, {3.8144408661596127`*^9, 3.8144408761595097`*^9}, { 3.8144461606147833`*^9, 3.8144462168850517`*^9}, {3.814446574311119*^9, 3.8144465849437847`*^9}, {3.814447076865971*^9, 3.8144470781463213`*^9}, { 3.814459239255537*^9, 3.8144592660637407`*^9}, {3.8144630190626287`*^9, 3.814463020284399*^9}, {3.814463365290266*^9, 3.8144633662889767`*^9}, { 3.8145164597988443`*^9, 3.814516466812581*^9}, {3.8145167844751053`*^9, 3.814516786031784*^9}, {3.814516950779499*^9, 3.814516981551613*^9}, { 3.814517344000019*^9, 3.8145173461806*^9}, 3.814518289619438*^9, { 3.8145189025237417`*^9, 3.814518919251699*^9}, {3.814525216178446*^9, 3.814525217526616*^9}, {3.814528877865128*^9, 3.814528878844674*^9}, { 3.814528924664089*^9, 3.8145289255991383`*^9}, {3.814726938987646*^9, 3.8147269400625668`*^9}, {3.814727068080759*^9, 3.8147270682188797`*^9}, { 3.814727143168509*^9, 3.814727149034535*^9}, {3.814731177252082*^9, 3.814731178283383*^9}, {3.8147314544340982`*^9, 3.814731455435177*^9}, 3.8147352586591463`*^9, {3.8149556573652763`*^9, 3.8149556963111563`*^9}, { 3.814955737818625*^9, 3.814955752204784*^9}, {3.814956303378141*^9, 3.814956357139718*^9}, {3.81495641092752*^9, 3.814956469997801*^9}, { 3.814956516455523*^9, 3.814956556701309*^9}, {3.8149566221787233`*^9, 3.81495662373589*^9}, {3.81495787756146*^9, 3.81495788505088*^9}, { 3.8149579985935497`*^9, 3.814958002138307*^9}, {3.814958223008951*^9, 3.8149582299329967`*^9}, {3.814963342703611*^9, 3.81496334412704*^9}, { 3.814963539482267*^9, 3.814963544253888*^9}, {3.814963702506652*^9, 3.814963703441136*^9}, 3.814971089264968*^9, {3.81497191507521*^9, 3.814971916025351*^9}, {3.814972326824293*^9, 3.814972327910605*^9}, { 3.814977236483583*^9, 3.814977237556244*^9}, {3.8149773360228577`*^9, 3.814977337177644*^9}, {3.814980972159051*^9, 3.8149809733054323`*^9}, { 3.8149822666393557`*^9, 3.81498226759524*^9}, {3.815068350307578*^9, 3.815068352466969*^9}, {3.8150781448340893`*^9, 3.815078192248948*^9}, { 3.815120825439077*^9, 3.815120827027547*^9}, 3.815120888623952*^9, { 3.815216275244216*^9, 3.815216302953288*^9}, {3.815217241078273*^9, 3.815217245413155*^9}, {3.815223203163213*^9, 3.8152232042759333`*^9}, { 3.8152234951189337`*^9, 3.8152234982547197`*^9}, {3.815232833599586*^9, 3.815232837930019*^9}, {3.815233418056861*^9, 3.815233428385111*^9}, { 3.8154724469330883`*^9, 3.815472504160824*^9}, {3.8154728254509373`*^9, 3.815472904226865*^9}, 3.815473753406251*^9, {3.815474458687583*^9, 3.8154744597431593`*^9}, {3.8154744926850653`*^9, 3.815474494185527*^9}, { 3.815474736555834*^9, 3.8154747368955517`*^9}, {3.8154752183959627`*^9, 3.815475220021121*^9}, {3.815475655349874*^9, 3.815475659309602*^9}, { 3.8154759834971647`*^9, 3.8154759875332336`*^9}, {3.81547898605225*^9, 3.815479006334119*^9}, {3.815481282898404*^9, 3.815481287150217*^9}, { 3.815483429593411*^9, 3.815483432238656*^9}, {3.8154837645962048`*^9, 3.815483766654634*^9}, {3.8155042741524057`*^9, 3.8155043000687113`*^9}, { 3.81550437404268*^9, 3.8155043797117243`*^9}, {3.816195870862892*^9, 3.816195874558445*^9}, {3.816196632829485*^9, 3.816196633912835*^9}, { 3.8161966697386913`*^9, 3.816196672804494*^9}, {3.816196740777363*^9, 3.816196742144246*^9}, {3.8162519276847887`*^9, 3.816252028298189*^9}, { 3.816252092376397*^9, 3.816252095031682*^9}, {3.8162522523946333`*^9, 3.81625225374443*^9}, {3.816256769046628*^9, 3.8162567726791353`*^9}, { 3.8199905469144*^9, 3.819990659503131*^9}, {3.819990693162197*^9, 3.819990710431182*^9}, {3.819990764818686*^9, 3.819990789064493*^9}, { 3.819991613266933*^9, 3.81999166597361*^9}, {3.819995260428939*^9, 3.819995261092719*^9}, {3.819995629160625*^9, 3.819995650175201*^9}, { 3.819996070955283*^9, 3.819996083588334*^9}, {3.819996161592619*^9, 3.819996164187941*^9}, {3.819996512548031*^9, 3.819996513264737*^9}, { 3.819996544496152*^9, 3.819996561321454*^9}, {3.8199965936964407`*^9, 3.819996594700782*^9}, {3.819996702170404*^9, 3.819996710057496*^9}, { 3.819999811494219*^9, 3.819999811917914*^9}, {3.8200047189526253`*^9, 3.820004719055056*^9}, {3.8200075666096163`*^9, 3.820007575678179*^9}, { 3.8200077570507393`*^9, 3.82000775749076*^9}, {3.820008013304978*^9, 3.820008022675972*^9}, {3.82000811930355*^9, 3.820008125084148*^9}, { 3.820008191780325*^9, 3.8200081965767117`*^9}, {3.8200110323926973`*^9, 3.820011037580586*^9}, {3.820014200995708*^9, 3.8200142022331533`*^9}, { 3.8200167787491693`*^9, 3.820016778830351*^9}, {3.820055009971817*^9, 3.820055010262191*^9}, {3.820057756222055*^9, 3.820057757461206*^9}, { 3.8200606696251097`*^9, 3.8200606701136017`*^9}, {3.8200651601927423`*^9, 3.8200651607692738`*^9}, {3.8200682929450817`*^9, 3.820068294182096*^9}, { 3.820078496939695*^9, 3.820078502533617*^9}, 3.8200827713255987`*^9, 3.820082924822163*^9, {3.820086812016876*^9, 3.820086812224516*^9}, { 3.820090451511016*^9, 3.82009045325327*^9}, {3.8200943355464077`*^9, 3.8200943358866997`*^9}, {3.820231811027454*^9, 3.820231852237211*^9}, { 3.820231883505843*^9, 3.82023189634659*^9}, {3.820233399250886*^9, 3.82023340147466*^9}, {3.820237669960801*^9, 3.8202376704874*^9}, { 3.8202430848503733`*^9, 3.820243089860195*^9}, {3.820247673739616*^9, 3.8202476739711847`*^9}, {3.8206678416278477`*^9, 3.8206679312371063`*^9}, {3.820667989552594*^9, 3.820668153707406*^9}, { 3.820668193908114*^9, 3.820668206480982*^9}, {3.82067048699671*^9, 3.8206704923642282`*^9}, {3.820683462032484*^9, 3.820683517273179*^9}, { 3.820683620210383*^9, 3.82068363514878*^9}, {3.820683877992989*^9, 3.820683889469738*^9}, {3.820684161625609*^9, 3.820684162972384*^9}, { 3.8206842222371273`*^9, 3.820684227629136*^9}, {3.8206843025502357`*^9, 3.820684304025134*^9}, {3.820684395397388*^9, 3.8206844112518806`*^9}, { 3.820687699241098*^9, 3.820687727403624*^9}, {3.8206897624353228`*^9, 3.820689783610249*^9}, {3.820693028720272*^9, 3.8206930289799137`*^9}, { 3.820696709931493*^9, 3.82069671209208*^9}, {3.820699093993175*^9, 3.82069909643738*^9}, {3.8207023043123627`*^9, 3.820702305829411*^9}, { 3.820704633009659*^9, 3.8207046374588833`*^9}, {3.8207454516244097`*^9, 3.820745456904602*^9}, {3.8207455272091227`*^9, 3.820745528383918*^9}, { 3.820745575019672*^9, 3.820745593872501*^9}, {3.82075227489986*^9, 3.8207522759134693`*^9}, {3.820753748955565*^9, 3.8207537495708313`*^9}, { 3.8209270595434103`*^9, 3.820927100690844*^9}, {3.820927170571104*^9, 3.8209271762454042`*^9}, {3.82092720687237*^9, 3.820927208477861*^9}, { 3.820927239445756*^9, 3.820927375967472*^9}, {3.8209274154525127`*^9, 3.820927539838739*^9}, {3.8209275801237173`*^9, 3.820927583789163*^9}, { 3.8209276618477297`*^9, 3.820927930940485*^9}, {3.820928038941081*^9, 3.8209280805245132`*^9}, {3.8209281373605213`*^9, 3.820928160099781*^9}, { 3.820928287691238*^9, 3.8209283320248003`*^9}, 3.8209283976822777`*^9, { 3.8209338510499*^9, 3.8209338727358418`*^9}, {3.820933921221362*^9, 3.820934438346562*^9}, {3.820934495143713*^9, 3.820934586667233*^9}, { 3.8209346405596867`*^9, 3.8209346782457323`*^9}, 3.8209347328067303`*^9, { 3.820934878618486*^9, 3.820934880351369*^9}, {3.820934961317541*^9, 3.8209349647827044`*^9}, {3.8209351012511663`*^9, 3.82093513069151*^9}, { 3.8209364613963013`*^9, 3.8209364736547813`*^9}, {3.820936709511107*^9, 3.820936717446322*^9}, {3.820936878769454*^9, 3.820936903772637*^9}, { 3.820937013832736*^9, 3.820937016215063*^9}, {3.820937101491418*^9, 3.820937104467169*^9}, {3.820940208266562*^9, 3.820940218228506*^9}, { 3.821206234365787*^9, 3.821206236859734*^9}, {3.821207168183267*^9, 3.821207174301566*^9}, {3.821207232531516*^9, 3.821207274537119*^9}, { 3.821207311831002*^9, 3.821207324211009*^9}, {3.821207459218523*^9, 3.821207460507162*^9}, {3.82121180336897*^9, 3.821211804613461*^9}, { 3.8212119802242737`*^9, 3.821211980999337*^9}, {3.821274132788734*^9, 3.821274132920595*^9}, {3.8212758742988777`*^9, 3.821276156048799*^9}, { 3.821292273239749*^9, 3.821292297847073*^9}, {3.8212923998137207`*^9, 3.8212925182306557`*^9}, {3.821293577391362*^9, 3.821293616427363*^9}, 3.8212955748290377`*^9, {3.8212962289058313`*^9, 3.821296229781703*^9}, { 3.821296389186611*^9, 3.821296391169788*^9}, {3.821303876336019*^9, 3.8213038798918257`*^9}, {3.821349833221972*^9, 3.821349833539125*^9}, { 3.82135177874067*^9, 3.8213517799287167`*^9}, {3.8213622581676083`*^9, 3.821362262528818*^9}, {3.821363388567175*^9, 3.821363390600565*^9}, { 3.821436970209441*^9, 3.8214369723818827`*^9}, {3.821440589932901*^9, 3.821440590862204*^9}, {3.8223953816411753`*^9, 3.822395382000898*^9}, { 3.822395618687676*^9, 3.822395619930676*^9}, {3.822407945212399*^9, 3.822407946482131*^9}, {3.82256977011626*^9, 3.82256977799228*^9}, { 3.822569951040628*^9, 3.822569957393498*^9}, {3.82257818146796*^9, 3.822578183164061*^9}, {3.8225814303460693`*^9, 3.822581436393128*^9}, { 3.822581552311409*^9, 3.8225815531233177`*^9}, {3.822581739202653*^9, 3.822581739537462*^9}, {3.822581872131124*^9, 3.822581872533125*^9}, { 3.8225820681084747`*^9, 3.822582068536951*^9}, {3.822582164965891*^9, 3.8225821703715487`*^9}, {3.822583056951191*^9, 3.822583058962035*^9}, { 3.8225844976723633`*^9, 3.822584497885995*^9}, {3.8225857576126537`*^9, 3.822585759205029*^9}, {3.8225877771658382`*^9, 3.822587778236129*^9}, { 3.822589234387109*^9, 3.822589235985677*^9}, {3.822589483452086*^9, 3.822589534668488*^9}, {3.8225895703644876`*^9, 3.822589573015369*^9}, { 3.822590031223529*^9, 3.822590049087798*^9}, {3.822591938667309*^9, 3.8225919413912067`*^9}, {3.823198955655016*^9, 3.823198956947657*^9}, { 3.823200251539215*^9, 3.823200252240623*^9}, {3.8232003306123943`*^9, 3.823200336032585*^9}, {3.823291420423072*^9, 3.823291424348111*^9}, { 3.8232919997215652`*^9, 3.823292081423339*^9}, {3.8232922953381433`*^9, 3.823292306986072*^9}, {3.823292396752125*^9, 3.823292397777382*^9}, { 3.823555221616041*^9, 3.8235552419239388`*^9}, {3.823555332568521*^9, 3.823555348668706*^9}, {3.8235570324125147`*^9, 3.8235570342026167`*^9}, { 3.8236115666192207`*^9, 3.8236115670515423`*^9}, {3.82429809668614*^9, 3.8242981090672293`*^9}, {3.824298177696061*^9, 3.82429817793465*^9}, { 3.8413250349134293`*^9, 3.8413250425952272`*^9}, 3.841557218314356*^9, { 3.841557548236739*^9, 3.841557651133493*^9}, {3.842333714829956*^9, 3.842333761239587*^9}, {3.8423338104883223`*^9, 3.842333824868032*^9}, { 3.842333862289102*^9, 3.8423340758878393`*^9}, {3.842334106437619*^9, 3.842334113988263*^9}, {3.842334160508748*^9, 3.84233417388593*^9}, { 3.842334252378809*^9, 3.8423342750946083`*^9}, {3.8423343127342377`*^9, 3.842334312931706*^9}, {3.842336616299777*^9, 3.8423366216217937`*^9}, 3.8423367104369507`*^9, 3.842341933733156*^9, {3.8423419942162952`*^9, 3.842342039600032*^9}, {3.843032795713039*^9, 3.843032798346848*^9}, 3.8430342413168488`*^9, {3.843635586512763*^9, 3.843635627067279*^9}, { 3.844176108806386*^9, 3.8441761143109818`*^9}, {3.8441761531922293`*^9, 3.844176197423791*^9}, {3.8441763056642103`*^9, 3.844176318081545*^9}, { 3.844513638067279*^9, 3.844513638824205*^9}, {3.8448605143002167`*^9, 3.844860533936371*^9}, {3.848757145749354*^9, 3.848757148127302*^9}, { 3.867750152575047*^9, 3.8677501780746202`*^9}, {3.8679199387516956`*^9, 3.8679199578571486`*^9}, {3.870171182010078*^9, 3.870171417090356*^9}, { 3.870172108165368*^9, 3.870172137353655*^9}, {3.87017618241955*^9, 3.870176182832943*^9}, {3.87017935194624*^9, 3.870179352203678*^9}, { 3.87024679582944*^9, 3.87024679752265*^9}, {3.87024698029574*^9, 3.8702469830294113`*^9}, {3.870339441584299*^9, 3.8703394557236156`*^9}, { 3.870339567645289*^9, 3.870339580001545*^9}, {3.870518202495804*^9, 3.870518225482174*^9}, {3.870518928091772*^9, 3.8705189454228573`*^9}, { 3.8705189860346203`*^9, 3.870519012693815*^9}, {3.870519490690826*^9, 3.870519513764144*^9}, {3.870519663236252*^9, 3.870519686200189*^9}},ExpressionUUID->"5119977a-4ed2-4fd0-a67b-\ 988e94b05652"], Cell[BoxData["\<\"ccrit,hs,ps =\"\>"], "Output", CellChangeTimes->{3.8430342766048603`*^9, 3.843035951903359*^9, 3.843043506834787*^9, 3.843635649578258*^9, 3.843636480753022*^9, 3.8436381621089487`*^9, 3.844176225072103*^9, 3.844261512105216*^9, 3.8448603892640467`*^9, 3.844860539188417*^9, 3.8679199913447905`*^9, 3.8679234648041754`*^9, 3.870172270428204*^9, 3.870176334654551*^9, 3.8703387619371033`*^9, 3.87033945629567*^9, 3.870520018901885*^9, 3.870524371633156*^9, 3.870775048372064*^9, 3.871046236896221*^9},ExpressionUUID->"f8977842-8d30-4a0f-8692-\ ce04a162e09b"], Cell[BoxData[ RowBox[{"{", RowBox[{"91.408825`", ",", GraphicsBox[RasterBox[CompressedData[" 1:eJzt0dEJgDAUQ9GCK7iAk7iDIwh+u7ojiOCPIFjRvBflFgpyoCaQbpyHqSml tPvdvs/P0uM4juM4/le/+07dB3/X1Tmq/+JHz86/Om79vubZ+bXu0uNrnp3/ 1F16uHp2/lN36eHq2fkqd+mR7dn50e7Sg31j3aUH+2rcpQf7xrpLD/bVuEsP do91lx7sq3GXHuwb6y492Ffj9e9XfJjbVw== "], {{0, 0}, {20, 20}}, {0, 1}], AxesLabel->FormBox["Automatic", TraditionalForm], Frame->Automatic, FrameLabel->{ FormBox["E1", TraditionalForm], FormBox["B1", TraditionalForm]}, FrameTicks->{{{1., FormBox["1.`", TraditionalForm]}, {10, FormBox["10", TraditionalForm]}}, {{0.55, FormBox["0.55`", TraditionalForm]}, {5, FormBox["5", TraditionalForm]}}}, GridLinesStyle->Directive[ GrayLevel[0.5, 0.4]], Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultPlotStyle" -> Automatic}]}], "}"}]], "Output", CellChangeTimes->{3.8430342766048603`*^9, 3.843035951903359*^9, 3.843043506834787*^9, 3.843635649578258*^9, 3.843636480753022*^9, 3.8436381621089487`*^9, 3.844176225072103*^9, 3.844261512105216*^9, 3.8448603892640467`*^9, 3.844860539188417*^9, 3.8679199913447905`*^9, 3.8679234648041754`*^9, 3.870172270428204*^9, 3.870176334654551*^9, 3.8703387619371033`*^9, 3.87033945629567*^9, 3.870520018901885*^9, 3.870524371633156*^9, 3.870775048372064*^9, 3.8710463283719053`*^9},ExpressionUUID->"3e04cccb-6e83-4cbf-8ce2-\ bb55241fd8af"], Cell[BoxData["\<\"E1B1Plots_Minx=0.55,Maxx=5,Miny=1.,Maxy=10,beta=1.,eta=7..\ txt\"\>"], "Output", CellChangeTimes->{3.8430342766048603`*^9, 3.843035951903359*^9, 3.843043506834787*^9, 3.843635649578258*^9, 3.843636480753022*^9, 3.8436381621089487`*^9, 3.844176225072103*^9, 3.844261512105216*^9, 3.8448603892640467`*^9, 3.844860539188417*^9, 3.8679199913447905`*^9, 3.8679234648041754`*^9, 3.870172270428204*^9, 3.870176334654551*^9, 3.8703387619371033`*^9, 3.87033945629567*^9, 3.870520018901885*^9, 3.870524371633156*^9, 3.870775048372064*^9, 3.871046328405675*^9},ExpressionUUID->"a46e903f-00b5-415c-92d8-\ 1b943aa30f80"] }, Open ]], Cell[BoxData[""], "Input", CellChangeTimes->{{3.8442583064582*^9, 3.844258323492936*^9}, { 3.870179485876191*^9, 3.870179533547607*^9}, {3.870180970648944*^9, 3.8701810003018084`*^9}, {3.8701813683511457`*^9, 3.870181466261743*^9}, { 3.8702458671878967`*^9, 3.870245875218087*^9}, 3.870246791087891*^9},ExpressionUUID->"6b37896f-2ee7-4eae-8aaf-\ 56ad297becab"] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1301, 744}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"11.2 for Linux x86 (64-bit) (September 10, 2017)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 711, 14, 346, "Title",ExpressionUUID->"553ed30a-6730-473d-8682-f67e849b952c"], Cell[CellGroupData[{ Cell[1316, 40, 697, 16, 341, "Section",ExpressionUUID->"a747abc8-3b73-4236-b76d-9f7e5a389046"], Cell[2016, 58, 3965, 96, 239, "Input",ExpressionUUID->"64d0e723-f54c-4c8f-907e-2dd284af16b7"] }, Open ]], Cell[CellGroupData[{ Cell[6018, 159, 350, 7, 68, "Section",ExpressionUUID->"8a8cbf01-963c-4b63-a4d5-d3c0812c0320"], Cell[CellGroupData[{ Cell[6393, 170, 4465, 117, 262, "Input",ExpressionUUID->"44499988-2afd-4ab3-bec9-71c5368b14bd"], Cell[10861, 289, 4933, 85, 72, "Output",ExpressionUUID->"4c270777-89bc-4627-805d-42019592b15e"], Cell[15797, 376, 6781, 145, 75, "Output",ExpressionUUID->"f3a2ff0b-d3d3-44b9-a6ad-c625cd4ac009"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[22627, 527, 508, 7, 68, "Section",ExpressionUUID->"9a1e7251-ec19-4baa-86eb-cecbf78ffcc2"], Cell[CellGroupData[{ Cell[23160, 538, 11263, 206, 285, "Input",ExpressionUUID->"c218540f-34cc-4e4f-bfe4-3d27be67ba86"], Cell[34426, 746, 17223, 469, 389, "Output",ExpressionUUID->"e230c138-949c-4822-a4a0-95eb5f08a269"], Cell[51652, 1217, 850, 17, 55, "Output",ExpressionUUID->"194cbf59-492e-415c-9c27-26abcbcad0f7"] }, Open ]], Cell[52517, 1237, 2407, 56, 239, "Input",ExpressionUUID->"287e0f4c-0953-4b8b-810f-69f8de8594ad"] }, Open ]], Cell[CellGroupData[{ Cell[54961, 1298, 292, 5, 65, "Subchapter",ExpressionUUID->"6ff83796-23ad-4212-8639-61206155d192"], Cell[CellGroupData[{ Cell[55278, 1307, 14683, 223, 354, "Input",ExpressionUUID->"4cb1ff4a-a751-482f-b501-ee2009ee5be2"], Cell[69964, 1532, 463, 7, 35, "Output",ExpressionUUID->"76f778b7-ad19-4460-9570-4122290b7c2f"], Cell[70430, 1541, 446, 7, 35, "Output",ExpressionUUID->"168eb232-a29e-437c-a1a2-915461a092a1"], Cell[70879, 1550, 460, 7, 35, "Output",ExpressionUUID->"8178d135-110f-48d2-8d24-bb2b83c38d68"], Cell[71342, 1559, 445, 7, 35, "Output",ExpressionUUID->"eeacdccc-f38a-49f5-bda9-4a4ccea9eb1b"], Cell[71790, 1568, 473, 7, 35, "Output",ExpressionUUID->"f19230cd-6fff-41d6-9914-8ac5d1b7e7a9"], Cell[72266, 1577, 447, 7, 35, "Output",ExpressionUUID->"133a27c2-1a89-447b-b1a3-55ee45f38e5e"], Cell[72716, 1586, 2456, 68, 92, "Output",ExpressionUUID->"e36be28a-ae22-4597-b095-8cce1c19edc1"], Cell[75175, 1656, 470, 7, 35, "Output",ExpressionUUID->"e77b46f3-a040-42aa-86d0-4a825e3c2005"], Cell[75648, 1665, 450, 7, 35, "Output",ExpressionUUID->"a01c009f-05c5-48ea-a6b9-28e15682a134"] }, Open ]], Cell[CellGroupData[{ Cell[76135, 1677, 352, 5, 68, "Section",ExpressionUUID->"eca8fcbe-e1d7-4292-acc1-708515933e00"], Cell[CellGroupData[{ Cell[76512, 1686, 17742, 341, 446, "Input",ExpressionUUID->"0f5a2bd7-6372-4d3a-a9e5-51ea3d601cab"], Cell[94257, 2029, 481, 9, 35, "Output",ExpressionUUID->"949a2f6a-e00e-4b68-9303-24e8d157f9a4"], Cell[94741, 2040, 481, 9, 35, "Output",ExpressionUUID->"ed50d2d5-ea87-41e9-a58b-5587d38d352f"], Cell[95225, 2051, 483, 9, 35, "Output",ExpressionUUID->"5b8e8da4-2463-40ca-9e58-466494018b78"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[95757, 2066, 775, 13, 68, "Section",ExpressionUUID->"21293870-4c70-42fd-883f-aec614437a87"], Cell[CellGroupData[{ Cell[96557, 2083, 41534, 943, 1942, "Input",ExpressionUUID->"5119977a-4ed2-4fd0-a67b-988e94b05652"], Cell[138094, 3028, 597, 9, 35, "Output",ExpressionUUID->"f8977842-8d30-4a0f-8692-ce04a162e09b"], Cell[138694, 3039, 1470, 32, 199, "Output",ExpressionUUID->"3e04cccb-6e83-4cbf-8ce2-bb55241fd8af"], Cell[140167, 3073, 647, 10, 35, "Output",ExpressionUUID->"a46e903f-00b5-415c-92d8-1b943aa30f80"] }, Open ]], Cell[140829, 3086, 375, 6, 31, "Input",ExpressionUUID->"6b37896f-2ee7-4eae-8aaf-56ad297becab"] }, Open ]] }, Open ]] }, Open ]] } ] *)