(* 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[ 115246, 2943] NotebookOptionsPosition[ 111625, 2823] NotebookOutlinePosition[ 111968, 2838] CellTagsIndexPosition[ 111925, 2835] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ StyleBox[ RowBox[{"Method", " ", "of", " ", "Least", " ", "Squares"}], "Title"], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{ RowBox[{"Kepler", "'"}], "s", " ", "Third", " ", "Law"}], ",", " ", RowBox[{ RowBox[{"Zipf", "'"}], "s", " ", "Law", " ", RowBox[{"(", RowBox[{"word", " ", "frequencies"}], ")"}]}]}], "Title"]}], "Input", CellChangeTimes->{{3.479591395375*^9, 3.479591400796875*^9}, { 3.479653177953125*^9, 3.479653182984375*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"periodT", " ", "=", " ", RowBox[{"{", RowBox[{ "0.2408467", ",", "0.61519726", ",", "1.0000174", ",", "1.8808476", ",", "11.862615", ",", "29.447498", ",", "84.016846", ",", " ", "164.79132"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"semimajoraxisL", " ", "=", " ", RowBox[{"{", RowBox[{ "0.387", ",", " ", "0.723", ",", " ", "1.000", ",", " ", "1.524", ",", " ", "5.203", ",", "9.539", ",", "19.182", ",", " ", "30.06"}], "}"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.47959141165625*^9, 3.479591476125*^9}, { 3.479591678265625*^9, 3.47959170821875*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "0.387`", ",", "0.723`", ",", "1.`", ",", "1.524`", ",", "5.203`", ",", "9.539`", ",", "19.182`", ",", "30.06`"}], "}"}]], "Output", CellChangeTimes->{3.479591477890625*^9, 3.4795917085*^9, 3.479678130078125*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"LogT", " ", "=", " ", RowBox[{"Log", "[", "periodT", "]"}]}], "\[IndentingNewLine]", RowBox[{"LogL", " ", "=", " ", RowBox[{"Log", "[", "semimajoraxisL", "]"}]}]}], "Input", CellChangeTimes->{{3.47959148053125*^9, 3.479591502828125*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "1.4235946474737198`"}], ",", RowBox[{"-", "0.4858123146367721`"}], ",", "0.000017399848621701225`", ",", "0.6317225263028842`", ",", "2.4733918583042294`", ",", "3.3826089489333198`", ",", "4.431017326355375`", ",", "5.1046799461853505`"}], "}"}]], "Output", CellChangeTimes->{{3.479591499421875*^9, 3.479591503453125*^9}, 3.479591709859375*^9, 3.479678131953125*^9, 3.479678274296875*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "0.9493305859523552`"}], ",", RowBox[{"-", "0.3243460568233724`"}], ",", "0.`", ",", "0.42133845726445457`", ",", "1.6492353823081665`", ",", "2.255388658163071`", ",", "2.9539723393057593`", ",", "3.4031953843248286`"}], "}"}]], "Output", CellChangeTimes->{{3.479591499421875*^9, 3.479591503453125*^9}, 3.479591709859375*^9, 3.479678131953125*^9, 3.479678274296875*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"x", " ", "=", " ", "LogL"}], ";", " ", RowBox[{"y", " ", "=", " ", "LogT"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xbar", " ", "=", " ", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}], " ", "/", " ", "8"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sigmaxsq", " ", "=", " ", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", "xbar"}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}], "/", "8"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"8", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], " ", "-", " ", RowBox[{"8", " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"8", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{"b", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "8"}], " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], " ", "+", " ", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"8", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.479591547625*^9, 3.47959159640625*^9}, { 3.47959163471875*^9, 3.47959164140625*^9}, {3.47959171378125*^9, 3.47959175015625*^9}, {3.47959178459375*^9, 3.47959188296875*^9}, { 3.4796543119375*^9, 3.479654314140625*^9}}], Cell[BoxData["1.4998576225279827`"], "Output", CellChangeTimes->{3.4795915969375*^9, 3.4795916370625*^9, 3.479591715609375*^9, 3.47959175075*^9, 3.47959188565625*^9, 3.47965431465625*^9, 3.479678134140625*^9, 3.479678275640625*^9}], Cell[BoxData["0.000148796268381272`"], "Output", CellChangeTimes->{3.4795915969375*^9, 3.4795916370625*^9, 3.479591715609375*^9, 3.47959175075*^9, 3.47959188565625*^9, 3.47965431465625*^9, 3.479678134140625*^9, 3.47967827565625*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "we", " ", "now", " ", "use", " ", "the", " ", "alternate", " ", "formulas", " ", "for", " ", "a", " ", "and", " ", "b"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"x2", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x1", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x1y1", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y1", " ", "=", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"i", ",", "1", ",", "8"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"8", " ", "x1y1"}], " ", "-", " ", RowBox[{"x1", " ", "y1"}]}], ")"}], " ", "/", " ", RowBox[{"(", RowBox[{ RowBox[{"8", " ", "x2"}], " ", "-", " ", RowBox[{"x1", " ", "x1"}]}], ")"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x1y1"}], " ", "-", " ", RowBox[{"x2", " ", "y1"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"x1", " ", "x1"}], " ", "-", " ", RowBox[{"x2", " ", "8"}]}], ")"}]}]}]}]], "Input", CellChangeTimes->{{3.47967813609375*^9, 3.479678268765625*^9}, { 3.47967831009375*^9, 3.479678350609375*^9}, {3.479678434109375*^9, 3.4796784375*^9}}], Cell[BoxData["1.4998576225279823`"], "Output", CellChangeTimes->{{3.479678269859375*^9, 3.479678277046875*^9}, 3.479678313015625*^9, 3.4796783513125*^9, 3.479678438203125*^9}], Cell[BoxData["0.00014879626838127196`"], "Output", CellChangeTimes->{{3.479678269859375*^9, 3.479678277046875*^9}, 3.479678313015625*^9, 3.4796783513125*^9, 3.479678438203125*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", "x", "]"}], ";", " ", RowBox[{"Clear", "[", "y", "]"}], ";", " ", RowBox[{"Clear", "[", "z", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"v", " ", "=", " ", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.479593742390625*^9, 3.47959376715625*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Cross", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", " ", RowBox[{"p", "[", "1", "]"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Cross", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", " ", RowBox[{"p", "[", "1", "]"}]}], "]"}], " ", ".", " ", RowBox[{"Cross", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", " ", RowBox[{"p", "[", "1", "]"}]}], "]"}], " "}]}], "Input", CellChangeTimes->{{3.47959381396875*^9, 3.479593825140625*^9}, { 3.479593976265625*^9, 3.47959399684375*^9}, {3.479594060359375*^9, 3.479594061203125*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"12", " ", "z"}], "13"]}], ",", FractionBox[ RowBox[{"5", " ", "z"}], "13"], ",", RowBox[{ FractionBox[ RowBox[{"12", " ", "x"}], "13"], "-", FractionBox[ RowBox[{"5", " ", "y"}], "13"]}]}], "}"}]], "Output", CellChangeTimes->{ 3.47959382575*^9, {3.479593979453125*^9, 3.479593997484375*^9}, 3.479594061796875*^9}], Cell[BoxData[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ FractionBox[ RowBox[{"12", " ", "x"}], "13"], "-", FractionBox[ RowBox[{"5", " ", "y"}], "13"]}], ")"}], "2"], "+", SuperscriptBox["z", "2"]}]], "Output", CellChangeTimes->{ 3.47959382575*^9, {3.479593979453125*^9, 3.479593997484375*^9}, 3.479594061796875*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"D", "[", RowBox[{ RowBox[{"ArcCos", "[", "t", "]"}], ",", "t"}], "]"}]], "Input", CellChangeTimes->{{3.479594447765625*^9, 3.4795944611875*^9}}], Cell[BoxData[ RowBox[{"-", FractionBox["1", SqrtBox[ RowBox[{"1", "-", SuperscriptBox["t", "2"]}]]]}]], "Output", CellChangeTimes->{{3.47959445246875*^9, 3.47959446190625*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "40"}], "+", RowBox[{"78", " ", "x"}]}], ",", RowBox[{ RowBox[{"-", "42"}], "+", RowBox[{"78", " ", "y"}]}], ",", RowBox[{ RowBox[{"-", "24"}], "+", RowBox[{"78", " ", "z"}]}]}], ")"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"40", "/", "78"}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"42", "/", "78"}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"24", "/", "78"}], ")"}], "^", "2"}]}]], "Input", CellChangeTimes->{{3.4795954535625*^9, 3.47959547134375*^9}}], Cell[BoxData[ FractionBox["985", "1521"]], "Output", CellChangeTimes->{3.479595472015625*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"1", " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"42", "/", "40"}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"24", "/", "40"}], ")"}], "^", "2"}]}]], "Input", CellChangeTimes->{{3.4795961959375*^9, 3.47959620446875*^9}}], Cell[BoxData[ FractionBox["197", "80"]], "Output", CellChangeTimes->{3.479596205015625*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"197", " ", "5"}]], "Input", CellChangeTimes->{{3.47959621328125*^9, 3.479596215046875*^9}, 3.479596291421875*^9}], Cell[BoxData["985"], "Output", CellChangeTimes->{3.479596216125*^9, 3.479596291796875*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"80", " ", "5"}]], "Input", CellChangeTimes->{{3.479596294953125*^9, 3.479596295875*^9}}], Cell[BoxData["400"], "Output", CellChangeTimes->{3.479596296234375*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Sqrt", "[", RowBox[{"80.", "/", "197"}], "]"}]], "Input", CellChangeTimes->{{3.47959635115625*^9, 3.47959635496875*^9}}], Cell[BoxData["0.6372529878771661`"], "Output", CellChangeTimes->{3.479596355390625*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(", RowBox[{"400", "/", "985"}], ")"}], " ", RowBox[{"(", RowBox[{"24", "/", "40"}], ")"}]}]], "Input", CellChangeTimes->{{3.479596629390625*^9, 3.479596663578125*^9}}], Cell[BoxData[ FractionBox["48", "197"]], "Output", CellChangeTimes->{{3.47959664296875*^9, 3.479596664078125*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"48", " ", "5"}]], "Input", CellChangeTimes->{{3.479596646875*^9, 3.479596668859375*^9}}], Cell[BoxData["240"], "Output", CellChangeTimes->{{3.479596648140625*^9, 3.479596672671875*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{ RowBox[{"Zipf", "'"}], "s", " ", "Law", " ", RowBox[{"(", RowBox[{"word", " ", "frequencies"}], ")"}]}], "Title"]], "Input", CellChangeTimes->{{3.479601042484375*^9, 3.479601048859375*^9}, { 3.479653187796875*^9, 3.479653191015625*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"citypop", " ", "=", " ", RowBox[{"{", RowBox[{ "8363710", ",", "3833995", ",", "2853114", ",", "2242193", ",", "1567924", ",", "1540351", ",", "1351305", ",", "1279910", ",", "1279329", ",", "948279", ",", "912062", ",", "808976", ",", "807815", ",", "798382", ",", "757688", ",", "754885", ",", "703073", ",", "687456", ",", "669651", ",", "636919", ",", "620535", ",", "613190", ",", "604477", ",", "598707", ",", "598541"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"ranks", " ", "=", " ", RowBox[{"{", RowBox[{ "1", ",", "2", ",", "3", ",", "4", ",", "5", ",", "6", ",", "7", ",", "8", ",", "9", ",", "10", ",", "11", ",", "12", ",", "13", ",", "14", ",", "15", ",", "16", ",", "17", ",", "18", ",", "19", ",", "20", ",", "21", ",", "22", ",", "23", ",", "24", ",", "25"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"both", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "25"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"both", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"both", ",", RowBox[{"{", RowBox[{ RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "both", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"logboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "25"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"logboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"logboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "logboth", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"loglogboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "25"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"loglogboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"loglogboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", RowBox[{"Log", "[", "n", "]"}]}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "loglogboth", "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.479601053234375*^9, 3.479601059359375*^9}, { 3.479601102234375*^9, 3.479601200125*^9}, {3.479601308875*^9, 3.479601325578125*^9}, {3.47960143678125*^9, 3.4796014558125*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJw1zjsLQQEUB/CbiUUGg1d5P0LyfkzuLbcsFjOJYnQlk6i7SOQzGBkZZJZ8 BZ/BIu5swf3fc06d8z+/znL8PaXZNwmC4Py3nv9SD4+OZKxaVbcsNsiCqI/5 TSabdat1t0S24X56lsh2+KuyHfBqmyd74N0yR/bB3QU7AA+1JDkEz+4JcgT2 HWPkGGzds+Nwq81OwiNvlJyC1xZ2Gr5OwuQMfPmEyFnjX4Wdg19KkJw3/A6Q C/DmzC7C0zG7BA9kdhl21dgV8QcTczOi "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960117028125*^9, 3.479601200703125*^9, 3.479601326109375*^9, 3.4796014571875*^9}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztW71vHEUUH9/ufdgXIDIpSISUU0zOwXITxzFQoCQ2GIiECAVKkwIrQrIL SDhSICJBjQQ1QktJID0SSNtQ0EBFswUpIgTl/gNukI7d+dqZ29/um7Vvby3F Try3njfv9+b93nzsm517e+fe7gcf7tzbu70zeGO0c3d37/Yng+07o6TIm2Ns bpcx9s+Apffj5FZexsn/8V56SUvo+4/5R49FARM/o5DFHCViQfIvSu7jcNTi sgUWRFw5+41DNtIqyW+KE0ReIepYl43C+DwvOJFVkz/CStaACTWB3pmK8XnS eOI9q8f9Yc79EQtjhV6z+z0z/qlOzMKRsNeVZVGQ2W/xqy+MBpH69CWJprrR FmlUuHVeQqSEJg3gLobKnTZU8Wqz2qnHx2WTzkAQX7ODXTqQcRjojuWpawol u0aRuZQz3eWT+il/Q40Z6nrKQ1Tfr8Fepwa/sshlA69Ov0rCptzjnWGiKboH GIPfKzArqmI6uZuyJylsVL9Ou+0Z+JlSL6d1GZ56/SwJq2+OddGmtNyaziUz em1R0/lTTn1KlSuHl/Qkly0srZI+3WqiKU0QkkKFET151UiIV9xR1LiIwmzl PsGvT+dMoRXBam0iT5pyVjupZtF0ljPXILw+9esz60/fy4HyJI6zh9gRvQof 0kvRS9hd1Z05iu7O7Vwv50ALuHZLzjXGk5ppSve+ICk/qyFSVnV93X0MiH7N 1upwb5BDrtM9GcQrAl0mU6LKO0bXUWrqcfLqAWXXSmQuTTGCIJ3NzFWUXSuR lTTlBstWWEXspmPZFijzik29a1RXcdx0LNsCZSVe8euym1w0Ob1RuTdjN9PL M0LeBfr9KZctzMDGUW/LUbK7eBz6J9XuDEM/PITu+hFqCyp7rphaNCFvTrms xFSVtcFJ7om1Y8qoR13+RHj9nyj79MH3f46vvPrtN19/9RvLkSAr/X3zvR/S Svuvb33HaMw3//j931/vfPTT+7/8/FdZ9aYpoOQPAUX3Jyl6CCi6T2JaFJVU nyUFG/zDzzYBGN89ahNqrwi5tVHtAZ3rgM1HNkRuu9RPEd4CDEvFl1VVY4PO L7Ztsf7I9npij4fy+iXltbE54KaTy0RRg+uO+zYIxuPJrr0NiH9MYlokl1Rv enRT8lVA0f4kRauAon0S06KopHrTFFDyc4AiLjHzV65yDvCUq4nRLbIonaYZ OWasTL7GP8g3v9Q8KmDyrzsefPk5Y8tfcFLVbig1uRZh3bh6nS2vSdrNN2q+ fjgcliP8+Bljr92SQbDfyVVuTcn7vcpYaM+YwrrIPwq2L+2dTip4h+1DF8EI Qlka0LHGUl+PNFDTGkEUetOjipIvAcZOu+lYjJ2GjC0Bxij0phk5ZmzacrSu nXHTsRg7AxlD6xqF3jQjYvpt59a67LSWZ53WQimTifc8YPhZrTn5Ep9OxwC2 FYkybJyn8b+HFqIVMYSoD8UQCVye2R7LLWzVMrg10By9oFXI6krXxpkne6ib DA4Q/gEciCisFHrTA5GSnwKMXXDTsRi7ABk7BRij0JtmhJIvAsZW3HQsxlYg Y4uAMQq9aUYo+UnA2KqbjsXYKmTsJGCMQm+aETFrFh1uHldKKfqA3Rck/sSp 2gNlP30QCXFqt2ed+gRZGk4Z+yBgCtA8fwjO71FtFS/Putq3g2SOIjYL8FAU 4M89ebTO5TSRPM6DnnLJTceK/iU4DudBWCn0pschJe8BxtALWqBjMbYOGesB xij0phmh5F3A2GU3HYuxy5CxLmCMQm+aEUreAYxtuOlYjG1AxjqAMQq9aUbE lNmBq6P+slO1/BGR/CIvK/jCS7UUEsVDPrWBb2K4pJAobiu6frZIqnWkSgq5 rvk1zswKhl1SRxGeHl4gk79l+thxgilYGGeeOfJ7+aU7Nvc/ddvgQA==\ \>"]], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQBGIQnR1WsOnXA30HBjD4YL9kLofl/EA9KJ/BgfOa+rvg XbpQPofDGdtq5692ML6AQ8Ol2b0X2nWgfBGHqF3mfLvrYHwJh6jWPz0brGF8 GYctpx9efioP4ys4vA07GbAKzldyOKY/+VJkmzaUr+JgEpN2LLMIxldzOLv/ uBajKYyv4fD+8II3LiYwvpaDU9+jWEc9GF/H4UZKimuFMIyv51Cc3ls9SRDG N3A4Zv1rqdZbLSjf0OHSNeWUVQ9hfCOHvYqJL62uwPjGDkU35CMm7YLxTRxa Qtq7rdbC+KYO/0SMbiouh/HNHORvZPGcmw/jmzvk2Ys8ezILxrdwaH0V0XsA zrd0AAAyTWUT "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960117028125*^9, 3.479601200703125*^9, 3.479601326109375*^9, 3.479601457234375*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQBGIQnR1WsOnXA30HBihYMpfDcn6gnsN7y1//nPSe2XNe U38XvEvXgZt5bdWryR/tz9hWO3+104XKf7NvuDS790K7joNx79vvjgd+2kft MufbXafjwKSg9Yt79R/7qNY/PRusdRzCo4zWFSn9t99y+uHlp/I6DsWrre03 LWJweBt2MmAVkA8xn9HhmP7kS5Ft2g5ioVt3b8xmcjCJSTuWWaTt4GNev/2x HrPD2f3HtRhNtR1q513eI/aQ2eH94QVvXEy0HQJ7DLVZWlkcnPoexTrqaTss 331d7JQMq8ONlBTXCmFtB3kPX0uxVawOxem91ZMEtaHuZ3M4Zv1rqdZbLYeq +l8Hs1exOVy6ppyy6qGWQ9vEbKtqeXaHvYqJL62uaDnobGZhM+tjdyi6IR8x aZeWw8THxdVO39kdWkLau63WajlsXJc/RS+aw+GfiNFNxeVaDieO2paX7OJw kL+RxXNuvpbDh80X/9aKcjrk2Ys8ezJLy+GnzqSa5TmcDq2vInoPAPmQ8ON0 AACMCa3Y "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960117028125*^9, 3.479601200703125*^9, 3.479601326109375*^9, 3.47960145728125*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"x", " ", "=", RowBox[{"Log", "[", RowBox[{"1.0", " ", "citypop"}], "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", " ", "=", " ", RowBox[{"Log", "[", RowBox[{"1.0", " ", "ranks"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numdo", " ", "=", " ", "5"}], ";"}], "\[IndentingNewLine]", RowBox[{"xbar", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], "\n", RowBox[{"sigmaxsq", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", "xbar"}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], "\n", RowBox[{"a", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"numdo", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "-", RowBox[{"numdo", " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\n", RowBox[{"b", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "numdo"}], " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "+", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\n"}], "Input", CellChangeTimes->{{3.479601339046875*^9, 3.47960135096875*^9}, { 3.479601468546875*^9, 3.47960148471875*^9}, {3.47960155009375*^9, 3.479601567109375*^9}, {3.479601614953125*^9, 3.479601619453125*^9}, { 3.479601766390625*^9, 3.479601798765625*^9}}], Cell[BoxData["14.970196024883842`"], "Output", CellChangeTimes->{{3.47960155296875*^9, 3.479601568203125*^9}, 3.47960161996875*^9, {3.479601773859375*^9, 3.479601799609375*^9}}], Cell[BoxData["0.3207960077739069`"], "Output", CellChangeTimes->{{3.47960155296875*^9, 3.479601568203125*^9}, 3.47960161996875*^9, {3.479601773859375*^9, 3.479601799609375*^9}}], Cell[BoxData[ RowBox[{"-", "0.9984576743228499`"}]], "Output", CellChangeTimes->{{3.47960155296875*^9, 3.479601568203125*^9}, 3.47960161996875*^9, {3.479601773859375*^9, 3.479601799609375*^9}}], Cell[BoxData["15.904605455719153`"], "Output", CellChangeTimes->{{3.47960155296875*^9, 3.479601568203125*^9}, 3.47960161996875*^9, {3.479601773859375*^9, 3.479601799625*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Log", "[", "135.", "]"}]], "Input", CellChangeTimes->{{3.479601410453125*^9, 3.4796014275*^9}}], Cell[BoxData["4.90527477843843`"], "Output", CellChangeTimes->{{3.479601413984375*^9, 3.47960142796875*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"citypop", " ", "=", " ", RowBox[{"{", RowBox[{ "8085742", ",", "3819951", ",", "2869121", ",", "2009690", ",", "1479339", ",", "1388416", ",", "1266753", ",", "1214725", ",", "1208318", ",", "911402", ",", "898349", ",", "783438", ",", "773781", ",", "751682", ",", "728432", ",", "672011", ",", "645978", ",", "628670", ",", "586941", ",", "585122", ",", "584658", ",", "584113", ",", "581616", ",", "569101", ",", "563384", ",", "557478", ",", "544765", ",", "538544", ",", "523303", ",", "517017", ",", "507658", ",", "475460", ",", "471856", ",", "469032", ",", "461324", ",", "451455", ",", "445335", ",", "442768", ",", "439467", ",", "432376", ",", "423019", ",", "404267", ",", "398844", ",", "387807", ",", "380149", ",", "376815", ",", "373188", ",", "370448", ",", "355007", ",", "354617", ",", "342510", ",", "332361", ",", "332223", ",", "325337", ",", "317647", ",", "317361", ",", "316802", ",", "308973", ",", "290418", ",", "285018", ",", "281514", ",", "280404", ",", "279208", ",", "277911", ",", "271466", ",", "271035", ",", "270951", ",", "266798", ",", "248762", ",", "247610", ",", "241991", ",", "241727", ",", "239097", ",", "236620", ",", "235594", ",", "232838", ",", "229110", ",", "226401", ",", "225090", ",", "219495", ",", "218432", ",", "218027", ",", "217989", ",", "215093", ",", "214852", ",", "212215", ",", "211299", ",", "210834", ",", "206872", ",", "206481", ",", "204525", ",", "200499", ",", "200123", ",", "199336", ",", "199060", ",", "198376", ",", "198364", ",", "197488", ",", "197388", ",", "196790", ",", "196624", ",", "196093", ",", "195601", ",", "195357", ",", "194729", ",", "194455", ",", "194248", ",", "193882", ",", "193464", ",", "193316", ",", "190299", ",", "190117", ",", "187873", ",", "185702", ",", "184053", ",", "181647", ",", "180872", ",", "179894", ",", "179599", ",", "178612", ",", "176365", ",", "175706", ",", "173278", ",", "170561", ",", "167402", ",", "167082", ",", "167029", ",", "164237", ",", "162917", ",", "162742", ",", "162184", ",", "161696", ",", "160368", ",", "158880", ",", "157063", ",", "156178", ",", "154887", ",", "154147", ",", "153938", ",", "153386", ",", "152157", ",", "151903", ",", "151725", ",", "151654", ",", "151640", ",", "150867", ",", "150782", ",", "148927", ",", "147840", ",", "146878", ",", "145757", ",", "145250", ",", "144502", ",", "144413", ",", "144001", ",", "143408", ",", "142914", ",", "142621", ",", "142474", ",", "142454", ",", "142185", ",", "141336", ",", "141114", ",", "139664", ",", "137894", ",", "136671", ",", "136093", ",", "136016", ",", "133834", ",", "132197", ",", "131249", ",", "129270", ",", "128923", ",", "128549", ",", "127759", ",", "127580", ",", "127573", ",", "127005", ",", "126182", ",", "125896", ",", "125740", ",", "124977", ",", "124512", ",", "124387", ",", "124372", ",", "124192", ",", "123570", ",", "123215", ",", "122542", ",", "122008", ",", "121740", ",", "120292", ",", "120107", ",", "119708", ",", "118737", ",", "118379", ",", "117881", ",", "117357", ",", "117115", ",", "116887"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ranks", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", " ", "=", " ", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "200"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"ranks", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"ranks", ",", "n"}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"both", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "200"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"both", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"both", ",", RowBox[{"{", RowBox[{ RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "both", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"logboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "200"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"logboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"logboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "logboth", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"loglogboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "200"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"loglogboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"loglogboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", RowBox[{"Log", "[", "n", "]"}]}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "loglogboth", "]"}]}], "Input", CellChangeTimes->{{3.47960251296875*^9, 3.47960259621875*^9}, { 3.479602671953125*^9, 3.4796026731875*^9}, {3.4796029646875*^9, 3.479602968015625*^9}, {3.479603016265625*^9, 3.479603024515625*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJw10H9olHUcB/Cv2uyyFVctW/6oM9caurbTublN0+9Kbdasa81pul3Xfugl xcR/UvGPr1q0DaplQRYbHSgxNokjRCYt+hKFphanFQ5G8JCRVxjeP+WIqbnP +/05OJ578f6+7/08z4K2rsbO6caY07e+U9dbH7diLF6Hn7nVU75YtIE2dsr+ 0lo6NGWzZWgVHRYPH6qiC8SuvJIuFNc1VdDzxHk9S+mI+EaX+hFxU2YxXYS9 HepicX/2MbpEvGhAvUj8Zq6YLhUHteoy8ewDRXRUfPzEQnqJeF+5eql4758L 6ArxmnH1MvHRi+pKcdtP6irsf61eLp4TV1eLdy5R14gbZqhrxf99GqFXiBNd 6pV4H7mH6SfEHe+qV4lLJh+iV4vDefNpK/bn5sEOTqVoUyeOlmgOT74zl/mT 4pYy2sF2puZPiY9dmMMcDtppswb5+QeZ08cKma8V5zbRDt43/ADzdeKfK2gH j96czfxpPO8p2sG2X/N6PM/Q/czh34/SZr145HIBc7hoPW2eQX+d5nDH+H3M nxWvbKcd/G1c8wZxqFlz+Pj4vcw3iE+fvIc5nOykzXPiwTLN4dh8zZ8Xp2Zp Drf8G2YeQ/992sJXemkH9/fQHh6p1/4L4kzr3ezTNbSDq3vvYh+e2E+bRjxf OW3hxG/57MOxb2gPR1+jzYvi9LU72YdDR2gHh/fSHq4+O4v9JrEZoi28eYB2 cEk/7eHBOu1vFEdqtQ+PXL6DfbjwLO3hvlHaNIu9pS3cvUz7cPJaiH348B7a bBKPvU5b2G6kHZxq0D7z5drfLA6qtA+7udqHE4Xah9N52n8J54324ezV29mH uwPaMx+jzRZx/Xe0hUdGtQ/bE9qHM0Pa34r3N6B9OPKx9uGxWu3DJZXabxGH rs5kH44O0w4efIP28O5S2rTi/gq0D6f/yWMf3n2F9nDqB9rExcm3aQtndmkf zhVrHw5duI39l8WxLtrSnbSDB9tpDyfztZ/A+7g0A47A6V9pC9tzdAI2X9EO dh/QKeadtIerS+kA7gvT5hXk2enchwvP0xbuO0Mn4NhJ2sHZQ3QKnuilPZx8 iw7g6EHatIn9Ad2HD+/QfTi3Xffhicd1Hzb5us/zf0/jPpw9Qwdw+EvatIvT R+gInE3RFg59SCdgt592cHgXnWI/qftwrEP3udeu+x3i+q26D3ur+3CmRvdh l6/7cPQXw3144nPaw9l+OmD/E9p04v6a6QhcspC2cHcenYAT79202IfTe+gU /6+N9nAuSgewf5Q223BeHYEz02gLJ368ASdgN0o72H9Bp+DwQdrDZhsdwNE4 bbbL1bbSEdg00haOLNZ9OFeg+3Dw13Xuw+572sN+hA64F6dNUq7pRjoC+xra wuHrk9zn+T9oB/edolNw5jPa8/xHdADv7KHNq/Z/7ePcZg== "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960259690625*^9, 3.479602674453125*^9, 3.4796029699375*^9, 3.4796030294375*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJw9jQlUzfkCx+8knotMKK9jScUgISpkeXzbJNpu+3br3ttdKiWvQ+ha7lCh UrbmMJgztlRmisjStNtaRkYqW5FtxjL1LKEIb96c//f9zvmd3/mc7/n8PuaK eB+Vnkgkqvnr/u/1u3tT+9M1a4j+Pi8X1K28bl/uNVVgEaKupsybXjZF4P6Q iy3dn/YjG0L/t2O14pjJAhuhPUOek+VGNsHS6jztvhHkUVA+ehtr0Zdshv5d O53i9MgWuHEx4fGABCuBx+HJ1+IxuTHk8XCOV5fsHkeeiJE3k8RzTcmTkPzr 5Imqr8mTsf7Hnz97fpkk8FRE9Gns17+FPA29H3tWXa0iT0faxbEnnpwi2+Bl dbPsl+1kWzSet04bl0W2Q3GnYVZpJnkGGkqNfBz/zzPhlmfdYJhBnoWk+K7c Dh3ZHoE7K5ZUa8mzIXth95U0kTwHX26USLfFkOfCxvviu+NK8jyMk64f1uNP /hfmnm8zrvUiz0emavnJZQvJC3DteJjVw3FkoOlI951eM4F1QOYb7aDno7k7 oOR9jcUyY+4O6Jza52zoAO6OEI+aVhLeh7sj8p+N1Y8WcXeCxxqbSbt7LYXd CWXz9+uhS2CRM4rNtCGDnnF3xjYUBzVf5+6C1DNXp7+o4+4Cz+yZHzJKuS/E g/Tl1wKKuC+EY0ytXdPP3F1x4cTwNe253F3xxqh654ij3BfBPX7NSM8s7ouQ VukjLs/k7obYCZ1/dq3l7oZXac+VmXHcFyPQ+PbbQLJuMe7Kf//pBwX3JSho apNEBHD/izua8cyfuzs+mj+bokfWuaM661yIzSLuHjjpqju3diJ3DzhkXDrc bsrdE70xdo7pJtw90ZF5JjZlOHcvlFpesPI34u6FnOb6h6+HcPdG+bF+0f/p KzC8oac1/2YMWecN387Zye/0Ba70hrf33i6DTxMFXwI7M2n6rNsCQ4Logm/G 32wRWCdBUUfG6uH1AldKcHhGhHpzHX0f9JbYJHdcpO+DS9Pitu2poO8Dp4Lm D+Fl9H0Ql2qSE3SWvi8SoucNNi6k74udLk/7fs6l7wuLVSlN+jn0fZGh73XI fC99P2Rb7Zcc/Y6+H7ZG1Q4pyqbvB6mH4ex95Eo/TM1XoTqTvj8eb9V4WZPh D5/iU4daN9P3x9DgQ8eOptD3x4bEzkfJyfQDMFZ2J7AhiX4A7p3sL2tZQz8A L6ruN7xbQT8A++0PdNfE0A9Ejavdy8Zo+oE4tiNldYOGfiAKt79vNSFXBqLq ddhpfxX9IJg+2HbZg4wglL1vyxNH0g+CcvTrhVcU9IPQoeeaPVtOPxh7h93r rpPRD8bnjXkpURH0gxFsE7RgTDj9YPQUWh7RSumHIH1uzIHmUPohEN10GTya rAtB2g2TFTtC6IfAwCjL4ftg+qForr8o3RBEPxSu+HXwSLIuFIsKagMNvOmH 4l2LpCPLi34Y9hTWH9W60Q/D+4om2y1O9MNgcKp25Pv59MOw69581Rx7+lJM Toj9IJ5JXwq/8Xl+Trb0pdDaj5wVa0NfCsV3+57ttaYfDufduSmmlvTD0dQl u2U7gX44NBNLPg0zox+O81mpG3r/ST8CScmNpjGD6Ucg4tUL33gD+hEwLU4o dyFXRmD1CDeXI/r0ZYi2sNqxUSSwmQzuQyqK3n6ZIPwnw6r0Q5tsPwssk6FQ k7Ult1dgnQwfbv/j+/hugX+U4dHAdZV5bwSulOHpoByLyx0Ct8uQL4vb8OC5 wCI5fMpKUg/+LrCZHE/vSBqmPGZfjrqB57q1j9iX46se0ZijD9iXo3aJpiC1 jX05kvLnFh1vZV+OM5lG+kPJ7XKszE3Pb7rLvgLraq9tryGbKRBh7nF60m32 Fbhw37Cl/Rb7CnTPiZSkNbKvwPrHqeZ1DewroHlityKhnn0FFtd4IfUK+wqY Jsq2ll1iPxIBB9BgcoH9SIzPiDOoqWY/Em1/ztgtqWI/EvsSWk/bVbAfCbtN 9UuXl7EfiV3OS2MjS9mPREP6Hw6tv7Afiatin7u1ZJESr6Q/jG4sYV+JtkyF WdVZ9pUId7cc/uYM+0r47v7kkXeSfSUu5DzdZX6cfSX6tuQ/2ZHDvhI9K7ap 1UfYV+Lh4WX2o8giFV6Wmqxt2se+CvpbPsg3Z7OvwvXcIMNbO9hXod/lhALD NPZVsDg1IOLEFvZVKE47lG27mX0VklKN9gxMZl+F8t7GRt9N7KsxfvP5g+5k MzVWNX+s6PMt+2oYi04fLFzHvhrlVaKZqWvZV6PR0qfokZZ9NUadsHpxaTX7 amSJH/i1JbKvxilP/YJgskiDwsuq+25kMw1c2lvFGSvZ10CSuKqzJ4F9DZKN w1p0/2Zfg/NedW+LlrGvwZia1WmJsexrcMVZ5x2wlH0N8vff6Z6hYj8Kpr+9 mRerZD8Kk6qWDD+nYD8KVn3v7cwIZz8Kk29a14qk7Edhe1bZ7Fsh7EfBYaN6 8fYg9qPQOvSPkm8D2Y+CvrnEa2MA+9H4L/Fra44= "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960259690625*^9, 3.479602674453125*^9, 3.4796029699375*^9, 3.479603029484375*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJwtlgk0VWsbx0+cyHyGfeQminL33uecuIak6fpTNMlYGSJDSaGSlKlBMpQx Q91K3e+mQhQhFRkiyvDRTYYGmqiolCJT6NvfWvZa79rrt9733c+7/8/wPhoe e2w9JVgsVi0z/v/e8KI99NojHbCmnvr9j43KrbTxdfHYpKn2e+MdjZHLdMsW QE7yesjHlG/G7jK0RY/Ugqn5IWP2v5l1Mt5i6Mf3DZvcGzV+HeeekbhGDIm5 wjG5nHFjn6qroWmzxLB30sv11/xlvK3rh6/mdDH25Sw1LkhnYcZg8opdEuKp 70/Dk2r/bll/EZQ3Ft3N95HAOyWZOVneIlgsOnK7S1sSK/dsL0mdL8Khv5+U Kr+RhGp7iMxSdRFs4nRF7Eg2Iv4rpjyVRMi8265cP3s6Dv9zfdLylxBzVq9f rJw9Ha6SzVIz2oRT55fC+M/RwMZKIUKOjFX6ZEshpnrejXeFQkQl+SwJnSON /qpWt7snhRAXsqUME6TRXKwTMz9RiKSufaGmw9Io+sJJLE0QIj93T6r25hlo KiVsTRmurVkeGFAyA2uu6jRx4oToL3w8cUggg5A9g1l9YUKMipMPZvrKwD65 Yl1VqHBKPxm4fTKY5nJAiMZZ/R9TObL49aTEJd5biJMRkndeuchCz7p6KGeb EPPWnm5zypLFfJfD/NGNzP9NNjYY9ctiaXGnoM5KiLQ2udGMhXJI8PTL320u hJMoaXVTkBwe5TiL3s4Xwqf2TvmTYjm0XB55Pj5XiMd3WgefjjDrB0LlP6oJ MaRQfHR8oTxKhms1dwuEqBlQsnPaK48v2pK3N8sKp/wrD5nZf5RskRSi7OGS L1SXPLJ757F3soTIMrlf8+A3BawP1hOmjtNoUzi+Y9RSAWV/npfAIA0Ff72l T8IVUDQ31Em+l4bltQQzzyIFxKPIofUxjQXLtTdUvFdA1K1G3U/1NLifCeV+ ZUVYnjIciyulsWLDuOF0c0W8ifV7tKmAxuW6SwfkAhRh6l1n0HKdhq/pyTLF i4q4f0M5+HUWjVv9J3VmNSpigKhKnnWFxv671suMRhRhsSdY1TKRnopPJcTc s5UpT6CRp3rtdvM6JfiSXz4PHqTxi2Xzyy1ACd9iPm5L2EUzfs3N1jqvBHvB sx/2DOcs09xB3VfCC/f31/72oCHPrwr061VCbkunjesmGgOTJY7h0hzk9rWi dyMNa19hw4/fOPip0btAgmHHc9cvCsQcVCXecdJbTWO9+YOa7uUc5K8Ku3OQ ohFMPxz1tOLAJK7m0mt1GrF6nh1pbhyMexuYxqowevryPh3fy0Ffwi3fSGUa nzrICJ1wDkrp+6KNBA1Vla6epGQOMlob3n7n0lP5x0F5ptTOr9MZPejUlOQC DiRCNbTmMJwxb+1zgyoO7L4sjhhiM/7gUt/OPebA2vrsoMIEhcVNd/958JoD g7kusYueUbg7K/VC8VcOduZq/d7eRkHDxOvE/kkOCvrigpQbKNS1zR/8Kc/F pYWu26PrKRjbDzpbqnIxXqIX0VdN4b2Opv5umouaP3bFn6mgoJYz3ua8iIsV ua1jW8ooOEunn1M142JXlEqGw20KvJ57VTm2XPjvXKYoyKOwTM8iWtGNi2Sz numTWRRoy3mE2S4uNAMjW9gZFPbl552wDuEijm2VrnGWwsqwoAndaC5Oic7b XDlNTdUfLk7sqOMWnKLgeSFDIvQfLlzWcxanMRze7RHffY0L7WxPVCVQaPEc kqGKueg+4WWlw7DspZEt5jVc2BYVpndEUzAv6/XGYy54jumZVyIpNBsns2d2 cnHkwJeuiAgKbntdftb3cDHP7bl9UwiFoZjcmS6DXLzMn+HWFkyBI9ZVbfrF xafKV01DARTKnbUq1eV4OG90YaTWm0LAjsxyK2UealcZ9DfvpODXNlDlrsFD ZlJkUJMXhbMVmYEbxTzknRzuUGG4q3QiTLSIh8rvzjc3elIYuSx9uNuEB/U3 8Q/W/5/XvZtx1IKHsuHOqzJbKXgHjNyaZs/DNrXv5g89KIjUlH/zcuehT2LV qcXu1FR95uEs/+VIvRsF3RJV8u1+HibDr0bucKWYehD9auQID456DsZztlD4 zHf0HDrBw2gefTnUhcLbNMOFL1J4iF3qfaF1M4UrS/RrMy/wwGo3U1RjOMWz P8U5k4eYJyoBSU4UNq+qkB27wYMCkWhyzpFCtfW9tLASHlobql2OOFCQ2jit sP8+D6vwX0VVhnXNzxWvbeRhdW6dvYI1hYuXawWJbTwMtdn0JVpRiC+UEZW/ 4uFMXsOV0DVMfB3Td3raw8NwRYv+8RUU/noTN+3lN8ZeYZ3q8J8UzJZ3H3w0 xkPKyz89lxhRiI2/onVdkg+xv++YjCGF3FVfbQPl+djw+9UNK/Qp3PbLcNQW 8BFqpLrIV4+CXduJuGY1PjxOp/We1WHyIXjlsm2/87EyNStSnaam7jc+Wgbd nuqTFNYVfltrs4gPL6pkgj+XQu9pf5PrxnwUJ0YdGZ/JxPd86e7hVXyERDSr eytS+EB/D9a15sP12ye7PQoUHF2rVjs58KFe5F9uxrCptluqvxsfQbPWmF1m U8j/Pus/ITv42KkpSgpnUWi4VXTV348PC25FwY9fJAr99g04BfERGJt+TH+S hEetZol+GB95XonHs8ZJdHvTSyai+Rh7Jn1uzwgJcX9C7u1EPrrkDt27OkDC x0HaZutffPTIZ2g+6CNxSmvCkfUfPrLddh1585FEoCCNlZjBh21ZSdTF9yQu hLWGKOUy65/bNC3oJjHqvm/iaBEf9XJ3RkK7SPDkZpa+L+Vj2ihrzpU3JLYc N+9dXs1H3Tqv3KhOEurF+25FNzB6ZC8tyOkgIa3O3lDdzMetBILNY7jJ3eHb wDM+9mfFZre8IKf6Bz4O1T06Wcuws8fbLnEPo5/G+pvCZyTmB05LMfzKx/1X nLbXT0l82VX9Q3+Ij5ElW21imkk4fX9hqDXBx+HuKI36JhILPqiGybIJeL0z CPBvIHHO8uhYlyyBtbVWiHpIwo5l1JTPJaB+wO1EWQ0Ja51/NANUCGy6gCaV +yRMK+3URXMI/B63S6G2ikRIsOhDmxaBzs8LU20qSYyNrEgPFBNI8++4aVBB wuuURICCPgGDYw0+fmXM+WfMCjuzmEDKSh/fraUk8kcSelVAoCn2g0nHXRIq rqXNCeYEGmVsX9QxHP7Q3mXMgsA3l7/VmktIJPjtz95sx9hL8JhbeZvE6rVx zwscCWyxoJUHbpEIbQ1WYLkRsEudWH81n8Sl7kiPFdsJ3M/oSdHIIVFkvH3k oC+B6W3Z75IySNwTmfZd8ycwGhC/fftlEq7tRVueBBF4e2m30WyGeddbd/Yf JtBfqnKwJY2c6t8IsI+PuUefYux/bC9XjCXwOMuB8zSJxOwXDuAkEZB64J/L iSERwa6rn/EXAc1CWdcbx0lEh1QHjZwnUBSTfko/mgQneu+2V+kEQqKIM3IR JP5wOZZZlkWgfLy52e4YE293lm1OyWX0ji6+aMHwUm//M243CQS2/qyQPErC oDk+UKuEgIB182LeIRIOqxM5byuY/ZUsw6iDJCz3cUJP1xBopm0LukJJiBLN m0wbCMy+IfpUE0Ri/BJb48O/BBJl3mzoPEDi32tlJ461ESi0ZOc6MlyQF6Wu 0kEg74HnqzUMix6q/Lj8hoDZ6w6ZuP0kctWeaAo/ELA5EPhl1J/E1helpZmf CUQInNvC9pJol4puV/tOoNiq/kfBbhIvHRKi44YJzKkNijngS0KD3fl6YJzA w5Vh1pt8SJT4iYftJATIPv98ZKEnCcPP7zuvSQug/u/AMt9tTH5b0lmT8gII K9cp3/Fg5lOqd6zhCSCa/jI5bgsJI29jvfiZAojbdepYLuRUfy3AycSyxU+d SFTOi6BYmgKYhG9fe9KBxM5tP4J0SAE6eB9KjtqTiGr86zcHsQBsDRur8E1M PpxPVQ3RFeB/InVLYA== "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.47960259690625*^9, 3.479602674453125*^9, 3.4796029699375*^9, 3.479603029515625*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"x", " ", "=", RowBox[{"Log", "[", RowBox[{"1.0", " ", "citypop"}], "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", " ", "=", " ", RowBox[{"Log", "[", RowBox[{"1.0", " ", "ranks"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numdo", " ", "=", " ", "200"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xbar", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], ";"}], "\n", RowBox[{ RowBox[{"sigmaxsq", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", "xbar"}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], ";"}], "\n", RowBox[{"a", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"numdo", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "-", RowBox[{"numdo", " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\n", RowBox[{"b", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "numdo"}], " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "+", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\n"}], "Input", CellChangeTimes->{ 3.4796026895625*^9, {3.47960297784375*^9, 3.4796029935*^9}, { 3.479603038546875*^9, 3.479603061625*^9}}], Cell[BoxData[ RowBox[{"-", "1.3507630916011883`"}]], "Output", CellChangeTimes->{ 3.479602690265625*^9, {3.479602979203125*^9, 3.479602994234375*^9}, { 3.47960303296875*^9, 3.47960306215625*^9}}], Cell[BoxData["21.074884255436405`"], "Output", CellChangeTimes->{ 3.479602690265625*^9, {3.479602979203125*^9, 3.479602994234375*^9}, { 3.47960303296875*^9, 3.47960306215625*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{"Word", " ", "Frequencies"}], "Title"]], "Input", CellChangeTimes->{ 3.4796030665*^9, {3.479604878203125*^9, 3.479604882890625*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Wordfreq", " ", "=", " ", RowBox[{"{", RowBox[{ "61847", ",", "29391", ",", "26817", ",", "21626", ",", "18214", ",", "16284", ",", "10875", ",", "9982", ",", "9343", ",", "9236", ",", "8875", ",", "8412", ",", "7308", ",", "6954", ",", "6810", ",", "6644", ",", "6575", ",", "6475", ",", "5096", ",", "4790", ",", "4735", ",", "4707", ",", "4626", ",", "4623", ",", "4599", ",", "4577", ",", "4452", ",", "4332", ",", "4285", ",", "4134", ",", "3801", ",", "3792", ",", "3719", ",", "3707", ",", "3578", ",", "3490", ",", "3430", ",", "3328", ",", "3227", ",", "3006", ",", "2802", ",", "2686", ",", "2608", ",", "2593", ",", "2551", ",", "2532", ",", "2493", ",", "2470", ",", "2436", ",", "2369", ",", "2354", ",", "2183", ",", "2087", ",", "2055", ",", "1962", ",", "1893", ",", "1795", ",", "1774", ",", "1733", ",", "1712", ",", "1712", ",", "1683", ",", "1649", ",", "1634", ",", "1632", ",", "1595", ",", "1561", ",", "1542", ",", "1542", ",", "1525", ",", "1524", ",", "1434", ",", "1383", ",", "1382", ",", "1364", ",", "1343", ",", "1336", ",", "1298", ",", "1277", ",", "1275", ",", "1254", ",", "1248", ",", "1241", ",", "1233", ",", "1220", ",", "1193", ",", "1186", ",", "1165", ",", "1145", ",", "1135", ",", "1119", ",", "1112", ",", "1085", ",", "1064", ",", "1033", ",", "1016", ",", "995", ",", "958", ",", "953", ",", "950", ",", "943", ",", "932", ",", "927", ",", "916", ",", "903", ",", "902", ",", "902", ",", "896", ",", "891", ",", "888", ",", "881", ",", "862"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Length", "[", "Wordfreq", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"citypop", " ", "=", " ", "Wordfreq"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ranks", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", " ", "=", " ", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", "112"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"ranks", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"ranks", ",", "n"}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"both", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "112"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"both", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"both", ",", RowBox[{"{", RowBox[{ RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "both", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"logboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "112"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"logboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"logboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", "n"}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "logboth", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"loglogboth", " ", "=", " ", RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"n", "=", "1"}], ",", " ", RowBox[{"n", "\[LessEqual]", " ", "112"}], ",", " ", RowBox[{"n", "++"}], ",", " ", RowBox[{"loglogboth", " ", "=", " ", RowBox[{"AppendTo", "[", RowBox[{"loglogboth", ",", RowBox[{"{", RowBox[{ RowBox[{"Log", "[", RowBox[{"citypop", "[", RowBox[{"[", "n", "]"}], "]"}], "]"}], ",", RowBox[{"Log", "[", "n", "]"}]}], "}"}]}], "]"}]}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"ListPlot", "[", "loglogboth", "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.47960350359375*^9, 3.479603558328125*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxdz01I03Ecx/EfnSRiSXjYIcJCJCRimZk92H7ZKstlZllZHcQWdKqOEUI/ Oo0yerIys1hl9EymED0Y/HoCA+kBCiI6GHjoMGSHUR081Py8fwcbbO+9/h++ jM1tP9i8b5ox5tC/d6GF12jluNW3XLzw6Qe/YzNZG/uGiybr5nzBxfK7T7hE ff8BR7XXvsWz1aOvcan2+Cs8T50eXKa95SUuVw94PF/tHMIVasVzvEDteYYX qvmnOKY+DF6kNgRXqj8H8WL18QCu4v8EL1H3Blers4KXqjODa9Tfj/Ay9Ufw cjUVvELN9uOV6ovgWrU8eJU68ADH1UywVSPY4V/32Ver+XvsOIFNnTpylx0b bNao1XfYccct9oQ6dpMdT/Sxr1XbsMN1YV+njt9gx8PYrFePhB0nw16vFoUd H7vOvkFNYYcj19g3qo0ZdhzBpkG9fZUdd11hT6ozetnxk8vsm9QS7LDvYW+c aoeTYd+svrnEjtPYNPEcW5zt5jlOYY+j2GyZaos/X+Qef8Qe77/AfbPae557 3I0dbsIeT3Rxv1UdxRbXY4fz57jHWWy28XvYYocdToV7nAj3LWpxuMd9Z7nH aexxFTbb1bEz3OMh7HAH9rgNmx38v9Pc4wx2uAx77E9xv1ONYotzJ7nHw9jj fmxa1cPY4tJwj/90co8z2OxSa7DFMez+s8cm3O9Wv57gHnvscBp7nDvO/R77 Fx6B+BU= "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{{3.47960353596875*^9, 3.479603559078125*^9}, 3.479604086515625*^9, 3.47960493525*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxd0H8s3Hccx/ErlrWrqvqxmVI/g6kYehxW7uXn+XnO3TnHYdphaBVR1aDt maY3P7p1VbSIMcIsqn5MV8uSM6tujo6sHV309MdM3XqWVUekP3S2fN9Lum/y yTePvD7Pfz42+3OF6TosFitv4/zzZ389w75g6ADWv9+f3CN2kgc9NfaMWVC0 md78No28GaNVa5sbV+0YG2LqXiWnr5lsAqmmerc0nWyG2+xAbnudLWMLfNPB rx/NJltDrtbsWwgj26J/crGFFUS2h1S+9as6V7IDas0+fOFoTHbCDWGwwmrc hrEzfuqIOTLQR3ZB6fqQ6lkH2RW+prPsxAayG7ZJq9WbasnuUKQKtr91muwB xbxjaKUReQ/05ctnf162ZsxGSM176odasie+U6PkvobshWyXEf6Lu2QO/jJN 4X3xn72xqV3HgDVH9oFShOqqWbIvTEd0rVKuk9+BZwCnP1VJ3gv1ufTbqiGy H7he5u0HO8n+CNzD1c8pJXNR/6pW5FVCBowWtxUE5jOWAwtV49MNubQH4HB5 VubRZNoDsFJ2/sAVAe2BcOAgb4ZHe+DG+4RX7vemPQj6EolpoxPtQdgxtXIs TI/2YNgaj4sbhq2YPRiHvKtV/FrGrBAULK9W2xTRHoJE68kin3zaQ6EoqWiT p9IeCp+A1+ZuSWnn4Yw48tLjMNp5CFcVrhZyaQ/DfJo6N9id9g0vBaQNmtAe Dq5zSijbgPZwxJ+ez3syvIvZI/Bo6DfXpbOM5RHoti+tNy+nPRJjlhxBj4j2 SHywvOr8ytu0R8H3dz2nujVLZo9CT0u7Z4WWMSsa8v46Rf1N2qMhSttbufMH 2vkvW87HYpLxk6bLtMdgsOq+5nIr7TFQv67NkZ2nXYDBeb13zesZQwCd9VZl 8wm6L8BEby3/oyzGwwLsMPtlXTeZ+tiXjVh0emx50BVLfSzKno8tPxdQH4uG z7kGXUbUCxEqGbc8t2TB9EJc5e+ST2sZy4WYVF5TZt1hPCyE26Pu7FOTjFki mF0UHXdRUS9Ckn6GZ0s39SIM5j5O1jZSL8KbvE9LYhuoF8PvWJQ6poJ6MWwO tNRknKRejLZcK5P249SLYavjckXvKPVxmE6ecHLOoT4O5cVlBVwJ9XFYL8q8 5i2gPg7dF7xPLvlRL4Fh9x8WHi7US7B2WFPnZUe9BE8VJ+50m1AvwcTdL2e5 BtTH48zHU9EmKzuZPh4Loy4/Gt5jLI8HJ27y4tgI4+F4PNNdHxjvY8yS4oZb TX1mE/VSqJo0nQ2F1EuRnN71q38e9VIcev9Wil4O9QmwcuttzU+nPgGzV5+O 5iRQn4D82kufPRRRnwD/mb6Ufh71idgXHRi33Yf6RGz1yZjTcqj/n4cTwavo jHjDnXoZepsPSpJ2Uy+DrDjDc9GRehlUyvVPrltRL8P3QQNexVuoT8LfA73O /g== "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{{3.47960353596875*^9, 3.479603559078125*^9}, 3.479604086515625*^9, 3.479604935296875*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxFlAs0lekax91KbHvb3+dS+b5tf3trxnUUuU9m/6nQTUok4ajTLrmEQTpU LnVGbo27iWqIoSKFo1wqUuSSVBp1jqRiDMWMkEEu46xlrd61nvWu33rWetfz /J/n/Qv2++8US0lISAQsxP9vo+oXRue4X0Ni8RzVcP69JHUV/jSfnrPW7xfF 5Kk8v3dgFVjS18Lep34UNcRPLsue0FjMT4ievI0zLb2ogbWJw39Z1U2JXAYT dF3EGpBidKZZRTOiV0bWovwMIXa7GpZ8L5wX3S6wz2zwFiKo6FtR2SUJRHYP 7uu3Ey6+L4my9oEcifVCqDpV1JT6SMElknUrQ1+IraYRt3r1pZG+4sy8ppIQ Jy523FZ9K42OnRti+K0C7Egw0JX5twyeFWw/Wl4qQGHNC9UWegmOz1W1zBQI wLfbZq56dQksVLqMXLMEi/UvBdsloVsyXYCwiOl7PleXIsbTQVE7UYAfkn0s wvmyiOnTtIkjBdArl1lqclYWCpGjKb+OMkjuDQq3/ksWG1P/2f1hiEFpiX+a /t5luN+N8HeDDJoaLEODq5fBW6/efv4Ng5Hyp7MnVOQwruJhe3WBp/RSjhf6 ykEyX4oj8ZpZ1E8OtY5IiO9i0KY28j6NKw+Vemm+RxuDpNPSlT3u8jC2Mi3z rGWgsTmj0/WyPLrTxK9aqhjw59pazUbkITJRy/ctZJDdyZoqMGbBeq1Iwe84 A1fdZLvHx1jIlB1yNAln4NNUebejigVygB1kHcjgaeWv4y8nWeiPb+3M8mcw wa6KmjFWQPCpw17H3Bk0jCk6ugYq4FPUTz6VDszifBXwtSkCXtgyuPPQ4g+t XoUFfTbF7TdjcNnqfkPjSjYUnJ1VsrUYdLLPeE3Zs0E8+XTCToYB+3vDbzui 2RAqte7KquPDvvjsRnEFG0fMElrs0/n4xlJ/V20/G0GjEwmCUD6IIWXVEVUO XJn2UPNAPtbvmjFZYsNBTHhsXqQnH/nNeUdZwRyYW8m/funCh6910h1OLgdJ u7ZcH7Pj4+ZI0mq1Ng42tYRMhIj4CKlxWGc2yUHfgW7/DQb8xf1URN+w1YEK ZT6uU8W3nm1RhEjHw8aIw8e8xI55z2BF7E7sC5iuU8emKyVXvzqviI9Vv+kP p6ijaJ3QS+u+IopXHc9UO6UOBaX60IBBRTTzTB1KHNUxNle9J1qWi+jRCZ0l q9Xh4KvT+mklFxbvZbQyJnnYk3UtV0WPi5KcfOPYIR622TQ29FlyEVmWEZP5 nId/aT+cEm/nwvHAujiqiYd4Q/GrbM8vzPYlP5wJ5GLATWn6/E0ePrzSPL06 mouK+HeDN3N5oFb0DiSncNGtOuS39yfe4v9byPfJ/EMtk4cQ7bTUlDIupOZy ay9G8FCgsfl/RvVcPLqRbn/2MA/fEFofs55yQaz475y0Ow/mj2tyGt984Rq1 tAtVf3JRaCj3+5UdPAisDsWGzHERNds8OuvAQ3PnqvHPCgSyLos4V0geRLvH 3ewpAjbOrby0YRr9q4Vrj2gTeGCvHtk5RINXNNPpZkqgvbax9nAPDTfZS1nU RgJrPhZ7/9BOL+xtXX3RTgIrrjme1Guhsc5wawzHk4CbwkHjnGIa2vYayhv9 CFT4j7kPZdMIKr0e6xBGYKXtz+E7smhsiDw2axBDwPLE1u7tsfSi/xAQ+OSk HjxNQ3yhQCo8h0CeP185/ySN6L79iX3FBIRSepUyx2g8F0/IaVUR6HR/pKXj R0M+b9LDpoHAqbCoIJEzDZs7g954SmAu1KvRzIHGM1GKzPJuAsXnzE4PW9Lw DHT/3DJAgFv8B22oR2MirmS5+ziByeDBDBMNGlw9A+rxPIHPMRE9xco07rp9 dU+dReLRm/90iTg0gr0K725XJZH045Ntyp8oBHSO1e8TkOhv0HvMfUvhXG1h qJMeCVOn9mvN9RR6b89G6pqSmJGeK28tpTCZL3uyz4pEx5rUTK/zC7zlt2VR W0m0nB8szAqh4B08eVNyNwl38ZXe7wIo6PJUVx7aR+LIoZceMn7Uoj+T4K+5 kRsopmBQTWm+CyHR9eBzg98easEPYnomI0gEpl+/9MGRwpDSHvFELInvXpR6 lNlSeJdtYtyVSmLfNmsnRXMKv1isbSq8QIJlfvD1kCmFVPFIqlvhF95rWys/ fYOEbWzh5uUGFB441GVHVpO4cdHX2U2XwlInyfKR+yT2hh00HtBcqMcmq2pz 20I/tXPJbXwKuflNKj92kni4vtwkTI5CYrmc7t0eEn8DszntBg== "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{{3.47960353596875*^9, 3.479603559078125*^9}, 3.479604086515625*^9, 3.47960493534375*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"x", " ", "=", RowBox[{"Log", "[", RowBox[{"1.0", " ", "citypop"}], "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", " ", "=", " ", RowBox[{"Log", "[", RowBox[{"1.0", " ", "ranks"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Length", "[", "x", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"numdo", " ", "=", " ", "100"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xbar", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], ";"}], "\n", RowBox[{ RowBox[{"sigmaxsq", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", "xbar"}], ")"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], "/", "numdo"}]}], ";"}], "\n", RowBox[{"a", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"numdo", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "-", RowBox[{"numdo", " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}], "\n", RowBox[{"b", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "numdo"}], " ", "xbar", " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], " ", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}], "+", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}], " ", RowBox[{"Sum", "[", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "numdo"}], "}"}]}], "]"}]}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"numdo", "^", "2"}], " ", "sigmaxsq"}], ")"}]}]}]}], "Input", CellChangeTimes->{{3.479603583234375*^9, 3.479603615*^9}, { 3.479604917953125*^9, 3.47960495140625*^9}}], Cell[BoxData["112"], "Output", CellChangeTimes->{{3.47960358453125*^9, 3.479603616078125*^9}, { 3.479604922328125*^9, 3.479604952125*^9}}], Cell[BoxData[ RowBox[{"-", "1.056013932612408`"}]], "Output", CellChangeTimes->{{3.47960358453125*^9, 3.479603616078125*^9}, { 3.479604922328125*^9, 3.479604952125*^9}}], Cell[BoxData["12.000789924272064`"], "Output", CellChangeTimes->{{3.47960358453125*^9, 3.479603616078125*^9}, { 3.479604922328125*^9, 3.479604952125*^9}}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", "bestfitword", "]"}], ";", " ", RowBox[{"Clear", "[", "loglogboth", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"bestfitword", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"11.03241887", ",", "0.350455347"}], "}"}], ",", RowBox[{"{", RowBox[{"11", ",", "0.38469"}], "}"}], ",", RowBox[{"{", RowBox[{"10.99", ",", "0.3952501"}], "}"}], ",", RowBox[{"{", RowBox[{"10.98", ",", "0.4058102"}], "}"}], ",", RowBox[{"{", RowBox[{"10.97", ",", "0.4163703"}], "}"}], ",", RowBox[{"{", RowBox[{"10.96", ",", "0.4269304"}], "}"}], ",", RowBox[{"{", RowBox[{"10.95", ",", "0.4374905"}], "}"}], ",", RowBox[{"{", RowBox[{"10.94", ",", "0.4480506"}], "}"}], ",", RowBox[{"{", RowBox[{"10.93", ",", "0.4586107"}], "}"}], ",", RowBox[{"{", RowBox[{"10.92", ",", "0.4691708"}], "}"}], ",", RowBox[{"{", RowBox[{"10.91", ",", "0.4797309"}], "}"}], ",", RowBox[{"{", RowBox[{"10.9", ",", "0.490291"}], "}"}], ",", RowBox[{"{", RowBox[{"10.89", ",", "0.5008511"}], "}"}], ",", RowBox[{"{", RowBox[{"10.88", ",", "0.5114112"}], "}"}], ",", RowBox[{"{", RowBox[{"10.87", ",", "0.5219713"}], "}"}], ",", RowBox[{"{", RowBox[{"10.86", ",", "0.5325314"}], "}"}], ",", RowBox[{"{", RowBox[{"10.85", ",", "0.5430915"}], "}"}], ",", RowBox[{"{", RowBox[{"10.84", ",", "0.5536516"}], "}"}], ",", RowBox[{"{", RowBox[{"10.83", ",", "0.5642117"}], "}"}], ",", RowBox[{"{", RowBox[{"10.82", ",", "0.5747718"}], "}"}], ",", RowBox[{"{", RowBox[{"10.81", ",", "0.5853319"}], "}"}], ",", RowBox[{"{", RowBox[{"10.8", ",", "0.595892"}], "}"}], ",", RowBox[{"{", RowBox[{"10.79", ",", "0.6064521"}], "}"}], ",", RowBox[{"{", RowBox[{"10.78", ",", "0.6170122"}], "}"}], ",", RowBox[{"{", RowBox[{"10.77", ",", "0.6275723"}], "}"}], ",", RowBox[{"{", RowBox[{"10.76", ",", "0.6381324"}], "}"}], ",", RowBox[{"{", RowBox[{"10.75", ",", "0.6486925"}], "}"}], ",", RowBox[{"{", RowBox[{"10.74", ",", "0.6592526"}], "}"}], ",", RowBox[{"{", RowBox[{"10.73", ",", "0.6698127"}], "}"}], ",", RowBox[{"{", RowBox[{"10.72", ",", "0.6803728"}], "}"}], ",", RowBox[{"{", RowBox[{"10.71", ",", "0.6909329"}], "}"}], ",", RowBox[{"{", RowBox[{"10.7", ",", "0.701493"}], "}"}], ",", RowBox[{"{", RowBox[{"10.69", ",", "0.7120531"}], "}"}], ",", RowBox[{"{", RowBox[{"10.68", ",", "0.7226132"}], "}"}], ",", RowBox[{"{", RowBox[{"10.67", ",", "0.7331733"}], "}"}], ",", RowBox[{"{", RowBox[{"10.66", ",", "0.7437334"}], "}"}], ",", RowBox[{"{", RowBox[{"10.65", ",", "0.7542935"}], "}"}], ",", RowBox[{"{", RowBox[{"10.64", ",", "0.7648536"}], "}"}], ",", RowBox[{"{", RowBox[{"10.63", ",", "0.7754137"}], "}"}], ",", RowBox[{"{", RowBox[{"10.62", ",", "0.7859738"}], "}"}], ",", RowBox[{"{", RowBox[{"10.61", ",", "0.7965339"}], "}"}], ",", RowBox[{"{", RowBox[{"10.6", ",", "0.807094"}], "}"}], ",", RowBox[{"{", RowBox[{"10.59", ",", "0.8176541"}], "}"}], ",", RowBox[{"{", RowBox[{"10.58", ",", "0.8282142"}], "}"}], ",", RowBox[{"{", RowBox[{"10.57", ",", "0.8387743"}], "}"}], ",", RowBox[{"{", RowBox[{"10.56", ",", "0.8493344"}], "}"}], ",", RowBox[{"{", RowBox[{"10.55", ",", "0.8598945"}], "}"}], ",", RowBox[{"{", RowBox[{"10.54", ",", "0.8704546"}], "}"}], ",", RowBox[{"{", RowBox[{"10.53", ",", "0.8810147"}], "}"}], ",", RowBox[{"{", RowBox[{"10.52", ",", "0.8915748"}], "}"}], ",", RowBox[{"{", RowBox[{"10.51", ",", "0.9021349"}], "}"}], ",", RowBox[{"{", RowBox[{"10.5", ",", "0.912695"}], "}"}], ",", RowBox[{"{", RowBox[{"10.49", ",", "0.9232551"}], "}"}], ",", RowBox[{"{", RowBox[{"10.48", ",", "0.9338152"}], "}"}], ",", RowBox[{"{", RowBox[{"10.47", ",", "0.9443753"}], "}"}], ",", RowBox[{"{", RowBox[{"10.46", ",", "0.9549354"}], "}"}], ",", RowBox[{"{", RowBox[{"10.45", ",", "0.9654955"}], "}"}], ",", RowBox[{"{", RowBox[{"10.44", ",", "0.9760556"}], "}"}], ",", RowBox[{"{", RowBox[{"10.43", ",", "0.9866157"}], "}"}], ",", RowBox[{"{", RowBox[{"10.42", ",", "0.9971758"}], "}"}], ",", RowBox[{"{", RowBox[{"10.41", ",", "1.0077359"}], "}"}], ",", RowBox[{"{", RowBox[{"10.4", ",", "1.018296"}], "}"}], ",", RowBox[{"{", RowBox[{"10.39", ",", "1.0288561"}], "}"}], ",", RowBox[{"{", RowBox[{"10.38", ",", "1.0394162"}], "}"}], ",", RowBox[{"{", RowBox[{"10.37", ",", "1.0499763"}], "}"}], ",", RowBox[{"{", RowBox[{"10.36", ",", "1.0605364"}], "}"}], ",", RowBox[{"{", RowBox[{"10.35", ",", "1.0710965"}], "}"}], ",", RowBox[{"{", RowBox[{"10.34", ",", "1.0816566"}], "}"}], ",", RowBox[{"{", RowBox[{"10.33", ",", "1.0922167"}], "}"}], ",", RowBox[{"{", RowBox[{"10.32", ",", "1.1027768"}], "}"}], ",", RowBox[{"{", RowBox[{"10.31", ",", "1.1133369"}], "}"}], ",", RowBox[{"{", RowBox[{"10.3", ",", "1.123897"}], "}"}], ",", RowBox[{"{", RowBox[{"10.29", ",", "1.1344571"}], "}"}], ",", RowBox[{"{", RowBox[{"10.28", ",", "1.1450172"}], "}"}], ",", RowBox[{"{", RowBox[{"10.27", ",", "1.1555773"}], "}"}], ",", RowBox[{"{", RowBox[{"10.26", ",", "1.1661374"}], "}"}], ",", RowBox[{"{", RowBox[{"10.25", ",", "1.1766975"}], "}"}], ",", RowBox[{"{", RowBox[{"10.24", ",", "1.1872576"}], "}"}], ",", RowBox[{"{", RowBox[{"10.23", ",", "1.1978177"}], "}"}], ",", RowBox[{"{", RowBox[{"10.22", ",", "1.2083778"}], "}"}], ",", RowBox[{"{", RowBox[{"10.21", ",", "1.2189379"}], "}"}], ",", RowBox[{"{", RowBox[{"10.2", ",", "1.229498"}], "}"}], ",", RowBox[{"{", RowBox[{"10.19", ",", "1.2400581"}], "}"}], ",", RowBox[{"{", RowBox[{"10.18", ",", "1.2506182"}], "}"}], ",", RowBox[{"{", RowBox[{"10.17", ",", "1.2611783"}], "}"}], ",", RowBox[{"{", RowBox[{"10.16", ",", "1.2717384"}], "}"}], ",", RowBox[{"{", RowBox[{"10.15", ",", "1.2822985"}], "}"}], ",", RowBox[{"{", RowBox[{"10.14", ",", "1.2928586"}], "}"}], ",", RowBox[{"{", RowBox[{"10.13", ",", "1.3034187"}], "}"}], ",", RowBox[{"{", RowBox[{"10.12", ",", "1.3139788"}], "}"}], ",", RowBox[{"{", RowBox[{"10.11", ",", "1.3245389"}], "}"}], ",", RowBox[{"{", RowBox[{"10.1", ",", "1.335099"}], "}"}], ",", RowBox[{"{", RowBox[{"10.09", ",", "1.3456591"}], "}"}], ",", RowBox[{"{", RowBox[{"10.08", ",", "1.3562192"}], "}"}], ",", RowBox[{"{", RowBox[{"10.07", ",", "1.3667793"}], "}"}], ",", RowBox[{"{", RowBox[{"10.06", ",", "1.3773394"}], "}"}], ",", RowBox[{"{", RowBox[{"10.05", ",", "1.3878995"}], "}"}], ",", RowBox[{"{", RowBox[{"10.04", ",", "1.3984596"}], "}"}], ",", RowBox[{"{", RowBox[{"10.03", ",", "1.4090197"}], "}"}], ",", RowBox[{"{", RowBox[{"10.02", ",", "1.4195798"}], "}"}], ",", RowBox[{"{", RowBox[{"10.01", ",", "1.4301399"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "1.4407"}], "}"}], ",", RowBox[{"{", RowBox[{"9.99", ",", "1.4512601"}], "}"}], ",", RowBox[{"{", RowBox[{"9.98", ",", "1.4618202"}], "}"}], ",", RowBox[{"{", RowBox[{"9.97", ",", "1.4723803"}], "}"}], ",", RowBox[{"{", RowBox[{"9.96", ",", "1.4829404"}], "}"}], ",", RowBox[{"{", RowBox[{"9.95", ",", "1.4935005"}], "}"}], ",", RowBox[{"{", RowBox[{"9.94", ",", "1.5040606"}], "}"}], ",", RowBox[{"{", RowBox[{"9.93", ",", "1.5146207"}], "}"}], ",", RowBox[{"{", RowBox[{"9.92", ",", "1.5251808"}], "}"}], ",", RowBox[{"{", RowBox[{"9.91", ",", "1.5357409"}], "}"}], ",", RowBox[{"{", RowBox[{"9.9", ",", "1.546301"}], "}"}], ",", RowBox[{"{", RowBox[{"9.89", ",", "1.5568611"}], "}"}], ",", RowBox[{"{", RowBox[{"9.88", ",", "1.5674212"}], "}"}], ",", RowBox[{"{", RowBox[{"9.87", ",", "1.5779813"}], "}"}], ",", RowBox[{"{", RowBox[{"9.86", ",", "1.5885414"}], "}"}], ",", RowBox[{"{", RowBox[{"9.85", ",", "1.5991015"}], "}"}], ",", RowBox[{"{", RowBox[{"9.84", ",", "1.6096616"}], "}"}], ",", RowBox[{"{", RowBox[{"9.83", ",", "1.6202217"}], "}"}], ",", RowBox[{"{", RowBox[{"9.82", ",", "1.6307818"}], "}"}], ",", RowBox[{"{", RowBox[{"9.81", ",", "1.6413419"}], "}"}], ",", RowBox[{"{", RowBox[{"9.8", ",", "1.651902"}], "}"}], ",", RowBox[{"{", RowBox[{"9.79", ",", "1.6624621"}], "}"}], ",", RowBox[{"{", RowBox[{"9.78", ",", "1.6730222"}], "}"}], ",", RowBox[{"{", RowBox[{"9.77", ",", "1.6835823"}], "}"}], ",", RowBox[{"{", RowBox[{"9.76", ",", "1.6941424"}], "}"}], ",", RowBox[{"{", RowBox[{"9.75", ",", "1.7047025"}], "}"}], ",", RowBox[{"{", RowBox[{"9.74", ",", "1.7152626"}], "}"}], ",", RowBox[{"{", RowBox[{"9.73", ",", "1.7258227"}], "}"}], ",", RowBox[{"{", RowBox[{"9.72", ",", "1.7363828"}], "}"}], ",", RowBox[{"{", RowBox[{"9.71", ",", "1.7469429"}], "}"}], ",", RowBox[{"{", RowBox[{"9.7", ",", "1.757503"}], "}"}], ",", RowBox[{"{", RowBox[{"9.69", ",", "1.7680631"}], "}"}], ",", RowBox[{"{", RowBox[{"9.68", ",", "1.7786232"}], "}"}], ",", RowBox[{"{", RowBox[{"9.67", ",", "1.7891833"}], "}"}], ",", RowBox[{"{", RowBox[{"9.66", ",", "1.7997434"}], "}"}], ",", RowBox[{"{", RowBox[{"9.65", ",", "1.8103035"}], "}"}], ",", RowBox[{"{", RowBox[{"9.64", ",", "1.8208636"}], "}"}], ",", RowBox[{"{", RowBox[{"9.63", ",", "1.8314237"}], "}"}], ",", RowBox[{"{", RowBox[{"9.62", ",", "1.8419838"}], "}"}], ",", RowBox[{"{", RowBox[{"9.61", ",", "1.8525439"}], "}"}], ",", RowBox[{"{", RowBox[{"9.6", ",", "1.863104"}], "}"}], ",", RowBox[{"{", RowBox[{"9.59", ",", "1.8736641"}], "}"}], ",", RowBox[{"{", RowBox[{"9.58", ",", "1.8842242"}], "}"}], ",", RowBox[{"{", RowBox[{"9.57", ",", "1.8947843"}], "}"}], ",", RowBox[{"{", RowBox[{"9.56", ",", "1.9053444"}], "}"}], ",", RowBox[{"{", RowBox[{"9.55", ",", "1.9159045"}], "}"}], ",", RowBox[{"{", RowBox[{"9.54", ",", "1.9264646"}], "}"}], ",", RowBox[{"{", RowBox[{"9.53", ",", "1.9370247"}], "}"}], ",", RowBox[{"{", RowBox[{"9.52", ",", "1.9475848"}], "}"}], ",", RowBox[{"{", RowBox[{"9.51", ",", "1.9581449"}], "}"}], ",", RowBox[{"{", RowBox[{"9.5", ",", "1.968705"}], "}"}], ",", RowBox[{"{", RowBox[{"9.49", ",", "1.9792651"}], "}"}], ",", RowBox[{"{", RowBox[{"9.48", ",", "1.9898252"}], "}"}], ",", RowBox[{"{", RowBox[{"9.47", ",", "2.0003853"}], "}"}], ",", RowBox[{"{", RowBox[{"9.46", ",", "2.0109454"}], "}"}], ",", RowBox[{"{", RowBox[{"9.45", ",", "2.0215055"}], "}"}], ",", RowBox[{"{", RowBox[{"9.44", ",", "2.0320656"}], "}"}], ",", RowBox[{"{", RowBox[{"9.43", ",", "2.0426257"}], "}"}], ",", RowBox[{"{", RowBox[{"9.42", ",", "2.0531858"}], "}"}], ",", RowBox[{"{", RowBox[{"9.41", ",", "2.0637459"}], "}"}], ",", RowBox[{"{", RowBox[{"9.4", ",", "2.074306"}], "}"}], ",", RowBox[{"{", RowBox[{"9.39", ",", "2.0848661"}], "}"}], ",", RowBox[{"{", RowBox[{"9.38", ",", "2.0954262"}], "}"}], ",", RowBox[{"{", RowBox[{"9.37", ",", "2.1059863"}], "}"}], ",", RowBox[{"{", RowBox[{"9.36", ",", "2.1165464"}], "}"}], ",", RowBox[{"{", RowBox[{"9.35", ",", "2.1271065"}], "}"}], ",", RowBox[{"{", RowBox[{"9.34", ",", "2.1376666"}], "}"}], ",", RowBox[{"{", RowBox[{"9.33", ",", "2.1482267"}], "}"}], ",", RowBox[{"{", RowBox[{"9.32", ",", "2.1587868"}], "}"}], ",", RowBox[{"{", RowBox[{"9.31", ",", "2.1693469"}], "}"}], ",", RowBox[{"{", RowBox[{"9.3", ",", "2.179907"}], "}"}], ",", RowBox[{"{", RowBox[{"9.29", ",", "2.1904671"}], "}"}], ",", RowBox[{"{", RowBox[{"9.28", ",", "2.2010272"}], "}"}], ",", RowBox[{"{", RowBox[{"9.27", ",", "2.2115873"}], "}"}], ",", RowBox[{"{", RowBox[{"9.26", ",", "2.2221474"}], "}"}], ",", RowBox[{"{", RowBox[{"9.25", ",", "2.2327075"}], "}"}], ",", RowBox[{"{", RowBox[{"9.24", ",", "2.2432676"}], "}"}], ",", RowBox[{"{", RowBox[{"9.23", ",", "2.2538277"}], "}"}], ",", RowBox[{"{", RowBox[{"9.22", ",", "2.2643878"}], "}"}], ",", RowBox[{"{", RowBox[{"9.21", ",", "2.2749479"}], "}"}], ",", RowBox[{"{", RowBox[{"9.2", ",", "2.285508"}], "}"}], ",", RowBox[{"{", RowBox[{"9.19", ",", "2.2960681"}], "}"}], ",", RowBox[{"{", RowBox[{"9.18", ",", "2.3066282"}], "}"}], ",", RowBox[{"{", RowBox[{"9.17", ",", "2.3171883"}], "}"}], ",", RowBox[{"{", RowBox[{"9.16", ",", "2.3277484"}], "}"}], ",", RowBox[{"{", RowBox[{"9.15", ",", "2.3383085"}], "}"}], ",", RowBox[{"{", RowBox[{"9.14", ",", "2.3488686"}], "}"}], ",", RowBox[{"{", RowBox[{"9.13", ",", "2.3594287"}], "}"}], ",", RowBox[{"{", RowBox[{"9.12", ",", "2.3699888"}], "}"}], ",", RowBox[{"{", RowBox[{"9.11", ",", "2.3805489"}], "}"}], ",", RowBox[{"{", RowBox[{"9.1", ",", "2.391109"}], "}"}], ",", RowBox[{"{", RowBox[{"9.09", ",", "2.4016691"}], "}"}], ",", RowBox[{"{", RowBox[{"9.08", ",", "2.4122292"}], "}"}], ",", RowBox[{"{", RowBox[{"9.07", ",", "2.4227893"}], "}"}], ",", RowBox[{"{", RowBox[{"9.06", ",", "2.4333494"}], "}"}], ",", RowBox[{"{", RowBox[{"9.05", ",", "2.4439095"}], "}"}], ",", RowBox[{"{", RowBox[{"9.04", ",", "2.4544696"}], "}"}], ",", RowBox[{"{", RowBox[{"9.03", ",", "2.4650297"}], "}"}], ",", RowBox[{"{", RowBox[{"9.02", ",", "2.4755898"}], "}"}], ",", RowBox[{"{", RowBox[{"9.01", ",", "2.4861499"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "2.49671"}], "}"}], ",", RowBox[{"{", RowBox[{"8.99", ",", "2.5072701"}], "}"}], ",", RowBox[{"{", RowBox[{"8.98", ",", "2.5178302"}], "}"}], ",", RowBox[{"{", RowBox[{"8.97", ",", "2.5283903"}], "}"}], ",", RowBox[{"{", RowBox[{"8.96", ",", "2.5389504"}], "}"}], ",", RowBox[{"{", RowBox[{"8.95", ",", "2.5495105"}], "}"}], ",", RowBox[{"{", RowBox[{"8.94", ",", "2.5600706"}], "}"}], ",", RowBox[{"{", RowBox[{"8.93", ",", "2.5706307"}], "}"}], ",", RowBox[{"{", RowBox[{"8.92", ",", "2.5811908"}], "}"}], ",", RowBox[{"{", RowBox[{"8.91", ",", "2.5917509"}], "}"}], ",", RowBox[{"{", RowBox[{"8.9", ",", "2.602311"}], "}"}], ",", RowBox[{"{", RowBox[{"8.89", ",", "2.6128711"}], "}"}], ",", RowBox[{"{", RowBox[{"8.88", ",", "2.6234312"}], "}"}], ",", RowBox[{"{", RowBox[{"8.87", ",", "2.6339913"}], "}"}], ",", RowBox[{"{", RowBox[{"8.86", ",", "2.6445514"}], "}"}], ",", RowBox[{"{", RowBox[{"8.85", ",", "2.6551115"}], "}"}], ",", RowBox[{"{", RowBox[{"8.84", ",", "2.6656716"}], "}"}], ",", RowBox[{"{", RowBox[{"8.83", ",", "2.6762317"}], "}"}], ",", RowBox[{"{", RowBox[{"8.82", ",", "2.6867918"}], "}"}], ",", RowBox[{"{", RowBox[{"8.81", ",", "2.6973519"}], "}"}], ",", RowBox[{"{", RowBox[{"8.8", ",", "2.707912"}], "}"}], ",", RowBox[{"{", RowBox[{"8.79", ",", "2.7184721"}], "}"}], ",", RowBox[{"{", RowBox[{"8.78", ",", "2.7290322"}], "}"}], ",", RowBox[{"{", RowBox[{"8.77", ",", "2.7395923"}], "}"}], ",", RowBox[{"{", RowBox[{"8.76", ",", "2.7501524"}], "}"}], ",", RowBox[{"{", RowBox[{"8.75", ",", "2.7607125"}], "}"}], ",", RowBox[{"{", RowBox[{"8.74", ",", "2.7712726"}], "}"}], ",", RowBox[{"{", RowBox[{"8.73", ",", "2.7818327"}], "}"}], ",", RowBox[{"{", RowBox[{"8.72", ",", "2.7923928"}], "}"}], ",", RowBox[{"{", RowBox[{"8.71", ",", "2.8029529"}], "}"}], ",", RowBox[{"{", RowBox[{"8.7", ",", "2.813513"}], "}"}], ",", RowBox[{"{", RowBox[{"8.69", ",", "2.8240731"}], "}"}], ",", RowBox[{"{", RowBox[{"8.68", ",", "2.8346332"}], "}"}], ",", RowBox[{"{", RowBox[{"8.67", ",", "2.8451933"}], "}"}], ",", RowBox[{"{", RowBox[{"8.66", ",", "2.8557534"}], "}"}], ",", RowBox[{"{", RowBox[{"8.65", ",", "2.8663135"}], "}"}], ",", RowBox[{"{", RowBox[{"8.64", ",", "2.8768736"}], "}"}], ",", RowBox[{"{", RowBox[{"8.63", ",", "2.8874337"}], "}"}], ",", RowBox[{"{", RowBox[{"8.62", ",", "2.8979938"}], "}"}], ",", RowBox[{"{", RowBox[{"8.61", ",", "2.9085539"}], "}"}], ",", RowBox[{"{", RowBox[{"8.6", ",", "2.919114"}], "}"}], ",", RowBox[{"{", RowBox[{"8.59", ",", "2.9296741"}], "}"}], ",", RowBox[{"{", RowBox[{"8.58", ",", "2.9402342"}], "}"}], ",", RowBox[{"{", RowBox[{"8.57", ",", "2.9507943"}], "}"}], ",", RowBox[{"{", RowBox[{"8.56", ",", "2.9613544"}], "}"}], ",", RowBox[{"{", RowBox[{"8.55", ",", "2.9719145"}], "}"}], ",", RowBox[{"{", RowBox[{"8.54", ",", "2.9824746"}], "}"}], ",", RowBox[{"{", RowBox[{"8.53", ",", "2.9930347"}], "}"}], ",", RowBox[{"{", RowBox[{"8.52", ",", "3.0035948"}], "}"}], ",", RowBox[{"{", RowBox[{"8.51", ",", "3.0141549"}], "}"}], ",", RowBox[{"{", RowBox[{"8.5", ",", "3.024715"}], "}"}], ",", RowBox[{"{", RowBox[{"8.49", ",", "3.0352751"}], "}"}], ",", RowBox[{"{", RowBox[{"8.48", ",", "3.0458352"}], "}"}], ",", RowBox[{"{", RowBox[{"8.47", ",", "3.0563953"}], "}"}], ",", RowBox[{"{", RowBox[{"8.46", ",", "3.0669554"}], "}"}], ",", RowBox[{"{", RowBox[{"8.45", ",", "3.0775155"}], "}"}], ",", RowBox[{"{", RowBox[{"8.44", ",", "3.0880756"}], "}"}], ",", RowBox[{"{", RowBox[{"8.43", ",", "3.0986357"}], "}"}], ",", RowBox[{"{", RowBox[{"8.42", ",", "3.1091958"}], "}"}], ",", RowBox[{"{", RowBox[{"8.41", ",", "3.1197559"}], "}"}], ",", RowBox[{"{", RowBox[{"8.4", ",", "3.130316"}], "}"}], ",", RowBox[{"{", RowBox[{"8.39", ",", "3.1408761"}], "}"}], ",", RowBox[{"{", RowBox[{"8.38", ",", "3.1514362"}], "}"}], ",", RowBox[{"{", RowBox[{"8.37", ",", "3.1619963"}], "}"}], ",", RowBox[{"{", RowBox[{"8.36", ",", "3.1725564"}], "}"}], ",", RowBox[{"{", RowBox[{"8.35", ",", "3.1831165"}], "}"}], ",", RowBox[{"{", RowBox[{"8.34", ",", "3.1936766"}], "}"}], ",", RowBox[{"{", RowBox[{"8.33", ",", "3.2042367"}], "}"}], ",", RowBox[{"{", RowBox[{"8.32", ",", "3.2147968"}], "}"}], ",", RowBox[{"{", RowBox[{"8.31", ",", "3.2253569"}], "}"}], ",", RowBox[{"{", RowBox[{"8.3", ",", "3.235917"}], "}"}], ",", RowBox[{"{", RowBox[{"8.29", ",", "3.2464771"}], "}"}], ",", RowBox[{"{", RowBox[{"8.28", ",", "3.2570372"}], "}"}], ",", RowBox[{"{", RowBox[{"8.27", ",", "3.2675973"}], "}"}], ",", RowBox[{"{", RowBox[{"8.26", ",", "3.2781574"}], "}"}], ",", RowBox[{"{", RowBox[{"8.25", ",", "3.2887175"}], "}"}], ",", RowBox[{"{", RowBox[{"8.24", ",", "3.2992776"}], "}"}], ",", RowBox[{"{", RowBox[{"8.23", ",", "3.3098377"}], "}"}], ",", RowBox[{"{", RowBox[{"8.22", ",", "3.3203978"}], "}"}], ",", RowBox[{"{", RowBox[{"8.21", ",", "3.3309579"}], "}"}], ",", RowBox[{"{", RowBox[{"8.2", ",", "3.341518"}], "}"}], ",", RowBox[{"{", RowBox[{"8.19", ",", "3.3520781"}], "}"}], ",", RowBox[{"{", RowBox[{"8.18", ",", "3.3626382"}], "}"}], ",", RowBox[{"{", RowBox[{"8.17", ",", "3.3731983"}], "}"}], ",", RowBox[{"{", RowBox[{"8.16", ",", "3.3837584"}], "}"}], ",", RowBox[{"{", RowBox[{"8.15", ",", "3.3943185"}], "}"}], ",", RowBox[{"{", RowBox[{"8.14", ",", "3.4048786"}], "}"}], ",", RowBox[{"{", RowBox[{"8.13", ",", "3.4154387"}], "}"}], ",", RowBox[{"{", RowBox[{"8.12", ",", "3.4259988"}], "}"}], ",", RowBox[{"{", RowBox[{"8.11", ",", "3.4365589"}], "}"}], ",", RowBox[{"{", RowBox[{"8.1", ",", "3.447119"}], "}"}], ",", RowBox[{"{", RowBox[{"8.09", ",", "3.4576791"}], "}"}], ",", RowBox[{"{", RowBox[{"8.08", ",", "3.4682392"}], "}"}], ",", RowBox[{"{", RowBox[{"8.07", ",", "3.4787993"}], "}"}], ",", RowBox[{"{", RowBox[{"8.06", ",", "3.4893594"}], "}"}], ",", RowBox[{"{", RowBox[{"8.05", ",", "3.4999195"}], "}"}], ",", RowBox[{"{", RowBox[{"8.04", ",", "3.5104796"}], "}"}], ",", RowBox[{"{", RowBox[{"8.03", ",", "3.5210397"}], "}"}], ",", RowBox[{"{", RowBox[{"8.02", ",", "3.5315998"}], "}"}], ",", RowBox[{"{", RowBox[{"8.01", ",", "3.5421599"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "3.55272"}], "}"}], ",", RowBox[{"{", RowBox[{"7.99", ",", "3.5632801"}], "}"}], ",", RowBox[{"{", RowBox[{"7.98", ",", "3.5738402"}], "}"}], ",", RowBox[{"{", RowBox[{"7.97", ",", "3.5844003"}], "}"}], ",", RowBox[{"{", RowBox[{"7.96", ",", "3.5949604"}], "}"}], ",", RowBox[{"{", RowBox[{"7.95", ",", "3.6055205"}], "}"}], ",", RowBox[{"{", RowBox[{"7.94", ",", "3.6160806"}], "}"}], ",", RowBox[{"{", RowBox[{"7.93", ",", "3.6266407"}], "}"}], ",", RowBox[{"{", RowBox[{"7.92", ",", "3.6372008"}], "}"}], ",", RowBox[{"{", RowBox[{"7.91", ",", "3.6477609"}], "}"}], ",", RowBox[{"{", RowBox[{"7.9", ",", "3.658321"}], "}"}], ",", RowBox[{"{", RowBox[{"7.89", ",", "3.6688811"}], "}"}], ",", RowBox[{"{", RowBox[{"7.88", ",", "3.6794412"}], "}"}], ",", RowBox[{"{", RowBox[{"7.87", ",", "3.6900013"}], "}"}], ",", RowBox[{"{", RowBox[{"7.86", ",", "3.7005614"}], "}"}], ",", RowBox[{"{", RowBox[{"7.85", ",", "3.7111215"}], "}"}], ",", RowBox[{"{", RowBox[{"7.84", ",", "3.7216816"}], "}"}], ",", RowBox[{"{", RowBox[{"7.83", ",", "3.7322417"}], "}"}], ",", RowBox[{"{", RowBox[{"7.82", ",", "3.7428018"}], "}"}], ",", RowBox[{"{", RowBox[{"7.81", ",", "3.7533619"}], "}"}], ",", RowBox[{"{", RowBox[{"7.8", ",", "3.763922"}], "}"}], ",", RowBox[{"{", RowBox[{"7.79", ",", "3.7744821"}], "}"}], ",", RowBox[{"{", RowBox[{"7.78", ",", "3.7850422"}], "}"}], ",", RowBox[{"{", RowBox[{"7.77", ",", "3.7956023"}], "}"}], ",", RowBox[{"{", RowBox[{"7.76", ",", "3.8061624"}], "}"}], ",", RowBox[{"{", RowBox[{"7.75", ",", "3.8167225"}], "}"}], ",", RowBox[{"{", RowBox[{"7.74", ",", "3.8272826"}], "}"}], ",", RowBox[{"{", RowBox[{"7.73", ",", "3.8378427"}], "}"}], ",", RowBox[{"{", RowBox[{"7.72", ",", "3.8484028"}], "}"}], ",", RowBox[{"{", RowBox[{"7.71", ",", "3.8589629"}], "}"}], ",", RowBox[{"{", RowBox[{"7.7", ",", "3.869523"}], "}"}], ",", RowBox[{"{", RowBox[{"7.69", ",", "3.8800831"}], "}"}], ",", RowBox[{"{", RowBox[{"7.68", ",", "3.8906432"}], "}"}], ",", RowBox[{"{", RowBox[{"7.67", ",", "3.9012033"}], "}"}], ",", RowBox[{"{", RowBox[{"7.66", ",", "3.9117634"}], "}"}], ",", RowBox[{"{", RowBox[{"7.65", ",", "3.9223235"}], "}"}], ",", RowBox[{"{", RowBox[{"7.64", ",", "3.9328836"}], "}"}], ",", RowBox[{"{", RowBox[{"7.63", ",", "3.9434437"}], "}"}], ",", RowBox[{"{", RowBox[{"7.62", ",", "3.9540038"}], "}"}], ",", RowBox[{"{", RowBox[{"7.61", ",", "3.9645639"}], "}"}], ",", RowBox[{"{", RowBox[{"7.6", ",", "3.975124"}], "}"}], ",", RowBox[{"{", RowBox[{"7.59", ",", "3.9856841"}], "}"}], ",", RowBox[{"{", RowBox[{"7.58", ",", "3.9962442"}], "}"}], ",", RowBox[{"{", RowBox[{"7.57", ",", "4.0068043"}], "}"}], ",", RowBox[{"{", RowBox[{"7.56", ",", "4.0173644"}], "}"}], ",", RowBox[{"{", RowBox[{"7.55", ",", "4.0279245"}], "}"}], ",", RowBox[{"{", RowBox[{"7.54", ",", "4.0384846"}], "}"}], ",", RowBox[{"{", RowBox[{"7.53", ",", "4.0490447"}], "}"}], ",", RowBox[{"{", RowBox[{"7.52", ",", "4.0596048"}], "}"}], ",", RowBox[{"{", RowBox[{"7.51", ",", "4.0701649"}], "}"}], ",", RowBox[{"{", RowBox[{"7.5", ",", "4.080725"}], "}"}], ",", RowBox[{"{", RowBox[{"7.49", ",", "4.0912851"}], "}"}], ",", RowBox[{"{", RowBox[{"7.48", ",", "4.1018452"}], "}"}], ",", RowBox[{"{", RowBox[{"7.47", ",", "4.1124053"}], "}"}], ",", RowBox[{"{", RowBox[{"7.46", ",", "4.1229654"}], "}"}], ",", RowBox[{"{", RowBox[{"7.45", ",", "4.1335255"}], "}"}], ",", RowBox[{"{", RowBox[{"7.44", ",", "4.1440856"}], "}"}], ",", RowBox[{"{", RowBox[{"7.43", ",", "4.1546457"}], "}"}], ",", RowBox[{"{", RowBox[{"7.42", ",", "4.1652058"}], "}"}], ",", RowBox[{"{", RowBox[{"7.41", ",", "4.1757659"}], "}"}], ",", RowBox[{"{", RowBox[{"7.4", ",", "4.186326"}], "}"}], ",", RowBox[{"{", RowBox[{"7.39", ",", "4.1968861"}], "}"}], ",", RowBox[{"{", RowBox[{"7.38", ",", "4.2074462"}], "}"}], ",", RowBox[{"{", RowBox[{"7.37", ",", "4.2180063"}], "}"}], ",", RowBox[{"{", RowBox[{"7.36", ",", "4.2285664"}], "}"}], ",", RowBox[{"{", RowBox[{"7.35", ",", "4.2391265"}], "}"}], ",", RowBox[{"{", RowBox[{"7.34", ",", "4.2496866"}], "}"}], ",", RowBox[{"{", RowBox[{"7.33", ",", "4.2602467"}], "}"}], ",", RowBox[{"{", RowBox[{"7.32", ",", "4.2708068"}], "}"}], ",", RowBox[{"{", RowBox[{"7.31", ",", "4.2813669"}], "}"}], ",", RowBox[{"{", RowBox[{"7.3", ",", "4.291927"}], "}"}], ",", RowBox[{"{", RowBox[{"7.29", ",", "4.3024871"}], "}"}], ",", RowBox[{"{", RowBox[{"7.28", ",", "4.3130472"}], "}"}], ",", RowBox[{"{", RowBox[{"7.27", ",", "4.3236073"}], "}"}], ",", RowBox[{"{", RowBox[{"7.26", ",", "4.3341674"}], "}"}], ",", RowBox[{"{", RowBox[{"7.25", ",", "4.3447275"}], "}"}], ",", RowBox[{"{", RowBox[{"7.24", ",", "4.3552876"}], "}"}], ",", RowBox[{"{", RowBox[{"7.23", ",", "4.3658477"}], "}"}], ",", RowBox[{"{", RowBox[{"7.22", ",", "4.3764078"}], "}"}], ",", RowBox[{"{", RowBox[{"7.21", ",", "4.3869679"}], "}"}], ",", RowBox[{"{", RowBox[{"7.2", ",", "4.397528"}], "}"}], ",", RowBox[{"{", RowBox[{"7.19", ",", "4.4080881"}], "}"}], ",", RowBox[{"{", RowBox[{"7.18", ",", "4.4186482"}], "}"}], ",", RowBox[{"{", RowBox[{"7.17", ",", "4.4292083"}], "}"}], ",", RowBox[{"{", RowBox[{"7.16", ",", "4.4397684"}], "}"}], ",", RowBox[{"{", RowBox[{"7.15", ",", "4.4503285"}], "}"}], ",", RowBox[{"{", RowBox[{"7.14", ",", "4.4608886"}], "}"}], ",", RowBox[{"{", RowBox[{"7.13", ",", "4.4714487"}], "}"}], ",", RowBox[{"{", RowBox[{"7.12", ",", "4.4820088"}], "}"}], ",", RowBox[{"{", RowBox[{"7.11", ",", "4.4925689"}], "}"}], ",", RowBox[{"{", RowBox[{"7.1", ",", "4.503129"}], "}"}], ",", RowBox[{"{", RowBox[{"7.09", ",", "4.5136891"}], "}"}], ",", RowBox[{"{", RowBox[{"7.08", ",", "4.5242492"}], "}"}], ",", RowBox[{"{", RowBox[{"7.07", ",", "4.5348093"}], "}"}], ",", RowBox[{"{", RowBox[{"7.06", ",", "4.5453694"}], "}"}], ",", RowBox[{"{", RowBox[{"7.05", ",", "4.5559295"}], "}"}], ",", RowBox[{"{", RowBox[{"7.04", ",", "4.5664896"}], "}"}], ",", RowBox[{"{", RowBox[{"7.03", ",", "4.5770497"}], "}"}], ",", RowBox[{"{", RowBox[{"7.02", ",", "4.5876098"}], "}"}], ",", RowBox[{"{", RowBox[{"7.01", ",", "4.5981699"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "4.60873"}], "}"}], ",", RowBox[{"{", RowBox[{"6.99", ",", "4.6192901"}], "}"}], ",", RowBox[{"{", RowBox[{"6.98", ",", "4.6298502"}], "}"}], ",", RowBox[{"{", RowBox[{"6.97", ",", "4.6404103"}], "}"}], ",", RowBox[{"{", RowBox[{"6.96", ",", "4.6509704"}], "}"}], ",", RowBox[{"{", RowBox[{"6.95", ",", "4.6615305"}], "}"}], ",", RowBox[{"{", RowBox[{"6.94", ",", "4.6720906"}], "}"}], ",", RowBox[{"{", RowBox[{"6.93", ",", "4.6826507"}], "}"}], ",", RowBox[{"{", RowBox[{"6.92", ",", "4.6932108"}], "}"}], ",", RowBox[{"{", RowBox[{"6.91", ",", "4.7037709"}], "}"}], ",", RowBox[{"{", RowBox[{"6.9", ",", "4.714331"}], "}"}], ",", RowBox[{"{", RowBox[{"6.89", ",", "4.7248911"}], "}"}], ",", RowBox[{"{", RowBox[{"6.88", ",", "4.7354512"}], "}"}], ",", RowBox[{"{", RowBox[{"6.87", ",", "4.7460113"}], "}"}], ",", RowBox[{"{", RowBox[{"6.86", ",", "4.7565714"}], "}"}], ",", RowBox[{"{", RowBox[{"6.85", ",", "4.7671315"}], "}"}], ",", RowBox[{"{", RowBox[{"6.84", ",", "4.7776916"}], "}"}], ",", RowBox[{"{", RowBox[{"6.83", ",", "4.7882517"}], "}"}], ",", RowBox[{"{", RowBox[{"6.82", ",", "4.7988118"}], "}"}], ",", RowBox[{"{", RowBox[{"6.81", ",", "4.8093719"}], "}"}], ",", RowBox[{"{", RowBox[{"6.8", ",", "4.819932"}], "}"}], ",", RowBox[{"{", RowBox[{"6.79", ",", "4.8304921"}], "}"}], ",", RowBox[{"{", RowBox[{"6.78", ",", "4.8410522"}], "}"}], ",", RowBox[{"{", RowBox[{"6.77", ",", "4.8516123"}], "}"}], ",", RowBox[{"{", RowBox[{"6.76", ",", "4.8621724"}], "}"}], ",", RowBox[{"{", RowBox[{"6.75", ",", "4.8727325"}], "}"}]}], "\n", "}"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"loglogboth", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"11.03241887", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"10.28844378", ",", "0.693147181"}], "}"}], ",", RowBox[{"{", RowBox[{"10.19679129", ",", "1.098612289"}], "}"}], ",", RowBox[{"{", RowBox[{"9.981651574", ",", "1.386294361"}], "}"}], ",", RowBox[{"{", RowBox[{"9.809945808", ",", "1.609437912"}], "}"}], ",", RowBox[{"{", RowBox[{"9.69793831", ",", "1.791759469"}], "}"}], ",", RowBox[{"{", RowBox[{"9.294221856", ",", "1.945910149"}], "}"}], ",", RowBox[{"{", RowBox[{"9.20853875", ",", "2.079441542"}], "}"}], ",", RowBox[{"{", RowBox[{"9.142382679", ",", "2.197224577"}], "}"}], ",", RowBox[{"{", RowBox[{"9.13086417", ",", "2.302585093"}], "}"}], ",", RowBox[{"{", RowBox[{"9.090993614", ",", "2.397895273"}], "}"}], ",", RowBox[{"{", RowBox[{"9.037414537", ",", "2.48490665"}], "}"}], ",", RowBox[{"{", RowBox[{"8.896724917", ",", "2.564949357"}], "}"}], ",", RowBox[{"{", RowBox[{"8.847072313", ",", "2.63905733"}], "}"}], ",", RowBox[{"{", RowBox[{"8.826147399", ",", "2.708050201"}], "}"}], ",", RowBox[{"{", RowBox[{"8.801469471", ",", "2.772588722"}], "}"}], ",", RowBox[{"{", RowBox[{"8.791029857", ",", "2.833213344"}], "}"}], ",", RowBox[{"{", RowBox[{"8.775703887", ",", "2.890371758"}], "}"}], ",", RowBox[{"{", RowBox[{"8.536211197", ",", "2.944438979"}], "}"}], ",", RowBox[{"{", RowBox[{"8.47428569", ",", "2.995732274"}], "}"}], ",", RowBox[{"{", RowBox[{"8.462737006", ",", "3.044522438"}], "}"}], ",", RowBox[{"{", RowBox[{"8.456806041", ",", "3.091042453"}], "}"}], ",", RowBox[{"{", RowBox[{"8.439447843", ",", "3.135494216"}], "}"}], ",", RowBox[{"{", RowBox[{"8.438799124", ",", "3.17805383"}], "}"}], ",", RowBox[{"{", RowBox[{"8.433594168", ",", "3.218875825"}], "}"}], ",", RowBox[{"{", RowBox[{"8.428799041", ",", "3.258096538"}], "}"}], ",", RowBox[{"{", RowBox[{"8.401108712", ",", "3.295836866"}], "}"}], ",", RowBox[{"{", RowBox[{"8.373784608", ",", "3.33220451"}], "}"}], ",", RowBox[{"{", RowBox[{"8.362875831", ",", "3.36729583"}], "}"}], ",", RowBox[{"{", RowBox[{"8.32700074", ",", "3.401197382"}], "}"}], ",", RowBox[{"{", RowBox[{"8.243019469", ",", "3.433987204"}], "}"}], ",", RowBox[{"{", RowBox[{"8.240648863", ",", "3.465735903"}], "}"}], ",", RowBox[{"{", RowBox[{"8.221210094", ",", "3.496507561"}], "}"}], ",", RowBox[{"{", RowBox[{"8.217978203", ",", "3.526360525"}], "}"}], ",", RowBox[{"{", RowBox[{"8.182559264", ",", "3.555348061"}], "}"}], ",", RowBox[{"{", RowBox[{"8.157657015", ",", "3.583518938"}], "}"}], ",", RowBox[{"{", RowBox[{"8.14031554", ",", "3.610917913"}], "}"}], ",", RowBox[{"{", RowBox[{"8.110126802", ",", "3.63758616"}], "}"}], ",", RowBox[{"{", RowBox[{"8.079308192", ",", "3.663561646"}], "}"}], ",", RowBox[{"{", RowBox[{"8.00836557", ",", "3.688879454"}], "}"}], ",", RowBox[{"{", RowBox[{"7.938088727", ",", "3.713572067"}], "}"}], ",", RowBox[{"{", RowBox[{"7.895808377", ",", "3.737669618"}], "}"}], ",", RowBox[{"{", RowBox[{"7.866338923", ",", "3.761200116"}], "}"}], ",", RowBox[{"{", RowBox[{"7.860570786", ",", "3.784189634"}], "}"}], ",", RowBox[{"{", RowBox[{"7.844240718", ",", "3.80666249"}], "}"}], ",", RowBox[{"{", RowBox[{"7.836764783", ",", "3.828641396"}], "}"}], ",", RowBox[{"{", RowBox[{"7.821242084", ",", "3.850147602"}], "}"}], ",", RowBox[{"{", RowBox[{"7.81197343", ",", "3.871201011"}], "}"}], ",", RowBox[{"{", RowBox[{"7.798112629", ",", "3.891820298"}], "}"}], ",", RowBox[{"{", RowBox[{"7.770223204", ",", "3.912023005"}], "}"}], ",", RowBox[{"{", RowBox[{"7.763871288", ",", "3.931825633"}], "}"}], ",", RowBox[{"{", RowBox[{"7.688455357", ",", "3.951243719"}], "}"}], ",", RowBox[{"{", RowBox[{"7.643482907", ",", "3.970291914"}], "}"}], ",", RowBox[{"{", RowBox[{"7.628031127", ",", "3.988984047"}], "}"}], ",", RowBox[{"{", RowBox[{"7.58171964", ",", "4.007333185"}], "}"}], ",", RowBox[{"{", RowBox[{"7.545918151", ",", "4.025351691"}], "}"}], ",", RowBox[{"{", RowBox[{"7.492760301", ",", "4.043051268"}], "}"}], ",", RowBox[{"{", RowBox[{"7.480992163", ",", "4.060443011"}], "}"}], ",", RowBox[{"{", RowBox[{"7.45760929", ",", "4.077537444"}], "}"}], ",", RowBox[{"{", RowBox[{"7.445417557", ",", "4.094344562"}], "}"}], ",", RowBox[{"{", RowBox[{"7.445417557", ",", "4.110873864"}], "}"}], ",", RowBox[{"{", RowBox[{"7.428333194", ",", "4.127134385"}], "}"}], ",", RowBox[{"{", RowBox[{"7.407924323", ",", "4.143134726"}], "}"}], ",", RowBox[{"{", RowBox[{"7.398786275", ",", "4.158883083"}], "}"}], ",", RowBox[{"{", RowBox[{"7.397561536", ",", "4.17438727"}], "}"}], ",", RowBox[{"{", RowBox[{"7.374629015", ",", "4.189654742"}], "}"}], ",", RowBox[{"{", RowBox[{"7.353081921", ",", "4.204692619"}], "}"}], ",", RowBox[{"{", RowBox[{"7.340835554", ",", "4.219507705"}], "}"}], ",", RowBox[{"{", RowBox[{"7.340835554", ",", "4.234106505"}], "}"}], ",", RowBox[{"{", RowBox[{"7.329749689", ",", "4.248495242"}], "}"}], ",", RowBox[{"{", RowBox[{"7.329093736", ",", "4.262679877"}], "}"}], ",", RowBox[{"{", RowBox[{"7.268223021", ",", "4.276666119"}], "}"}], ",", RowBox[{"{", RowBox[{"7.232010332", ",", "4.290459441"}], "}"}], ",", RowBox[{"{", RowBox[{"7.231287004", ",", "4.304065093"}], "}"}], ",", RowBox[{"{", RowBox[{"7.218176838", ",", "4.317488114"}], "}"}], ",", RowBox[{"{", RowBox[{"7.202661197", ",", "4.33073334"}], "}"}], ",", RowBox[{"{", RowBox[{"7.197435354", ",", "4.343805422"}], "}"}], ",", RowBox[{"{", RowBox[{"7.168579897", ",", "4.356708827"}], "}"}], ",", RowBox[{"{", RowBox[{"7.152268856", ",", "4.369447852"}], "}"}], ",", RowBox[{"{", RowBox[{"7.150701458", ",", "4.382026635"}], "}"}], ",", RowBox[{"{", RowBox[{"7.134093721", ",", "4.394449155"}], "}"}], ",", RowBox[{"{", RowBox[{"7.129297549", ",", "4.406719247"}], "}"}], ",", RowBox[{"{", RowBox[{"7.123672785", ",", "4.418840608"}], "}"}], ",", RowBox[{"{", RowBox[{"7.117205503", ",", "4.430816799"}], "}"}], ",", RowBox[{"{", RowBox[{"7.106606138", ",", "4.442651256"}], "}"}], ",", RowBox[{"{", RowBox[{"7.084226422", ",", "4.454347296"}], "}"}], ",", RowBox[{"{", RowBox[{"7.07834158", ",", "4.465908119"}], "}"}], ",", RowBox[{"{", RowBox[{"7.060476366", ",", "4.477336814"}], "}"}], ",", RowBox[{"{", RowBox[{"7.043159916", ",", "4.48863637"}], "}"}], ",", RowBox[{"{", RowBox[{"7.03438793", ",", "4.49980967"}], "}"}], ",", RowBox[{"{", RowBox[{"7.020190708", ",", "4.510859507"}], "}"}], ",", RowBox[{"{", RowBox[{"7.013915475", ",", "4.521788577"}], "}"}], ",", RowBox[{"{", RowBox[{"6.989335266", ",", "4.532599493"}], "}"}], ",", RowBox[{"{", RowBox[{"6.96979067", ",", "4.543294782"}], "}"}], ",", RowBox[{"{", RowBox[{"6.940222469", ",", "4.553876892"}], "}"}], ",", RowBox[{"{", RowBox[{"6.923628628", ",", "4.564348191"}], "}"}], ",", RowBox[{"{", RowBox[{"6.902742737", ",", "4.574710979"}], "}"}], ",", RowBox[{"{", RowBox[{"6.864847778", ",", "4.584967479"}], "}"}], ",", RowBox[{"{", RowBox[{"6.859614904", ",", "4.59511985"}], "}"}], ",", RowBox[{"{", RowBox[{"6.856461985", ",", "4.605170186"}], "}"}], ",", RowBox[{"{", RowBox[{"6.849066283", ",", "4.615120517"}], "}"}], ",", RowBox[{"{", RowBox[{"6.837332815", ",", "4.624972813"}], "}"}], ",", RowBox[{"{", RowBox[{"6.831953566", ",", "4.634728988"}], "}"}], ",", RowBox[{"{", RowBox[{"6.820016365", ",", "4.644390899"}], "}"}], ",", RowBox[{"{", RowBox[{"6.805722553", ",", "4.65396035"}], "}"}], ",", RowBox[{"{", RowBox[{"6.80461452", ",", "4.663439094"}], "}"}], ",", RowBox[{"{", RowBox[{"6.80461452", ",", "4.672828834"}], "}"}], ",", RowBox[{"{", RowBox[{"6.797940413", ",", "4.682131227"}], "}"}], ",", RowBox[{"{", RowBox[{"6.792344427", ",", "4.691347882"}], "}"}], ",", RowBox[{"{", RowBox[{"6.788971743", ",", "4.700480366"}], "}"}], ",", RowBox[{"{", RowBox[{"6.781057626", ",", "4.709530201"}], "}"}], ",", RowBox[{"{", RowBox[{"6.759255271", ",", "4.718498871"}], "}"}]}], "\n", "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.479603999015625*^9, 3.47960400609375*^9}, { 3.47960416175*^9, 3.47960419965625*^9}, {3.479604275671875*^9, 3.479604296078125*^9}, 3.479604346109375*^9, {3.47960440078125*^9, 3.4796044120625*^9}, {3.47960447553125*^9, 3.479604518375*^9}, { 3.4796046665625*^9, 3.47960467896875*^9}, {3.47960474809375*^9, 3.479604750734375*^9}, {3.47960506415625*^9, 3.47960506553125*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{"{", RowBox[{"loglogboth", ",", "bestfitword"}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.47960401621875*^9, 3.479604060234375*^9}, { 3.479604371078125*^9, 3.479604387921875*^9}, {3.479604524375*^9, 3.479604528921875*^9}}], Cell[BoxData[ GraphicsBox[{ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxFlAs0lekax7fcbbb9fbmU79u+vbeUom1Mbl20/5njUmdGSqgmhqZtVGNQ pEKuJaU5JkW26hCHRBcaNt0Uw3HcOqaSSdpmkty6oOSWffZZq7V613rWu/7r We9v/d/3ed5HsD10o2QOi8UKU8b/96rO32yzuQvB+rS8Spa9vJKxANO/2iuc RX1ig4D3D+7tWICq9rqDgxkj4hXNM1o542ZomRmZdRaNi7saUx3KzpshpOr+ hzV3J8XFC9MtN0vMoBFhMsUumRGXZjiLCzKFcD4hurJHqBDXSDyy6ncJ8a8S V3H5BRZaCoYC+9yFyI3MU/JVkHSrN5f1lRCRqRU3y3bPQXWalixTJMQfC47I notUUXk+VbForpL38sktoz9V8ZaNFKZZgOxrAku1w2pY7ue173qZAD1jvUZN tDpaaquaZgoFIOvXLje6pI7x77ptt0oFn/xrYHL0eLfKaQFm0yfu7b6kAUPj TfqLTwggF4auiGY0EXDD3PUYKcDC6Vl1+581kaL+5uSjUT5+tU+Kdv6giYSO oO6hYT5irQ+cEn2rhYLdTtF/DfDB2mcTFXFDC0s+1HkoevgwN3v4MdZQG39Y +7tdUuqdjkdiin7URnWWBof1jI/DBqPK99OG5Juv0o538dHp92bwFFcHidtY jH8rHzH+qlVyPx1872BfHlDDh6ZLesfWizqo9ZA8barmY35kW7PjWx3IhPMK fiziY8Me48lCOzb61ol1Q2L4GI+Ld2/bz8bC+/1e9tF8sC/V3HlQzUZnGmev czgf58rq33VOsBGKpg5pKB9OqxsSZux0ETkcHLzfT+lHqNi4NVwXA+2Zu6s8 +Z/qqwufx6vDHrvxsVTF7bXFc12cVFl3bLsjH5KLLfUN8/XQ+NzbMMeCj4js pOBJDz2MMWOx7mp8dKWIVj5I1MPmzpZN0rsMnB6dcZFU6CGz6ViTx2kGBfEL NtX06aHGfDxNEMXg8F7G6K0RB5v0HkctD2cgNlPYq7ty0MFOyY8PYPD6csk+ dgQHnTL1Z52bGXwdEXebk8fBVJLX1TF3BtyXUmuTVg4K3WPGI8UMLr7wWuU4 wcGhDT2hf7NhPvWnPo4vE++oMGBQKjkn+/3v+hhCoKsth8GTA1sVARH6qHAf Dpu6awoj9u1L5mf1sTOoV/TqpCnEuTbBFnX64PrGZpkkmeImuz0qbEAfl8dt Pa94meIoS7YlUZOLvi1jS9StTWE2KWp+P5+LHLGGReYEDyp7ruUZWnGR/KzQ LnWYByuirb7XiQu30zkpWQ95YBobJyXrudiQvOYY1cjDrcDApzkBn/Xv4fpD R8O5uGXInTpbyYPUaFGydSIXZ+nBgco8HgT+T/t/OcnFgLw35NszPGTH2in/ HxdjhN53Jlk8RLw7nXGynAuDledrzsfx8NPVtU9sa5V51fMeP+/koWyLYETa rvQ/2TGr6sdD+b+rcht6PmuWRs656jfK8zFaL4s38MAODkqNnOVi+m796EdP HgwEC95N6xJ4MLiKU0zycIH7bpsHReDKvHbeqVc0XPKZZT8tJtB3eV58xzCN JYsnO7Y5EJC219TslNO4NlsspVwIqAordh25T6Nh5mZtyUYC7mobD1k10bhd 4ZHCCSDAJAXa5ZbS2JpvYeASQiC7Z9RvOIfG/pKiVM+DBL5wK4zeIKVxwuTQ R5sUAipL13avT6VRml+unD8E7ublZQQl0/hy7oU50bkEfhuiDAoO0TD4PuhE bykBX+vFVWr7aYwYTWpbVBOw0mu3WBJCI8f+rb9rPYGe5ri9Yh8aPttf7EI7 gQAirMHRk8aU6IyacTeBNKlj8isnGstWeE839St58tf0l1Y0NMyvG/u9I5DV P5hpb0bjeJwl1aZQ3sfqoLzUgEZx8dJ7pmwSlLCqS8yhUZ9eeGe9EYm9iv9+ Y/Cewk7bkdpAAYnKfaI27p8UaDo/ytuKhI5b6+X/1FK4ka2WYOlAYmn+9PXm Mgo/5Chie9eQ8IpKzwo+S4H38JVWwtckRunBImkkhTW82UoVXxLLfK4+Xx1G YVXr3Pk/BJKI9Jf7q4VQn+YzCZuYyrxwCQXfCNNFf0WSKLBUaQjZQiHhRaJ8 Io7EIx3ZhSEvCpMCb8l4KonQMJl/uRsFRmRr15Wh5Jmt9NZfTuHJ4S8ai86R +Ke55NmwA4Xr+19nbCv6rOW+Mp2payRSPxStM7ahEC66kxN/g4Tu0R0+2ywp OEk+lr+tI9F8NsiufxGFRY/PVa9rJXGgU/FLK0NhwrXO8B8dJOQ61fYHtSnI 8tQs78hJ/A8C0vXU "]]}, {Hue[0.9060679774997897, 0.6, 0.6], PointBox[CompressedData[" 1:eJwt2GdQU4nXx/HQIUBIgj3tVmEVu6wV88OGhcUFCxZ0wbao2MCydkUUCyqu FQsrKjbEhgp/UXHtBV1sCBYQUFFARXqVB5577pu8yMyZk/mez8ydsJPn+kwz lclkcSYyWdNnUsbt7lHKtlh8PmzA29npRtn/P20ROXTwk5vRGcZlzc545NSI 0AWOcLw3PNNY4ey37Wa5iFPrx459VJVpnO8mf3WkWESP2ICD/8W+NhZ5J+nX FYq4dWtm7nOfN8bpTc8nEb/nhjhlNLwx5ixtGijinWxF0NvTb41N05zfiphh WH/+/bh3xldH5jVOFFHhtq3ig2WW0adp3FMRa/329vmSkGV8nNr0iHBYFrP6 q3+2cUjTwHsinni+zK2yfG+8Vd60oIhmb8dZLJv43ti03rxrIibMynKqT3hv bBrnnSQipmbysFXyHGO3pidBRP7G/CCTgBxj43KNA0V0bB20bW1ijvGXpgVP iFhwovi8hSLXeLRpvSMiknssfBE+NddoaBoXLcLkXnWFTXKucV9U0yPCY8zK 1hGqPGPTr126U8TWj6Z9FYF5xqaf67dNxIsF4ZMir+cZm9Zz2ySijbndGnXz D8amcfp1IgJ2RB7ZOeuD1GO1iONc87stbn6QeiwV8fV81Oe9rT5KPRaK6O6u t9XM/Sj1mCdiadrhDgfvfJR6zBJx4w+n3w3aT1KP6SKsvscFxwR/knoEiPht Zedd/INPUg8/ETvsLyXGGvKlHr4iXh/o9dppUb7Uw1sE43K97mRqvtTDU8T0 5P4GF/6z1MNDRPywe+5nlnyWevQXUZo5fGrntM9SDzcRvWakrb/Q9ovUo6eI 1VWjTnZf8UXq0U3E3fDMR5eff5F6dBRh33LSt57tCqQeziJGHstVJq8ukHrw IqJc/+zm9qpA6qEX8f524eiUDoVSj9Yi2o6a95d7WKHUo5mIoLyyfbdeF0o9 HEQkBC+5NqhLkdRDLqLa5Gf2vfAiqYeFCOP2UNNhWUXkQ8R6xkpM7f6VfAhI PbvZw2vzV/IhQG1UzkzL+Uo+BIx7sjPCp+c38iHgn4mtz77Y+o18CPhYdPDp mI/fyIcAl+VcWUaf7+RDQLDt8RYT/v5OPgQk7Wvf693n7+RDQEOzjPizhmLy IeCIX5uU4mHF5EPA4KN+aV0WFZMPAV8Ko3OCY4rJh4CIbjklCanF5ENAp2W8 eXllMfkQ8OzmtOa/8j/Ih4BF8hNtF3v9IB8CWvsU9Eha8oN8CLga5TK0+ugP 8iHgj5w543un/SAfAkx/OT9rWe0P8iHg2LzS5VfblpAPAUOTXLfWe5eQDwFF DYv/6beihHwI2OZx5dyqEyXkQ0DXbXX/3nheQj4EvEzv99ykoYR8CFiiX/Oh f7tS8iFAO/1W+drRpeRDQEq8hdWd1aXkQ8Dkco9WlqdLyYcAC7dNv3i8KiUf Ak6EpfbeYFpGPgR4pio8H3QoIx8Cvjt6T5SPKyMfAv6esGPO8LAy8iHA9cjL VRFny8iHgIyCltsfvy4jHwKWdR1/WGFZTj4EGJYeSBjRpZx8CLj5b9btSL9y 8iFgmg2b/jS8nHwIsPaekq9OKCcfAuL2xlaNzConHwJGvM+32WVTQT4ElDi1 06R3ryAfAnbNDXJp6V9BPgT0TDzjNnZzBfkQ8OZnsVfU5QryIWDl4G7+r3Mq yIcAbuvC+Rr7SvLB487LxFC/npXkg0egrmbHwSmV5IOH7bS+sVlbK8kHjzOn V142XKkkHzy8y27c8/9YST54lPcxy4xRVpEPHnvXDirI7VNFPnj0eRRey/9Z RT54ZKkf2k37u4p88Fgz3k5/7FoV+eAhHvbqlP+5inzwuP8lEs7NqskHj1ld nnvPMFaTDx6KJc2nnJpZTT54nL/hu6BwVzX54DHKet86l3+ryQePqhFvd88u qiYfPPbv0Z8407KGfPDol+3/v+/9a8gHj5y2Rx52nlNDPniEzfn4Zn5UDfng 4XTZ6euF2zXkg8ej+hk/S7/XkA8ecwaddnDV1JIPHqot35hFg2vJB4+LLzp3 TZxfSz54+GpDBlQdqCUfPGqmXBrV634t+eARHVc5bWlpLfng4V7aa3Gyvo58 8PjQe/mGuqF15INHeOj1KLeFdeSDR7uHsriVh+rIB48nqgFXUx7VkQ8e88et eyyrrCMfPJrF3Mty5+rJB4/EzzbFob/Vkw8e4zt7mtz+q5588KhfvFVtcbSe fPCISUnjB/9XTz54DLRydA2vqScfPPK9Rg++L/4kHzw27d7ja+P9k3zw6JCV GThs+U/ywSNN1C7dfPwn+eCxYPakzanPfpIPHi0vHTpg//Mn+eBxpS433uuX BvLBY+JAMWXbqAbywUMW8Wda2qoG8sHh6POTOaq4BvLBwUNTVOKT3kA+OCzs 0tL0oUwGyUfj96OTNs1sIYPkg0OrJeMcbV1kkHxwKDhQsz/OXQbJB4fkG/sF T18ZJB8ctnzoG18UJIPkg8Mk6yzXLaEySD44dHZZdb3DXhkkHxxMf2c8nsTL IPng8Dzk3//m3JJB8sEhds/ksQ6ZMkg+OCxKNss5+00GyQeHIdlHZ/xubgLJ B4c2ZoNLilubQPLBobBt/tLtnUwg+eBwbdgGs66DTCD54LB1zi8Rz8abQPLB wf/vh81C5plA8sGhy+VZBx3Xm0DywcHstV3bi/tNIPng8LI+/syo8yaQfHA4 xo7oUX7XBJIPDn8NKk7Z9dYEkg8OQ2dsH/JriQkkHxw0W7o+TbcyheSDw9dz z8ct1plKPWZxuP5iQW7LbqZSj+kcIquaz0oaYir1COAQoE0sHTvJVOrhx6Eb xi6vDjGVevhysJhabb5vo6nUw5tDevi+Lb3/MZV6eHI4EdenxZuLplIPDw5L /nsbveyhqdSjP4fhpSuctO9NpR5uHLQtDeeulptKPXpy+Nb7Rs+JtmZSj24c bkwK+LeeMZN6dOTQ9HoT/auZ1MOZw5RjR5718zSTevAcuj8cOCE7wEzqoedg +e1j3qrFZlKP1hwyVOFBzBYzqUczDiddnctvHDaTejhwWDbuwYqAJDOph5yD 54qZlqZPzKQeFhz0MbbbDueZST1kHIpvn245oNpM6lHD4t/Pvx3KU5hLPcpZ 7LD77hwmmEs9illM7Rx5XuhtTj5Y/DqqS+87I8zJBwurvxrfYKaZkw8WmftD hlsuMycfLOJSmr04FmlOPlgsz7vk53HMnHyw8LLy/ZifbE4+WBjaV83e8NSc fLD44RVV4ZxvTj5Y3AruvepBnTn5YLFz9xurmWoL8sFi+pXlkXJnC/LBokeW rnWcmwX5YGFjmhIzfKQF+WDxRvRvVxRoQT5YnB5qkhCx0oJ8sFg5+3CfDjst yAeLEdsH3H580oJ8sGAvffCck2JBPliUZKx7qXhpQT5Y3K5rO+lsgQX5YLGb uf9phMySfLD4c+CMucXNLckHi16B8qrI9pbkg4U8Im51F3dL8sHi7VlPm2dj LMkHizPPv24PDrIkHyxWVW5t4xhqST5YeGs6H0nYY0k+WHDGp+1HxVuSDxZl k4Mvlt20JB8s7q53dNuVYUk+WOw5dfGO6zdL8sFixpPRXulmVuSDRe+SivRF ra3IBwu7Fnv/aNnJinyweNer1+fEgVbkg8XZia/njR1vRT5YrFmzrLpqrhX5 YOETqw2NWmdFPlgID67Je++3Ih8syosm7Xh9zop8sLinlGmX3bUiHyyiuscc 1by1Ih8sZo7t3+HqDyvywaJv0wFaWZMPFvaHwvrVa63JB4vsW+K9g12tyQeD c/l3R/QbYk0+GITaBmZkTbQmHwxGdbIJWBViTT4YiCNPfTFstCYfDCoXDQ++ EW1NPhjc31dU43/Rmnww2Hd9y1qTh9bkg0FQbke7w9nW5IOBm2Xazv7l1uSD gUO7+bo8uQ35YPD+N/WxtYwN+WBwYX5CR+FXG/LBIGzXqMTbw23IB4PR/ys3 TguwIR8MnN7tvm+x2IZ8MKiS9fQ+FmFDPhg8FDIzBx+2IR8M9g9ZOjk/0YZ8 MJgdpCkMf2xDPhgYI6+GOOfZkA8GyosT6+5X2ZAPBrmvfobNUMjJB4OE2n/s 5YKcfDBYZ3DffaqXnHww8B2Qqx8+Qk4+GDj/ufZ44VQ5+WBQs0noHLFUTj4Y PDpzJ8klUk4+GBx8Nt39caycfDCYW2H1cHaynHwwQJuTPoqncvLBQN1v2Jsz n+Tkg0FeQOGUEXVy8sHg0rqIou8qW/LBYP3JDgsjnWzJB4Oxj5/Ud3azJR8M 2v2Yu/6pjy35YFDbTOUQHGhLPhg87nlhj3qlLflgEO03kknYYUs+GMxbXXZi 5Elb8sGg/9FdXcqu25IPBo73f72y84Ut+WDwsfBVf9cCW/LB4LLDkkcvG2zJ B4MN3dqMWtTcjnwwGO+b/LZFezvywaD9Mr9pibAjHwzqo+u/+o6xIx8MntyM XlQ1y448GHDok7Fh7xo7un8D5stzwnvtsaN7N2BAx1Dl69N2dN8GNPfho5be tKN7NuDTwtusJsOO7teApKhpp5K/2tG9GrDxmmU3PzN7uk8DJuQcT65rZU/3 aEAHi6EDD3a0p/sz4KdzQarbQHu6NwPSPDePzhpnT/dlQMw8l6yVc+3pngwI 2fl4umGdPd2PAQOT5nxP2WdP92JAi7cOf/mfs6f7MOBzwzmZyV17ugcD/sf7 bIx5Y0/9DdjsUarq/8OeehvgN2vnvlxLBfU1oOM2V36tVkE9DZAlpMfxXRXU z4Cn6Yu73/ZQUC8DjtS0ujZ1ooL6GLBAf2WQRYiCehgwuP+EJ7EbFNRDj1bT 68YMjlZQDz2+bDyY/SlBQT30SI7vFxj+QEE99Ih4ml3slK2gHnpMKl+95H6Z gnro0ak1ZzpD7kA99DBxu7XJhnGgHno895/qeMrVgXrocTTM4sCw4Q7UQ49F J44Jhf4O1EMPj1SP+M2LHKiHHq2LP7u6RDhQDz0KHTddT41xoB56XO3R3mN2 ogP10GPrhNT/7B87UA89/lg1e+yZXAfqoUeXI4ocryoH6qHHmrxXFystldRD j5kLnavvqZTUQ4+RVkvc9mqV1EOPvlEP1gQ6KamHHkL7Nnd7dlVSDz3sr82U 27gpqYceFV7JXpkeSuqhR/Z72x0nfZTUQ4f7wX6vlkxUUg8dzpnHa4YFKqmH DlG76/9oE6KkHjqEOnsdLVihpB46zLoS/fnKBiX10GGU53eXzTuU1EOHvlnG +ROildRDB3Fe5KX2J5XUQweFaU51bYKSeuhQuaNLv9TrSuqhw3sxNPTAAyX1 aNw38dndoBdK6qHD+aG8rVu2knrosO9NyAj7AiX10GHt7Ns73pUpqYcOQU1/ mDUoqUfjvtunaVfKVdRDBzf+sr9XcxX10KHtJctYPaOiHjo4ePh++dZORT0a 98043iHFVUU9GvedWTV/G1TUQ4cHdUMu/zFcRT10uLA1qqbTGBX10GE/U9BP FqCiHlqsvdB7bdosFfXQImjg5nuHFqmohxaj09/Yzl+joh5a9At0+d09QkU9 tHCqWb5TtUdFPbRwiHickROjoh5aVOn0ugunVdRDi5yzcwJCE1XUQ4uH7imx PjdV1EOLhOcOBdxjFfXQYv80/46lr1TUQ4uwynPBt3JV1EOL2RtNEnd8VVEP LcZofGqnVKmohxbG+MPG7mZq6tG4r7F0rblCTT20UD4dcP9FKzX10KJ68k67 WF5NPbTILfvw+8KOauqhxaP1rrsG9VJTj8Z9W63PbD5QTT20OHAqXffJS009 tFjX12ny5XFq6qHFnCeLj62fqqYeWvj63y8YM1dNPTQwlrTq5LRUTT00cA6b EVIZpqYeGqhaXEm8t01NPTSoOS6v27NPTT00yO01AYGxauqhwaNHcWE9z6mp hwZNrx/WyWrqocHB7572mXfU1EOD9WsOep9MU1MPDeY4ftu15I2aemjgG9vv 9dBPauqhAXps07f5oaYeGvzyIHtyQa2aemigntD5+BVLR+rRuG/R6sJNKkfq oUHeyqedJmgdqYcGqUpuQXsnR+qhwaXDwUm1XRyphwbR3W/VPerrSD0a973r 6H7Aw5F6aDB37NR1QT6O1EODsQUXH/Sd6Eg9NHBfbqGwD3SkHhq0U4zxeRfs SD0a9z10bHf8Ckf8H7tabiQ= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, ImageSize->{805., Automatic}, PlotRange->Automatic, PlotRangeClipping->True]], "Output", CellChangeTimes->{{3.47960403590625*^9, 3.479604063015625*^9}, 3.4796041661875*^9, 3.47960431809375*^9, {3.479604371859375*^9, 3.47960441353125*^9}, {3.479604520921875*^9, 3.4796045296875*^9}, { 3.4796046721875*^9, 3.479604680859375*^9}, 3.479604753390625*^9, 3.479605069890625*^9}] }, Open ]], Cell[BoxData["\[AliasDelimiter]"], "Input", CellChangeTimes->{3.479604434734375*^9}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxdz01I03Ecx/EfnSRiSXjYIcJCJCRimZk92H7ZKstlZllZHcQWdKqOEUI/ Oo0yerIys1hl9EymED0Y/HoCA+kBCiI6GHjoMGSHUR081Py8fwcbbO+9/h++ jM1tP9i8b5ox5tC/d6GF12jluNW3XLzw6Qe/YzNZG/uGiybr5nzBxfK7T7hE ff8BR7XXvsWz1aOvcan2+Cs8T50eXKa95SUuVw94PF/tHMIVasVzvEDteYYX qvmnOKY+DF6kNgRXqj8H8WL18QCu4v8EL1H3Blers4KXqjODa9Tfj/Ay9Ufw cjUVvELN9uOV6ovgWrU8eJU68ADH1UywVSPY4V/32Ver+XvsOIFNnTpylx0b bNao1XfYccct9oQ6dpMdT/Sxr1XbsMN1YV+njt9gx8PYrFePhB0nw16vFoUd H7vOvkFNYYcj19g3qo0ZdhzBpkG9fZUdd11hT6ozetnxk8vsm9QS7LDvYW+c aoeTYd+svrnEjtPYNPEcW5zt5jlOYY+j2GyZaos/X+Qef8Qe77/AfbPae557 3I0dbsIeT3Rxv1UdxRbXY4fz57jHWWy28XvYYocdToV7nAj3LWpxuMd9Z7nH aexxFTbb1bEz3OMh7HAH9rgNmx38v9Pc4wx2uAx77E9xv1ONYotzJ7nHw9jj fmxa1cPY4tJwj/90co8z2OxSa7DFMez+s8cm3O9Wv57gHnvscBp7nDvO/R77 Fx6B+BU= "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, ImageSize->{791., Automatic}, PlotRange->Automatic, PlotRangeClipping->True]], "Input"] }, WindowSize->{942, 694}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 499, 13, 167, "Input"], Cell[CellGroupData[{ Cell[1069, 37, 664, 15, 72, "Input"], Cell[1736, 54, 262, 6, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2035, 65, 271, 5, 52, "Input"], Cell[2309, 72, 467, 9, 30, "Output"], Cell[2779, 83, 447, 9, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3263, 97, 2939, 89, 132, "Input"], Cell[6205, 188, 238, 3, 30, "Output"], Cell[6446, 193, 239, 3, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6722, 201, 2227, 68, 152, "Input"], Cell[8952, 271, 180, 2, 30, "Output"], Cell[9135, 275, 184, 2, 30, "Output"] }, Open ]], Cell[9334, 280, 400, 10, 72, "Input"], Cell[CellGroupData[{ Cell[9759, 294, 707, 19, 52, "Input"], Cell[10469, 315, 439, 15, 47, "Output"], Cell[10911, 332, 365, 12, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11313, 349, 175, 4, 31, "Input"], Cell[11491, 355, 195, 6, 56, "Output"] }, Open ]], Cell[11701, 364, 295, 11, 31, "Input"], Cell[CellGroupData[{ Cell[12021, 379, 364, 11, 31, "Input"], Cell[12388, 392, 95, 2, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12520, 399, 290, 8, 31, "Input"], Cell[12813, 409, 93, 2, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12943, 416, 142, 3, 31, "Input"], Cell[13088, 421, 91, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13216, 427, 113, 2, 31, "Input"], Cell[13332, 431, 72, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13441, 437, 148, 3, 31, "Input"], Cell[13592, 442, 88, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13717, 448, 215, 6, 31, "Input"], Cell[13935, 456, 116, 2, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14088, 463, 113, 2, 31, "Input"], Cell[14204, 467, 96, 1, 30, "Output"] }, Open ]], Cell[14315, 471, 287, 7, 58, "Input"], Cell[CellGroupData[{ Cell[14627, 482, 3177, 80, 292, "Input"], Cell[17807, 564, 2618, 48, 240, 504, 12, "CachedBoxData", "BoxData", "Output"], Cell[20428, 614, 716, 16, 246, "Output"], Cell[21147, 632, 894, 19, 246, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22078, 656, 3038, 97, 192, "Input"], Cell[25119, 755, 182, 2, 30, "Output"], Cell[25304, 759, 182, 2, 30, "Output"], Cell[25489, 763, 199, 3, 30, "Output"], Cell[25691, 768, 179, 2, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25907, 775, 120, 2, 31, "Input"], Cell[26030, 779, 109, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26176, 785, 6075, 122, 572, "Input"], Cell[32254, 909, 1785, 33, 239, "Output"], Cell[34042, 944, 2877, 51, 243, "Output"], Cell[36922, 997, 4382, 76, 243, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41341, 1078, 2970, 98, 192, "Input"], Cell[44314, 1178, 202, 4, 30, "Output"], Cell[44519, 1184, 185, 3, 30, "Output"] }, Open ]], Cell[44719, 1190, 170, 4, 58, "Input"], Cell[CellGroupData[{ Cell[44914, 1198, 4431, 102, 432, "Input"], Cell[49348, 1302, 1077, 22, 239, "Output"], Cell[50428, 1326, 1809, 34, 246, "Output"], Cell[52240, 1362, 2649, 47, 246, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[54926, 1414, 2997, 97, 192, "Input"], Cell[57926, 1513, 141, 2, 30, "Output"], Cell[58070, 1517, 173, 3, 30, "Output"], Cell[58246, 1522, 157, 2, 30, "Output"] }, Open ]], Cell[58418, 1527, 42563, 1098, 2352, "Input"], Cell[CellGroupData[{ Cell[101006, 2629, 287, 6, 31, "Input"], Cell[101296, 2637, 9223, 157, 516, "Output"] }, Open ]], Cell[110534, 2797, 85, 1, 31, "Input"], Cell[110622, 2800, 999, 21, 503, "Input"] } ] *) (* End of internal cache information *)