(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 34552, 978] NotebookOptionsPosition[ 33061, 931] NotebookOutlinePosition[ 33402, 946] CellTagsIndexPosition[ 33359, 943] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"numtries", "=", "3000"}], RowBox[{"(*", " ", RowBox[{"times", " ", "to", " ", "try", " ", "random", " ", "numbers"}], " ", "*)"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ StyleBox[" ", "Subsubtitle"], StyleBox[ RowBox[{ RowBox[{"Define", " ", "lines"}], ",", " ", RowBox[{"as", " ", "vectors", " ", "in", " ", RowBox[{"I", "^", RowBox[{"{", RowBox[{"1", ",", "6"}], "}"}]}]}]}], "Subsubtitle"], StyleBox[" ", "Subsubtitle"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"line", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "27"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "7"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5500779031465406`*^9, 3.5500779631439724`*^9}, { 3.5500782898726606`*^9, 3.550078290438693*^9}, {3.5500787502069902`*^9, 3.550078759880543*^9}, {3.5501505487388945`*^9, 3.5501505804987106`*^9}, { 3.5501508091077867`*^9, 3.550150814530097*^9}, 3.5501516970925765`*^9, { 3.5501519120878735`*^9, 3.5501519161151037`*^9}, {3.550152285227216*^9, 3.550152292802649*^9}, {3.5747944072657905`*^9, 3.5747944111220107`*^9}, { 3.5747959474738855`*^9, 3.574795950370051*^9}}], Cell[BoxData["3000"], "Output", CellChangeTimes->{3.5501505836388903`*^9, 3.5501508153521433`*^9, 3.5501517050650325`*^9, 3.5501517561159525`*^9, 3.5501518360165224`*^9, 3.550151883955264*^9, 3.5501519219314365`*^9, 3.5501521953870773`*^9, 3.574795098052301*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"s", "=", "1"}], ";"}]], "Input", CellChangeTimes->{{3.550077989450477*^9, 3.5500779967988977`*^9}, 3.5500788774792695`*^9}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "i"}], "]"}], "]"}], "=", "1"}], ";", RowBox[{"s", "++"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.550078009468622*^9, 3.550078117775817*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "7"}], "]"}], "]"}], "=", "2"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", "1"}], ",", RowBox[{"k", "\[LessEqual]", "6"}], ",", RowBox[{"k", "++"}], ",", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"k", "\[NotEqual]", "i"}], ",", RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "k"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}]}], "]"}]}], "]"}], ";", "\[IndentingNewLine]", RowBox[{"s", "++"}]}]}], "\[IndentingNewLine]", "]"}]], "Input", CellChangeTimes->{{3.5500781367879047`*^9, 3.550078161567322*^9}, { 3.5500782050068064`*^9, 3.5500782788010273`*^9}, {3.5500783334581532`*^9, 3.550078367801118*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "\[LessEqual]", "6"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "7"}], "]"}], "]"}], "=", "1"}], ";", " ", RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "i"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", " ", RowBox[{ RowBox[{"line", "[", RowBox[{"[", RowBox[{"s", ",", "j"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", " ", RowBox[{"s", "++"}]}]}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.550078410662569*^9, 3.550078481401615*^9}, { 3.550078527434248*^9, 3.5500785518096423`*^9}, {3.550078602080518*^9, 3.5500786036896095`*^9}, {3.5500786632580166`*^9, 3.550078668945342*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"line", "//", "MatrixForm"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.550078582677408*^9, 3.550078588463739*^9}, { 3.550078779831684*^9, 3.550078785231993*^9}}], Cell[BoxData[{ RowBox[{ StyleBox[ RowBox[{ RowBox[{"Define", " ", "roots"}], ",", " ", RowBox[{"as", " ", "vectors", " ", "in", " ", RowBox[{"I", "^", RowBox[{"{", RowBox[{"1", ",", "6"}], "}"}]}]}]}], "Subsubtitle"], StyleBox[" ", "Subsubtitle"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"root", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "36"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "7"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5500787924634066`*^9, 3.550078838296028*^9}, 3.5747944214656024`*^9, {3.5747959784576573`*^9, 3.5747959815698357`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"r", "=", "1"}], ";"}]], "Input", CellChangeTimes->{{3.550078849416664*^9, 3.5500788733110313`*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "\[LessEqual]", "6"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "i"}], "]"}], "]"}], "=", "1"}], ";", RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "j"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", " ", RowBox[{"r", "++"}]}]}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.5500789124982724`*^9, 3.5500789358396077`*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "\[LessEqual]", "6"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"k", "=", RowBox[{"j", "+", "1"}]}], ",", RowBox[{"k", "\[LessEqual]", "6"}], ",", RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "7"}], "]"}], "]"}], "=", "1"}], ";", RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "i"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "j"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", " ", RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"r", ",", "k"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}], ";", RowBox[{"r", "++"}]}]}], "]"}]}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.5500790145281086`*^9, 3.550079069923277*^9}, { 3.5500791185050554`*^9, 3.550079119601118*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"36", ",", "7"}], "]"}], "]"}], "=", "2"}], ";", " ", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{"root", "[", RowBox[{"[", RowBox[{"36", ",", "i"}], "]"}], "]"}], "=", RowBox[{"-", "1"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5500791964065113`*^9, 3.5500792411750717`*^9}, 3.5500792812793655`*^9, 3.550079355710623*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"root", "//", "MatrixForm"}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"Lorenzian", " ", "dot", " ", "product", " ", "in", " ", RowBox[{"I", "^", RowBox[{"{", RowBox[{"1", ",", "6"}], "}"}], " "}]}], "Subsubtitle"]}]], "Input", CellChangeTimes->{{3.550079422063418*^9, 3.550079428358778*^9}, { 3.5500795086713715`*^9, 3.5500795139676743`*^9}, {3.5747944367524767`*^9, 3.5747944375125203`*^9}, {3.5747959954656305`*^9, 3.57479599896883*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"mdot", "[", RowBox[{"x_", ",", "y_"}], "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "7", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "7", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "1", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "1", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "2", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "3", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "3", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "4", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "4", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "5", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "5", "]"}], "]"}]}], "-", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "6", "]"}], "]"}], "*", RowBox[{"y", "[", RowBox[{"[", "6", "]"}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.5500796091281176`*^9, 3.5500796152634687`*^9}, { 3.55007966284019*^9, 3.5500796819282813`*^9}, {3.5500797309680862`*^9, 3.550079760879797*^9}, {3.5500797966718445`*^9, 3.550079867850916*^9}, { 3.550079913120505*^9, 3.550079937534901*^9}, {3.5500799782252283`*^9, 3.550080011912155*^9}, {3.5500802388661366`*^9, 3.550080246608579*^9}, { 3.5500805758644114`*^9, 3.550080591553309*^9}, 3.550080793225844*^9, { 3.5500808540573235`*^9, 3.550080944841516*^9}, 3.550333542541541*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Prods", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "36"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "27"}], "}"}]}], "]"}]}], ";"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.550081756549943*^9, 3.550081787823732*^9}, { 3.550082107071992*^9, 3.550082124992017*^9}, 3.5500821570328493`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", RowBox[{"n", "\[LessEqual]", "36"}], ",", RowBox[{"n", "++"}], ",", " ", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"m", "=", "1"}], ",", RowBox[{"m", "\[LessEqual]", "27"}], ",", RowBox[{"m", "++"}], ",", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"Prods", "[", RowBox[{"[", RowBox[{"n", ",", "m"}], "]"}], "]"}], "=", RowBox[{"mdot", "[", " ", RowBox[{ RowBox[{"root", "[", RowBox[{"[", "n", "]"}], "]"}], ",", RowBox[{"line", "[", RowBox[{"[", "m", "]"}], "]"}]}], " ", "]"}]}]}], " ", "]"}]}], "]"}], " ", "*)"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.550079765753076*^9, 3.5500797855202065`*^9}, { 3.5500806927370963`*^9, 3.5500806928151007`*^9}, {3.5500807287521563`*^9, 3.5500807726556673`*^9}, {3.550080810202815*^9, 3.550080829904942*^9}, { 3.550080970705995*^9, 3.5500810261361656`*^9}, {3.5500814410328965`*^9, 3.5500815355063*^9}, {3.5500817251551476`*^9, 3.550081750361589*^9}, { 3.5500817967222404`*^9, 3.5500818139142237`*^9}, {3.550081891377655*^9, 3.550081909008663*^9}, {3.550082026410378*^9, 3.550082041962268*^9}, { 3.5500821311783705`*^9, 3.5500821371447115`*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Prods", "//", "MatrixForm"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.550080343182103*^9, 3.5500803553668003`*^9}, { 3.550081854855566*^9, 3.5500818611829276`*^9}, {3.5500819351521587`*^9, 3.550082009759426*^9}, {3.550082050791773*^9, 3.5500820547910013`*^9}, { 3.5500821177356014`*^9, 3.550082141257947*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"broot", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "7"}], "}"}]}], "]"}]}], " ", ";"}], RowBox[{"(*", " ", RowBox[{ RowBox[{ "a", " ", "basis", " ", "of", " ", "basis", " ", "roots", " ", "of", " ", "E6"}], ",", " ", RowBox[{"the", " ", "first", " ", "6", " ", "roots", " ", "chosen"}]}], " ", "*)"}]}]], "Input", CellChangeTimes->{{3.550083151273717*^9, 3.550083220704688*^9}, { 3.5747945330729856`*^9, 3.574794552865118*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "1", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "21", "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5500832250879383`*^9, 3.550083244785065*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "2", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.550083248461275*^9, 3.550083262376071*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "3", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "6", "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5500832640111647`*^9, 3.55008327808797*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "4", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "10", "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.550083279793068*^9, 3.550083289848643*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "5", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "13", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "6", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", "15", "]"}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.550083291461735*^9, 3.5500833221684914`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"ep", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "21"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "12"}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"edges", "-", "to", "-", "endpoints"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ StyleBox["Paths", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["from", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["the", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["center", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["of", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["the", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["tree", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["T", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["to", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], StyleBox["the", "Subsubtitle"], StyleBox[" ", "Subsubtitle"], RowBox[{ StyleBox["endpoints", "Subsubtitle"], StyleBox[".", "Subsubtitle"], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"A", ",", "B", ",", RowBox[{"C", " ", "central", " ", "branches"}], ",", " ", RowBox[{ RowBox[{"A0", "..."}], "C1", " ", "the", " ", "next", " ", "6"}], ",", " ", RowBox[{ RowBox[{"A00", "..."}], "C11", " ", "the", " ", "next", " ", "12"}]}], " ", "*)"}]}]}]}], "Input", CellChangeTimes->{{3.5500833700442295`*^9, 3.5500834421523542`*^9}, { 3.5747944763607426`*^9, 3.5747944773457985`*^9}, {3.574796021618126*^9, 3.574796049089697*^9}}], Cell[BoxData[ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "0"}], ",", RowBox[{"i", "<", "3"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "1"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "1"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "1"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "2"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "1"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "3"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "1"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "4"}]}], "]"}], "]"}], "=", "1"}]}]}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"A", " ", "B", " ", "C"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "2"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "1"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "2"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "2"}]}], "]"}], "]"}], "=", "1"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"A0", " ", "B0", " ", "C0"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "3"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "3"}]}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "3"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "4"}]}], "]"}], "]"}], "=", "1"}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"A1", " ", "B1", " ", "C1"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "4"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "1"}]}], "]"}], "]"}], "=", "1"}], ";", RowBox[{"(*", " ", RowBox[{"A00", " ", "B00", " ", "C00"}], " ", "*)"}], "\n", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "5"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "2"}]}], "]"}], "]"}], "=", "1"}], ";", RowBox[{"(*", " ", RowBox[{"A01", " ", "B01", " ", "C01"}], " ", "*)"}], "\n", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "6"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "3"}]}], "]"}], "]"}], "=", "1"}], ";", RowBox[{"(*", " ", RowBox[{"A10", " ", "B10", " ", "C10"}], " ", "*)"}], "\n", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{ RowBox[{ RowBox[{"7", "*", "i"}], "+", "7"}], ",", RowBox[{ RowBox[{"4", "*", "i"}], "+", "4"}]}], "]"}], "]"}], "=", "1"}]}]}], " ", RowBox[{"(*", " ", RowBox[{"A11", " ", "B11", " ", "C11"}], " ", "*)"}], "]"}]], "Input", CellChangeTimes->{{3.550083449492774*^9, 3.550083477816394*^9}, { 3.550083621616619*^9, 3.5500837919813633`*^9}, {3.5500838240621977`*^9, 3.550083836792926*^9}, {3.5500838852316966`*^9, 3.5500838860487432`*^9}, { 3.5500839381607237`*^9, 3.5500839642312155`*^9}, {3.5500840274618316`*^9, 3.5500840604687195`*^9}, {3.5500842867826643`*^9, 3.550084295754177*^9}, 3.5501412507390795`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"weight", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}], ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "1", ",", "1", ",", "1", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "1", ",", RowBox[{"-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0"}], "}"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.5500840842770815`*^9, 3.550084280280292*^9}, { 3.55014131001847*^9, 3.550141494181003*^9}, {3.5501415308451004`*^9, 3.5501415742745843`*^9}, {3.550141627572633*^9, 3.550141657331335*^9}, { 3.5501416925493493`*^9, 3.550141744948346*^9}}], Cell[BoxData[""], "Input", CellChangeTimes->{ 3.55008422900836*^9, 3.550084282487418*^9, {3.550141339481155*^9, 3.550141349952754*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"the", " ", "next", " ", "6", " ", "roots"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5747945158590016`*^9, 3.5747945223283715`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"fnum", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"num", "=", "1"}], ",", RowBox[{"num", "\[LessEqual]", "numtries"}], ",", RowBox[{"num", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{"fnum", "[", RowBox[{"[", "i", "]"}], "]"}], "=", RowBox[{"RandomInteger", "[", RowBox[{"{", RowBox[{"1", ",", "36"}], "}"}], "]"}]}]}], " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{ RowBox[{"fnum", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "3", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "4", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "5", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "6", "]"}], "]"}]}], "]"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"froot", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "7"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{"froot", "[", RowBox[{"[", "i", "]"}], "]"}], "=", RowBox[{"root", "[", RowBox[{"[", RowBox[{"fnum", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}]}]}], " ", "]"}], ";", "\[IndentingNewLine]", RowBox[{"l", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"e", ",", "21"}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "27"}], "}"}], ",", RowBox[{"{", RowBox[{"f", ",", "6"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"e", "=", "1"}], ",", RowBox[{"e", "\[LessEqual]", "21"}], ",", RowBox[{"e", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"s", "=", "1"}], ",", RowBox[{"s", "\[LessEqual]", "27"}], ",", RowBox[{"s", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"f", "=", "1"}], ",", RowBox[{"f", "\[LessEqual]", "6"}], ",", RowBox[{"f", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"l", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "f"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{"e", ",", RowBox[{"6", "+", "f"}]}], "]"}], "]"}], "*", RowBox[{"mdot", "[", RowBox[{ RowBox[{"froot", "[", RowBox[{"[", "f", "]"}], "]"}], ",", RowBox[{"line", "[", RowBox[{"[", "s", "]"}], "]"}]}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "6"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"l", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "f"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"l", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "f"}], "]"}], "]"}], "-", RowBox[{ RowBox[{"ep", "[", RowBox[{"[", RowBox[{"e", ",", "i"}], "]"}], "]"}], "*", RowBox[{"mdot", "[", RowBox[{ RowBox[{"froot", "[", RowBox[{"[", "f", "]"}], "]"}], ",", RowBox[{"weight", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "*", RowBox[{"mdot", "[", RowBox[{ RowBox[{"broot", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"line", "[", RowBox[{"[", "s", "]"}], "]"}]}], "]"}]}]}]}], ";"}]}], " ", "\[IndentingNewLine]", "]"}]}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"q", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"e", ",", "21"}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "27"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "21"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"e", "=", "1"}], ",", RowBox[{"e", "\[LessEqual]", "21"}], ",", RowBox[{"e", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"s", "=", "1"}], ",", RowBox[{"s", "\[LessEqual]", "27"}], ",", RowBox[{"s", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<=", "6"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"k", "=", RowBox[{ RowBox[{"i", "*", RowBox[{ RowBox[{"(", RowBox[{"i", "-", "1"}], ")"}], "/", "2"}]}], "+", "j"}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"q", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "k"}], "]"}], "]"}], "=", RowBox[{ RowBox[{"l", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "i"}], "]"}], "]"}], "*", RowBox[{"l", "[", RowBox[{"[", RowBox[{"e", ",", "s", ",", "j"}], "]"}], "]"}]}]}], ";"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"Q", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"e", ",", "21"}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "21"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"e", "=", "1"}], ",", RowBox[{"e", "\[LessEqual]", "21"}], ",", RowBox[{"e", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"s", "=", "1"}], ",", RowBox[{"s", "\[LessEqual]", "27"}], ",", RowBox[{"s", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"Q", "[", RowBox[{"[", "e", "]"}], "]"}], "+=", RowBox[{ RowBox[{"q", "[", RowBox[{"[", RowBox[{"e", ",", "s"}], "]"}], "]"}], "/", "6"}]}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", RowBox[{"mrank", "=", RowBox[{"MatrixRank", "[", "Q", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"mrank", "\[Equal]", "21"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"Det", "[", "Q", "]"}], ",", "\"\< \>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "1", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "2", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "3", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "4", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "5", "]"}], "]"}], ",", "\"\<-\>\"", ",", RowBox[{"fnum", "[", RowBox[{"[", "6", "]"}], "]"}]}], "]"}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], "Input", CellChangeTimes->{{3.5501437283367896`*^9, 3.5501437592335567`*^9}, { 3.550150780779166*^9, 3.5501508277788544`*^9}, {3.5501508840430727`*^9, 3.550150884996127*^9}, {3.55015104047602*^9, 3.550151071966821*^9}, { 3.5501511055307407`*^9, 3.5501511311502066`*^9}, {3.550151163991085*^9, 3.5501511660862045`*^9}, 3.550151250789049*^9, 3.5501513686337895`*^9, { 3.550151402279714*^9, 3.5501514686395097`*^9}, {3.5501515224435873`*^9, 3.5501515406876307`*^9}, {3.5501516086625185`*^9, 3.5501516506469197`*^9}, {3.550151726950284*^9, 3.550151743755245*^9}, { 3.550151790974946*^9, 3.5501518293191395`*^9}, {3.5501518670662985`*^9, 3.5501518981610765`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TimeUsed", "[", "]"}]], "Input", CellChangeTimes->{{3.5501520080783634`*^9, 3.5501520117455735`*^9}, { 3.550152091355127*^9, 3.550152097305467*^9}, {3.5501521488104134`*^9, 3.5501521742498684`*^9}}], Cell[BoxData["4336.7339999999995`"], "Output", CellChangeTimes->{ 3.550152016346837*^9, 3.5501520979225025`*^9, {3.55015214962346*^9, 3.550152178038085*^9}, 3.550152273800562*^9, 3.5747976441629305`*^9}] }, Open ]] }, WindowSize->{905, 954}, WindowMargins->{{Automatic, 0}, {0, Automatic}}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 1296, 30, 92, "Input"], Cell[1878, 54, 272, 4, 30, "Output"] }, Open ]], Cell[2165, 61, 163, 4, 31, "Input"], Cell[2331, 67, 408, 12, 31, "Input"], Cell[2742, 81, 1047, 28, 112, "Input"], Cell[3792, 111, 1129, 30, 72, "Input"], Cell[4924, 143, 212, 4, 31, "Input"], Cell[5139, 149, 691, 19, 52, "Input"], Cell[5833, 170, 135, 3, 31, "Input"], Cell[5971, 175, 818, 23, 72, "Input"], Cell[6792, 200, 1442, 40, 92, "Input"], Cell[8237, 242, 578, 17, 31, "Input"], Cell[8818, 261, 566, 12, 72, "Input"], Cell[9387, 275, 1651, 46, 52, "Input"], Cell[11041, 323, 432, 12, 31, "Input"], Cell[11476, 337, 1412, 31, 92, "Input"], Cell[12891, 370, 367, 6, 31, "Input"], Cell[13261, 378, 631, 18, 31, "Input"], Cell[13895, 398, 251, 7, 31, "Input"], Cell[14149, 407, 248, 7, 31, "Input"], Cell[14400, 416, 249, 7, 31, "Input"], Cell[14652, 425, 249, 7, 31, "Input"], Cell[14904, 434, 437, 13, 52, "Input"], Cell[15344, 449, 1697, 48, 72, "Input"], Cell[17044, 499, 4193, 126, 172, "Input"], Cell[21240, 627, 1226, 31, 52, "Input"], Cell[22469, 660, 142, 3, 31, "Input"], Cell[22614, 665, 183, 3, 31, "Input"], Cell[22800, 670, 9779, 245, 792, "Input"], Cell[CellGroupData[{ Cell[32604, 919, 228, 4, 31, "Input"], Cell[32835, 925, 210, 3, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)