(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 261590, 5881] NotebookOptionsPosition[ 260616, 5844] NotebookOutlinePosition[ 261050, 5863] CellTagsIndexPosition[ 261007, 5860] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Female reinforcement with learning - program \ \>", "Title", CellChangeTimes->{{3.514027121633904*^9, 3.5140271567415752`*^9}, { 3.521290889531502*^9, 3.521290897122304*^9}, {3.521473152463529*^9, 3.521473153868655*^9}, {3.523700380890628*^9, 3.5237003819779997`*^9}, { 3.573392098669653*^9, 3.573392103029467*^9}}], Cell["Maria Servedio 2012", "Subsubtitle", CellChangeTimes->{{3.514027138292103*^9, 3.514027153691615*^9}, { 3.5733920891750803`*^9, 3.573392093813245*^9}, {3.573397735171698*^9, 3.5733977364578753`*^9}}], Cell[CellGroupData[{ Cell["program", "Section", CellChangeTimes->{{3.480178905886136*^9, 3.480178906494956*^9}}], Cell[BoxData[ RowBox[{"FemReinfLearnLoopProg", " ", ":=", " ", RowBox[{"(", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "arrays", " ", "of", " ", "plots", " ", "in", " ", "case", " ", "I", " ", "want", " ", "to", " ", "look", " ", "at", " ", "them"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a12plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"a22plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"m12plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"m22plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"n12plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"n22plot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D1amplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D1anplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D1nmplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D1amnplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D2amplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D2anplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D2nmplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"D2amnplot", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtaba12", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "14"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabn12", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabn22", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "14"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabn2m2p2", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabn1m1p2", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabhybp2", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabD1nm", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabD2nm", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "14"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabkappa", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"grandtabgamma", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "7"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"kappaparams", " ", "=", RowBox[{"Table", "[", RowBox[{"{", RowBox[{"1.0", ",", " ", "1.1"}], "}"}], "]"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"gammaparams", " ", "=", " ", RowBox[{"Table", "[", RowBox[{"{", RowBox[{ "1.2", ",", " ", "1.3", ",", " ", "1.4", ",", " ", "1.5", ",", " ", "1.6", ",", "1.7", ",", " ", "1.8"}], "}"}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", " ", RowBox[{"taba2p1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"(*", RowBox[{ "setting", " ", "up", " ", "tables", " ", "for", " ", "figures"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"tabn2p1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabm2p1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDamp1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDanp1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDnmp1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDamnp1", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"taba2p2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabn2p2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabm2p2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDamp2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDanp2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDnmp2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"tabDamnp2", "=", " ", RowBox[{"Table", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"i", ",", "100000"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"\[Gamma]1", " ", "=", " ", RowBox[{"gammaparams", "[", RowBox[{"[", "k1", "]"}], "]"}]}], ";", RowBox[{"(*", " ", RowBox[{"learning", ",", " ", RowBox[{"pop", " ", "1"}], ",", " ", RowBox[{"\[Gamma]", " ", ">", " ", "1"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"\[Gamma]2", " ", "=", " ", "\[Gamma]1"}], ";", " ", "\[IndentingNewLine]", RowBox[{"\[Kappa]", " ", "=", RowBox[{"kappaparams", "[", RowBox[{"[", "k2", "]"}], "]"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{ "multiplier", " ", "of", " ", "female", " ", "preference", " ", "strength", " ", "for", " ", "A2"}], ",", " ", RowBox[{"\[Kappa]", " ", ">", " ", "1"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"count", " ", "=", " ", "0"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabkappa", "[", RowBox[{"[", RowBox[{"k1", ",", "k2"}], "]"}], "]"}], " ", "=", " ", "\[Kappa]"}], " ", ";", RowBox[{"(*", " ", RowBox[{ "storing", " ", "values", " ", "for", " ", "later", " ", "table"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabgamma", "[", RowBox[{"[", RowBox[{"k1", ",", "k2"}], "]"}], "]"}], " ", "=", " ", "\[Gamma]1"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"starting", " ", "values"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"a12", " ", "=", " ", "0.01"}], ";", "\[IndentingNewLine]", RowBox[{"n12", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"m12", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"a22", " ", "=", " ", "0.01"}], ";", "\[IndentingNewLine]", RowBox[{"n22", " ", "=", "1.0"}], ";", "\[IndentingNewLine]", RowBox[{"m22", " ", "=", " ", "1.0"}], ";", "\[IndentingNewLine]", RowBox[{"D1an", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D1am", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D1nm", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D1anm", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D2an", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D2am", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D2nm", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", RowBox[{"D2anm", " ", "=", " ", "0.0"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "genotypes", " ", "in", " ", "terms", " ", "of", " ", "allele", " ", "frequencies"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}]}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], "D1nm"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], " ", "D1am"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}], "D1an"}], " ", "-", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], "m12"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], "D1nm"}], "-", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], "D1am"}], " ", "+", " ", RowBox[{"m12", " ", "D1an"}], " ", "+", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], " ", "n12", " ", RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], "D1nm"}], " ", "+", " ", RowBox[{"n12", " ", "D1am"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}], "D1an"}], " ", "+", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], " ", "n12", " ", "m12"}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a12"}], ")"}], "D1nm"}], "-", RowBox[{"n12", " ", "D1am"}], " ", "-", " ", RowBox[{"m12", " ", "D1an"}], " ", "-", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "5", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a12", RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}]}], " ", "+", RowBox[{"a12", " ", "D1nm"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], "D1am"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}], "D1an"}], " ", "+", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "6", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a12", RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], "m12"}], " ", "-", " ", RowBox[{"a12", " ", "D1nm"}], "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n12"}], ")"}], "D1am"}], " ", "-", " ", RowBox[{"m12", " ", "D1an"}], " ", "-", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "7", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a12", " ", "n12", " ", RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}]}], "-", " ", RowBox[{"a12", " ", "D1nm"}], " ", "-", " ", RowBox[{"n12", " ", "D1am"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m12"}], ")"}], " ", "D1an"}], " ", "-", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "8", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a12", " ", "n12", " ", "m12"}], " ", "+", " ", RowBox[{"a12", " ", "D1nm"}], " ", "+", " ", RowBox[{"n12", " ", "D1am"}], " ", "+", " ", RowBox[{"m12", " ", "D1an"}], " ", "+", " ", "D1anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}]}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], "D2nm"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], " ", "D2am"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}], "D2an"}], " ", "-", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], "m22"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], "D2nm"}], "-", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], "D2am"}], " ", "+", " ", RowBox[{"m22", " ", "D2an"}], " ", "+", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], " ", "n22", " ", RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], "D2nm"}], " ", "+", " ", RowBox[{"n22", " ", "D2am"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}], "D2an"}], " ", "+", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], " ", "n22", " ", "m22"}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "a22"}], ")"}], "D2nm"}], "-", RowBox[{"n22", " ", "D2am"}], " ", "-", " ", RowBox[{"m22", " ", "D2an"}], " ", "-", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "5", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a22", RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}]}], " ", "+", RowBox[{"a22", " ", "D2nm"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], "D2am"}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}], "D2an"}], " ", "+", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "6", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a22", RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], "m22"}], " ", "-", " ", RowBox[{"a22", " ", "D2nm"}], "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "n22"}], ")"}], "D2am"}], " ", "-", " ", RowBox[{"m22", " ", "D2an"}], " ", "-", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "7", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a22", " ", "n22", " ", RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}]}], "-", " ", RowBox[{"a22", " ", "D2nm"}], " ", "-", " ", RowBox[{"n22", " ", "D2am"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "m22"}], ")"}], " ", "D2an"}], " ", "-", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "8", "]"}], " ", "=", " ", RowBox[{ RowBox[{"a22", " ", "n22", " ", "m22"}], " ", "+", " ", RowBox[{"a22", " ", "D2nm"}], " ", "+", " ", RowBox[{"n22", " ", "D2am"}], " ", "+", " ", RowBox[{"m22", " ", "D2an"}], " ", "+", " ", "D2anm"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"r1", " ", "=", " ", RowBox[{"(", "0.5", ")"}]}], ";", RowBox[{"(*", RowBox[{ RowBox[{"setting", " ", "recombination", " ", "rates"}], " ", "=", " ", RowBox[{"1", "/", "2"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"r2", " ", "=", " ", RowBox[{"(", "0.5", ")"}]}], ";", "\[IndentingNewLine]", RowBox[{"r3", " ", "=", " ", RowBox[{"(", "0.5", ")"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"starting", " ", "next", " ", "gen", " ", "freqs"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"x1t1", "[", "i", "]"}], " ", "=", " ", "1.0"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"x2t1", "[", "i", "]"}], " ", "=", " ", "1.0"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"gen", " ", "=", " ", "1"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"While", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"gen", " ", "<", " ", "10"}], ")"}], "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "1", "]"}], " ", "-", " ", RowBox[{"x1", "[", "1", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", " ", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "2", "]"}], " ", "-", " ", RowBox[{"x1", "[", "2", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "3", "]"}], " ", "-", " ", RowBox[{"x1", "[", "3", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "4", "]"}], " ", "-", " ", RowBox[{"x1", "[", "4", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "5", "]"}], " ", "-", " ", RowBox[{"x1", "[", "5", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "6", "]"}], " ", "-", " ", RowBox[{"x1", "[", "6", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "7", "]"}], " ", "-", " ", RowBox[{"x1", "[", "7", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x1t1", "[", "8", "]"}], " ", "-", " ", RowBox[{"x1", "[", "8", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], " ", "-", " ", RowBox[{"x2", "[", "1", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], " ", "-", " ", RowBox[{"x2", "[", "2", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "3", "]"}], " ", "-", " ", RowBox[{"x2", "[", "3", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "4", "]"}], " ", "-", " ", RowBox[{"x2", "[", "4", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "5", "]"}], " ", "-", " ", RowBox[{"x2", "[", "5", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "6", "]"}], " ", "-", " ", RowBox[{"x2", "[", "6", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "7", "]"}], " ", "-", " ", RowBox[{"x2", "[", "7", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}], " ", "||", RowBox[{"(", RowBox[{ RowBox[{"Abs", "[", RowBox[{ RowBox[{"x2t1", "[", "8", "]"}], " ", "-", " ", RowBox[{"x2", "[", "8", "]"}]}], "]"}], " ", ">", " ", "0.0000000001"}], ")"}]}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", "recursions", " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"gen", " ", ">", " ", "1"}], ",", RowBox[{ RowBox[{ RowBox[{"x1", "[", "1", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "1", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "2", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "2", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "3", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "3", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "4", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "4", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "5", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "5", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "6", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "6", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "7", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "7", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", "[", "8", "]"}], " ", "=", " ", RowBox[{"x1t1", "[", "8", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "1", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "1", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "2", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "2", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "3", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "3", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "4", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "4", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "5", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "5", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "6", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "6", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "7", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "7", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2", "[", "8", "]"}], " ", "=", " ", RowBox[{"x2t1", "[", "8", "]"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"gen", "++"}], ";", "\[IndentingNewLine]", RowBox[{"count", "++"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"taba2p1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "a12"}], ";", " ", RowBox[{"(*", RowBox[{"for", " ", "plotting"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tabn2p1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "n12"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabm2p1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", "m12"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDamp1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D1am"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDanp1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D1an"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDnmp1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D1nm"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDamnp1", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D1anm"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"taba2p2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "a22"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabn2p2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "n22"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabm2p2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", "m22"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDamp2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D2am"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDanp2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D2an"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDnmp2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D2nm"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tabDamnp2", "[", RowBox[{"[", "count", "]"}], "]"}], " ", "=", " ", "D2anm"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"perturbation", " ", "-", " ", RowBox[{ "currently", " ", "set", " ", "to", " ", "work", " ", "only", " ", "when", " ", "A2", " ", "is", " ", "high"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"count", "\[Equal]", " ", "perturb"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tempx1", "[", "5", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "5", "]"}], "*", "0.95"}]}], ";", "\[IndentingNewLine]", RowBox[{"perttot1", " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "5", "]"}], "-", RowBox[{"tempx1", "[", "5", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"t1", "[", "k", "]"}], " ", "=", " ", RowBox[{"RandomReal", "[", "]"}]}], ",", " ", RowBox[{"{", RowBox[{"k", ",", "7"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"sumrand1", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"t1", "[", "k", "]"}], ",", " ", RowBox[{"{", RowBox[{"k", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "1", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "1", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "2", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "2", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "3", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "3", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "4", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "4", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "6", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "6", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "5", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "7", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "7", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "6", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx1", "[", "8", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x1", "[", "8", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t1", "[", "7", "]"}], "*", RowBox[{"perttot1", "/", "sumrand1"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"x1", "[", "k", "]"}], " ", "=", " ", RowBox[{"tempx1", "[", "k", "]"}]}], ",", " ", RowBox[{"{", RowBox[{"k", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"tempx2", "[", "8", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x2", "[", "8", "]"}], "*", "0.95"}]}], ";", "\[IndentingNewLine]", RowBox[{"perttot2", " ", "=", " ", RowBox[{ RowBox[{"x2", "[", "8", "]"}], "-", RowBox[{"tempx2", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"t2", "[", "k", "]"}], " ", "=", " ", RowBox[{"RandomReal", "[", "]"}]}], ",", " ", RowBox[{"{", RowBox[{"k", ",", "7"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"sumrand2", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"t2", "[", "k", "]"}], ",", " ", RowBox[{"{", RowBox[{"k", ",", " ", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"tempx2", "[", "k", "]"}], " ", "=", " ", RowBox[{ RowBox[{"x2", "[", "k", "]"}], " ", "+", " ", RowBox[{ RowBox[{"t2", "[", "k", "]"}], "*", RowBox[{"perttot2", "/", "sumrand2"}]}]}]}], ",", " ", RowBox[{"{", RowBox[{"k", ",", "7"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"x2", "[", "k", "]"}], " ", "=", " ", RowBox[{"tempx2", "[", "k", "]"}]}], ",", " ", RowBox[{"{", RowBox[{"k", ",", " ", "8"}], "}"}]}], "]"}], ";"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", "migration", " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"xm1", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{"x1", "[", "i", "]"}], RowBox[{"(", RowBox[{"1.0", "-", "m"}], ")"}]}], "+", RowBox[{ RowBox[{"x2", "[", "i", "]"}], " ", "m"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}]}], "]"}]}], ";", "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"xm2", "[", "i", "]"}], "=", RowBox[{ RowBox[{ RowBox[{"x2", "[", "i", "]"}], RowBox[{"(", RowBox[{"1.0", "-", "m"}], ")"}]}], "+", RowBox[{ RowBox[{"x1", "[", "i", "]"}], " ", "m"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "mean", " ", "fitnesses", " ", "for", " ", "natural", " ", "selection"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"wbar1", "=", RowBox[{ RowBox[{"xm1", "[", "1", "]"}], "+", RowBox[{"xm1", "[", "4", "]"}], "+", RowBox[{"xm1", "[", "5", "]"}], "+", RowBox[{"xm1", "[", "8", "]"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"xm1", "[", "2", "]"}], "+", RowBox[{"xm1", "[", "3", "]"}], "+", RowBox[{"xm1", "[", "6", "]"}], "+", RowBox[{"xm1", "[", "7", "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"wbar2", "=", RowBox[{ RowBox[{"xm2", "[", "1", "]"}], "+", RowBox[{"xm2", "[", "4", "]"}], "+", RowBox[{"xm2", "[", "5", "]"}], "+", RowBox[{"xm2", "[", "8", "]"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"xm2", "[", "2", "]"}], "+", RowBox[{"xm2", "[", "3", "]"}], "+", RowBox[{"xm2", "[", "6", "]"}], "+", RowBox[{"xm2", "[", "7", "]"}]}], ")"}], "*", RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}]}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "natural", " ", "selection", " ", "in", " ", "each", " ", "population"}], " ", "*)"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"xns1", "[", "1", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "1", "]"}], "/", "wbar1"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "2", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "2", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar1"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "3", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "3", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar1"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "4", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "4", "]"}], "/", "wbar1"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "5", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "5", "]"}], "/", "wbar1"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "6", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "6", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar1"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "7", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "7", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar1"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns1", "[", "8", "]"}], "=", RowBox[{ RowBox[{"xm1", "[", "8", "]"}], "/", "wbar1"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"xns2", "[", "1", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "1", "]"}], "/", "wbar2"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "2", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "2", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar2"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "3", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "3", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar2"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "4", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "4", "]"}], "/", "wbar2"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "5", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "5", "]"}], "/", "wbar2"}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "6", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "6", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar2"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "7", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "7", "]"}], " ", RowBox[{ RowBox[{"(", RowBox[{"1.0", "-", "s"}], ")"}], "/", "wbar2"}]}]}], ";", "\n", " ", RowBox[{ RowBox[{"xns2", "[", "8", "]"}], "=", RowBox[{ RowBox[{"xm2", "[", "8", "]"}], "/", "wbar2"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "frequencies", " ", "of", " ", "males", " ", "at", " ", "the", " ", "M", " ", "and", " ", "N", " ", "loci"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y1", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns1", "[", "1", "]"}], " ", "+", " ", RowBox[{"xns1", "[", "5", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y1", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns1", "[", "2", "]"}], " ", "+", " ", RowBox[{"xns1", "[", "6", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y1", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns1", "[", "3", "]"}], " ", "+", " ", RowBox[{"xns1", "[", "7", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y1", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns1", "[", "4", "]"}], " ", "+", " ", RowBox[{"xns1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"y2", "[", "1", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns2", "[", "1", "]"}], " ", "+", " ", RowBox[{"xns2", "[", "5", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y2", "[", "2", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns2", "[", "2", "]"}], " ", "+", " ", RowBox[{"xns2", "[", "6", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y2", "[", "3", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns2", "[", "3", "]"}], " ", "+", " ", RowBox[{"xns2", "[", "7", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"y2", "[", "4", "]"}], " ", "=", " ", RowBox[{ RowBox[{"xns2", "[", "4", "]"}], " ", "+", " ", RowBox[{"xns2", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"frequencies", " ", "of", " ", "encounter"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"mate1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"xns1", "[", "i", "]"}], " ", RowBox[{"xns1", "[", "j", "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"mate2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"xns2", "[", "i", "]"}], " ", RowBox[{"xns2", "[", "j", "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ "G", " ", "is", " ", "a", " ", "matrix", " ", "of", " ", "coefficients"}], ",", " ", RowBox[{ "F", " ", "is", " ", "the", " ", "mating", " ", "table"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"G1", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\n", " ", RowBox[{"F1", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"G2", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\n", " ", RowBox[{"F2", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Population", " ", "1"}], " ", "*)"}], "\n", RowBox[{"(*", " ", RowBox[{"A1N1M1", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\[IndentingNewLine]", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"1", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y1", "[", "1", "]"}], "\[Gamma]1", " ", "\[Alpha]"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}], " ", "+", " ", RowBox[{"y1", "[", "4", "]"}]}], ")"}], "\[Alpha]"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N2M2", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"4", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "1", "]"}], "+", RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}]}], ")"}], "\[Alpha]"}], "+", " ", RowBox[{ RowBox[{"y1", "[", "4", "]"}], "\[Gamma]1", " ", "\[Alpha]"}]}]}]}], " ", "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N1M1", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "1"}], " ", "&&", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}]}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "1", "]"}], "+", RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}]}], ")"}], "+", " ", RowBox[{ RowBox[{"y1", "[", "4", "]"}], "/", "\[Gamma]2"}]}]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N2M2", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "0"}], " ", "&&", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}]}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y1", "[", "1", "]"}], "/", "\[Gamma]2"}], "+", RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}], "+", " ", RowBox[{"y1", "[", "4", "]"}]}], ")"}]}]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"purebred", " ", RowBox[{"w", "/", "hybrid"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "1"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", "0"}]}], ")"}], " ", "&&", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "2"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], " ", "\[Equal]", " ", "3"}]}], ")"}]}], ",", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "1", "]"}], "+", " ", RowBox[{"y1", "[", "4", "]"}]}], ")"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}]}], " ", ")"}], "/", "\[Gamma]h"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"hybrid", " ", "females"}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "2"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", "3"}]}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A2N1M1", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"5", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y1", "[", "1", "]"}], "\[Gamma]1", " ", "\[Kappa]", " ", "\[Alpha]"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}], " ", "+", " ", RowBox[{"y1", "[", "4", "]"}]}], ")"}], "\[Kappa]", " ", "\[Alpha]"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A2N2M2", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"8", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y1", "[", "1", "]"}], "+", RowBox[{"y1", "[", "2", "]"}], " ", "+", " ", RowBox[{"y1", "[", "3", "]"}]}], ")"}], "\[Kappa]", " ", "\[Alpha]"}], "+", " ", RowBox[{ RowBox[{"y1", "[", "4", "]"}], "\[Gamma]1", " ", "\[Kappa]", " ", "\[Alpha]"}]}]}]}], " ", "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Population", " ", "2"}], " ", "*)"}], "\n", RowBox[{"(*", " ", RowBox[{"A1N1M1", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\[IndentingNewLine]", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"1", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y2", "[", "1", "]"}], "\[Gamma]1", " ", "\[Alpha]"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}], " ", "+", " ", RowBox[{"y2", "[", "4", "]"}]}], ")"}], "\[Alpha]"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N2M2", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"4", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "1", "]"}], "+", RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}]}], ")"}], "\[Alpha]"}], "+", " ", RowBox[{ RowBox[{"y2", "[", "4", "]"}], "\[Gamma]1", " ", "\[Alpha]"}]}]}]}], " ", "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N1M1", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "1"}], " ", "&&", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}]}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "1", "]"}], "+", RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}]}], ")"}], "+", " ", RowBox[{ RowBox[{"y2", "[", "4", "]"}], "/", "\[Gamma]2"}]}]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A1N2M2", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "0"}], " ", "&&", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}]}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y2", "[", "1", "]"}], "/", "\[Gamma]2"}], "+", RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}], "+", " ", RowBox[{"y2", "[", "4", "]"}]}], ")"}]}]}]}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"purebred", " ", RowBox[{"w", "/", "hybrid"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "1"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", "0"}]}], ")"}], " ", "&&", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "2"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], " ", "\[Equal]", " ", "3"}]}], ")"}]}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "1", "]"}], "+", " ", RowBox[{"y2", "[", "4", "]"}]}], ")"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}]}], " ", ")"}], "/", "\[Gamma]h"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"hybrid", " ", "females"}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], " ", "\[Equal]", " ", "2"}], " ", "||", " ", RowBox[{ RowBox[{"Mod", "[", RowBox[{"i", ",", "4"}], "]"}], "\[Equal]", "3"}]}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", "1"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A2N1M1", " ", RowBox[{"w", "/", "M1N1"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "1"}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"5", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"y2", "[", "1", "]"}], "\[Gamma]1", " ", "\[Kappa]", " ", "\[Alpha]"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}], " ", "+", " ", RowBox[{"y2", "[", "4", "]"}]}], ")"}], "\[Kappa]", " ", "\[Alpha]"}]}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\n", " ", RowBox[{"(*", " ", RowBox[{"A2N2M2", " ", RowBox[{"w", "/", "M2N2"}]}], " ", "*)"}], "\n", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Mod", "[", RowBox[{"j", ",", "4"}], "]"}], "\[Equal]", " ", "0"}], ",", " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"8", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y2", "[", "1", "]"}], "+", RowBox[{"y2", "[", "2", "]"}], " ", "+", " ", RowBox[{"y2", "[", "3", "]"}]}], ")"}], "\[Kappa]", " ", "\[Alpha]"}], "+", " ", RowBox[{ RowBox[{"y2", "[", "4", "]"}], "\[Gamma]1", " ", "\[Kappa]", " ", "\[Alpha]"}]}]}]}], " ", "]"}], ",", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", "normalizations", " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"z1", "[", "i", "]"}], " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"xns1", "[", "j", "]"}], " ", RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"z2", "[", "i", "]"}], " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"xns2", "[", "j", "]"}], " ", RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"j", ",", " ", "8"}], "}"}]}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", " ", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"mating", " ", "table"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"F1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"mate1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", RowBox[{ RowBox[{"G1", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "/", RowBox[{"z1", "[", "i", "]"}]}]}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"F2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"mate2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", RowBox[{ RowBox[{"G2", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "/", RowBox[{"z2", "[", "i", "]"}]}]}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "8"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"j", ",", "8"}], "}"}]}], "]"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"recombination", " ", "pop", " ", "1"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "1", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "2", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "3", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "4", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "4"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "5", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "6", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "6"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "7", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "7"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x1t1", "[", "8", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"4", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"6", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"7", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F1", "[", RowBox[{"[", RowBox[{"8", ",", "8"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"recombination", " ", "pop", " ", "2"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "3", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "4", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "4"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "5", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "6", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "6"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "7", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "7"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"x2t1", "[", "8", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "8"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"1", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"2", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"3", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r2", "-", RowBox[{"r2", "*", "r1"}], "+", "r3", "+", RowBox[{"r2", "*", "r3"}], "-", RowBox[{"r1", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"4", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"5", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r3"}], "+", "r2", "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"r3", "*", "r2"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r1", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"6", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"r1", "-", RowBox[{"r1", "*", "r2"}], "+", "r3", "+", RowBox[{"r1", "*", "r3"}], "-", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r2", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", "r3", ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"7", ",", "8"}], "]"}], "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "3"}], ")"}], "*", RowBox[{"(", RowBox[{"3.0", "-", RowBox[{"2.0", "*", "r1"}], "-", RowBox[{"2.0", "*", "r2"}], "+", RowBox[{"r1", "*", "r2"}], "-", RowBox[{"2.0", "*", "r3"}], "+", RowBox[{"r1", "*", "r3"}], "+", RowBox[{"r2", "*", "r3"}]}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "1"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r1"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "2"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r2"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "4"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"(", RowBox[{"1.0", "-", "r3"}], ")"}], "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "5"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "6"}], "]"}], "]"}]}], "+", RowBox[{"4", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "7"}], "]"}], "]"}]}], "+", RowBox[{"8", "*", RowBox[{"F2", "[", RowBox[{"[", RowBox[{"8", ",", "8"}], "]"}], "]"}]}], "+", "0"}], ")"}]}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"allele", " ", "freqs", " ", "for", " ", "plotting"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"a12", " ", "=", " ", RowBox[{ RowBox[{"x1t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"n12", " ", "=", " ", RowBox[{ RowBox[{"x1t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"m12", " ", "=", " ", RowBox[{ RowBox[{"x1t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D1nm", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "5", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "6", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "7", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D1am", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "3", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "4", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "7", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D1an", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "2", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "4", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x1t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x1t1", "[", "6", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D1anm", " ", "=", " ", RowBox[{ RowBox[{"x1t1", "[", "8", "]"}], " ", "-", " ", RowBox[{"a12", " ", "n12", " ", "m12"}], " ", "-", " ", RowBox[{"a12", " ", "D1nm"}], " ", "-", " ", RowBox[{"n12", " ", "D1am"}], " ", "-", " ", RowBox[{"m12", " ", "D1an"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"a22", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"n22", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"m22", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D2nm", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "5", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "6", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "7", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D2am", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "3", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "4", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "7", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D2an", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "2", "]"}]}], ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "7", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}], ")"}]}], "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "4", "]"}]}], " ", ")"}], RowBox[{"(", RowBox[{ RowBox[{"x2t1", "[", "5", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "6", "]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"D2anm", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "8", "]"}], " ", "-", " ", RowBox[{"a22", " ", "n22", " ", "m22"}], " ", "-", " ", RowBox[{"a22", " ", "D2nm"}], " ", "-", " ", RowBox[{"n22", " ", "D2am"}], " ", "-", " ", RowBox[{"m22", " ", "D2an"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"n2m2p2", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "4", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "8", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"n1m1p2", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "1", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "5", "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"hybp2", " ", "=", " ", RowBox[{ RowBox[{"x2t1", "[", "2", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "3", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "6", "]"}], " ", "+", " ", RowBox[{"x2t1", "[", "7", "]"}]}]}], ";"}]}], "]"}], ";", " ", RowBox[{"(*", " ", "loop", " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtaba12", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "a12"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabn12", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "n12"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabn22", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "n22"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabD1nm", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "D1nm"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabD2nm", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "D2nm"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabn2m2p2", "[", RowBox[{"[", RowBox[{"k1", ",", "k2"}], "]"}], "]"}], " ", "=", " ", "n2m2p2"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabn1m1p2", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", "n1m1p2"}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"grandtabhybp2", "[", RowBox[{"[", RowBox[{"k1", ",", "k2"}], "]"}], "]"}], " ", "=", " ", "hybp2"}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Print", "[", RowBox[{ "\"\<\[Gamma]1=\>\"", ",", " ", "\[Gamma]1", ",", " ", "\"\< \[Kappa]=\>\"", ",", " ", "\[Kappa]", ",", " ", "\"\< run: \>\"", ",", " ", "k1", ",", " ", "\"\<,\>\"", ",", " ", "k2", ",", "\"\< gen = \>\"", ",", " ", "gen"}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"plots", " ", "pop", " ", "1"}], "*)"}], " ", "\[IndentingNewLine]", RowBox[{"a12plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}]}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"taba2p1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Red"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"m12plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabm2p1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Blue"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"n12plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabn2p1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Green"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D1amplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDamp1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Red"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D1anplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDanp1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Blue"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D1nmplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDnmp1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Green"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D1amnplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDamnp1", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Orange"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"plots", " ", "pop", " ", "2"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a22plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"taba2p2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Dashed", ",", "Red"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"m22plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabm2p2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Dashed", ",", " ", "Blue"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"n22plot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabn2p2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Dashed", ",", " ", "Green"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{"0", ",", " ", "1"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D2amplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDamp2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Dashed", ",", " ", "Red"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D2anplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDanp2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Dashed", ",", " ", "Blue"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D2nmplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDnmp2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", "Dashed", ",", " ", "Green"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"D2amnplot", "[", RowBox[{"[", RowBox[{"k1", ",", " ", "k2"}], "]"}], "]"}], " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"tabDamnp2", ",", " ", RowBox[{"Joined", " ", "\[Rule]", " ", "True"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"AbsoluteThickness", "[", "2", "]"}], ",", " ", "Dashed", ",", "Orange"}], "}"}]}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "0.25"}], ",", " ", "0.25"}], "}"}]}]}], "]"}]}]}], ",", " ", RowBox[{"{", RowBox[{"k1", ",", "7"}], "}"}]}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"k2", ",", " ", "2"}], "}"}]}], "]"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", ")"}]}]], "Input", CellChangeTimes->{{3.443832124313095*^9, 3.443832128792643*^9}, { 3.4801762399053926`*^9, 3.4801766988017073`*^9}, {3.480176733681239*^9, 3.480176836535383*^9}, {3.480176880648589*^9, 3.480176954327414*^9}, 3.4801771204805937`*^9, {3.480177175115868*^9, 3.4801772435409813`*^9}, { 3.4801773122992897`*^9, 3.48017742972188*^9}, {3.480177473734274*^9, 3.480177486937045*^9}, {3.481302347204659*^9, 3.481302438723559*^9}, { 3.491753927626013*^9, 3.491754131319603*^9}, {3.491754375578094*^9, 3.491754438374323*^9}, {3.4917545111907454`*^9, 3.491754703170113*^9}, { 3.491754739144246*^9, 3.4917553107234097`*^9}, {3.491755351219783*^9, 3.491755475312735*^9}, {3.491755512768173*^9, 3.491755649907276*^9}, { 3.491756188961295*^9, 3.4917561896611357`*^9}, {3.491756226155265*^9, 3.491756244530486*^9}, 3.491756374462903*^9, {3.491757147432742*^9, 3.491757159891451*^9}, {3.491757195391528*^9, 3.491757229280201*^9}, { 3.491757788084031*^9, 3.4917578447286053`*^9}, 3.4917579126690483`*^9, { 3.491757951207048*^9, 3.491758007780802*^9}, {3.491758097164484*^9, 3.4917581211704082`*^9}, {3.491758175522084*^9, 3.49175820199441*^9}, { 3.503422328017503*^9, 3.503422567056016*^9}, {3.503429994205681*^9, 3.503430148265085*^9}, 3.503443363264385*^9, 3.503443410972497*^9, { 3.50344344273806*^9, 3.5034435269650917`*^9}, 3.503443561671504*^9, { 3.503443737867126*^9, 3.503443760428548*^9}, {3.5034438355150013`*^9, 3.503443970988378*^9}, {3.503444078390135*^9, 3.5034441190723257`*^9}, { 3.5034458693812838`*^9, 3.5034460112188807`*^9}, 3.503446044789296*^9, { 3.5034788118946943`*^9, 3.503479110331469*^9}, {3.50347914225695*^9, 3.503479231800033*^9}, {3.5034999261406717`*^9, 3.503499959805025*^9}, { 3.514027260915021*^9, 3.514027602100378*^9}, {3.514027646461671*^9, 3.51402771055396*^9}, {3.514027794068149*^9, 3.514027967812831*^9}, { 3.5140280026558723`*^9, 3.514028560305862*^9}, {3.514029232894362*^9, 3.514029252289242*^9}, {3.5140295976064577`*^9, 3.5140296827111053`*^9}, { 3.514029769243894*^9, 3.514029781657205*^9}, {3.514030596708047*^9, 3.514030674858756*^9}, {3.514030717696323*^9, 3.5140308286544847`*^9}, { 3.5140311429500303`*^9, 3.514031192115156*^9}, {3.514031258952015*^9, 3.5140313503881493`*^9}, {3.514031380554289*^9, 3.514031406059251*^9}, { 3.5140314484121943`*^9, 3.514031451680059*^9}, 3.5140315138529882`*^9, { 3.5140316090697517`*^9, 3.514031949128314*^9}, {3.514032007102929*^9, 3.514032169620919*^9}, {3.5140326978567677`*^9, 3.514032703435986*^9}, { 3.5140327678914757`*^9, 3.51403281256848*^9}, {3.514033334423712*^9, 3.514033383665491*^9}, {3.514033512631987*^9, 3.514033544846099*^9}, { 3.514140155232731*^9, 3.514140169541915*^9}, {3.514140458580222*^9, 3.5141404721013002`*^9}, {3.514221305610467*^9, 3.514221312579563*^9}, { 3.514221349292633*^9, 3.514221355124131*^9}, {3.515945491446329*^9, 3.515945523860302*^9}, {3.515945563856113*^9, 3.515945761560706*^9}, { 3.515945816534878*^9, 3.515945937531435*^9}, {3.515945974155365*^9, 3.515946234351449*^9}, {3.515946334314558*^9, 3.5159463371631317`*^9}, { 3.5159463932519407`*^9, 3.515946393906062*^9}, {3.51594647576606*^9, 3.515946486726145*^9}, {3.515946661391719*^9, 3.515946816093258*^9}, { 3.521214152339319*^9, 3.521214231138836*^9}, {3.52121458713118*^9, 3.52121458736968*^9}, {3.521214618079221*^9, 3.5212149770248327`*^9}, { 3.521215015809986*^9, 3.5212151486385317`*^9}, {3.521217593779735*^9, 3.521217601748266*^9}, {3.521288718027377*^9, 3.52128872206948*^9}, 3.52146030652293*^9, {3.521460445705998*^9, 3.521460528162426*^9}, { 3.521460574815021*^9, 3.521460580326273*^9}, {3.521472561358705*^9, 3.521472583700378*^9}, {3.523714343506624*^9, 3.523714668946105*^9}, { 3.523714769480163*^9, 3.523714841274226*^9}, {3.523714872517419*^9, 3.523714877325761*^9}, {3.523716583862556*^9, 3.523716604232429*^9}, 3.523716999988557*^9, {3.523717317116232*^9, 3.523717342593493*^9}, { 3.523718095123952*^9, 3.523718207655731*^9}, {3.5237196142687597`*^9, 3.5237196457296743`*^9}, {3.523720201039669*^9, 3.523720217341668*^9}, { 3.523722202767844*^9, 3.523722240791367*^9}, {3.5237224054979677`*^9, 3.5237224639736147`*^9}, {3.523722800250311*^9, 3.523722833698604*^9}, { 3.523722951042984*^9, 3.523722981665699*^9}, {3.523796015957466*^9, 3.523796021194807*^9}, {3.52379605645716*^9, 3.5237961055605392`*^9}, { 3.523796200468729*^9, 3.523796201386476*^9}, {3.523796258184675*^9, 3.5237963385145807`*^9}, {3.523796469824966*^9, 3.523796471769545*^9}, { 3.523796506636196*^9, 3.5237965264918756`*^9}, {3.529769813080151*^9, 3.529769848531612*^9}, 3.5303750771922827`*^9, 3.530531383621875*^9, { 3.536048638637637*^9, 3.536048709006999*^9}, 3.5360512374321127`*^9, 3.536051326220503*^9, 3.536053002786614*^9, {3.5363005990967493`*^9, 3.536300609730329*^9}, {3.536300643713402*^9, 3.536300644435194*^9}, 3.536300675675618*^9, {3.536301001518881*^9, 3.536301008228723*^9}, 3.5363010433662167`*^9, {3.536477622150187*^9, 3.53647778758324*^9}, { 3.5364778177299013`*^9, 3.536477829388431*^9}, {3.5364778673466187`*^9, 3.53647792963899*^9}, {3.536572833312173*^9, 3.5365728488476963`*^9}, { 3.5365729130984087`*^9, 3.536572916782699*^9}, 3.5365729512758007`*^9, { 3.536653261309209*^9, 3.536653315145152*^9}, {3.5369110353816633`*^9, 3.536911128076499*^9}, 3.5369111862429047`*^9, {3.5369112414145823`*^9, 3.536911247887134*^9}, {3.536911285603506*^9, 3.536911316887877*^9}, { 3.536911739659046*^9, 3.536911750660429*^9}, {3.536916713636198*^9, 3.53691672018804*^9}, {3.5369172685215807`*^9, 3.5369172728273067`*^9}, { 3.537166144759726*^9, 3.537166174251789*^9}, {3.537166226918206*^9, 3.537166396528126*^9}, 3.537166441397664*^9, {3.537168898108651*^9, 3.537168899998633*^9}, {3.53716933245805*^9, 3.537169339148738*^9}, { 3.537179859927849*^9, 3.537179906677608*^9}, {3.5371844635813513`*^9, 3.537184483023384*^9}, {3.537184648434805*^9, 3.537184654175001*^9}, 3.540015962967949*^9, 3.5400161886195393`*^9, 3.540017828394025*^9, { 3.547968682181552*^9, 3.547968683096402*^9}, 3.547968726038128*^9, { 3.548006013834094*^9, 3.5480060168993073`*^9}, {3.5480653275401363`*^9, 3.548065328600849*^9}, {3.564926314730556*^9, 3.5649263458345127`*^9}, { 3.564926380854205*^9, 3.564926390609222*^9}, {3.564926500471252*^9, 3.564926500894775*^9}, {3.5649265982737827`*^9, 3.5649266025266037`*^9}, { 3.564926841001355*^9, 3.564926851031899*^9}, {3.568646707368765*^9, 3.5686467090613213`*^9}, 3.568646744955523*^9, {3.570899657496279*^9, 3.57089967746838*^9}, {3.570899888506546*^9, 3.570899899448162*^9}, { 3.570899930868803*^9, 3.5708999797637987`*^9}, {3.570900045444524*^9, 3.57090004567065*^9}, 3.5709001116548433`*^9, {3.5709001749606953`*^9, 3.570900236846974*^9}, 3.570900337146482*^9, {3.570900367307336*^9, 3.57090036753386*^9}, {3.573392113080944*^9, 3.573392147200727*^9}, { 3.57339218670393*^9, 3.5733921968181257`*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["sample parameters for a run", "Section", CellChangeTimes->{{3.564926352915195*^9, 3.5649263616751823`*^9}, { 3.5733922230015507`*^9, 3.573392228673739*^9}, {3.573392378583311*^9, 3.573392380807152*^9}}], Cell[CellGroupData[{ Cell["\<\ m = 0.01, alpha = 1.5, s = 0.5 \[Gamma]1 = \[Gamma]2\ \>", "Subsection", CellChangeTimes->{{3.516026547065349*^9, 3.516026584717346*^9}, 3.516100148216528*^9, {3.516468142607476*^9, 3.516468144197953*^9}, 3.516468196896242*^9, 3.516636727514839*^9, {3.523196711560729*^9, 3.523196711831583*^9}, {3.564926521038846*^9, 3.564926543009675*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"m", " ", "=", " ", "0.01"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"migration", " ", "rate"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Gamma]h", " ", "=", " ", "1"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"learning", " ", "hybrid", " ", RowBox[{"males", ":", " ", "1"}]}], " ", "=", " ", RowBox[{"no", " ", "learning"}]}], ",", " ", RowBox[{"\[Gamma]", " ", "=", " ", RowBox[{"as", " ", "much", " ", "as", " ", "wrong", " ", "type"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Alpha]", " ", "=", " ", "1.5"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ "female", " ", "preference", " ", "strength", " ", "for", " ", "A1"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s", " ", "=", " ", "0.5"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"selection", " ", "against", " ", "hybrids"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"perturb", " ", "=", " ", "30000"}], ";", " ", RowBox[{"(*", " ", RowBox[{"perturbation", " ", "-", " ", RowBox[{ "set", " ", "to", " ", "work", " ", "only", " ", "when", " ", "A2", " ", "is", " ", "high"}]}], " ", "*)"}], "\[IndentingNewLine]", "FemReinfLearnLoopProg"}]}], "Input", CellChangeTimes->{{3.514032222202791*^9, 3.5140322343509073`*^9}, { 3.5140323878894997`*^9, 3.514032468174387*^9}, 3.514032647605874*^9, { 3.51403278290648*^9, 3.514032815756278*^9}, {3.5140328815729103`*^9, 3.514032975487468*^9}, {3.514033016671938*^9, 3.5140333200079803`*^9}, { 3.514033419636463*^9, 3.5140334318486967`*^9}, {3.51413402623825*^9, 3.514134026836852*^9}, {3.514134134058956*^9, 3.514134160503242*^9}, 3.5141343794817877`*^9, 3.514134436632209*^9, {3.514134524831188*^9, 3.514134566931458*^9}, {3.514135297589431*^9, 3.514135315556625*^9}, 3.514135362845748*^9, 3.514135420572118*^9, {3.514140185767137*^9, 3.514140194350168*^9}, {3.514140257439637*^9, 3.5141402588372707`*^9}, { 3.5141403346529293`*^9, 3.514140445712631*^9}, {3.514215112669895*^9, 3.51421511316461*^9}, {3.5142151607687283`*^9, 3.514215165487211*^9}, { 3.51421533810504*^9, 3.514215338974161*^9}, 3.514215599994832*^9, 3.514215953442871*^9, {3.5142159918068323`*^9, 3.514215996933469*^9}, { 3.514216121595604*^9, 3.514216128152278*^9}, 3.514216444300692*^9, 3.514216833026417*^9, 3.514217080467054*^9, {3.514217239897188*^9, 3.514217240862001*^9}, 3.514217676026309*^9, 3.5142178800641317`*^9, 3.514217993637972*^9, {3.514218142872869*^9, 3.514218144005908*^9}, 3.514218435730377*^9, 3.514218558158923*^9, 3.514218713812604*^9, { 3.514220809838345*^9, 3.51422081404235*^9}, {3.514298354135023*^9, 3.514298361860251*^9}, 3.514298395752501*^9, {3.5142988579775267`*^9, 3.514298859189425*^9}, {3.514298977102168*^9, 3.514298977188586*^9}, { 3.514299092784359*^9, 3.514299094048225*^9}, {3.514299360050045*^9, 3.5142993765826807`*^9}, {3.5142995055046883`*^9, 3.51429952026224*^9}, 3.514299617476015*^9, {3.514299695851323*^9, 3.5142996964501867`*^9}, 3.51429978095313*^9, 3.514299864107677*^9, {3.5143000796838303`*^9, 3.51430008162372*^9}, {3.514300199710643*^9, 3.514300201582322*^9}, { 3.514300405471455*^9, 3.514300408827907*^9}, {3.514300528948205*^9, 3.514300529114991*^9}, {3.5143006406146927`*^9, 3.51430064118293*^9}, { 3.514300801375379*^9, 3.514300841804954*^9}, {3.514300960843096*^9, 3.5143009613776817`*^9}, {3.514301079655219*^9, 3.514301080949728*^9}, { 3.514302325708337*^9, 3.514302343464243*^9}, {3.514302585229569*^9, 3.514302588220537*^9}, {3.5143026986811666`*^9, 3.5143027024572363`*^9}, { 3.514302935263352*^9, 3.514302936368393*^9}, {3.514303303030798*^9, 3.514303306517708*^9}, {3.5143033958882227`*^9, 3.514303397870742*^9}, { 3.514303605713953*^9, 3.514303611327189*^9}, {3.514303705365494*^9, 3.514303709584547*^9}, {3.514303865619038*^9, 3.51430387115067*^9}, { 3.5143040923470087`*^9, 3.514304096982119*^9}, {3.5143047208930798`*^9, 3.5143047322567377`*^9}, 3.514304881463066*^9, 3.514305031272566*^9, 3.5143053369387083`*^9, {3.515946215075808*^9, 3.5159462766217747`*^9}, { 3.515946501743834*^9, 3.515946503574054*^9}, {3.51594685039323*^9, 3.515946851015648*^9}, {3.515947210237327*^9, 3.515947250262351*^9}, { 3.5160265138198557`*^9, 3.5160265400815887`*^9}, 3.5160265890794086`*^9, { 3.516037661297078*^9, 3.516037669831325*^9}, {3.5161001544267473`*^9, 3.5161001644260273`*^9}, {3.5163828653473797`*^9, 3.5163828694686947`*^9}, {3.516468153947225*^9, 3.516468154223599*^9}, { 3.516468199410553*^9, 3.516468206362729*^9}, 3.516707680942259*^9, 3.547967932886571*^9, {3.573392277431497*^9, 3.573392279174079*^9}}], Cell["\<\ There are many potential things you can look at as output. Here is a sample.\ \ \>", "Text", CellChangeTimes->{{3.57339235502435*^9, 3.573392368110811*^9}}], Cell[BoxData["grandtabn12"], "Input", CellChangeTimes->{{3.5159465810210876`*^9, 3.5159465835486193`*^9}}], Cell[BoxData["grandtabn22"], "Input", CellChangeTimes->{{3.515946862937921*^9, 3.515946869143119*^9}, 3.573392343327477*^9}], Cell[BoxData["grandtabD2nm"], "Input"], Cell[BoxData["grandtabkappa"], "Input"], Cell[BoxData["grandtabgamma"], "Input"] }, Open ]] }, Closed]] }, Open ]] }, Evaluator->"Kernel 6", WindowSize->{1082, 1119}, WindowMargins->{{78, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{1, Automatic}, FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 6, \ 2010)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 333, 6, 76, "Title"], Cell[903, 30, 210, 3, 26, "Subsubtitle"], Cell[CellGroupData[{ Cell[1138, 37, 92, 1, 67, "Section"], Cell[1233, 40, 253260, 5674, 7123, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[254530, 5719, 214, 3, 37, "Section"], Cell[CellGroupData[{ Cell[254769, 5726, 366, 6, 34, "Subsection"], Cell[255138, 5734, 4900, 86, 103, "Input"], Cell[260041, 5822, 168, 4, 26, "Text"], Cell[260212, 5828, 107, 1, 27, "Input"], Cell[260322, 5831, 129, 2, 27, "Input"], Cell[260454, 5835, 38, 0, 27, "Input"], Cell[260495, 5837, 39, 0, 27, "Input"], Cell[260537, 5839, 39, 0, 27, "Input"] }, Open ]] }, Closed]] }, Open ]] } ] *) (* End of internal cache information *)