FixedPreferences.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    148972,       3444]
NotebookOptionsPosition[    139857,       3294]
NotebookOutlinePosition[    140306,       3311]
CellTagsIndexPosition[    140263,       3308]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[StyleBox["Model definition for fixed preferences", "Title"]], \
"Text",
 CellChangeTimes->{{3.858861697550729*^9, 3.858861700203829*^9}, {
  3.858862625159627*^9, 3.858862626175812*^9}, {3.869289107748027*^9, 
  3.869289113022397*^9}},ExpressionUUID->"bca58f23-7d8f-4b00-a95a-\
4a9c3983f9f4"],

Cell["\<\
Pablo Rosillo Rodes, 2022
Final Master Thesis\
\>", "Text",
 CellChangeTimes->{{3.8692889628899393`*^9, 
  3.869288976228018*^9}},ExpressionUUID->"eb1e4875-dc77-47c3-ae5c-\
03672dcc3d66"],

Cell["Rates definition:", "Text",
 CellChangeTimes->{{3.869289144557233*^9, 
  3.869289162386273*^9}},ExpressionUUID->"80d6cbf0-a7d4-4343-8444-\
c7375179ad8b"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"PX1Y1", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", "s1"}], ")"}], 
      RowBox[{"(", 
       RowBox[{"y1", "+", "y2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PX1Y2", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PX1X2", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PX2Y1", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PX2Y2", "=", 
     RowBox[{"s2", 
      RowBox[{"(", 
       RowBox[{"y1", "+", "y2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PX2X1", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY1X1", "=", " ", 
     RowBox[{"s1", 
      RowBox[{"(", 
       RowBox[{"x1", "+", "x2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY1Y2", "=", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY1X2", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY2X1", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY2X2", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", "s2"}], ")"}], 
      RowBox[{"(", 
       RowBox[{"x1", "+", "x2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PY2Y1", "=", "0"}], ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.860152513534335*^9, 3.8601525240418367`*^9}, 
   3.860152577446747*^9, {3.8601526083390827`*^9, 3.8601526128924093`*^9}, {
   3.8601527804483137`*^9, 3.860152788167676*^9}, 3.8601532350946493`*^9, {
   3.861081484011465*^9, 3.861081529463724*^9}, {3.8612507634575377`*^9, 
   3.861250830599481*^9}, {3.861251520782751*^9, 3.861251530824423*^9}, {
   3.8612515771898823`*^9, 3.861251607010626*^9}, {3.861254600399337*^9, 
   3.8612546062120533`*^9}, {3.861255525052267*^9, 3.861255552469657*^9}, {
   3.861692700523755*^9, 3.8616927216708612`*^9}, {3.8616953784001493`*^9, 
   3.8616954101472*^9}, {3.861698250929429*^9, 3.8616982906458607`*^9}, 
   3.861784002947122*^9, 3.86928912043865*^9},
 CellLabel->
  "In[397]:=",ExpressionUUID->"af4b99ea-e871-4906-81b7-b78ed789b70b"],

Cell["Dynamical equations construction:", "Text",
 CellChangeTimes->{{3.869289165066123*^9, 
  3.869289171465263*^9}},ExpressionUUID->"85e8d96f-828b-43a9-b6a2-\
7bddc2a8570c"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"x1eq", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "x1"}], " ", 
    RowBox[{"(", 
     RowBox[{"PX1Y1", "+", "PX1Y2", "+", "PX1X2"}], ")"}]}], "+", 
   RowBox[{"x2", " ", "PX2X1"}], "+", 
   RowBox[{"y1", " ", "PY1X1"}], "+", 
   RowBox[{"y2", " ", "PY2X1"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"x2eq", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "x2"}], " ", 
    RowBox[{"(", 
     RowBox[{"PX2Y1", "+", "PX2Y2", "+", "PX2X1"}], ")"}]}], "+", 
   RowBox[{"x1", " ", "PX1X2"}], "+", 
   RowBox[{"y1", " ", "PY1X2"}], "+", 
   RowBox[{"y2", " ", "PY2X2"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"y1eq", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "y1"}], 
    RowBox[{"(", 
     RowBox[{"PY1X1", "+", "PY1Y2", "+", "PY1X2"}], ")"}]}], "+", 
   RowBox[{"x1", " ", "PX1Y1"}], "+", 
   RowBox[{"x2", " ", "PX2Y1"}], "+", " ", 
   RowBox[{"y2", " ", "PY2Y1"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"y2eq", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "y2"}], 
    RowBox[{"(", 
     RowBox[{"PY2X1", "+", "PY2X2", "+", "PY2Y1"}], ")"}]}], "+", 
   RowBox[{"x1", " ", "PX1Y2"}], "+", 
   RowBox[{"x2", " ", "PX2Y2"}], "+", 
   RowBox[{"y1", " ", "PY1Y2"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"x1eq", "+", "x2eq", "+", "y1eq", "+", "y2eq"}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"x1eq", "+", "x2eq"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"y1eq", "+", "y2eq"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"x1eq", "+", "y1eq"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"x2eq", "+", "y2eq"}], "]"}]}], "Input",
 CellChangeTimes->{{3.858744968604349*^9, 3.858745129063846*^9}, {
   3.858745257287199*^9, 3.85874537299953*^9}, 3.8588619245278254`*^9, {
   3.858862125331993*^9, 3.8588621443545647`*^9}, {3.858862375574929*^9, 
   3.8588623764170923`*^9}, {3.8588626820069637`*^9, 3.858862688995494*^9}, {
   3.85886295707443*^9, 3.858863218022505*^9}, 3.858867819295282*^9, 
   3.858867994243079*^9, {3.8588680550918617`*^9, 3.858868067098027*^9}, {
   3.8588681021166477`*^9, 3.8588681375656977`*^9}, {3.858868473848609*^9, 
   3.858868476040607*^9}, {3.8588685233856983`*^9, 3.858868525419894*^9}, {
   3.858868567718862*^9, 3.858868655706195*^9}, {3.861698297941472*^9, 
   3.861698298713286*^9}, {3.8616986734586067`*^9, 3.861698674121381*^9}, {
   3.863496646508053*^9, 3.863496667640366*^9}, {3.863511186784711*^9, 
   3.86351118993181*^9}, {3.86351140493859*^9, 3.863511415673964*^9}, {
   3.8700579113098583`*^9, 3.8700579350679073`*^9}},
 CellLabel->
  "In[409]:=",ExpressionUUID->"551d2428-46e7-4cf5-ba2c-362ae31e2494"],

Cell[BoxData[
 RowBox[{
  RowBox[{"s1", " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.8716998667718973`*^9},
 CellLabel->
  "Out[409]=",ExpressionUUID->"5a60df6f-acfe-4c9c-ae6b-b8c0191d2e8c"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "s2"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], "-", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.871699866772752*^9},
 CellLabel->
  "Out[410]=",ExpressionUUID->"08069168-1955-4769-b46b-798dfcfe30cd"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "s1"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.8716998667735023`*^9},
 CellLabel->
  "Out[411]=",ExpressionUUID->"1d4366de-e75d-4326-bae4-c774f93fc4ac"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", "s2"}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], ")"}]}], "+", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.871699866774247*^9},
 CellLabel->
  "Out[412]=",ExpressionUUID->"e2d9cdad-d17b-4e39-bb09-9a0588625511"],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.8716998667749367`*^9},
 CellLabel->
  "Out[413]=",ExpressionUUID->"4ca77c4f-682a-417d-ab7b-7ac4d530869a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"s1", " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s2"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}], "-", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.871699866775674*^9},
 CellLabel->
  "Out[414]=",ExpressionUUID->"dbfdf09d-6bdd-4001-9805-fe1e3c65b2d8"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "s1"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s2"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}], "+", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.8716998667783337`*^9},
 CellLabel->
  "Out[415]=",ExpressionUUID->"20613bae-7fa8-437e-862b-4979a5f7489d"],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.8716998667790728`*^9},
 CellLabel->
  "Out[416]=",ExpressionUUID->"a2d353e1-2759-4b04-a11e-f2337db561f9"],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.863496657604657*^9, 3.863496668023798*^9}, 
   3.863499224802866*^9, 3.86350010810042*^9, 3.86351119308646*^9, 
   3.863511410488317*^9, 3.8643543618500223`*^9, 3.864358254374349*^9, 
   3.8643674157817307`*^9, 3.864369228064681*^9, 3.8648789160124693`*^9, 
   3.864878962550952*^9, 3.864879144418295*^9, 3.8649752567365513`*^9, 
   3.865059056483795*^9, 3.865141740033313*^9, 3.8659118620690117`*^9, 
   3.865912553460565*^9, 3.866023524068924*^9, 3.866023596474749*^9, {
   3.866024632255033*^9, 3.866024660804166*^9}, {3.866446427980405*^9, 
   3.866446455546529*^9}, 3.866625197273931*^9, 3.867404804635585*^9, {
   3.867404849751561*^9, 3.867404875386568*^9}, 3.867467477470193*^9, 
   3.867467517387738*^9, 3.867579887644443*^9, 3.867752506042671*^9, 
   3.867764725625094*^9, 3.867924151760211*^9, 3.8685118115111837`*^9, 
   3.8687891577952957`*^9, 3.869289266190576*^9, 3.869369125365221*^9, {
   3.87005793769543*^9, 3.8700579436532707`*^9}, 3.871210179614257*^9, 
   3.871210543970446*^9, 3.871699866779707*^9},
 CellLabel->
  "Out[417]=",ExpressionUUID->"1d114ae8-ab81-4ce3-977e-0edfef173971"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["x1eq"], "Input",
 CellChangeTimes->{{3.8588684884444838`*^9, 3.85886848935054*^9}},
 CellLabel->
  "In[418]:=",ExpressionUUID->"2983a36d-74b9-41f8-85c8-4ccf5f2daead"],

Cell[BoxData[
 RowBox[{
  RowBox[{"s1", " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "-", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.858868490171268*^9, {3.858868660180188*^9, 3.858868674469836*^9}, 
   3.858868735863825*^9, 3.858868925409233*^9, 3.858869020944479*^9, 
   3.858931968336874*^9, 3.860152466932103*^9, 3.860152535196206*^9, 
   3.860152625537751*^9, 3.860152800223133*^9, 3.860153107905122*^9, 
   3.860153172237775*^9, 3.8601532430635653`*^9, 3.860153441777686*^9, 
   3.8601535689111443`*^9, {3.8601537759527493`*^9, 3.860153805424054*^9}, 
   3.860173353789002*^9, 3.860387951820673*^9, 3.860405542383161*^9, 
   3.860406029071403*^9, 3.860645193806717*^9, 3.861081667779497*^9, 
   3.861082653849979*^9, 3.861250865190516*^9, 3.861251616697472*^9, 
   3.861253705544664*^9, 3.8612546534491673`*^9, 3.861255594815064*^9, 
   3.861255709431615*^9, 3.861692738736107*^9, 3.8616954219168463`*^9, 
   3.861698513301406*^9, 3.861698590949222*^9, 3.8616986528481483`*^9, {
   3.8616986877138853`*^9, 3.86169871702894*^9}, 3.8616989783109827`*^9, 
   3.8616997125570517`*^9, 3.861699808810865*^9, 3.861700015520158*^9, 
   3.861771545690524*^9, 3.861771679224846*^9, 3.861780313814085*^9, 
   3.861780407816905*^9, 3.861780589177126*^9, 3.861781539507964*^9, 
   3.861782265946273*^9, 3.861791543656928*^9, 3.862057671384601*^9, 
   3.862416478419366*^9, 3.862477653596377*^9, 3.863065105439657*^9, {
   3.863086979689714*^9, 3.8630870051919737`*^9}, 3.863087366949091*^9, 
   3.863145274757678*^9, 3.863146310721583*^9, 3.8634177191888523`*^9, 
   3.863495570631689*^9, 3.863499224815494*^9, 3.863500108102892*^9, 
   3.863511193092516*^9, 3.864354361861142*^9, 3.86435825439554*^9, 
   3.864367415791524*^9, 3.864369228087412*^9, 3.8648789160322104`*^9, 
   3.864878962556999*^9, 3.8648791444345407`*^9, 3.864975256744266*^9, 
   3.865059056491323*^9, 3.865141740039858*^9, 3.8659118620754633`*^9, 
   3.865912553469159*^9, 3.866023524077485*^9, 3.866023596484761*^9, {
   3.866024632263612*^9, 3.8660246608124647`*^9}, {3.866446427994298*^9, 
   3.866446455562202*^9}, 3.8666251972800694`*^9, 3.867404804643094*^9, {
   3.867404849758863*^9, 3.86740487539393*^9}, 3.8674674774782467`*^9, 
   3.867467517395804*^9, 3.8675798876492233`*^9, 3.867752506047956*^9, 
   3.867764725630232*^9, 3.867924151766988*^9, 3.868511811523356*^9, 
   3.868789157800782*^9, 3.869289266196362*^9, 3.869369125379932*^9, 
   3.871210179623445*^9, 3.871210543995043*^9, 3.8716998667898083`*^9},
 CellLabel->
  "Out[418]=",ExpressionUUID->"e48602e6-f464-45d2-9b32-bd98475db257"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Xsimp", "[", 
   RowBox[{"x1_", ",", " ", "x2_", ",", " ", "y1_", ",", " ", "y2_"}], "]"}], 
  " ", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"x1eq", "+", "x2eq"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.85874516208361*^9, 3.8587452345144043`*^9}, {
   3.858745376663151*^9, 3.858745388263069*^9}, {3.858745446766433*^9, 
   3.858745473106357*^9}, {3.858745521363174*^9, 3.858745545673217*^9}, {
   3.858745589778235*^9, 3.858745591286046*^9}, {3.858745639593884*^9, 
   3.858745650271687*^9}, {3.85875400221306*^9, 3.858754007948798*^9}, {
   3.858861519121765*^9, 3.858861546269418*^9}, {3.858861581456911*^9, 
   3.858861630993597*^9}, {3.8588616617497463`*^9, 3.8588616629543324`*^9}, {
   3.858861741829862*^9, 3.858861798551693*^9}, {3.858861834780342*^9, 
   3.85886184028736*^9}, {3.858861932525635*^9, 3.85886195035944*^9}, {
   3.858861985438102*^9, 3.858861995133304*^9}, {3.858862029926497*^9, 
   3.8588620322074013`*^9}, 3.8588622080859127`*^9, 3.85886229597012*^9, 
   3.858862533087508*^9, 3.8588626010009127`*^9, 3.858863232327042*^9, 
   3.858866407604615*^9},
 CellLabel->
  "In[419]:=",ExpressionUUID->"cf407a73-21af-4f5b-a822-837fcf26ebef"],

Cell[BoxData[
 RowBox[{
  RowBox[{"s1", " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s2"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}], "-", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.858866410923758*^9, {3.858866698953274*^9, 3.858866725010623*^9}, 
   3.858867204421648*^9, {3.8588672573256702`*^9, 3.858867305024558*^9}, 
   3.858867393664032*^9, {3.8588674268460827`*^9, 3.858867522201447*^9}, 
   3.8588677702039127`*^9, {3.8588678065875196`*^9, 3.85886784878653*^9}, {
   3.858867896720846*^9, 3.858867915774563*^9}, 3.85886799829993*^9, {
   3.858868070280911*^9, 3.8588681412123327`*^9}, 3.858868481189891*^9, {
   3.8588686601847563`*^9, 3.858868674484004*^9}, 3.858868735875943*^9, 
   3.8588689254237556`*^9, 3.8588690209499283`*^9, 3.858931968355131*^9, 
   3.860152466939705*^9, 3.860152535203041*^9, 3.860152625547406*^9, 
   3.8601528002296333`*^9, 3.860153107914216*^9, 3.860153172250051*^9, 
   3.8601532430703297`*^9, 3.8601534417806787`*^9, 3.860153568922249*^9, {
   3.860153775968973*^9, 3.860153805433679*^9}, 3.860173353796253*^9, 
   3.860387951828699*^9, 3.860405542407424*^9, 3.860406029074707*^9, 
   3.860645193820654*^9, 3.8610816678026667`*^9, 3.861082653852519*^9, 
   3.8612508652046347`*^9, 3.861251616709269*^9, 3.861253705555697*^9, 
   3.8612546534622297`*^9, 3.861255594833562*^9, 3.86125570943652*^9, 
   3.8616927387513103`*^9, 3.861695421935451*^9, 3.861698513312909*^9, 
   3.861698590965651*^9, 3.8616986528527193`*^9, {3.861698687731934*^9, 
   3.8616987170368357`*^9}, 3.861698978322411*^9, 3.8616997125736017`*^9, 
   3.861699808825103*^9, 3.861700015536182*^9, 3.86177154571754*^9, 
   3.86177167922991*^9, 3.8617803138416023`*^9, 3.861780407835088*^9, 
   3.861780589200089*^9, 3.861781539523497*^9, 3.86178226597285*^9, 
   3.8617915436681623`*^9, 3.862057671401074*^9, 3.862416478432631*^9, 
   3.8624776536016293`*^9, 3.8630651054511414`*^9, {3.86308697970817*^9, 
   3.863087005206657*^9}, 3.863087366963395*^9, 3.863145274779375*^9, 
   3.863146310737061*^9, 3.863417719205171*^9, 3.863495570647187*^9, 
   3.863499224828629*^9, 3.863500108109395*^9, 3.86351119310186*^9, 
   3.8643543618733387`*^9, 3.864358254405018*^9, 3.8643674158069963`*^9, 
   3.8643692281052427`*^9, 3.864878916046434*^9, 3.864878962567697*^9, 
   3.864879144447341*^9, 3.864975256757934*^9, 3.865059056505265*^9, 
   3.8651417400582027`*^9, 3.86591186209156*^9, 3.8659125534799223`*^9, 
   3.866023524093981*^9, 3.8660235965016212`*^9, {3.866024632280525*^9, 
   3.866024660824651*^9}, {3.8664464280048723`*^9, 3.866446455566818*^9}, 
   3.866625197293044*^9, 3.867404804657714*^9, {3.867404849773816*^9, 
   3.8674048754077387`*^9}, 3.867467477492612*^9, 3.867467517402068*^9, 
   3.86757988766543*^9, 3.867752506061344*^9, 3.867764725643526*^9, 
   3.8679241517855883`*^9, 3.868511811538506*^9, 3.86878915781292*^9, 
   3.8692892662069063`*^9, 3.869369125392482*^9, 3.871210179626173*^9, 
   3.871210543998548*^9, 3.871699866792897*^9},
 CellLabel->
  "Out[419]=",ExpressionUUID->"6d3dde4e-80a4-44c6-b8eb-aa2011b4ae75"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Xeq", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Xsimp", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"X", "/", "2"}], "+", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{"X", "/", "2"}], "-", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.858862000538527*^9, 3.858862112738143*^9}, {
   3.858862159804182*^9, 3.858862170630928*^9}, {3.8588622031948433`*^9, 
   3.858862217217596*^9}, 3.8588622927070723`*^9, {3.858862335189515*^9, 
   3.8588623713715487`*^9}, {3.8588624030121393`*^9, 3.858862423920089*^9}, 
   3.858862534571671*^9, 3.85886259983747*^9, 3.858863233722219*^9, {
   3.858866456946732*^9, 3.8588664662466383`*^9}, 3.858866707258143*^9, {
   3.858867014532146*^9, 3.858867017095325*^9}, 3.858867273636587*^9, {
   3.858867362804555*^9, 3.858867366198247*^9}, {3.86015302495588*^9, 
   3.86015304226884*^9}, {3.8616982183250103`*^9, 3.861698233491702*^9}, {
   3.861698315387947*^9, 3.861698350597335*^9}, {3.863145139171679*^9, 
   3.863145154285096*^9}, {3.863500096777033*^9, 3.863500111311249*^9}},
 CellLabel->
  "In[420]:=",ExpressionUUID->"2a458930-9ac6-47d4-a410-4dde926116e4"],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", 
     RowBox[{"(", 
      RowBox[{"s1", "-", "s2"}], ")"}], " ", 
     SuperscriptBox["X", "2"]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", "\[Omega]"}], 
    "+", 
    RowBox[{"X", " ", 
     RowBox[{"(", 
      RowBox[{"1", "+", "s1", "-", 
       RowBox[{"3", " ", "s2"}], "-", 
       RowBox[{"2", " ", "\[Alpha]"}], "+", 
       RowBox[{"2", " ", "s1", " ", "\[Alpha]"}], "+", 
       RowBox[{"2", " ", "s2", " ", "\[Alpha]"}], "+", 
       RowBox[{"2", " ", "\[Omega]"}], "-", 
       RowBox[{"2", " ", "s1", " ", "\[Omega]"}], "-", 
       RowBox[{"2", " ", "s2", " ", "\[Omega]"}]}], ")"}]}]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{{3.858867282005418*^9, 3.858867305028775*^9}, 
   3.858867393679269*^9, {3.8588674268533707`*^9, 3.8588675222112093`*^9}, 
   3.858867770207775*^9, {3.858867806590564*^9, 3.858867848790469*^9}, {
   3.858867896735847*^9, 3.858867915779181*^9}, 3.858867998307207*^9, {
   3.858868070294032*^9, 3.858868141224883*^9}, 3.858868481203677*^9, {
   3.858868660196684*^9, 3.858868674495759*^9}, 3.858868735879496*^9, 
   3.858868925427186*^9, 3.8588690209631567`*^9, 3.858931968358777*^9, 
   3.860152466950375*^9, 3.8601525352164516`*^9, 3.8601526255500927`*^9, 
   3.860152800240898*^9, 3.86015310791783*^9, 3.8601531722531977`*^9, 
   3.860153243079257*^9, 3.860153441792322*^9, 3.860153568925741*^9, {
   3.8601537759738417`*^9, 3.860153805437243*^9}, 3.860173353808614*^9, 
   3.860387951832019*^9, 3.8604055424112997`*^9, 3.860406029085499*^9, 
   3.860645193824031*^9, 3.8610816678152313`*^9, 3.8610826538553877`*^9, 
   3.8612508652073812`*^9, 3.8612516167124243`*^9, 3.861253705565742*^9, 
   3.861254653465076*^9, 3.861255594836357*^9, 3.861255709450314*^9, 
   3.861692738767109*^9, 3.861695421949988*^9, 3.861698513317726*^9, 
   3.861698590970326*^9, 3.861698652858292*^9, {3.86169868774832*^9, 
   3.861698717050139*^9}, 3.861698978338728*^9, 3.861699712588023*^9, 
   3.861699808839128*^9, 3.861700015553525*^9, 3.8617715457429323`*^9, 
   3.861771679245329*^9, 3.8617803138613987`*^9, 3.86178040785034*^9, 
   3.86178058921723*^9, 3.861781539544114*^9, 3.861782265996117*^9, 
   3.861791543680668*^9, 3.862057671417923*^9, 3.86241647844715*^9, 
   3.8624776536190453`*^9, 3.863065105465067*^9, {3.863086979723831*^9, 
   3.8630870052229013`*^9}, 3.863087366975773*^9, 3.863145274796938*^9, 
   3.8631463107522793`*^9, 3.86341771922013*^9, 3.8634955706587152`*^9, 
   3.863499224836419*^9, {3.863500108124824*^9, 3.863500113108425*^9}, 
   3.863511193115471*^9, 3.864354361886752*^9, 3.8643582544218397`*^9, 
   3.864367415824747*^9, 3.864369228118043*^9, 3.864878916056055*^9, 
   3.864878962581849*^9, 3.8648791444558764`*^9, 3.864975256773211*^9, 
   3.865059056520219*^9, 3.865141740076857*^9, 3.865911862105064*^9, 
   3.865912553493766*^9, 3.866023524105116*^9, 3.866023596516403*^9, {
   3.866024632291827*^9, 3.8660246608299*^9}, {3.8664464280172377`*^9, 
   3.86644645557504*^9}, 3.8666251973022337`*^9, 3.8674048046739683`*^9, {
   3.867404849789221*^9, 3.867404875423011*^9}, 3.8674674775114193`*^9, 
   3.867467517414691*^9, 3.867579887681788*^9, 3.8677525060757837`*^9, 
   3.867764725657322*^9, 3.867924151801195*^9, 3.8685118115550823`*^9, 
   3.8687891578280153`*^9, 3.869289266219664*^9, 3.869369125407488*^9, 
   3.871210179637514*^9, 3.871210544008057*^9, 3.871699866804178*^9},
 CellLabel->
  "Out[420]=",ExpressionUUID->"4b5076b9-042e-4de8-ae0b-53d7115a5600"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"\[Omega]simp", "[", 
     RowBox[{"x1_", ",", " ", "x2_", ",", " ", "y1_", ",", " ", "y2_"}], 
     "]"}], " ", "=", " ", 
    RowBox[{"Simplify", "[", 
     RowBox[{"x1eq", "-", "x2eq"}], "]"}]}], ";"}], "\n"}], "\n", 
 RowBox[{"\[Omega]eq", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"\[Omega]simp", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"X", "/", "2"}], "+", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{"X", "/", "2"}], "-", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}]}], "]"}], "]"}]}], "\n"}], "Input",
 CellChangeTimes->{{3.858862497684667*^9, 3.858862550995639*^9}, {
   3.858862591790657*^9, 3.8588625929918613`*^9}, {3.85886324267227*^9, 
   3.8588632445659637`*^9}, {3.8588664885327263`*^9, 3.8588664967520523`*^9}, 
   3.858866715625643*^9, {3.8588669977183523`*^9, 3.85886700049039*^9}, 
   3.858867279443709*^9, {3.858867344290998*^9, 3.85886739066299*^9}, {
   3.8601530488352737`*^9, 3.8601530510385647`*^9}, {3.8616983842233953`*^9, 
   3.861698411323135*^9}, {3.861698649765191*^9, 3.861698677853909*^9}, {
   3.863145150896463*^9, 3.8631451558937817`*^9}},
 CellLabel->
  "In[421]:=",ExpressionUUID->"1523053d-5e99-4475-a0ae-b61faed2e050"],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
     SuperscriptBox["X", "2"]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "-", "s2"}], ")"}], " ", "\[Omega]"}], 
    "+", 
    RowBox[{"X", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3"}], "+", "s1", "+", 
       RowBox[{"3", " ", "s2"}], "+", 
       RowBox[{"2", " ", "\[Alpha]"}], "+", 
       RowBox[{"2", " ", "s1", " ", "\[Alpha]"}], "-", 
       RowBox[{"2", " ", "s2", " ", "\[Alpha]"}], "-", 
       RowBox[{"2", " ", "s1", " ", "\[Omega]"}], "+", 
       RowBox[{"2", " ", "s2", " ", "\[Omega]"}]}], ")"}]}]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{{3.858867282032443*^9, 3.8588673050540247`*^9}, 
   3.85886735043851*^9, 3.858867393717901*^9, {3.858867426884988*^9, 
   3.858867522243573*^9}, 3.8588677702708673`*^9, {3.8588678066214323`*^9, 
   3.8588678488538446`*^9}, {3.8588678967974577`*^9, 3.8588679158115253`*^9}, 
   3.858867998368438*^9, {3.858868070336669*^9, 3.858868141254753*^9}, 
   3.858868481232396*^9, {3.8588686602505493`*^9, 3.858868674561034*^9}, 
   3.858868735940387*^9, 3.858868925476754*^9, 3.85886902101335*^9, 
   3.8589319684361887`*^9, 3.860152467020193*^9, 3.860152535239037*^9, 
   3.860152625613163*^9, 3.860152800283431*^9, 3.860153107965496*^9, 
   3.860153172315708*^9, 3.860153243144227*^9, 3.860153441855175*^9, 
   3.860153568996172*^9, {3.860153776062327*^9, 3.8601538055109463`*^9}, 
   3.860173353856845*^9, 3.860387951929304*^9, 3.860405542513476*^9, 
   3.8604060291452017`*^9, 3.8606451939344473`*^9, 3.86108166792449*^9, 
   3.8610826538776093`*^9, 3.861250865282031*^9, 3.8612516167858353`*^9, 
   3.861253705646572*^9, 3.861254653543522*^9, 3.8612555949176083`*^9, 
   3.8612557095358877`*^9, 3.861692738870607*^9, 3.861695422019733*^9, 
   3.861698513325014*^9, 3.8616985909819613`*^9, 3.86169865287298*^9, {
   3.861698687766081*^9, 3.8616987170684967`*^9}, 3.861698978356106*^9, 
   3.861699712604026*^9, 3.861699808859638*^9, 3.861700015571217*^9, 
   3.8617715457886333`*^9, 3.861771679262348*^9, 3.86178031388442*^9, 
   3.8617804078687277`*^9, 3.86178058923534*^9, 3.861781539566641*^9, 
   3.861782266026038*^9, 3.861791543697299*^9, 3.862057671434417*^9, 
   3.8624164784646187`*^9, 3.862477653646536*^9, 3.863065105483674*^9, {
   3.86308697974091*^9, 3.863087005239463*^9}, 3.8630873669910793`*^9, 
   3.863145274813877*^9, 3.863146310770735*^9, 3.863417719236714*^9, 
   3.8634955706751747`*^9, 3.8634992248464537`*^9, 3.863500108139492*^9, 
   3.8635111931336718`*^9, 3.864354361901305*^9, 3.864358254439834*^9, 
   3.864367415838118*^9, 3.8643692281459723`*^9, 3.864878916069276*^9, 
   3.86487896259686*^9, 3.864879144467458*^9, 3.864975256790559*^9, 
   3.865059056536964*^9, 3.865141740093924*^9, 3.8659118621151533`*^9, 
   3.865912553509789*^9, 3.866023524120275*^9, 3.866023596534788*^9, {
   3.866024632308351*^9, 3.8660246608438377`*^9}, {3.866446428035275*^9, 
   3.8664464555879917`*^9}, 3.866625197313993*^9, 3.86740480469063*^9, {
   3.8674048498046103`*^9, 3.867404875439085*^9}, 3.86746747752775*^9, 
   3.867467517427215*^9, 3.867579887698251*^9, 3.867752506093811*^9, 
   3.867764725673753*^9, 3.8679241518200808`*^9, 3.868511811571783*^9, 
   3.868789157844467*^9, 3.869289266235806*^9, 3.869369125425091*^9, 
   3.871210179652671*^9, 3.871210544020751*^9, 3.871699866814583*^9},
 CellLabel->
  "Out[422]=",ExpressionUUID->"5e2019ba-8543-4539-b7de-1aa51350ddf5"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Ysimp", "[", 
   RowBox[{"x1_", ",", " ", "x2_", ",", " ", "y1_", ",", " ", "y2_"}], "]"}], 
  " ", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"y1eq", "+", "y2eq"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.85874516208361*^9, 3.8587452345144043`*^9}, {
   3.858745376663151*^9, 3.858745388263069*^9}, {3.858745446766433*^9, 
   3.858745473106357*^9}, {3.858745521363174*^9, 3.858745545673217*^9}, {
   3.858745589778235*^9, 3.858745591286046*^9}, {3.858745639593884*^9, 
   3.858745650271687*^9}, {3.85875400221306*^9, 3.858754007948798*^9}, {
   3.858861519121765*^9, 3.858861546269418*^9}, {3.858861581456911*^9, 
   3.858861630993597*^9}, {3.8588616617497463`*^9, 3.8588616629543324`*^9}, {
   3.858861741829862*^9, 3.858861798551693*^9}, {3.858861834780342*^9, 
   3.85886184028736*^9}, {3.858861932525635*^9, 3.85886195035944*^9}, {
   3.858861985438102*^9, 3.858861995133304*^9}, {3.858862029926497*^9, 
   3.8588620322074013`*^9}, 3.8588622080859127`*^9, 3.85886229597012*^9, 
   3.858862533087508*^9, 3.8588626010009127`*^9, 3.858863232327042*^9, 
   3.858866407604615*^9, {3.8659124863387337`*^9, 3.865912491496647*^9}},
 CellLabel->
  "In[423]:=",ExpressionUUID->"5e8ec10b-9b57-4a69-bef3-09c7a4fadedf"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "s1"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y1"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s2"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], " ", "y2"}], "-", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s1"}], ")"}], " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}], "+", 
  RowBox[{"s2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"y1", "+", "y2"}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.858866410923758*^9, {3.858866698953274*^9, 3.858866725010623*^9}, 
   3.858867204421648*^9, {3.8588672573256702`*^9, 3.858867305024558*^9}, 
   3.858867393664032*^9, {3.8588674268460827`*^9, 3.858867522201447*^9}, 
   3.8588677702039127`*^9, {3.8588678065875196`*^9, 3.85886784878653*^9}, {
   3.858867896720846*^9, 3.858867915774563*^9}, 3.85886799829993*^9, {
   3.858868070280911*^9, 3.8588681412123327`*^9}, 3.858868481189891*^9, {
   3.8588686601847563`*^9, 3.858868674484004*^9}, 3.858868735875943*^9, 
   3.8588689254237556`*^9, 3.8588690209499283`*^9, 3.858931968355131*^9, 
   3.860152466939705*^9, 3.860152535203041*^9, 3.860152625547406*^9, 
   3.8601528002296333`*^9, 3.860153107914216*^9, 3.860153172250051*^9, 
   3.8601532430703297`*^9, 3.8601534417806787`*^9, 3.860153568922249*^9, {
   3.860153775968973*^9, 3.860153805433679*^9}, 3.860173353796253*^9, 
   3.860387951828699*^9, 3.860405542407424*^9, 3.860406029074707*^9, 
   3.860645193820654*^9, 3.8610816678026667`*^9, 3.861082653852519*^9, 
   3.8612508652046347`*^9, 3.861251616709269*^9, 3.861253705555697*^9, 
   3.8612546534622297`*^9, 3.861255594833562*^9, 3.86125570943652*^9, 
   3.8616927387513103`*^9, 3.861695421935451*^9, 3.861698513312909*^9, 
   3.861698590965651*^9, 3.8616986528527193`*^9, {3.861698687731934*^9, 
   3.8616987170368357`*^9}, 3.861698978322411*^9, 3.8616997125736017`*^9, 
   3.861699808825103*^9, 3.861700015536182*^9, 3.86177154571754*^9, 
   3.86177167922991*^9, 3.8617803138416023`*^9, 3.861780407835088*^9, 
   3.861780589200089*^9, 3.861781539523497*^9, 3.86178226597285*^9, 
   3.8617915436681623`*^9, 3.862057671401074*^9, 3.862416478432631*^9, 
   3.8624776536016293`*^9, 3.8630651054511414`*^9, {3.86308697970817*^9, 
   3.863087005206657*^9}, 3.863087366963395*^9, 3.863145274779375*^9, 
   3.863146310737061*^9, 3.863417719205171*^9, 3.863495570647187*^9, 
   3.863499224828629*^9, 3.863500108109395*^9, 3.86351119310186*^9, 
   3.8643543618733387`*^9, 3.864358254405018*^9, 3.8643674158069963`*^9, 
   3.8643692281052427`*^9, 3.864878916046434*^9, 3.864878962567697*^9, 
   3.864879144447341*^9, 3.864975256757934*^9, 3.865059056505265*^9, 
   3.8651417400582027`*^9, 3.86591186209156*^9, 3.865912553513988*^9, 
   3.866023524125662*^9, 3.866023596540655*^9, {3.866024632323143*^9, 
   3.8660246608483057`*^9}, {3.866446428039379*^9, 3.866446455592537*^9}, 
   3.866625197317623*^9, 3.867404804694906*^9, {3.867404849809228*^9, 
   3.867404875443219*^9}, 3.8674674775324783`*^9, 3.8674675174316187`*^9, 
   3.867579887702406*^9, 3.8677525060981817`*^9, 3.867764725678021*^9, 
   3.8679241518339863`*^9, 3.868511811575932*^9, 3.868789157848484*^9, 
   3.8692892662399397`*^9, 3.869369125429104*^9, 3.871210179655549*^9, 
   3.871210544026617*^9, 3.871699866823131*^9},
 CellLabel->
  "Out[423]=",ExpressionUUID->"f180403d-b79d-411f-ada9-b395aa5e00c2"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Yeq", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Ysimp", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"X", "/", "2"}], "+", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{"X", "/", "2"}], "-", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.858862000538527*^9, 3.858862112738143*^9}, {
   3.858862159804182*^9, 3.858862170630928*^9}, {3.8588622031948433`*^9, 
   3.858862217217596*^9}, 3.8588622927070723`*^9, {3.858862335189515*^9, 
   3.8588623713715487`*^9}, {3.8588624030121393`*^9, 3.858862423920089*^9}, 
   3.858862534571671*^9, 3.85886259983747*^9, 3.858863233722219*^9, {
   3.858866456946732*^9, 3.8588664662466383`*^9}, 3.858866707258143*^9, {
   3.858867014532146*^9, 3.858867017095325*^9}, 3.858867273636587*^9, {
   3.858867362804555*^9, 3.858867366198247*^9}, {3.86015302495588*^9, 
   3.86015304226884*^9}, {3.8616982183250103`*^9, 3.861698233491702*^9}, {
   3.861698315387947*^9, 3.861698350597335*^9}, {3.863145139171679*^9, 
   3.863145154285096*^9}, {3.863500096777033*^9, 3.863500111311249*^9}, {
   3.8659124940351152`*^9, 3.8659124978813133`*^9}},
 CellLabel->
  "In[424]:=",ExpressionUUID->"469ee195-6aa8-4281-b976-bebe95fa1d0a"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"s1", "-", "s2"}], ")"}], " ", 
   SuperscriptBox["X", "2"]}], "-", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", "\[Omega]"}], 
  "-", 
  RowBox[{
   FractionBox["1", "2"], " ", "X", " ", 
   RowBox[{"(", 
    RowBox[{"1", "+", "s1", "-", 
     RowBox[{"3", " ", "s2"}], "-", 
     RowBox[{"2", " ", "\[Alpha]"}], "+", 
     RowBox[{"2", " ", "s1", " ", "\[Alpha]"}], "+", 
     RowBox[{"2", " ", "s2", " ", "\[Alpha]"}], "+", 
     RowBox[{"2", " ", "\[Omega]"}], "-", 
     RowBox[{"2", " ", "s1", " ", "\[Omega]"}], "-", 
     RowBox[{"2", " ", "s2", " ", "\[Omega]"}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.858867282005418*^9, 3.858867305028775*^9}, 
   3.858867393679269*^9, {3.8588674268533707`*^9, 3.8588675222112093`*^9}, 
   3.858867770207775*^9, {3.858867806590564*^9, 3.858867848790469*^9}, {
   3.858867896735847*^9, 3.858867915779181*^9}, 3.858867998307207*^9, {
   3.858868070294032*^9, 3.858868141224883*^9}, 3.858868481203677*^9, {
   3.858868660196684*^9, 3.858868674495759*^9}, 3.858868735879496*^9, 
   3.858868925427186*^9, 3.8588690209631567`*^9, 3.858931968358777*^9, 
   3.860152466950375*^9, 3.8601525352164516`*^9, 3.8601526255500927`*^9, 
   3.860152800240898*^9, 3.86015310791783*^9, 3.8601531722531977`*^9, 
   3.860153243079257*^9, 3.860153441792322*^9, 3.860153568925741*^9, {
   3.8601537759738417`*^9, 3.860153805437243*^9}, 3.860173353808614*^9, 
   3.860387951832019*^9, 3.8604055424112997`*^9, 3.860406029085499*^9, 
   3.860645193824031*^9, 3.8610816678152313`*^9, 3.8610826538553877`*^9, 
   3.8612508652073812`*^9, 3.8612516167124243`*^9, 3.861253705565742*^9, 
   3.861254653465076*^9, 3.861255594836357*^9, 3.861255709450314*^9, 
   3.861692738767109*^9, 3.861695421949988*^9, 3.861698513317726*^9, 
   3.861698590970326*^9, 3.861698652858292*^9, {3.86169868774832*^9, 
   3.861698717050139*^9}, 3.861698978338728*^9, 3.861699712588023*^9, 
   3.861699808839128*^9, 3.861700015553525*^9, 3.8617715457429323`*^9, 
   3.861771679245329*^9, 3.8617803138613987`*^9, 3.86178040785034*^9, 
   3.86178058921723*^9, 3.861781539544114*^9, 3.861782265996117*^9, 
   3.861791543680668*^9, 3.862057671417923*^9, 3.86241647844715*^9, 
   3.8624776536190453`*^9, 3.863065105465067*^9, {3.863086979723831*^9, 
   3.8630870052229013`*^9}, 3.863087366975773*^9, 3.863145274796938*^9, 
   3.8631463107522793`*^9, 3.86341771922013*^9, 3.8634955706587152`*^9, 
   3.863499224836419*^9, {3.863500108124824*^9, 3.863500113108425*^9}, 
   3.863511193115471*^9, 3.864354361886752*^9, 3.8643582544218397`*^9, 
   3.864367415824747*^9, 3.864369228118043*^9, 3.864878916056055*^9, 
   3.864878962581849*^9, 3.8648791444558764`*^9, 3.864975256773211*^9, 
   3.865059056520219*^9, 3.865141740076857*^9, 3.865911862105064*^9, 
   3.865912553527357*^9, 3.866023524141255*^9, 3.866023596553664*^9, {
   3.866024632338049*^9, 3.866024660864986*^9}, {3.86644642805057*^9, 
   3.8664464556048*^9}, 3.8666251973310223`*^9, 3.867404804706688*^9, {
   3.8674048498219433`*^9, 3.867404875455394*^9}, 3.867467477543931*^9, 
   3.867467517443364*^9, 3.8675798877146063`*^9, 3.867752506111412*^9, 
   3.86776472569168*^9, 3.867924151841566*^9, 3.868511811587631*^9, 
   3.868789157860358*^9, 3.869289266252439*^9, 3.869369125441017*^9, 
   3.871210179661607*^9, 3.8712105440359*^9, 3.871699866830064*^9},
 CellLabel->
  "Out[424]=",ExpressionUUID->"d2a3cb84-b01b-4553-9c7d-9101fb66b1d5"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"zsimp", "[", 
     RowBox[{"x1_", ",", " ", "x2_", ",", " ", "y1_", ",", " ", "y2_"}], 
     "]"}], " ", "=", " ", 
    RowBox[{"Simplify", "[", 
     RowBox[{"y2eq", "-", "y1eq"}], "]"}]}], ";"}], "\n"}], "\n", 
 RowBox[{"zeq", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"zsimp", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"X", "/", "2"}], "+", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{"X", "/", "2"}], "-", 
      RowBox[{"\[Omega]", "/", "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "-", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "X"}], ")"}], "/", "2"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Omega]", "-", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "\[Alpha]"}], "-", "1"}], ")"}]}], ")"}], "/", 
       "2"}]}]}], "]"}], "]"}]}], "\n"}], "Input",
 CellChangeTimes->{{3.858862497684667*^9, 3.858862550995639*^9}, {
   3.858862591790657*^9, 3.8588625929918613`*^9}, {3.85886324267227*^9, 
   3.8588632445659637`*^9}, {3.8588664885327263`*^9, 3.8588664967520523`*^9}, 
   3.858866715625643*^9, {3.8588669977183523`*^9, 3.85886700049039*^9}, 
   3.858867279443709*^9, {3.858867344290998*^9, 3.85886739066299*^9}, {
   3.8601530488352737`*^9, 3.8601530510385647`*^9}, {3.8616983842233953`*^9, 
   3.861698411323135*^9}, {3.861698649765191*^9, 3.861698677853909*^9}, {
   3.863145150896463*^9, 3.8631451558937817`*^9}, {3.8659125010516*^9, 
   3.865912512819278*^9}},
 CellLabel->
  "In[425]:=",ExpressionUUID->"2af68567-fc3a-4029-8aa7-eead961c2027"],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
     SuperscriptBox["X", "2"]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "-", "s2"}], ")"}], " ", "\[Omega]"}], 
    "+", 
    RowBox[{"X", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3"}], "+", "s1", "+", 
       RowBox[{"3", " ", "s2"}], "+", 
       RowBox[{"2", " ", "\[Alpha]"}], "+", 
       RowBox[{"2", " ", "s1", " ", "\[Alpha]"}], "-", 
       RowBox[{"2", " ", "s2", " ", "\[Alpha]"}], "-", 
       RowBox[{"2", " ", "s1", " ", "\[Omega]"}], "+", 
       RowBox[{"2", " ", "s2", " ", "\[Omega]"}]}], ")"}]}]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{{3.858867282032443*^9, 3.8588673050540247`*^9}, 
   3.85886735043851*^9, 3.858867393717901*^9, {3.858867426884988*^9, 
   3.858867522243573*^9}, 3.8588677702708673`*^9, {3.8588678066214323`*^9, 
   3.8588678488538446`*^9}, {3.8588678967974577`*^9, 3.8588679158115253`*^9}, 
   3.858867998368438*^9, {3.858868070336669*^9, 3.858868141254753*^9}, 
   3.858868481232396*^9, {3.8588686602505493`*^9, 3.858868674561034*^9}, 
   3.858868735940387*^9, 3.858868925476754*^9, 3.85886902101335*^9, 
   3.8589319684361887`*^9, 3.860152467020193*^9, 3.860152535239037*^9, 
   3.860152625613163*^9, 3.860152800283431*^9, 3.860153107965496*^9, 
   3.860153172315708*^9, 3.860153243144227*^9, 3.860153441855175*^9, 
   3.860153568996172*^9, {3.860153776062327*^9, 3.8601538055109463`*^9}, 
   3.860173353856845*^9, 3.860387951929304*^9, 3.860405542513476*^9, 
   3.8604060291452017`*^9, 3.8606451939344473`*^9, 3.86108166792449*^9, 
   3.8610826538776093`*^9, 3.861250865282031*^9, 3.8612516167858353`*^9, 
   3.861253705646572*^9, 3.861254653543522*^9, 3.8612555949176083`*^9, 
   3.8612557095358877`*^9, 3.861692738870607*^9, 3.861695422019733*^9, 
   3.861698513325014*^9, 3.8616985909819613`*^9, 3.86169865287298*^9, {
   3.861698687766081*^9, 3.8616987170684967`*^9}, 3.861698978356106*^9, 
   3.861699712604026*^9, 3.861699808859638*^9, 3.861700015571217*^9, 
   3.8617715457886333`*^9, 3.861771679262348*^9, 3.86178031388442*^9, 
   3.8617804078687277`*^9, 3.86178058923534*^9, 3.861781539566641*^9, 
   3.861782266026038*^9, 3.861791543697299*^9, 3.862057671434417*^9, 
   3.8624164784646187`*^9, 3.862477653646536*^9, 3.863065105483674*^9, {
   3.86308697974091*^9, 3.863087005239463*^9}, 3.8630873669910793`*^9, 
   3.863145274813877*^9, 3.863146310770735*^9, 3.863417719236714*^9, 
   3.8634955706751747`*^9, 3.8634992248464537`*^9, 3.863500108139492*^9, 
   3.8635111931336718`*^9, 3.864354361901305*^9, 3.864358254439834*^9, 
   3.864367415838118*^9, 3.8643692281459723`*^9, 3.864878916069276*^9, 
   3.86487896259686*^9, 3.864879144467458*^9, 3.864975256790559*^9, 
   3.865059056536964*^9, 3.865141740093924*^9, 3.8659118621151533`*^9, 
   3.86591255353032*^9, 3.8660235241449423`*^9, 3.86602359655763*^9, {
   3.866024632341991*^9, 3.8660246608691483`*^9}, {3.86644642805392*^9, 
   3.866446455609644*^9}, 3.866625197334734*^9, 3.867404804710738*^9, {
   3.867404849825692*^9, 3.867404875458535*^9}, 3.867467477547759*^9, 
   3.8674675174468193`*^9, 3.8675798877179823`*^9, 3.867752506115013*^9, 
   3.8677647256952868`*^9, 3.8679241518545237`*^9, 3.868511811590743*^9, 
   3.868789157863401*^9, 3.869289266256111*^9, 3.8693691254438887`*^9, 
   3.87121017967041*^9, 3.871210544039381*^9, 3.871699866839962*^9},
 CellLabel->
  "Out[426]=",ExpressionUUID->"ca2a849b-0c8d-45d0-9b90-98fe0447712e"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8660239728373117`*^9, 3.8660239744063797`*^9}, {
   3.8660244323847303`*^9, 3.866024457068246*^9}, {3.86602454086482*^9, 
   3.8660245694739523`*^9}, {3.8660246049894733`*^9, 3.866024621489273*^9}, 
   3.86602465436162*^9},
 CellLabel->
  "In[427]:=",ExpressionUUID->"94bd60e5-ec42-43be-9bae-a13d8a1753a9"],

Cell["Eigenvalues and fixed points:", "Text",
 CellChangeTimes->{{3.869289185266778*^9, 
  3.869289189895204*^9}},ExpressionUUID->"075ab1d6-c7e9-4aca-a604-\
e709613023c5"],

Cell[BoxData[
 RowBox[{
  RowBox[{"CharacteristicPolynomial", "[", 
   RowBox[{
    RowBox[{"{", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"Xeq", ",", " ", "X"}], "]"}], ",", 
        RowBox[{"D", "[", 
         RowBox[{"Xeq", ",", " ", "\[Omega]"}], "]"}]}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"\[Omega]eq", ",", " ", "X"}], "]"}], ",", " ", 
        RowBox[{"D", "[", 
         RowBox[{"\[Omega]eq", ",", " ", "\[Omega]"}], "]"}]}], "}"}]}], 
     "\[IndentingNewLine]", "}"}], ",", "\[Lambda]"}], "]"}], ";"}]], "Input",\

 CellChangeTimes->{{3.85886654133745*^9, 3.85886663066705*^9}, {
   3.8588666697272882`*^9, 3.8588666740290823`*^9}, {3.858866736483482*^9, 
   3.858866787328065*^9}, {3.858867026564919*^9, 3.858867094742923*^9}, {
   3.858867215150569*^9, 3.858867225474916*^9}, 3.858867301132814*^9, 
   3.858868088685663*^9, {3.8588687083956003`*^9, 3.858868712122464*^9}, {
   3.8601530551201*^9, 3.860153062700697*^9}, {3.8601531229512444`*^9, 
   3.860153123491501*^9}, {3.8616984259444532`*^9, 3.861698460616787*^9}},
 CellLabel->
  "In[428]:=",ExpressionUUID->"0be0300d-1738-49bd-acd7-a10a35141679"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"gEig", " ", "=", " ", 
   RowBox[{"Eigenvalues", "[", 
    RowBox[{"{", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"Xeq", ",", " ", "X"}], "]"}], ",", 
        RowBox[{"D", "[", 
         RowBox[{"Xeq", ",", " ", "\[Omega]"}], "]"}]}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"\[Omega]eq", ",", " ", "X"}], "]"}], ",", " ", 
        RowBox[{"D", "[", 
         RowBox[{"\[Omega]eq", ",", " ", "\[Omega]"}], "]"}]}], "}"}]}], 
     "\[IndentingNewLine]", "}"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"FullSimplify", "[", "gEig", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Simplify", "[", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Xeq", " ", "==", " ", "0"}], ",", " ", 
      RowBox[{"\[Omega]eq", " ", "==", " ", "0"}]}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"X", ",", " ", "\[Omega]"}], "}"}]}], "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.85886654133745*^9, 3.85886663066705*^9}, {
   3.8588666697272882`*^9, 3.8588666740290823`*^9}, {3.858866736483482*^9, 
   3.858866787328065*^9}, {3.858867026564919*^9, 3.858867094742923*^9}, {
   3.858867215150569*^9, 3.858867225474916*^9}, 3.858867301132814*^9, 
   3.8588680844444313`*^9, {3.858868714685268*^9, 3.858868723034511*^9}, {
   3.858868935150717*^9, 3.858868944892967*^9}, {3.860153065214823*^9, 
   3.860153072076612*^9}, {3.86015312782381*^9, 3.86015312823001*^9}, 
   3.861698468936082*^9, 3.8630869369289637`*^9, {3.863086997081547*^9, 
   3.863087008197392*^9}, 3.863087357964196*^9, {3.8648788856687117`*^9, 
   3.864878913261455*^9}, 3.864878959345837*^9, {3.864878994553845*^9, 
   3.864879009134941*^9}, {3.864879089173069*^9, 3.8648791322195663`*^9}, {
   3.864879166932878*^9, 3.864879171199074*^9}, {3.865912520613352*^9, 
   3.865912534588524*^9}, {3.865912564941828*^9, 3.865912566211576*^9}, {
   3.865912603427758*^9, 3.865912609633834*^9}, {3.86591710805732*^9, 
   3.865917123693993*^9}, 3.866026364989078*^9, 3.866446250504384*^9, 
   3.8666249832205763`*^9, {3.867467474078169*^9, 3.867467474746876*^9}, {
   3.867467509907022*^9, 3.8674675150432673`*^9}, {3.867467572757125*^9, 
   3.86746760515073*^9}, {3.8674680756707563`*^9, 3.8674681062599173`*^9}, {
   3.867468356711096*^9, 3.867468361080101*^9}, {3.867470227343932*^9, 
   3.867470229232098*^9}},
 CellLabel->
  "In[429]:=",ExpressionUUID->"e0b86100-1777-4f21-ac6f-fca2bb511093"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{"s1", "-", 
      RowBox[{"2", " ", "s2"}], "-", 
      RowBox[{"3", " ", "s1", " ", "X"}], "+", 
      RowBox[{"3", " ", "s2", " ", "X"}], "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", "\[Alpha]"}], "+", 
      RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]", "-", 
      RowBox[{"s1", " ", "\[Omega]"}], "-", 
      RowBox[{"s2", " ", "\[Omega]"}], "-", 
      RowBox[{
       FractionBox["1", "2"], " ", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Alpha]", "-", "\[Omega]", "+", 
              RowBox[{"s2", " ", 
               RowBox[{"(", 
                RowBox[{"2", "-", 
                 RowBox[{"3", " ", "X"}], "-", "\[Alpha]", "+", "\[Omega]"}], 
                ")"}]}], "+", 
              RowBox[{"s1", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "1"}], "+", 
                 RowBox[{"3", " ", "X"}], "-", "\[Alpha]", "+", "\[Omega]"}], 
                ")"}]}]}], ")"}], "2"]}], "+", 
          RowBox[{"8", " ", 
           RowBox[{"(", 
            RowBox[{"2", "-", 
             RowBox[{"2", " ", "\[Alpha]"}], "+", 
             RowBox[{"X", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "5"}], "+", 
                RowBox[{"4", " ", "X"}], "+", 
                RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
             RowBox[{"s1", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "2"}], "+", 
                RowBox[{"4", " ", "s2", " ", 
                 SuperscriptBox[
                  RowBox[{"(", 
                   RowBox[{"1", "-", 
                    RowBox[{"2", " ", "X"}]}], ")"}], "2"]}], "+", 
                RowBox[{"X", " ", 
                 RowBox[{"(", 
                  RowBox[{"7", "-", 
                   RowBox[{"8", " ", "X"}], "-", 
                   RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
                RowBox[{"2", " ", "\[Alpha]"}], "-", "\[Omega]"}], ")"}]}], 
             "+", "\[Omega]", "-", 
             RowBox[{"s2", " ", 
              RowBox[{"(", 
               RowBox[{"4", "-", 
                RowBox[{"2", " ", "\[Alpha]"}], "+", 
                RowBox[{"X", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "11"}], "+", 
                   RowBox[{"8", " ", "X"}], "+", 
                   RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
                "\[Omega]"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], ")"}]}], ",", 
   RowBox[{
    FractionBox["1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{"s1", "-", 
      RowBox[{"2", " ", "s2"}], "-", 
      RowBox[{"3", " ", "s1", " ", "X"}], "+", 
      RowBox[{"3", " ", "s2", " ", "X"}], "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", "\[Alpha]"}], "+", 
      RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]", "-", 
      RowBox[{"s1", " ", "\[Omega]"}], "-", 
      RowBox[{"s2", " ", "\[Omega]"}], "+", 
      RowBox[{
       FractionBox["1", "2"], " ", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"4", " ", 
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Alpha]", "-", "\[Omega]", "+", 
              RowBox[{"s2", " ", 
               RowBox[{"(", 
                RowBox[{"2", "-", 
                 RowBox[{"3", " ", "X"}], "-", "\[Alpha]", "+", "\[Omega]"}], 
                ")"}]}], "+", 
              RowBox[{"s1", " ", 
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"-", "1"}], "+", 
                 RowBox[{"3", " ", "X"}], "-", "\[Alpha]", "+", "\[Omega]"}], 
                ")"}]}]}], ")"}], "2"]}], "+", 
          RowBox[{"8", " ", 
           RowBox[{"(", 
            RowBox[{"2", "-", 
             RowBox[{"2", " ", "\[Alpha]"}], "+", 
             RowBox[{"X", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "5"}], "+", 
                RowBox[{"4", " ", "X"}], "+", 
                RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
             RowBox[{"s1", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "2"}], "+", 
                RowBox[{"4", " ", "s2", " ", 
                 SuperscriptBox[
                  RowBox[{"(", 
                   RowBox[{"1", "-", 
                    RowBox[{"2", " ", "X"}]}], ")"}], "2"]}], "+", 
                RowBox[{"X", " ", 
                 RowBox[{"(", 
                  RowBox[{"7", "-", 
                   RowBox[{"8", " ", "X"}], "-", 
                   RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
                RowBox[{"2", " ", "\[Alpha]"}], "-", "\[Omega]"}], ")"}]}], 
             "+", "\[Omega]", "-", 
             RowBox[{"s2", " ", 
              RowBox[{"(", 
               RowBox[{"4", "-", 
                RowBox[{"2", " ", "\[Alpha]"}], "+", 
                RowBox[{"X", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "11"}], "+", 
                   RowBox[{"8", " ", "X"}], "+", 
                   RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], "+", 
                "\[Omega]"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], ")"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.864878916096121*^9, 3.864878962604074*^9, {3.864879001112512*^9, 
   3.864879009469846*^9}, 3.864879099274914*^9, {3.864879133610364*^9, 
   3.864879144541222*^9}, 3.864975256929865*^9, 3.8650590566905537`*^9, 
   3.8651417402509403`*^9, 3.865911862268268*^9, {3.865912553655555*^9, 
   3.865912568530415*^9}, 3.865912611498682*^9, 3.865917129049665*^9, 
   3.8660235243858633`*^9, 3.866023596757762*^9, 3.866024661023293*^9, {
   3.8664464282206593`*^9, 3.8664464557569532`*^9}, 3.866625197487474*^9, 
   3.867404804858536*^9, {3.867404849954516*^9, 3.867404875588387*^9}, 
   3.867467479597342*^9, 3.8674675194760313`*^9, {3.86746759619562*^9, 
   3.867467608775434*^9}, 3.8674680805667467`*^9, 3.86746836164482*^9, 
   3.8675798896623583`*^9, 3.867752508074881*^9, 3.867764727616933*^9, 
   3.8679241539236097`*^9, 3.868511813541338*^9, 3.868789159824354*^9, 
   3.8692892681771173`*^9, 3.869369127371667*^9, 3.8712101815883713`*^9, 
   3.871210545947267*^9, 3.871699868819395*^9},
 CellLabel->
  "Out[430]=",ExpressionUUID->"101be212-c79a-4131-9dc8-8b46822fe084"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"X", "\[Rule]", "0"}], ",", 
     RowBox[{"\[Omega]", "\[Rule]", "0"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"X", "\[Rule]", "1"}], ",", 
     RowBox[{"\[Omega]", "\[Rule]", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       RowBox[{"2", " ", "\[Alpha]"}]}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"X", "\[Rule]", 
      FractionBox[
       RowBox[{"1", "+", 
        RowBox[{"s2", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
        "+", 
        RowBox[{"s1", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s2"}]}], ")"}]}]]}], ",", 
     RowBox[{"\[Omega]", "\[Rule]", 
      FractionBox[
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "s2", "+", 
          RowBox[{"s1", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "3"}], "+", 
             RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"s2", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
          "+", 
          RowBox[{"s1", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", 
             RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}], 
       
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s2"}]}], ")"}]}]]}]}], "}"}]}], "}"}]], "Output",\

 CellChangeTimes->{
  3.864878916096121*^9, 3.864878962604074*^9, {3.864879001112512*^9, 
   3.864879009469846*^9}, 3.864879099274914*^9, {3.864879133610364*^9, 
   3.864879144541222*^9}, 3.864975256929865*^9, 3.8650590566905537`*^9, 
   3.8651417402509403`*^9, 3.865911862268268*^9, {3.865912553655555*^9, 
   3.865912568530415*^9}, 3.865912611498682*^9, 3.865917129049665*^9, 
   3.8660235243858633`*^9, 3.866023596757762*^9, 3.866024661023293*^9, {
   3.8664464282206593`*^9, 3.8664464557569532`*^9}, 3.866625197487474*^9, 
   3.867404804858536*^9, {3.867404849954516*^9, 3.867404875588387*^9}, 
   3.867467479597342*^9, 3.8674675194760313`*^9, {3.86746759619562*^9, 
   3.867467608775434*^9}, 3.8674680805667467`*^9, 3.86746836164482*^9, 
   3.8675798896623583`*^9, 3.867752508074881*^9, 3.867764727616933*^9, 
   3.8679241539236097`*^9, 3.868511813541338*^9, 3.868789159824354*^9, 
   3.8692892681771173`*^9, 3.869369127371667*^9, 3.8712101815883713`*^9, 
   3.871210545947267*^9, 3.8716998688329782`*^9},
 CellLabel->
  "Out[431]=",ExpressionUUID->"d62712ab-5cd6-4437-a092-826494985c7c"]
}, Open  ]],

Cell["\<\
Conditions under which we have 2 or 3 fixed points inside the phase plane\
\>", "Text",
 CellChangeTimes->{{3.871210060154893*^9, 
  3.87121007417246*^9}},ExpressionUUID->"2f6a8f69-80f6-47c3-8f2b-\
acf7d2de24a3"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Reduce", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"s2", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
         "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}], 
       "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}], "<=", "1"}], ",", 
     " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"s2", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
         "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}], 
       "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}], ">=", "0"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "-", "s2", "+", 
           RowBox[{"s1", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "3"}], "+", 
              RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", 
           RowBox[{"s2", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]",
            "+", 
           RowBox[{"s1", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", 
              RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}],
         ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}], "<=", "1"}], ",", 
     " ", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "-", "s2", "+", 
           RowBox[{"s1", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "3"}], "+", 
              RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{"1", "+", 
           RowBox[{"s2", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]",
            "+", 
           RowBox[{"s1", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1"}], "+", 
              RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}],
         ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}], ">=", 
      RowBox[{"-", "1"}]}], ",", 
     RowBox[{"1", ">=", " ", "s1", ">", "0.5"}], ",", " ", 
     RowBox[{"0.5", "<", "s2", "<", "s1"}], ",", " ", 
     RowBox[{"0", "<=", "\[Alpha]", "<=", "1"}]}], " ", "}"}], ",", " ", 
   "\[Alpha]", ",", " ", "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.871210076276485*^9, 3.871210215705594*^9}, {
  3.871210292125964*^9, 3.871210346131948*^9}, {3.871210377505197*^9, 
  3.871210383772278*^9}, {3.871210418811241*^9, 3.871210437218626*^9}},
 CellLabel->
  "In[432]:=",ExpressionUUID->"8affef03-a1d7-4f62-8f18-391362e5765e"],

Cell[BoxData[
 TemplateBox[{
  "Reduce", "ratnz", 
   "\"Reduce was unable to solve the system with inexact coefficients. The \
answer was obtained by solving a corresponding exact system and numericizing \
the result.\"", 2, 432, 27, 34335648624018903052, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{{3.871210352572826*^9, 3.8712103842498617`*^9}, {
   3.8712104214321413`*^9, 3.87121043785049*^9}, 3.87121054602712*^9, 
   3.871699868905802*^9},
 CellLabel->
  "During evaluation of \
In[432]:=",ExpressionUUID->"f3fa63c5-d384-4004-826a-41e39eb3e18c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"0.5`", "<", "s2", "<", "1.`"}], "&&", 
  RowBox[{"s2", "<", "s1", "\[LessEqual]", "1.`"}], "&&", 
  RowBox[{
   FractionBox[
    RowBox[{
     RowBox[{"-", "1.`"}], "+", "s1", "+", 
     RowBox[{"2.`", " ", "s2"}], "-", 
     RowBox[{"2.`", " ", "s1", " ", "s2"}]}], 
    RowBox[{
     RowBox[{"-", "1.`"}], "+", "s1", "+", "s2"}]], "\[LessEqual]", 
   "\[Alpha]", "\[LessEqual]", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "1.`"}], " ", "s1"}], "+", 
     RowBox[{"2.`", " ", "s1", " ", "s2"}]}], 
    RowBox[{
     RowBox[{"-", "1.`"}], "+", "s1", "+", "s2"}]]}]}]], "Output",
 CellChangeTimes->{
  3.871210216609336*^9, {3.8712103420834637`*^9, 3.871210384254216*^9}, {
   3.8712104214599667`*^9, 3.87121043785385*^9}, 3.871210546030506*^9, 
   3.871699868909787*^9},
 CellLabel->
  "Out[432]=",ExpressionUUID->"a6e9ec48-d213-4bd7-869c-c309d1851c19"]
}, Open  ]],

Cell["Curve between extinction and coexistence", "Text",
 CellChangeTimes->{{3.8673057756967907`*^9, 3.8673057821270638`*^9}, {
  3.867305881332798*^9, 3.8673058931480417`*^9}, {3.869289196794005*^9, 
  3.869289208847541*^9}},ExpressionUUID->"38c6c49c-15ac-4399-8d8b-\
af23274b30fa"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"1", "+", 
      RowBox[{"s2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "s2"}]}], ")"}]}]], "==", "0"}], ",", " ", "s2"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.867305792022735*^9, 3.867305810415432*^9}, {
  3.867305857785925*^9, 3.8673058628693647`*^9}, {3.8673061687378674`*^9, 
  3.8673061688350983`*^9}},
 CellLabel->
  "In[433]:=",ExpressionUUID->"3596d2f7-3612-4668-86ec-1157a7b276a6"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"s2", "\[Rule]", 
    FractionBox[
     RowBox[{
      RowBox[{"-", "1"}], "+", "s1", "+", "\[Alpha]", "-", 
      RowBox[{"s1", " ", "\[Alpha]"}]}], 
     RowBox[{
      RowBox[{"-", "2"}], "+", 
      RowBox[{"2", " ", "s1"}], "+", "\[Alpha]"}]]}], "}"}], "}"}]], "Output",\

 CellChangeTimes->{{3.8673058075438347`*^9, 3.867305810936224*^9}, 
   3.867305863530056*^9, 3.8673061699814863`*^9, 3.867404804892837*^9, {
   3.867404849971945*^9, 3.8674048756067677`*^9}, 3.8674674796324997`*^9, 
   3.867467519516377*^9, 3.867579889719322*^9, 3.8677525081310596`*^9, 
   3.867764727652687*^9, 3.867924153985216*^9, 3.8685118136045094`*^9, 
   3.868789159880375*^9, 3.869289268220923*^9, 3.869369127433268*^9, 
   3.8712101817270517`*^9, 3.871210546039765*^9, 3.871699868918412*^9},
 CellLabel->
  "Out[433]=",ExpressionUUID->"f0822bb4-0e31-46b2-ac1a-023105d8d834"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{
      RowBox[{"-", "1"}], "+", "s1", "+", "\[Alpha]", "-", 
      RowBox[{"s1", " ", "\[Alpha]"}]}], 
     RowBox[{
      RowBox[{"-", "2"}], "+", 
      RowBox[{"2", " ", "s1"}], "+", "\[Alpha]"}]], "==", "s1"}], ",", " ", 
   "s1"}], "]"}]], "Input",
 CellChangeTimes->{{3.867314293914281*^9, 3.86731429836952*^9}},
 CellLabel->
  "In[434]:=",ExpressionUUID->"a21023b5-72a6-4efd-afa9-5fe2f1517edf"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"s1", "\[Rule]", 
     FractionBox["1", "2"]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"s1", "\[Rule]", 
     RowBox[{"1", "-", "\[Alpha]"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.867314299019885*^9, 3.867404804896612*^9, {3.867404849974648*^9, 
   3.8674048756095667`*^9}, 3.867467479644*^9, 3.867467519522245*^9, 
   3.8675798897232733`*^9, 3.867752508138156*^9, 3.867764727655726*^9, 
   3.867924153988883*^9, 3.868511813607272*^9, 3.868789159883391*^9, 
   3.869289268244998*^9, 3.8693691274363403`*^9, 3.871210181729822*^9, 
   3.8712105460428247`*^9, 3.871699868921752*^9},
 CellLabel->
  "Out[434]=",ExpressionUUID->"cf7a8405-70b2-40ce-9c32-1172b4d8601e"]
}, Open  ]],

Cell["Curve between dominance and coexistence", "Text",
 CellChangeTimes->{{3.8673057756967907`*^9, 3.8673057821270638`*^9}, {
  3.867305881332798*^9, 3.867305913926736*^9}, {3.869289212081176*^9, 
  3.869289217310479*^9}},ExpressionUUID->"3fa0a123-270f-4500-99c6-\
bee4299e23d2"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"1", "+", 
      RowBox[{"s2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "+", 
        RowBox[{"2", " ", "s2"}]}], ")"}]}]], "==", "1"}], ",", " ", "s2"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.867305792022735*^9, 3.867305810415432*^9}, {
  3.867305857785925*^9, 3.8673058628693647`*^9}, {3.867305929259783*^9, 
  3.867305929327148*^9}, {3.867306171690845*^9, 3.8673061717833023`*^9}},
 CellLabel->
  "In[435]:=",ExpressionUUID->"79740bb2-3be6-4188-b756-767a9a5122b5"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"s2", "\[Rule]", 
    FractionBox[
     RowBox[{"s1", "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", "\[Alpha]"}]}], 
     RowBox[{
      RowBox[{"2", " ", "s1"}], "-", "\[Alpha]"}]]}], "}"}], "}"}]], "Output",\

 CellChangeTimes->{{3.8673058075438347`*^9, 3.867305810936224*^9}, 
   3.867305863530056*^9, 3.8673059306174726`*^9, 3.867306172279272*^9, 
   3.867404804908681*^9, 3.8674048756140757`*^9, 3.8674674796477537`*^9, 
   3.86746751953225*^9, 3.867579889732944*^9, 3.867752508148191*^9, 
   3.867764727658998*^9, 3.867924154002013*^9, 3.868511813610512*^9, 
   3.8687891598966627`*^9, 3.869289268248629*^9, 3.8693691274469976`*^9, 
   3.8712101817376328`*^9, 3.871210546051762*^9, 3.871699868924973*^9},
 CellLabel->
  "Out[435]=",ExpressionUUID->"44709be9-46fc-49ab-b49f-310cd00437cd"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{"s1", "-", "\[Alpha]", "+", 
      RowBox[{"s1", " ", "\[Alpha]"}]}], 
     RowBox[{
      RowBox[{"2", " ", "s1"}], "-", "\[Alpha]"}]], "==", "s1"}], ",", " ", 
   "s1"}], "]"}]], "Input",
 CellChangeTimes->{{3.867314330382997*^9, 3.8673143378169193`*^9}},
 CellLabel->
  "In[436]:=",ExpressionUUID->"4eb24018-4e80-451f-88e1-c61faaacfdb3"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"s1", "\[Rule]", 
     FractionBox["1", "2"]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"s1", "\[Rule]", "\[Alpha]"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.867314351681118*^9, 3.86740480491186*^9, 
  3.867404875625968*^9, 3.867467479650589*^9, 3.8674675195372057`*^9, 
  3.8675798897358427`*^9, 3.8677525081522083`*^9, 3.867764727661541*^9, 
  3.867924154005829*^9, 3.868511813621097*^9, 3.86878915989944*^9, 
  3.8692892682514553`*^9, 3.869369127459489*^9, 3.871210181740325*^9, 
  3.8712105460547953`*^9, 3.871699868927804*^9},
 CellLabel->
  "Out[436]=",ExpressionUUID->"66691ece-03a9-43db-986b-86506fdd9994"]
}, Open  ]],

Cell["Expression for \[Alpha] < 0.5", "Text",
 CellChangeTimes->{{3.867354223417178*^9, 3.867354240891694*^9}, {
  3.867354374813612*^9, 3.867354375327463*^9}, {3.8692892238504477`*^9, 
  3.8692892253839912`*^9}},ExpressionUUID->"4db1f3ed-df6b-4b36-80c3-\
76e2f096c051"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      FractionBox[
       RowBox[{"s1", "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}]}], 
       RowBox[{
        RowBox[{"2", " ", "s1"}], "-", "\[Alpha]"}]], ",", " ", 
      RowBox[{"{", 
       RowBox[{"s1", ",", "0.5", ",", "1"}], "}"}], ",", " ", 
      RowBox[{"Assumptions", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"\[Alpha]", ">", "0"}], ",", " ", 
         RowBox[{"\[Alpha]", "<", "0.5"}]}], "}"}]}]}], "]"}]}], "+", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "+", "\[Alpha]", "-", 
       RowBox[{"s1", " ", "\[Alpha]"}]}], 
      RowBox[{
       RowBox[{"-", "2"}], "+", 
       RowBox[{"2", " ", "s1"}], "+", "\[Alpha]"}]], ",", " ", 
     RowBox[{"{", 
      RowBox[{"s1", ",", "0.5", ",", 
       RowBox[{"1", "-", "\[Alpha]"}]}], "}"}], ",", " ", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"\[Alpha]", ">", "0"}], ",", " ", 
        RowBox[{"\[Alpha]", "<", "0.5"}]}], "}"}]}]}], "]"}], "+", 
   RowBox[{"Integrate", "[", 
    RowBox[{"s1", ",", " ", 
     RowBox[{"{", 
      RowBox[{"s1", ",", 
       RowBox[{"1", "-", "\[Alpha]"}], ",", "1"}], "}"}], ",", " ", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"\[Alpha]", ">", "0"}], ",", " ", 
        RowBox[{"\[Alpha]", "<", "0.5"}]}], "}"}]}]}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.867354243621537*^9, 3.867354357516706*^9}, {
   3.867354512968048*^9, 3.867354570549488*^9}, 3.867354808190084*^9, {
   3.867354855174409*^9, 3.867354937384017*^9}, {3.867404816154417*^9, 
   3.8674048383142357`*^9}, {3.867404869367298*^9, 3.867404871938602*^9}, {
   3.867405986329959*^9, 3.867405987926032*^9}},
 CellLabel->
  "In[437]:=",ExpressionUUID->"4099a828-afe9-46fd-9a8b-bd91da51607b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{"0.`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], ")"}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"0.25`", "\[VeryThinSpace]", "-", 
     RowBox[{"0.25`", " ", "\[Alpha]"}]}], ")"}], " ", "\[Alpha]", " ", 
   RowBox[{"Log", "[", 
    RowBox[{"2.`", "\[VeryThinSpace]", "-", 
     RowBox[{"1.`", " ", "\[Alpha]"}]}], "]"}]}], "+", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "0.25`"}], "+", 
     RowBox[{"0.25`", " ", "\[Alpha]"}]}], ")"}], " ", "\[Alpha]", " ", 
   RowBox[{"Log", "[", "\[Alpha]", "]"}]}]}]], "Output",
 CellChangeTimes->{{3.867354279916294*^9, 3.8673543198796997`*^9}, 
   3.8673543591999903`*^9, 3.867354539606945*^9, 3.867354810301694*^9, 
   3.8673549473517933`*^9, {3.867404805247085*^9, 3.867404818479583*^9}, {
   3.867404857324539*^9, 3.867404876036429*^9}, 3.867405989400771*^9, 
   3.8674674800818453`*^9, 3.86746751995788*^9, 3.867579890188779*^9, 
   3.8677525085973673`*^9, 3.867764728072397*^9, 3.867924154531847*^9, 
   3.868511814071527*^9, 3.868789160365562*^9, 3.869289268689521*^9, 
   3.869369127895471*^9, 3.8712101821719313`*^9, 3.871210546463442*^9, 
   3.871699869356204*^9},
 CellLabel->
  "Out[437]=",ExpressionUUID->"a52f89fe-931a-4615-933b-6bc743c8d601"]
}, Open  ]],

Cell["Expression for \[Alpha] > 0.5. ", "Text",
 CellChangeTimes->{{3.867354223417178*^9, 3.867354240891694*^9}, {
   3.8673543778692513`*^9, 3.867354378494666*^9}, {3.867409795224506*^9, 
   3.867409801024151*^9}, 3.8674105520057287`*^9, {3.8674129342923193`*^9, 
   3.8674129510474052`*^9}, {3.86747962812358*^9, 3.867479631177524*^9}, {
   3.869289228584495*^9, 
   3.8692892301999893`*^9}},ExpressionUUID->"4072e238-bf53-43ce-ae39-\
397f909a4369"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"Integrate", "[", 
     RowBox[{
      FractionBox[
       RowBox[{"s1", "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}]}], 
       RowBox[{
        RowBox[{"2", " ", "s1"}], "-", "\[Alpha]"}]], ",", " ", 
      RowBox[{"{", 
       RowBox[{"s1", ",", "\[Alpha]", ",", "1"}], "}"}], ",", " ", 
      RowBox[{"Assumptions", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"\[Alpha]", "<=", "1"}], ",", " ", 
         RowBox[{"\[Alpha]", ">=", "0.5"}]}], "}"}]}]}], "]"}]}], "+", 
   RowBox[{"Integrate", "[", 
    RowBox[{"s1", ",", " ", 
     RowBox[{"{", 
      RowBox[{"s1", ",", "\[Alpha]", ",", "1"}], "}"}], ",", " ", 
     RowBox[{"Assumptions", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"\[Alpha]", "<=", "1"}], ",", " ", 
        RowBox[{"\[Alpha]", ">=", "0.5"}]}], "}"}]}]}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.867354243621537*^9, 3.867354357516706*^9}, {
   3.8673543984839067`*^9, 3.8673544008753366`*^9}, {3.867404796105361*^9, 
   3.867404800757278*^9}, {3.867409765024239*^9, 3.8674098047771273`*^9}, 
   3.867409840496789*^9, {3.867410171440239*^9, 3.8674101754520082`*^9}, {
   3.8674105186291656`*^9, 3.867410567008356*^9}, 3.867411579632271*^9, {
   3.8674117329904346`*^9, 3.8674117405806217`*^9}, {3.867411803067007*^9, 
   3.867411828741955*^9}, {3.867412628373221*^9, 3.86741263760183*^9}, {
   3.867479635639123*^9, 3.867479649780385*^9}, 3.867764718783648*^9},
 CellLabel->
  "In[438]:=",ExpressionUUID->"9dbffce9-575e-4d54-93a0-34e715dcc8fd"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox["1", "2"]}], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], " ", "\[Alpha]", " ", 
  RowBox[{"ArcTanh", "[", 
   RowBox[{"1", "-", "\[Alpha]"}], "]"}]}]], "Output",
 CellChangeTimes->{3.867412638855022*^9, 3.8674674804457417`*^9, 
  3.867467520318161*^9, 3.867479653000277*^9, 3.867579890503381*^9, 
  3.86775250890994*^9, 3.867764728419133*^9, 3.867924154907922*^9, 
  3.86851181438449*^9, 3.868789160676165*^9, 3.869289269000144*^9, 
  3.869369128211259*^9, 3.871210182484129*^9, 3.8712105468093567`*^9, 
  3.87169986971377*^9},
 CellLabel->
  "Out[438]=",ExpressionUUID->"71605782-9458-4a87-a1c7-89f72116fb0c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindMaximum", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"0.75", "-", 
        RowBox[{
         RowBox[{"1", "/", "2"}], " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"1", "-", "\[Alpha]"}], ")"}], "^", "2"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.75", "+", 
           RowBox[{"0.", " ", "I"}]}], ")"}], " ", "\[Alpha]"}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.5", "+", 
           RowBox[{"0.", " ", "I"}]}], ")"}], " ", 
         RowBox[{"\[Alpha]", "^", "2"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.25", "-", 
           RowBox[{"0.25", " ", "\[Alpha]"}]}], ")"}], " ", "\[Alpha]", " ", 
         RowBox[{"Log", "[", 
          RowBox[{"1.", "-", 
           RowBox[{"1.", " ", "\[Alpha]"}]}], "]"}]}], "-", 
        RowBox[{"0.25", " ", 
         RowBox[{"(", 
          RowBox[{"1.", "+", 
           RowBox[{"1.", " ", "\[Alpha]"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"\[Alpha]", "-", 
              RowBox[{"1.", " ", 
               RowBox[{"\[Alpha]", "^", "2"}]}]}], ")"}], " ", 
            RowBox[{"Log", "[", 
             RowBox[{"1.", "-", "\[Alpha]"}], "]"}]}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"-", "1."}], "+", "\[Alpha]"}], ")"}], " ", "\[Alpha]", 
            " ", 
            RowBox[{"Log", "[", 
             RowBox[{"2.", "-", 
              RowBox[{"1.", " ", "\[Alpha]"}]}], "]"}]}]}], ")"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "0.25"}], "+", 
           RowBox[{"0.25", " ", "\[Alpha]"}]}], ")"}], " ", "\[Alpha]", " ", 
         RowBox[{"Log", "[", "\[Alpha]", "]"}]}]}], ")"}], "/", 
      RowBox[{
       RowBox[{"(", "0.5", ")"}], "^", "3"}]}], ",", " ", 
     RowBox[{"0", "<=", "\[Alpha]", "<=", "0.5"}]}], "}"}], ",", " ", 
   "\[Alpha]"}], "]"}]], "Input",
 CellChangeTimes->{{3.86735460383138*^9, 3.867354677972377*^9}, {
   3.867354725172124*^9, 3.867354744041183*^9}, 3.8673548254495783`*^9, {
   3.867354959567803*^9, 3.867355012429397*^9}},
 CellLabel->
  "In[439]:=",ExpressionUUID->"514dce6d-9e07-4dc4-b66b-5ec1e2a520aa"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.7322307002434056`", ",", 
   RowBox[{"{", 
    RowBox[{"\[Alpha]", "\[Rule]", "0.2721034709307097`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.8673546396763973`*^9, 3.867354679236589*^9}, 
   3.867354746069661*^9, 3.8673548275462513`*^9, 3.867354961253039*^9, {
   3.867354993847177*^9, 3.8673550141751833`*^9}, 3.867404805397773*^9, {
   3.8674048575648193`*^9, 3.8674048762227383`*^9}, 3.867467480480133*^9, 
   3.8674675203402643`*^9, 3.8675798905814857`*^9, 3.867752508984221*^9, 
   3.8677647284438133`*^9, 3.867924155043201*^9, 3.868511814460773*^9, 
   3.868789160747621*^9, 3.869289269071865*^9, 3.86936912827927*^9, 
   3.871210182568894*^9, 3.871210546834855*^9, 3.871699869744131*^9},
 CellLabel->
  "Out[439]=",ExpressionUUID->"fabc8192-ede8-4a10-ba0c-a3e73232397e"]
}, Open  ]],

Cell["Curves for \[Omega]", "Text",
 CellChangeTimes->{{3.867470190576433*^9, 3.867470197279125*^9}, {
  3.869289236679665*^9, 
  3.869289239555356*^9}},ExpressionUUID->"f00b8832-c761-4a5c-b6bf-\
f73cd26b37ac"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "s2", "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "+", 
            RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"s2", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
         "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s2"}]}], ")"}]}]], "==", "0"}], ",", " ", "s2"}], 
   "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.8674702838268833`*^9, 3.8674702975401163`*^9}},
 CellLabel->
  "In[440]:=",ExpressionUUID->"d0a65125-dfed-4b76-8647-69c9cc6b1cbf"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"s2", "\[Rule]", 
     FractionBox[
      RowBox[{"1", "-", 
       RowBox[{"3", " ", "s1"}]}], 
      RowBox[{"1", "-", 
       RowBox[{"4", " ", "s1"}]}]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"s2", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1", "+", "\[Alpha]", "-", 
       RowBox[{"s1", " ", "\[Alpha]"}]}], 
      RowBox[{
       RowBox[{"-", "2"}], "+", 
       RowBox[{"2", " ", "s1"}], "+", "\[Alpha]"}]]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.867470298103097*^9, 3.86757989059334*^9, 
  3.867752508996087*^9, 3.8677647284565363`*^9, 3.86792415509274*^9, 
  3.868511814472603*^9, 3.8687891607638474`*^9, 3.869289269083809*^9, 
  3.869369128291863*^9, 3.871210182580894*^9, 3.871210546846487*^9, 
  3.871699869756214*^9},
 CellLabel->
  "Out[440]=",ExpressionUUID->"52c90303-d70f-48b4-a43a-22bca5262b4d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FullSimplify", "[", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", "s2", "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "3"}], "+", 
            RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{"1", "+", 
         RowBox[{"s2", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
         "+", 
         RowBox[{"s1", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1"}], "+", 
            RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         RowBox[{"2", " ", "s2"}]}], ")"}]}]], "==", 
     RowBox[{
      RowBox[{"2", "\[Alpha]"}], "-", "1"}]}], ",", " ", "s2"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8674703223680277`*^9, 3.867470341613029*^9}},
 CellLabel->
  "In[441]:=",ExpressionUUID->"06745f53-3eb3-46fa-9cad-5e93a8e75656"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"s2", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{"-", "1"}], "+", "s1"}], 
      RowBox[{
       RowBox[{"-", "3"}], "+", 
       RowBox[{"4", " ", "s1"}]}]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"s2", "\[Rule]", 
     FractionBox[
      RowBox[{"s1", "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "s1"}], ")"}], " ", "\[Alpha]"}]}], 
      RowBox[{
       RowBox[{"2", " ", "s1"}], "-", "\[Alpha]"}]]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.86747034205856*^9, 3.8675798906045923`*^9, 
  3.86775250901007*^9, 3.8677647284741993`*^9, 3.867924155104825*^9, 
  3.868511814487105*^9, 3.868789160778945*^9, 3.86928926910166*^9, 
  3.869369128307795*^9, 3.871210182592258*^9, 3.871210546857312*^9, 
  3.871699869771181*^9},
 CellLabel->
  "Out[441]=",ExpressionUUID->"d0c57490-ff45-4d04-9b90-0e8873c3e73d"]
}, Open  ]],

Cell["\<\
Obviously they are the same.
Some phase plane analysis:
\
\>", "Text",
 CellChangeTimes->{{3.867478646813808*^9, 3.8674786515579643`*^9}, 
   3.86747871096912*^9, {3.869289243181975*^9, 
   3.869289258407484*^9}},ExpressionUUID->"1b7dfab8-ba7a-4c13-8d2a-\
37106568fe44"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Row", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"params", "=", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s1", " ", "->", " ", "a"}], ",", 
           RowBox[{"s2", "->", "b"}], ",", " ", 
           RowBox[{"\[Alpha]", " ", "->", " ", "d"}]}], "}"}]}], ";", " ", 
        RowBox[{"sol", " ", "=", " ", 
         RowBox[{"NSolve", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Xeq", "==", "0"}], ",", " ", 
              RowBox[{"\[Omega]eq", " ", "==", " ", "0"}], ",", " ", 
              RowBox[{"X", " ", ">=", "0"}], ",", " ", 
              RowBox[{"X", " ", "<=", " ", "1"}], ",", " ", 
              RowBox[{"\[Omega]", ">=", " ", 
               RowBox[{"-", "1"}]}], ",", " ", 
              RowBox[{"\[Omega]", "<=", " ", "1"}]}], "}"}], "/.", "params"}],
            ",", " ", 
           RowBox[{"{", 
            RowBox[{"X", ",", "\[Omega]"}], "}"}], ",", "\[IndentingNewLine]",
            "Reals"}], "]"}]}]}], ",", 
       RowBox[{
        RowBox[{"params", "=", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s1", "->", "a"}], ",", " ", 
           RowBox[{"s2", "->", "b"}], ",", " ", 
           RowBox[{"\[Alpha]", "->", "d"}]}], "}"}]}], ";", " ", 
        RowBox[{"sol", " ", "=", " ", 
         RowBox[{"NSolve", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Xeq", "==", "0"}], ",", " ", 
              RowBox[{"\[Omega]eq", " ", "==", " ", "0"}], ",", " ", 
              RowBox[{"X", " ", ">=", "0"}], ",", " ", 
              RowBox[{"X", " ", "<=", " ", "1"}], ",", " ", 
              RowBox[{"\[Omega]", ">=", " ", 
               RowBox[{"-", "1"}]}], ",", " ", 
              RowBox[{"\[Omega]", "<=", " ", "1"}]}], "}"}], "/.", "params"}],
            ",", " ", 
           RowBox[{"{", 
            RowBox[{"X", ",", "\[Omega]"}], "}"}], ",", "\[IndentingNewLine]",
            "Reals"}], "]"}]}], ";", " ", 
        RowBox[{
         RowBox[{"gEig", " ", "/.", " ", "params"}], " ", "/.", " ", 
         "sol"}]}], " ", ",", " ", 
       RowBox[{
        RowBox[{"params", " ", "=", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"s1", "->", "a"}], ",", " ", 
           RowBox[{"s2", "->", "b"}], ",", " ", 
           RowBox[{"\[Alpha]", "->", "d"}]}], "}"}]}], ";", 
        RowBox[{"Show", "[", 
         RowBox[{"StreamPlot", "[", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{
            RowBox[{"{", 
             RowBox[{"Xeq", ",", " ", "\[Omega]eq"}], "}"}], "/.", "params"}],
            ",", " ", 
           RowBox[{"{", 
            RowBox[{"X", ",", "0", ",", "1"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"\[Omega]", ",", " ", 
             RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{"RegionFunction", "\[Rule]", 
            RowBox[{"Function", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{"X", ",", " ", "\[Omega]"}], "}"}], ",", 
              RowBox[{
               RowBox[{"X", "\[LessEqual]", "1"}], " ", "&&", " ", 
               RowBox[{"\[Omega]", " ", "<=", " ", "1"}], " ", "&&", 
               RowBox[{"X", "\[GreaterEqual]", "0"}], " ", "&&", " ", 
               RowBox[{"\[Omega]", "\[GreaterEqual]", " ", 
                RowBox[{"-", "1"}]}]}]}], "]"}]}], ",", "\[IndentingNewLine]", 
           RowBox[{"RegionFillingStyle", "\[Rule]", "None"}], ",", " ", 
           RowBox[{"RegionBoundaryStyle", "\[Rule]", 
            RowBox[{"Directive", "[", 
             RowBox[{"Black", ",", " ", 
              RowBox[{"Thickness", "[", "0.005", "]"}]}], "]"}]}], ",", " ", 
           RowBox[{"FrameLabel", "\[Rule]", 
            RowBox[{"{", 
             RowBox[{"\"\<X\>\"", ",", "\"\<\[Omega]\>\""}], "}"}]}], ",", 
           " ", 
           RowBox[{"ImageSize", "\[Rule]", "500"}], ",", " ", 
           RowBox[{"StreamColorFunction", "\[Rule]", "Blue"}]}], "]"}], 
         "]"}]}]}], " ", "}"}], ",", " ", "\"\<\\n ----- \\n\>\""}], "]"}], 
   ",", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"a", ",", "0.4", ",", " ", "\"\<s1 (0, 1)\>\""}], "}"}], ",", 
     " ", "0", ",", " ", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"b", ",", "0.4", ",", " ", "\"\<s2\>\""}], "}"}], ",", " ", "0",
      ",", " ", "1"}], "}"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"d", ",", "1", ",", " ", "\"\<\[Alpha]\>\""}], "}"}], ",", " ", 
     "0", ",", " ", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.858866788945181*^9, 3.858866879190977*^9}, {
   3.858867098741412*^9, 3.858867103392097*^9}, 3.8588691158813*^9, {
   3.858869221213497*^9, 3.8588692356209793`*^9}, {3.858869265645424*^9, 
   3.858869267424342*^9}, 3.858931989302166*^9, {3.860153083772064*^9, 
   3.860153084114306*^9}, {3.860153272437931*^9, 3.860153273710642*^9}, {
   3.860153678519973*^9, 3.860153683075061*^9}, {3.861081541488422*^9, 
   3.861081606948578*^9}, {3.861253641833795*^9, 3.861253677176689*^9}, {
   3.861254619924552*^9, 3.861254636424799*^9}, {3.8612555582055073`*^9, 
   3.861255571435094*^9}, {3.861255680424262*^9, 3.8612556819976263`*^9}, {
   3.861698475805439*^9, 3.86169848215394*^9}, {3.861698547100659*^9, 
   3.861698570217772*^9}, {3.861771527047102*^9, 3.86177155328891*^9}, {
   3.8617802894553747`*^9, 3.8617802984085703`*^9}, {3.861780389179859*^9, 
   3.861780417176011*^9}, {3.861782040760109*^9, 3.861782066259344*^9}, 
   3.861783094910935*^9, {3.862477660818738*^9, 3.862477706935204*^9}, {
   3.8624777489230137`*^9, 3.862477895491427*^9}, {3.863145163523374*^9, 
   3.863145184241681*^9}, {3.863146468335793*^9, 3.863146471589459*^9}, {
   3.871699856961631*^9, 3.871699863421898*^9}},
 CellLabel->
  "In[442]:=",ExpressionUUID->"cb5a4693-e76e-4990-a1ed-f3b4ee78c7d7"],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`a$$ = 1, $CellContext`b$$ = 
    0.85, $CellContext`d$$ = 0.27, Typeset`show$$ = True, 
    Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{{
       Hold[$CellContext`a$$], 0.4, "s1 (0, 1)"}, 0, 1}, {{
       Hold[$CellContext`b$$], 0.4, "s2"}, 0, 1}, {{
       Hold[$CellContext`d$$], 1, "\[Alpha]"}, 0, 1}}, Typeset`size$$ = {
    494., {385., 12.791915893554688`}}, Typeset`update$$ = 0, 
    Typeset`initDone$$, Typeset`skipInitDone$$ = True}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`a$$ = 0.4, $CellContext`b$$ = 
        0.4, $CellContext`d$$ = 1}, "ControllerVariables" :> {}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> 
      Row[{$CellContext`params = {$CellContext`s1 -> $CellContext`a$$, \
$CellContext`s2 -> $CellContext`b$$, $CellContext`\[Alpha] -> \
$CellContext`d$$}; $CellContext`sol = NSolve[
            
            ReplaceAll[{$CellContext`Xeq == 0, $CellContext`\[Omega]eq == 
              0, $CellContext`X >= 0, $CellContext`X <= 
              1, $CellContext`\[Omega] >= -1, $CellContext`\[Omega] <= 
              1}, $CellContext`params], {$CellContext`X, \
$CellContext`\[Omega]}, 
            Reals], $CellContext`params = {$CellContext`s1 -> \
$CellContext`a$$, $CellContext`s2 -> $CellContext`b$$, $CellContext`\[Alpha] -> \
$CellContext`d$$}; $CellContext`sol = NSolve[
            
            ReplaceAll[{$CellContext`Xeq == 0, $CellContext`\[Omega]eq == 
              0, $CellContext`X >= 0, $CellContext`X <= 
              1, $CellContext`\[Omega] >= -1, $CellContext`\[Omega] <= 
              1}, $CellContext`params], {$CellContext`X, \
$CellContext`\[Omega]}, Reals]; ReplaceAll[
           
           ReplaceAll[$CellContext`gEig, $CellContext`params], \
$CellContext`sol], $CellContext`params = {$CellContext`s1 -> \
$CellContext`a$$, $CellContext`s2 -> $CellContext`b$$, $CellContext`\[Alpha] -> \
$CellContext`d$$}; Show[
           StreamPlot[
            
            ReplaceAll[{$CellContext`Xeq, $CellContext`\[Omega]eq}, \
$CellContext`params], {$CellContext`X, 0, 1}, {$CellContext`\[Omega], -1, 1}, 
            RegionFunction -> 
            Function[{$CellContext`X, $CellContext`\[Omega]}, 
              
              And[$CellContext`X <= 1, $CellContext`\[Omega] <= 
               1, $CellContext`X >= 0, $CellContext`\[Omega] >= -1]], 
            RegionFillingStyle -> None, RegionBoundaryStyle -> 
            Directive[Black, 
              Thickness[0.005]], FrameLabel -> {"X", "\[Omega]"}, ImageSize -> 
            500, StreamColorFunction -> Blue]]}, "\n ----- \n"], 
      "Specifications" :> {{{$CellContext`a$$, 0.4, "s1 (0, 1)"}, 0, 
         1}, {{$CellContext`b$$, 0.4, "s2"}, 0, 
         1}, {{$CellContext`d$$, 1, "\[Alpha]"}, 0, 1}}, "Options" :> {}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{539., {310., 316.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$},
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.858931968573983*^9, 3.8589319912834377`*^9}, 
   3.860152467089087*^9, 3.8601525352777643`*^9, 3.86015262567493*^9, 
   3.860152800331601*^9, 3.860153108036916*^9, 3.860153172353202*^9, 
   3.860153243203228*^9, 3.8601532869515953`*^9, 3.860153441897812*^9, 
   3.860153569047038*^9, {3.860153776134721*^9, 3.860153805567068*^9}, 
   3.860173353937909*^9, 3.86038795198038*^9, 3.860405542564948*^9, 
   3.860406029215342*^9, {3.8604069766301947`*^9, 3.860406978354844*^9}, 
   3.860645194002349*^9, {3.860645309384914*^9, 3.8606453113029346`*^9}, {
   3.860645388610216*^9, 3.860645426344595*^9}, {3.86108166801042*^9, 
   3.8610816815738077`*^9}, {3.86108265393524*^9, 3.861082662633116*^9}, 
   3.861250865324832*^9, 3.8612516168187733`*^9, 3.8612537056829042`*^9, {
   3.861254163684724*^9, 3.861254164561536*^9}, 3.861254653575507*^9, 
   3.861255594954206*^9, 3.8612557095605717`*^9, 3.8616927389523582`*^9, 
   3.861695422077867*^9, 3.861698513350925*^9, 3.8616985910132723`*^9, 
   3.86169865288443*^9, {3.861698687805884*^9, 3.8616987277667017`*^9}, 
   3.861698978377988*^9, 3.861699712637854*^9, 3.861699808889062*^9, {
   3.861699899704608*^9, 3.861699901621787*^9}, 3.861700015604261*^9, 
   3.861771545842497*^9, 3.861771679299528*^9, 3.8617802285761766`*^9, {
   3.861780299895774*^9, 3.861780313937201*^9}, {3.8617804043066883`*^9, 
   3.861780418731707*^9}, 3.8617805892687483`*^9, {3.861781235368318*^9, 
   3.8617812362184687`*^9}, 3.861781539608645*^9, {3.861781988039761*^9, 
   3.8617819889079123`*^9}, {3.861782038067231*^9, 3.861782090397736*^9}, 
   3.8617822660624847`*^9, {3.861783033381675*^9, 3.861783034932373*^9}, 
   3.861783096520219*^9, {3.86178313726598*^9, 3.8617831384974327`*^9}, 
   3.861783392560945*^9, 3.861791543726005*^9, 3.862057671457205*^9, 
   3.862416478488386*^9, 3.862477653671631*^9, 3.862477709083613*^9, 
   3.862477795237444*^9, {3.862477825973917*^9, 3.862477830083207*^9}, {
   3.862477874484*^9, 3.862477897122025*^9}, 3.8630651055167837`*^9, {
   3.863065165663289*^9, 3.8630651672120657`*^9}, {3.863086979792445*^9, 
   3.863087005359519*^9}, 3.863087367020277*^9, 3.863145274855899*^9, 
   3.863145352444034*^9, 3.863146310807551*^9, {3.863146449415991*^9, 
   3.863146472949383*^9}, 3.863146544990814*^9, {3.863146989003093*^9, 
   3.8631469901838427`*^9}, {3.8631478354642467`*^9, 3.8631478534509277`*^9}, 
   3.863417719281996*^9, 3.8634955707237797`*^9, 3.86349922488586*^9, 
   3.8635001081804333`*^9, 3.86351119316546*^9, {3.864354361920147*^9, 
   3.864354368457181*^9}, 3.864358254477199*^9, {3.864367415883051*^9, 
   3.8643674222478228`*^9}, {3.864369228181562*^9, 3.864369232740456*^9}, 
   3.864878916112282*^9, 3.8648789626236*^9, 3.8648791445936213`*^9, 
   3.864975256981759*^9, {3.865059056768909*^9, 3.8650590818618727`*^9}, 
   3.865141740308852*^9, 3.865911862328021*^9, 3.865912553727116*^9, 
   3.866023524429022*^9, 3.8660235968028917`*^9, 3.866024661084649*^9, {
   3.866446428240185*^9, 3.866446455778544*^9}, 3.867404805414383*^9, {
   3.86740485757934*^9, 3.867404876230185*^9}, 3.8674674804968147`*^9, 
   3.867467520354206*^9, {3.86747082055307*^9, 3.867470823704432*^9}, {
   3.867579890632525*^9, 3.86757991161213*^9}, {3.867752509017624*^9, 
   3.867752514490576*^9}, 3.867764728483375*^9, 3.867924155123609*^9, {
   3.8679242381416187`*^9, 3.8679242393095922`*^9}, {3.86792703928331*^9, 
   3.867927041150558*^9}, 3.86851181449677*^9, 3.868789160798369*^9, {
   3.8692892623708*^9, 3.869289269108593*^9}, 3.8693691283415813`*^9, {
   3.869369173222622*^9, 3.869369175837977*^9}, 3.871210182607565*^9, 
   3.87121054687188*^9, {3.8716998697793427`*^9, 3.871699876254136*^9}},
 CellLabel->
  "Out[442]=",ExpressionUUID->"fca7857c-d951-46ee-adef-78fb839e906f"]
}, Open  ]],

Cell["\<\
Proof of the fact that there is always one and only one stable fixed point:\
\>", "Text",
 CellChangeTimes->{{3.865915705940113*^9, 3.865915710546555*^9}, {
  3.869289282369295*^9, 3.869289298817587*^9}},
 EmphasizeSyntaxErrors->
  True,ExpressionUUID->"449afb7b-93d2-434b-8eb0-50188e0944c9"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{
    RowBox[{"Xsol", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Omega]sol", " ", "=", " ", "0"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"eig1", "  ", "=", " ", 
     RowBox[{
      RowBox[{"1", "/", "2"}], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "-", 
        RowBox[{
         RowBox[{"1", "/", "2"}], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "^", "2"}]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                RowBox[{"Xsol", "^", "2"}]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "^", "2"}]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{"eig2", " ", "=", " ", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "2"]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["Xsol", "2"]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "2"]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Xsol", " ", "=", " ", "1"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Omega]sol", " ", "=", " ", 
     RowBox[{
      RowBox[{"2", "\[Alpha]"}], "-", "1"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"eig11", "  ", "=", " ", 
     RowBox[{
      RowBox[{"1", "/", "2"}], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "-", 
        RowBox[{
         RowBox[{"1", "/", "2"}], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "^", "2"}]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                RowBox[{"Xsol", "^", "2"}]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "^", "2"}]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{"eig21", " ", "=", " ", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "2"]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["Xsol", "2"]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "2"]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Xsol", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"s2", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
        "+", 
        RowBox[{"s1", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Omega]sol", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", "s2", "+", 
          RowBox[{"s1", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "3"}], "+", 
             RowBox[{"4", " ", "s2"}]}], ")"}]}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"1", "+", 
          RowBox[{"s2", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}]}], "-", "\[Alpha]", 
          "+", 
          RowBox[{"s1", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"-", "1"}], "+", 
             RowBox[{"2", " ", "s2"}], "+", "\[Alpha]"}], ")"}]}]}], ")"}]}], 
       ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s1"}]}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", "s1", "+", "s2"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "s2"}]}], ")"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"eig12", "  ", "=", " ", 
     RowBox[{
      RowBox[{"1", "/", "2"}], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "-", 
        RowBox[{
         RowBox[{"1", "/", "2"}], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "^", "2"}]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                RowBox[{"Xsol", "^", "2"}]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   RowBox[{
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "^", "2"}]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   RowBox[{"Xsol", "^", "2"}]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\n", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"eig22", " ", "=", " ", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      RowBox[{"(", 
       RowBox[{"s1", "-", 
        RowBox[{"2", " ", "s2"}], "-", 
        RowBox[{"3", " ", "s1", " ", "Xsol"}], "+", 
        RowBox[{"3", " ", "s2", " ", "Xsol"}], "-", "\[Alpha]", "+", 
        RowBox[{"s1", " ", "\[Alpha]"}], "+", 
        RowBox[{"s2", " ", "\[Alpha]"}], "+", "\[Omega]sol", "-", 
        RowBox[{"s1", " ", "\[Omega]sol"}], "-", 
        RowBox[{"s2", " ", "\[Omega]sol"}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         RowBox[{"\[Sqrt]", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"4", " ", 
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Alpha]", "-", "\[Omega]sol", "+", 
                RowBox[{"s2", " ", 
                 RowBox[{"(", 
                  RowBox[{"2", "-", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}], "+", 
                RowBox[{"s1", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"3", " ", "Xsol"}], "-", "\[Alpha]", "+", 
                   "\[Omega]sol"}], ")"}]}]}], ")"}], "2"]}], "+", 
            RowBox[{"8", " ", 
             RowBox[{"(", 
              RowBox[{"2", "-", 
               RowBox[{"5", " ", "Xsol"}], "+", 
               RowBox[{"4", " ", 
                SuperscriptBox["Xsol", "2"]}], "-", 
               RowBox[{"2", " ", "\[Alpha]"}], "+", 
               RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
               RowBox[{"s1", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "2"}], "+", 
                  RowBox[{"4", " ", "s2", " ", 
                   SuperscriptBox[
                    RowBox[{"(", 
                    RowBox[{"1", "-", 
                    RowBox[{"2", " ", "Xsol"}]}], ")"}], "2"]}], "+", 
                  RowBox[{"7", " ", "Xsol"}], "-", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "+", 
                  RowBox[{"2", " ", "\[Alpha]"}], "-", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "-", 
                  "\[Omega]sol"}], ")"}]}], "+", "\[Omega]sol", "-", 
               RowBox[{"s2", " ", 
                RowBox[{"(", 
                 RowBox[{"4", "-", 
                  RowBox[{"11", " ", "Xsol"}], "+", 
                  RowBox[{"8", " ", 
                   SuperscriptBox["Xsol", "2"]}], "-", 
                  RowBox[{"2", " ", "\[Alpha]"}], "+", 
                  RowBox[{"2", " ", "Xsol", " ", "\[Alpha]"}], "+", 
                  "\[Omega]sol"}], ")"}]}]}], ")"}]}]}], ")"}]}]}]}], 
       ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"FindInstance", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"eig12", " ", "<", " ", "0"}], " ", "&&", " ", 
      RowBox[{"eig22", "<", " ", "0"}], " ", "&&", " ", 
      RowBox[{"eig1", " ", "<", " ", "0"}], " ", "&&", " ", 
      RowBox[{"eig2", " ", "<", " ", "0"}], " ", "&&", " ", 
      RowBox[{"s1", " ", ">", " ", "0.5"}], " ", "&&", " ", 
      RowBox[{"s1", " ", "<=", " ", "1"}], " ", "&&", " ", 
      RowBox[{"s2", " ", "<=", "1"}], " ", "&&", " ", 
      RowBox[{"s2", " ", ">=", " ", "0.5"}], " ", "&&", " ", 
      RowBox[{"\[Alpha]", " ", ">=", " ", "0"}], " ", "&&", " ", 
      RowBox[{"\[Alpha]", " ", "<=", " ", "1"}]}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
     "Reals"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.8659144386221457`*^9, 3.86591444548584*^9}, {
  3.865914629867428*^9, 3.865914632947584*^9}, {3.8659151733460608`*^9, 
  3.865915251227811*^9}, {3.865915309973053*^9, 3.86591535731148*^9}, {
  3.865915420823588*^9, 3.865915457214643*^9}, {3.865915490542314*^9, 
  3.86591555079107*^9}, {3.865915592844904*^9, 3.865915689104127*^9}, {
  3.868789152779052*^9, 3.868789155529029*^9}},
 CellLabel->
  "In[443]:=",ExpressionUUID->"3ca43e1e-ddbd-447d-b579-354d277bf9e9"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.867404811543532*^9, {3.867404863670363*^9, 3.867404882308198*^9}, 
   3.867467486600728*^9, 3.867467526473866*^9, 3.867579896747492*^9, 
   3.8677525152484417`*^9, 3.8677647345382442`*^9, 3.867924161490707*^9, 
   3.868511820668153*^9, 3.868789165986928*^9, 3.869289274320198*^9, 
   3.869369133682674*^9, 3.8712101878608913`*^9, 3.871210551933676*^9, 
   3.871699874921714*^9},
 CellLabel->
  "Out[455]=",ExpressionUUID->"923bfe99-7bf0-4b8b-b390-34989e325b94"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"eig11", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig21", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig1", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig2", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"s1", " ", ">=", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{"s1", " ", "<=", " ", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", "<=", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", ">=", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{"\[Alpha]", " ", ">=", " ", "0"}], " ", "&&", " ", 
    RowBox[{"\[Alpha]", " ", "<=", " ", "1"}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.868789134689683*^9, 3.868789142230661*^9}},
 CellLabel->
  "In[456]:=",ExpressionUUID->"5aeff5f3-8666-4127-9e4d-0df7aa403bdd"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.867404814298839*^9, {3.8674048664373693`*^9, 3.867404885077054*^9}, 
   3.867467489363205*^9, 3.867467529228161*^9, 3.8675798995105667`*^9, 
   3.867752518015082*^9, 3.8677647373231487`*^9, 3.867924164367985*^9, 
   3.868511823435416*^9, 3.86878916728959*^9, 3.869289275611272*^9, 
   3.869369135002418*^9, 3.871210189171838*^9, 3.871210553216374*^9, 
   3.8716998762374067`*^9},
 CellLabel->
  "Out[456]=",ExpressionUUID->"518ff038-45b9-44f9-82d9-5951deedca63"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"eig11", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig21", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig12", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"eig22", " ", "<", " ", "0"}], " ", "&&", " ", 
    RowBox[{"s1", " ", ">", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{"s1", " ", "<=", " ", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", "<=", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", ">=", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{"\[Alpha]", " ", ">=", " ", "0"}], " ", "&&", " ", 
    RowBox[{"\[Alpha]", " ", "<=", " ", "1"}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.8664513550511093`*^9, 3.866451356056759*^9}, {
  3.868789147453623*^9, 3.868789149893615*^9}},
 CellLabel->
  "In[457]:=",ExpressionUUID->"86ab6656-9e08-452b-9796-b5c44ad52dda"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{{3.865915297766872*^9, 3.865915365892737*^9}, {
   3.8659155213300323`*^9, 3.8659155512830772`*^9}, {3.865915598281645*^9, 
   3.8659156963369427`*^9}, 3.866023533084003*^9, 3.866023605274208*^9, 
   3.866024667923985*^9, {3.866446434318404*^9, 3.86644646183081*^9}, 
   3.867404815729196*^9, {3.8674048678700943`*^9, 3.867404886521366*^9}, 
   3.8674674908117437`*^9, 3.867467530656106*^9, 3.867579900958168*^9, 
   3.867752519485716*^9, 3.867764738780074*^9, 3.867924165923925*^9, 
   3.868511824874218*^9, 3.8687891741621428`*^9, 3.869289282440097*^9, 
   3.869369141986496*^9, 3.871210196089405*^9, 3.87121056006831*^9, 
   3.871699883117909*^9},
 CellLabel->
  "Out[457]=",ExpressionUUID->"291b8789-7267-420d-893f-30a92f004ce0"]
}, Open  ]],

Cell["Proof that eigenvalues are always real:", "Text",
 CellChangeTimes->{{3.866446299046276*^9, 3.866446316666811*^9}, {
  3.869289310203562*^9, 3.869289317215094*^9}},
 EmphasizeSyntaxErrors->
  True,ExpressionUUID->"1ee1f362-280a-4221-be27-6c7181c0e5f6"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"Im", "[", "eig1", "]"}], " ", "!=", " ", "0"}], " ", "||", " ", 
    
    RowBox[{
     RowBox[{"Im", "[", "eig2", "]"}], " ", "!=", " ", "0"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.866446331717065*^9, 3.86644642437904*^9}},
 CellLabel->
  "In[458]:=",ExpressionUUID->"bf9198b7-995a-47e1-965e-cc8711d97e25"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{{3.866446434375643*^9, 3.866446461858843*^9}, 
   3.867404815760995*^9, {3.86740486789633*^9, 3.867404886553508*^9}, 
   3.8674674908434668`*^9, 3.867467530682776*^9, 3.867579900988529*^9, 
   3.867752519517687*^9, 3.8677647388136883`*^9, 3.8679241659565783`*^9, 
   3.868511824902636*^9, 3.8687891741914663`*^9, 3.869289282464579*^9, 
   3.869369142014612*^9, 3.871210196113387*^9, 3.871210560091757*^9, 
   3.8716998831445503`*^9},
 CellLabel->
  "Out[458]=",ExpressionUUID->"8211249e-f266-4cf1-8040-a1444155e4f7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"Im", "[", "eig11", "]"}], " ", "!=", " ", "0"}], " ", "||", " ", 
    RowBox[{
     RowBox[{"Im", "[", "eig21", "]"}], " ", "!=", " ", "0"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.866446331717065*^9, 3.86644644910677*^9}},
 CellLabel->
  "In[459]:=",ExpressionUUID->"cdd863ab-eea8-41cf-97b0-cac36fcfba6d"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{{3.866446434375643*^9, 3.866446461871352*^9}, 
   3.867404815773497*^9, {3.867404867908101*^9, 3.867404886562192*^9}, 
   3.867467490861518*^9, 3.867467530698145*^9, 3.8675799010064297`*^9, 
   3.8677525195345707`*^9, 3.867764738827264*^9, 3.867924165978918*^9, 
   3.8685118249216223`*^9, 3.86878917419934*^9, 3.8692892824725323`*^9, 
   3.869369142028595*^9, 3.87121019612295*^9, 3.871210560100862*^9, 
   3.8716998831592627`*^9},
 CellLabel->
  "Out[459]=",ExpressionUUID->"6c166ba2-0b87-46d3-9cbb-4cc95eb2f4b8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"Im", "[", "eig12", "]"}], " ", "!=", " ", "0"}], " ", "||", " ", 
    RowBox[{
     RowBox[{"Im", "[", "eig22", "]"}], " ", "!=", " ", "0"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", " ", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.866446331717065*^9, 3.8664464524981422`*^9}},
 CellLabel->
  "In[460]:=",ExpressionUUID->"c4f51ab6-71e2-463b-a5e5-80786403c817"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{{3.866446434375643*^9, 3.866446461957384*^9}, 
   3.867404815860435*^9, {3.8674048679954777`*^9, 3.86740488665279*^9}, 
   3.867467490950437*^9, 3.8674675307871113`*^9, 3.867579901094849*^9, 
   3.867752519633802*^9, 3.867764738915173*^9, 3.867924166072104*^9, 
   3.868511825015864*^9, 3.868789174284878*^9, 3.869289282558588*^9, 
   3.869369142116076*^9, 3.8712101962160463`*^9, 3.871210560189649*^9, 
   3.871699883248452*^9},
 CellLabel->
  "Out[460]=",ExpressionUUID->"27b605a0-4ae1-46a1-b526-fc3d8abc71d9"]
}, Open  ]],

Cell["Coexistence is always stable?", "Text",
 CellChangeTimes->{{3.871210550003965*^9, 3.871210558238459*^9}, {
  3.871210591399549*^9, 
  3.871210596323874*^9}},ExpressionUUID->"b3d10c9f-4cea-4525-aef4-\
c7cb57d0ecc9"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FindInstance", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"eig12", ">=", "0"}], "&&", 
    RowBox[{"eig22", ">=", "0"}], "&&", 
    RowBox[{"s1", " ", ">", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{"s1", " ", "<=", " ", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", "<=", "1"}], " ", "&&", " ", 
    RowBox[{"s2", " ", ">=", " ", "0.5"}], " ", "&&", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1."}], "+", "s1", "+", 
        RowBox[{"2.", " ", "s2"}], "-", 
        RowBox[{"2.", " ", "s1", " ", "s2"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1."}], "+", "s1", "+", "s2"}], ")"}]}], "<=", 
     "\[Alpha]", "<=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "1."}], " ", "s1"}], "+", 
        RowBox[{"2.", " ", "s1", " ", "s2"}]}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1."}], "+", "s1", "+", "s2"}], ")"}]}]}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"s1", ",", " ", "s2", ",", " ", "\[Alpha]"}], "}"}], ",", 
   "Reals"}], "]"}]], "Input",
 CellChangeTimes->{{3.871210600021102*^9, 3.871210676661789*^9}},
 CellLabel->
  "In[461]:=",ExpressionUUID->"6bfb16a3-980a-4483-9c72-67c4810eee35"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{3.8712106796063967`*^9, 3.8716998858861017`*^9},
 CellLabel->
  "Out[461]=",ExpressionUUID->"d4c02563-e772-42fc-b691-81fb12a83f06"]
}, Open  ]]
},
WindowSize->{668, 796},
WindowMargins->{{Automatic, 243}, {Automatic, 0}},
TaggingRules->Association["TryRealOnly" -> False],
FrontEndVersion->"12.3 for Mac OS X ARM (64-bit) (September 30, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"abf73b34-79e6-4863-9799-885faf629fd8"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 306, 5, 135, "Text",ExpressionUUID->"bca58f23-7d8f-4b00-a95a-4a9c3983f9f4"],
Cell[867, 27, 197, 6, 58, "Text",ExpressionUUID->"eb1e4875-dc77-47c3-ae5c-03672dcc3d66"],
Cell[1067, 35, 159, 3, 35, "Text",ExpressionUUID->"80d6cbf0-a7d4-4343-8444-c7375179ad8b"],
Cell[1229, 40, 2249, 54, 325, "Input",ExpressionUUID->"af4b99ea-e871-4906-81b7-b78ed789b70b"],
Cell[3481, 96, 175, 3, 35, "Text",ExpressionUUID->"85e8d96f-828b-43a9-b6a2-7bddc2a8570c"],
Cell[CellGroupData[{
Cell[3681, 103, 2731, 62, 199, "Input",ExpressionUUID->"551d2428-46e7-4cf5-ba2c-362ae31e2494"],
Cell[6415, 167, 1401, 26, 34, "Output",ExpressionUUID->"5a60df6f-acfe-4c9c-ae6b-b8c0191d2e8c"],
Cell[7819, 195, 1399, 26, 34, "Output",ExpressionUUID->"08069168-1955-4769-b46b-798dfcfe30cd"],
Cell[9221, 223, 1420, 27, 34, "Output",ExpressionUUID->"1d4366de-e75d-4326-bae4-c774f93fc4ac"],
Cell[10644, 252, 1451, 28, 34, "Output",ExpressionUUID->"e2d9cdad-d17b-4e39-bb09-9a0588625511"],
Cell[12098, 282, 1158, 17, 34, "Output",ExpressionUUID->"4ca77c4f-682a-417d-ab7b-7ac4d530869a"],
Cell[13259, 301, 1682, 36, 34, "Output",ExpressionUUID->"dbfdf09d-6bdd-4001-9805-fe1e3c65b2d8"],
Cell[14944, 339, 1703, 37, 34, "Output",ExpressionUUID->"20613bae-7fa8-437e-862b-4979a5f7489d"],
Cell[16650, 378, 1158, 17, 34, "Output",ExpressionUUID->"a2d353e1-2759-4b04-a11e-f2337db561f9"],
Cell[17811, 397, 1156, 17, 34, "Output",ExpressionUUID->"1d114ae8-ab81-4ce3-977e-0edfef173971"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19004, 419, 180, 3, 30, "Input",ExpressionUUID->"2983a36d-74b9-41f8-85c8-4ccf5f2daead"],
Cell[19187, 424, 2732, 45, 34, "Output",ExpressionUUID->"e48602e6-f464-45d2-9b32-bd98475db257"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21956, 474, 1206, 21, 30, "Input",ExpressionUUID->"cf407a73-21af-4f5b-a822-837fcf26ebef"],
Cell[23165, 497, 3443, 61, 34, "Output",ExpressionUUID->"6d3dde4e-80a4-44c6-b8eb-aa2011b4ae75"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26645, 563, 1848, 45, 52, "Input",ExpressionUUID->"2a458930-9ac6-47d4-a410-4dde926116e4"],
Cell[28496, 610, 3635, 65, 74, "Output",ExpressionUUID->"4b5076b9-042e-4de8-ae0b-53d7115a5600"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32168, 680, 1904, 50, 115, "Input",ExpressionUUID->"1523053d-5e99-4475-a0ae-b61faed2e050"],
Cell[34075, 732, 3649, 66, 74, "Output",ExpressionUUID->"5e2019ba-8543-4539-b7de-1aa51350ddf5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37761, 803, 1254, 21, 30, "Input",ExpressionUUID->"5e8ec10b-9b57-4a69-bef3-09c7a4fadedf"],
Cell[39018, 826, 3464, 62, 34, "Output",ExpressionUUID->"f180403d-b79d-411f-ada9-b395aa5e00c2"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42519, 893, 1902, 46, 52, "Input",ExpressionUUID->"469ee195-6aa8-4281-b976-bebe95fa1d0a"],
Cell[44424, 941, 3547, 62, 88, "Output",ExpressionUUID->"d2a3cb84-b01b-4553-9c7d-9101fb66b1d5"]
}, Open  ]],
Cell[CellGroupData[{
Cell[48008, 1008, 1931, 51, 115, "Input",ExpressionUUID->"2af68567-fc3a-4029-8aa7-eead961c2027"],
Cell[49942, 1061, 3653, 66, 74, "Output",ExpressionUUID->"ca2a849b-0c8d-45d0-9b90-98fe0447712e"]
}, Open  ]],
Cell[53610, 1130, 355, 6, 30, "Input",ExpressionUUID->"94bd60e5-ec42-43be-9bae-a13d8a1753a9"],
Cell[53968, 1138, 171, 3, 35, "Text",ExpressionUUID->"075ab1d6-c7e9-4aca-a604-e709613023c5"],
Cell[54142, 1143, 1280, 29, 94, "Input",ExpressionUUID->"0be0300d-1738-49bd-acd7-a10a35141679"],
Cell[CellGroupData[{
Cell[55447, 1176, 2574, 51, 136, "Input",ExpressionUUID->"e0b86100-1777-4f21-ac6f-fca2bb511093"],
Cell[58024, 1229, 6540, 156, 281, "Output",ExpressionUUID->"101be212-c79a-4131-9dc8-8b46822fe084"],
Cell[64567, 1387, 3337, 90, 122, "Output",ExpressionUUID->"d62712ab-5cd6-4437-a092-826494985c7c"]
}, Open  ]],
Cell[67919, 1480, 222, 5, 35, "Text",ExpressionUUID->"2f6a8f69-80f6-47c3-8f2b-acf7d2de24a3"],
Cell[CellGroupData[{
Cell[68166, 1489, 4586, 140, 157, "Input",ExpressionUUID->"8affef03-a1d7-4f62-8f18-391362e5765e"],
Cell[72755, 1631, 579, 12, 46, "Message",ExpressionUUID->"f3fa63c5-d384-4004-826a-41e39eb3e18c"],
Cell[73337, 1645, 908, 25, 51, "Output",ExpressionUUID->"a6e9ec48-d213-4bd7-869c-c309d1851c19"]
}, Open  ]],
Cell[74260, 1673, 283, 4, 35, "Text",ExpressionUUID->"38c6c49c-15ac-4399-8d8b-af23274b30fa"],
Cell[CellGroupData[{
Cell[74568, 1681, 963, 29, 49, "Input",ExpressionUUID->"3596d2f7-3612-4668-86ec-1157a7b276a6"],
Cell[75534, 1712, 918, 20, 51, "Output",ExpressionUUID->"f0822bb4-0e31-46b2-ac1a-023105d8d834"]
}, Open  ]],
Cell[CellGroupData[{
Cell[76489, 1737, 492, 14, 48, "Input",ExpressionUUID->"a21023b5-72a6-4efd-afa9-5fe2f1517edf"],
Cell[76984, 1753, 749, 17, 50, "Output",ExpressionUUID->"cf7a8405-70b2-40ce-9c32-1172b4d8601e"]
}, Open  ]],
Cell[77748, 1773, 280, 4, 35, "Text",ExpressionUUID->"3fa0a123-270f-4500-99c6-bee4299e23d2"],
Cell[CellGroupData[{
Cell[78053, 1781, 1007, 29, 49, "Input",ExpressionUUID->"79740bb2-3be6-4188-b756-767a9a5122b5"],
Cell[79063, 1812, 851, 18, 51, "Output",ExpressionUUID->"44709be9-46fc-49ab-b49f-310cd00437cd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[79951, 1835, 431, 12, 48, "Input",ExpressionUUID->"4eb24018-4e80-451f-88e1-c61faaacfdb3"],
Cell[80385, 1849, 694, 15, 50, "Output",ExpressionUUID->"66691ece-03a9-43db-986b-86506fdd9994"]
}, Open  ]],
Cell[81094, 1867, 270, 4, 35, "Text",ExpressionUUID->"4db1f3ed-df6b-4b36-80c3-76e2f096c051"],
Cell[CellGroupData[{
Cell[81389, 1875, 1996, 52, 173, "Input",ExpressionUUID->"4099a828-afe9-46fd-9a8b-bd91da51607b"],
Cell[83388, 1929, 1300, 28, 34, "Output",ExpressionUUID->"a52f89fe-931a-4615-933b-6bc743c8d601"]
}, Open  ]],
Cell[84703, 1960, 451, 7, 35, "Text",ExpressionUUID->"4072e238-bf53-43ce-ae39-397f909a4369"],
Cell[CellGroupData[{
Cell[85179, 1971, 1625, 37, 105, "Input",ExpressionUUID->"9dbffce9-575e-4d54-93a0-34e715dcc8fd"],
Cell[86807, 2010, 705, 16, 50, "Output",ExpressionUUID->"71605782-9458-4a87-a1c7-89f72116fb0c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[87549, 2031, 2311, 62, 115, "Input",ExpressionUUID->"514dce6d-9e07-4dc4-b66b-5ec1e2a520aa"],
Cell[89863, 2095, 832, 15, 34, "Output",ExpressionUUID->"fabc8192-ede8-4a10-ba0c-a3e73232397e"]
}, Open  ]],
Cell[90710, 2113, 210, 4, 35, "Text",ExpressionUUID->"f00b8832-c761-4a5c-b6bf-f73cd26b37ac"],
Cell[CellGroupData[{
Cell[90945, 2121, 1338, 41, 107, "Input",ExpressionUUID->"d0a65125-dfed-4b76-8647-69c9cc6b1cbf"],
Cell[92286, 2164, 939, 26, 51, "Output",ExpressionUUID->"52c90303-d70f-48b4-a43a-22bca5262b4d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[93262, 2195, 1390, 43, 107, "Input",ExpressionUUID->"06745f53-3eb3-46fa-9cad-5e93a8e75656"],
Cell[94655, 2240, 958, 28, 51, "Output",ExpressionUUID->"d0c57490-ff45-4d04-9b90-0e8873c3e73d"]
}, Open  ]],
Cell[95628, 2271, 280, 8, 81, "Text",ExpressionUUID->"1b7dfab8-ba7a-4c13-8d2a-37106568fe44"],
Cell[CellGroupData[{
Cell[95933, 2283, 6088, 135, 388, "Input",ExpressionUUID->"cb5a4693-e76e-4990-a1ed-f3b4ee78c7d7"],
Cell[102024, 2420, 7462, 125, 668, "Output",ExpressionUUID->"fca7857c-d951-46ee-adef-78fb839e906f"]
}, Open  ]],
Cell[109501, 2548, 302, 6, 35, "Text",ExpressionUUID->"449afb7b-93d2-434b-8eb0-50188e0944c9"],
Cell[CellGroupData[{
Cell[109828, 2558, 20752, 496, 1624, "Input",ExpressionUUID->"3ca43e1e-ddbd-447d-b579-354d277bf9e9"],
Cell[130583, 3056, 531, 10, 34, "Output",ExpressionUUID->"923bfe99-7bf0-4b8b-b390-34989e325b94"]
}, Open  ]],
Cell[CellGroupData[{
Cell[131151, 3071, 928, 19, 73, "Input",ExpressionUUID->"5aeff5f3-8666-4127-9e4d-0df7aa403bdd"],
Cell[132082, 3092, 532, 10, 34, "Output",ExpressionUUID->"518ff038-45b9-44f9-82d9-5951deedca63"]
}, Open  ]],
Cell[CellGroupData[{
Cell[132651, 3107, 980, 20, 73, "Input",ExpressionUUID->"86ab6656-9e08-452b-9796-b5c44ad52dda"],
Cell[133634, 3129, 798, 13, 34, "Output",ExpressionUUID->"291b8789-7267-420d-893f-30a92f004ce0"]
}, Open  ]],
Cell[134447, 3145, 258, 4, 35, "Text",ExpressionUUID->"1ee1f362-280a-4221-be27-6c7181c0e5f6"],
Cell[CellGroupData[{
Cell[134730, 3153, 519, 14, 30, "Input",ExpressionUUID->"bf9198b7-995a-47e1-965e-cc8711d97e25"],
Cell[135252, 3169, 581, 10, 34, "Output",ExpressionUUID->"8211249e-f266-4cf1-8040-a1444155e4f7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[135870, 3184, 516, 13, 52, "Input",ExpressionUUID->"cdd863ab-eea8-41cf-97b0-cac36fcfba6d"],
Cell[136389, 3199, 580, 10, 34, "Output",ExpressionUUID->"6c166ba2-0b87-46d3-9cbb-4cc95eb2f4b8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[137006, 3214, 519, 13, 52, "Input",ExpressionUUID->"c4f51ab6-71e2-463b-a5e5-80786403c817"],
Cell[137528, 3229, 577, 10, 34, "Output",ExpressionUUID->"27b605a0-4ae1-46a1-b526-fc3d8abc71d9"]
}, Open  ]],
Cell[138120, 3242, 220, 4, 35, "Text",ExpressionUUID->"b3d10c9f-4cea-4525-aef4-c7cb57d0ecc9"],
Cell[CellGroupData[{
Cell[138365, 3250, 1277, 35, 73, "Input",ExpressionUUID->"6bfb16a3-980a-4483-9c72-67c4810eee35"],
Cell[139645, 3287, 196, 4, 34, "Output",ExpressionUUID->"d4c02563-e772-42fc-b691-81fb12a83f06"]
}, Open  ]]
}
]
*)