Eden Model2.nb

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

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

(* CreatedBy='Mathematica 11.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     82160,       1685]
NotebookOptionsPosition[     78097,       1614]
NotebookOutlinePosition[     78449,       1630]
CellTagsIndexPosition[     78406,       1627]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 "Physics 234  \nPS#9(B)\n",
 StyleBox["Eden Model: Border versus Area",
  FontSize->36]
}], "Title",
 CellChangeTimes->{{3.731933741464285*^9, 3.731933811667729*^9}, {
   3.733086138314348*^9, 3.7330861788780413`*^9}, 3.733089827843464*^9, {
   3.733587756424402*^9, 
   3.733587756706798*^9}},ExpressionUUID->"4ff34a9b-c9e3-4e77-b2ad-\
1e53638bf927"],

Cell[TextData[StyleBox["Ziyang Gao\n4/18/2018",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.731933819521956*^9, 3.731933824538842*^9}, {
  3.731939225488948*^9, 3.7319392305278673`*^9}, {3.733088578727078*^9, 
  3.733088579284164*^9}},ExpressionUUID->"60176971-891a-498f-b05e-\
1c7b25d8f0be"],

Cell[CellGroupData[{

Cell["Section 1: What is the problem about", "Section",
 CellChangeTimes->{{3.731933877737628*^9, 3.731933884921731*^9}, {
  3.731934100769842*^9, 3.7319341152049837`*^9}, {3.7319341676901703`*^9, 
  3.731934206224195*^9}, {3.7319396024353*^9, 3.7319396027873163`*^9}, {
  3.7330864596273613`*^9, 
  3.733086466532071*^9}},ExpressionUUID->"9463184e-e0ea-4946-b109-\
adcbb5c7ef73"],

Cell[TextData[{
 StyleBox["This problem is about investigating the Eden Model and trying to \
determine how the number of sites b in the border of the cluster depends on \
the number of sites n in the cluster itself. In particular, you will assume \
that b = ",
  FontSize->18],
 Cell[BoxData[
  FormBox[
   SuperscriptBox["cn", "d"], TraditionalForm]],ExpressionUUID->
  "0151a163-4cdc-4be1-a99b-af7c4a6ef2f1"],
 StyleBox[", and then estimate the best values and uncertainties for log[c] \
and d.",
  FontSize->18]
}], "Text",
 CellChangeTimes->{{3.731934230992525*^9, 3.731934244038454*^9}, {
  3.731934291164131*^9, 3.731934299507243*^9}, {3.7319343516085052`*^9, 
  3.73193435360321*^9}, {3.73193440429664*^9, 3.731934488422429*^9}, {
  3.731935615701275*^9, 3.731935618086338*^9}, {3.7319391777147427`*^9, 
  3.7319392070884123`*^9}, {3.733086212566287*^9, 3.733086219183638*^9}, {
  3.733086323993503*^9, 3.7330864097420692`*^9}},
 FontSize->16,ExpressionUUID->"4fad08c4-ae96-4325-9755-20c424ac9425"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Section 2: Statistical functions", "Section",
 CellChangeTimes->{{3.731935483103862*^9, 3.731935498483976*^9}, {
  3.731935551408482*^9, 3.731935582598387*^9}, {3.73308646927219*^9, 
  3.733086473055293*^9}},ExpressionUUID->"79597f95-25b4-4b2b-9164-\
b6fa06ab1942"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"ave", "[", "myList_", "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"Total", "[", "myList", "]"}], "/", 
     RowBox[{"Length", "[", "myList", "]"}]}], "//", "N"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"aveSqr", "[", "myList_", "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"Total", "[", 
      RowBox[{"myList", "^", "2"}], "]"}], "/", 
     RowBox[{"Length", "[", "myList", "]"}]}], "//", "N"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"var", "[", "myList_", "]"}], ":=", 
   RowBox[{
    RowBox[{"aveSqr", "[", "myList", "]"}], "-", 
    RowBox[{
     RowBox[{"ave", "[", "myList", "]"}], "^", "2"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"unc", "[", "myList_", "]"}], ":=", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"var", "[", "myList", "]"}], "/", 
     RowBox[{"Length", "[", "myList", "]"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.731935592727912*^9, 3.731935609122292*^9}, {
   3.731935825100363*^9, 3.7319358338529053`*^9}, {3.7319361293425293`*^9, 
   3.731936130006646*^9}, 3.7319367009544573`*^9, {3.731937088503779*^9, 
   3.731937096008078*^9}, {3.731939458631688*^9, 3.731939463061267*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"06f8333f-1096-4432-b651-057b0df4e38c"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Section 3: Develop a module for generating sizeList", "Section",
 CellChangeTimes->{{3.731939470272645*^9, 3.731939478755679*^9}, {
  3.7319395843885317`*^9, 3.731939597387458*^9}, {3.7319400104232397`*^9, 
  3.731940018718193*^9}, {3.7319443041699667`*^9, 3.7319443055850058`*^9}, {
  3.731949509578136*^9, 3.731949509808877*^9}, {3.733088547548583*^9, 
  3.733088553676586*^9}},ExpressionUUID->"cb4ce077-be6b-4232-94f6-\
9320176a7592"],

Cell[TextData[StyleBox["Code for a single growth:",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.733089589318679*^9, 
  3.733089611617515*^9}},ExpressionUUID->"b3c68032-dd45-4fbf-8e7b-\
cd945415abf5"],

Cell[BoxData[{
 RowBox[{"Clear", "[", "growCells", "]"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"growCells", "[", "oldCells_", "]"}], ":=", "\[IndentingNewLine]", 
   
   RowBox[{"Module", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "choices", ",", "oldCluster", ",", "oldBorder", ",", "newSite", ",", 
       "nnNewSite", ",", "tempBorder", ",", "newBorder", ",", "newCluster", 
       ",", "newCells"}], "}"}], ",", "\n", "\t", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"choices", "=", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"1", ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"0", ",", "1"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"0", ",", 
           RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"oldCluster", "=", 
       RowBox[{"oldCells", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"oldBorder", "=", 
       RowBox[{"oldCells", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"newSite", "=", 
       RowBox[{"oldBorder", "[", 
        RowBox[{"[", 
         RowBox[{"RandomInteger", "[", 
          RowBox[{"{", 
           RowBox[{"1", ",", 
            RowBox[{"Length", "[", "oldBorder", "]"}]}], "}"}], "]"}], "]"}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"newCluster", "=", 
       RowBox[{"Join", "[", 
        RowBox[{"oldCluster", ",", 
         RowBox[{"{", "newSite", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"nnNewSite", "=", 
       RowBox[{"Map", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"#", "+", "newSite"}], ")"}], "&"}], ",", "choices"}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"tempBorder", "=", 
       RowBox[{"Union", "[", 
        RowBox[{"oldBorder", ",", "nnNewSite"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"newBorder", "=", 
       RowBox[{"Complement", "[", 
        RowBox[{"tempBorder", ",", "newCluster"}], "]"}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"newCells", "=", 
       RowBox[{"{", 
        RowBox[{"newCluster", ",", "newBorder"}], "}"}]}]}]}], "\t\t", 
    "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input",
 CellLabel->"In[9]:=",ExpressionUUID->"b7df351a-bfac-4d07-8ad5-7f5fafa20ca2"],

Cell[TextData[StyleBox["A for-loop generating sizeList of 1000 times of \
growth:",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.7330896244824953`*^9, 3.733089647764377*^9}, {
  3.733089682888604*^9, 3.7330896831314363`*^9}, {3.7331596246666307`*^9, 
  3.733159646776009*^9}},ExpressionUUID->"15d372f7-6c98-4b21-9ceb-\
2509e7a9fbb5"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"oldCells", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"1", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", 
         RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sizeList", " ", "=", " ", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"newPic", "=", 
      RowBox[{"growCells", "[", "oldCells", "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"clusterArea", " ", "=", " ", 
      RowBox[{"Length", "[", 
       RowBox[{"newPic", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"borderArea", " ", "=", " ", 
      RowBox[{"Length", "[", 
       RowBox[{"newPic", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"AppendTo", "[", 
      RowBox[{"sizeList", ",", 
       RowBox[{"{", 
        RowBox[{"clusterArea", ",", "borderArea"}], "}"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"oldCells", "=", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"newPic", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"newPic", "[", 
         RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{"sizeList", ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Blue", "[", 
         RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"PointSize", "[", "0.01", "]"}]}], "}"}]}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<n\>\"", ",", "\"\<b\>\""}], "}"}]}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<cluster size n  versus  bordering b\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "12"}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}]}], "]"}]}], ",", 
     RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"model", "=", 
   RowBox[{"LinearModelFit", "[", 
    RowBox[{"sizeList", ",", "x", ",", "x"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot2", " ", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"model", "[", "\"\<BestFit\>\"", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "2", ",", "1000"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", "plot2"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"model", "[", "\"\<ParameterTable\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"model", "[", "\"\<BestFit\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"model", "[", "\"\<ParameterErrors\>\"", "]"}]}], "Input",
 CellChangeTimes->{{3.733089487731042*^9, 3.7330895068450003`*^9}, {
   3.733089712326661*^9, 3.733089773033499*^9}, {3.73308984984695*^9, 
   3.733089857553776*^9}, {3.733089976835205*^9, 3.7330900839813128`*^9}, {
   3.733090191850224*^9, 3.73309034019281*^9}, {3.733090376325292*^9, 
   3.733090390667686*^9}, {3.733091018487739*^9, 3.733091056932004*^9}, {
   3.733091099089858*^9, 3.7330911255347548`*^9}, 3.733091206108118*^9, {
   3.733091274859486*^9, 3.733091276546442*^9}, {3.733091391954361*^9, 
   3.733091459545443*^9}, 3.73309199381398*^9, {3.733092070017365*^9, 
   3.7330921384838123`*^9}, {3.733095293151875*^9, 3.733095308791675*^9}, {
   3.733095434603629*^9, 3.733095461697607*^9}, {3.733095512707172*^9, 
   3.733095585486483*^9}, {3.733095644477758*^9, 3.733095689486764*^9}, {
   3.733095834288021*^9, 3.733095834955641*^9}, 3.733095865432475*^9, {
   3.733095908857945*^9, 3.7330959093796043`*^9}, {3.733095955014286*^9, 
   3.733095956491598*^9}, {3.733096059198093*^9, 3.733096063718175*^9}, {
   3.733096096564747*^9, 3.733096157665468*^9}, {3.733096190390738*^9, 
   3.733096200430606*^9}, {3.733096232670431*^9, 3.733096241114482*^9}, {
   3.733096273486993*^9, 3.733096275014914*^9}, {3.733097221174057*^9, 
   3.7330972563354*^9}, {3.7330986980681467`*^9, 3.733098705416177*^9}, {
   3.733143484269719*^9, 3.733143495534583*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"35ad2452-f27c-4c52-b62d-92a2a587d859"],

Cell[BoxData[
 FormBox[
  GraphicsBox[{{{}, {{}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.01], 
       AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJx12zuw61YVxvE90LhUQeGCQjAM48kwjHib94a8HJ4iEOIkQJQX8c0luUrC
TZwH3J0XcZILUXlKlS5VnlKlS5UuVbpU6ZJh+JaZ9Z8hjeY331p7b8n2key7
84Wnbz383GdCCKfPhvCfo/6L/z3MdZzp+HkdMx1zHT+n4xfR92X0LdB3D/rM
X9LxKxjHvIDvge34VRwLHL+G49f/z/EbqPsm8m+h7tvwd+Al+r6Lcb+H/PsY
9wcY54fwjzDujzFO9OMks43zEz9OgsNPkZttnHu9072Y5z7k9yG/368j3Y91
PID8AfQ/6OvTgzivFfIVzush7/QQ5v8Z+s3W/3Md7bxluw7hFzraect2HcIv
ZTtvOPwKORx+jRwOpXeU7boms65bb/V2HX8j67pFOMl2HXs4POwdH0Y/3KM+
/Fa59cMJ7uHwO/TDSbbXuUceHsH6H8H6Lbf1W27X7/feUb68by3X69HD4VFf
Hx/14yW4t3p7n6+x/rU/vwT3qA+PKbfXH06P+foeeXgc/Y/79SWz9cPhCfTD
6Qk/Xo88/MGPF+XL+Vtur7/ZPtd/VJ0cZftcJ8vl3nJ9jsOfvKN8+Zwj7+Hw
pK+PT/o8WW7zw6Hy9XmF8SzX35XKcvs7A7eo7zHeiDw85Z0/hfmRV3BCfWu5
XS/kIxyexvxP+/OJyCuz/R2FW7iHR4wfnpGV53CEq2cwP9zCvfVrvhF5eNY7
h+OzmN9yuw/ALdzLdl8Z4fCc6uRctvtKhCvUJ9nuOy3cy3YfGq3f7kN/1lH3
lRyOst2XKvlyX4Jb2e5TPTzC4XkddV/K4fi8r6+QJ9nuey3yHv3j8369YeOd
b3x93PjxKjhZv80P9xu/vhH94YafL7/h83jD91c3/PjJ+m1+uLfxdN8ebTy7
r7/g8/wF3x/hCk7ob+Ee9SMcbvrzyW/6843IKzjdxPw3/fg98tH67fz/4p3L
9pwTkVdwku05qIV72Z6LRji8qHHk/EWfR9memyrUJ7hFfw+PqA8v+Tx/yefx
JcwvX57brF9u4d7Gk0c43PLOb/n+eMvPV8EJbm/59fbIR4wfaj9/JttzZQ4X
sj1nRrhEfWXjyzXqE/JGtufU1ur13NlZvdzDAzzCk63Pnktf1jjKs5d9nsv2
XFvA0frlEq7g2sa3517ZnnMby/Wc2sKdzS/38GDjyaPNL082nj3HvuLz7BWf
53ABR7iU7bm5slzPvbVsz8kJeWP9cmu5nns7W6/yHvUD8hGeUB9e9c5e9fX5
q37+Qrbn9Gj9com8Ql4jT3CD+hZ5h7y39coDPNr69dw/2Xj2PeGv3hmcw4Vs
3yuibN8jStm+Z1TIazjJ9r2jgVu4Q3+PfIBH1E9YX7jtnd3255fLl/NHHtFf
Iq+Q17f9ehLc2Hx2/rf9+XSyfQ/rUT+gfkQ+oT+8pnXI2Ws+z1/z/QUc4RLj
VbJ9L6zly/dC2b4XNrJ9L2yRd3CP8QZ4RP2EPLzu15fB+ev+fAo4or5EXsE1
6hPyBm5lez062a53L1/e/7K9n0bkE/rD1o+fbX2eIy/gCJfor7Z+/nrr15dQ
36C+RX0H9+gf0D8in+Dwhq/P4PwNX1/AUb6cP1zBtXz53US290Mj2/uzRd6h
v8f4A/IR/RPGD296Z3AOF3CES7iCazi96dfXvOnX3yLv4B4e0D/Kdn0m5OEt
7wzO4QKOb/nxS7iCa/Qn5A3cvuXPr5Pt+vXIB9n+3o2y/T2dkIe3fZ7B+du+
voDj2349JVyhvkaekDeYv0XeIe/hQbbfEUd4gsPftE47f9l+V8xl+12xkO13
xSjb74qlbL8rVnANJ/Q3GL/F/B3W12P9AzzCExz+7p3J9rtrDhdwhEvZrm+F
vIYT5m+Qt8g75D08YP4RnuTL7853vDM4v+P7izt+vgiXcAXXcLrjz69B3mL+
DuvrZft8DPLl+Ue+/P2HQ/L1s+THy5Kfb578enLkC/QX8BLzRXgFl8mvd43x
KniD/hr9W+QJ/TvkDfIr5C28x3wd8muM18MHeMD1PSIfkZ+QT5j/jDy84/tn
cPaOr58jz+EFXKB/iTzCK/nyeYLXsn1eKuQb5DXyLZxkez/vsJ4G9VcYv5Xt
7/EeeYf+a7jH/AfMP6D+iPoR9SfkE/Iz8vCuH3/2rl9/9q4/vznyHP0L1Bey
3b+WyCO8Qn0p2/1wLdv9soI3cA1v4QTvZLs/N3IuX8GtXMl7uJOTfI28hw/w
ILfyER7hE+ab4DPGD+/5fAZn7/n6OZyjfgEX8BL9EfnqPX8+JfI1+it4I0e5
hrcYL6F/Bzfov4JbeA938DXG75EfkA/wER7hEzzBZzi87+efyZfXX7bXY448
hxfv+/ELeIn5omyfpxVcon+N+SqsbyOPcm3j6d+9tza//bu4jS/vkDfIr2x+
uUW+R97B17L9u3xv65cPyAfkR+Qj8hPyCfkZefhA48gzOIPnH/j+HF7Itq+g
gJcYLyJfwaVs+xLWsu1LqJBv0F9jvi3qE8bbIW/gK7hF/x55h/wa7mXbV3GQ
bV/FAB9l22cxwid7PeTJzl8+Iw//8PkMzuC5bPs2ctn2bSxk26dRwEvZ9nVE
eIXxSniN+Sv0b2TbF1LDW6wnoX+H+RrkV3CL8fZwh/prrKeXbZ/KAR7go52/
9pGMtl75BE/w2dZj+1w+1HrkmWz7XDJ4jvoc+eJDP36BfIk8Il/BJbyGK4y3
QV5jvVvkCf071DfwlWz7fFrZ9vnsZdvn08m2z+ca7uED+gfkR3iUbd/Qyc5H
nuCznb/tE9rp/OQZnMm2r2gO57LtM1ogL5Av4Yj6FVzCa6yvgjeor+EtnLCe
HfIGvkJ9i3yP9XSyXe9ruJft9TkgH5Af4VG+vP6yvT8m2d5PZ1u/ff4/8vkM
zj7y480/8vPlyBfoL5Av0R/hFVzCa9muTyVfXn+4Rv0WeUK+Q97Il9cfbmV7
P+zhzq6Hvf7Ie/iA+sGun/bNHeERPqF/suun/Iz68LH37GNfn8m2D3Au276/
XLZ9fwvZ9v0VyJdwxHgruET9Gq4w3wZ5jfG2su1jTMh3cIPxrpC3GG+P69Wh
/hr1PfIDPKD+iHzE+Z9k22c5IT8jD594zz7x9RnyOfJcvrz+cIH+Jfoj6lfI
S+Rr2c6/gjeyXa8a+RZ5Qr5D3sj2el4hb9G/x3o7nM813MMH9A/Ij8hHzH+C
J6z3jDzc9Z7d9fWZfPn8wzm8kO3vR4Hxlnf9+iO8uuvPt0S+xnoreIP6Gt6i
PiHfYf5GtvfvFdyifg93qL+262P7lOEDPMi2T/oIj6g/2eshT8jPcPinr5/J
tq86k21f9RzOZdtnvUBeIF/Kl33YyFdwifo18gr5Bq7hrWz7vJNs+8B3qG/g
K/S3dv3kPfIO419jvB75AR7gI8Yf4RM8ybZv/QyHf3nP4Ey2fe5z2fa55/AC
LtC/xPgRXqG+xHhr5BW8kW3ffY3+LeoT8h3yBvkV8hbe43w62d4f1/Ll/zOQ
L68/8gH5ER7hEzxhvDPy8KnGsc8/nH3q6+fIc3gBF/AS40Xkq//532Mecgs=

        "]]}, {}}, {}, {}, {}, {}}, {{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJwVjnk4lAkAh8cRsbWVW7alZN25kvTY/Cp30hijGFRkxvCV0S6SMzmS+L5m
h6lF1tVkdklNbTZHjrbEaFJblmW3C9OmRAoda9f+8T6/5/3n97yrI3gMtiKN
RtuxwP/rc7W2Xp1GA6s9Y1kFyx3qDknPx75RhCBbkagPd0dPg6fxjVFl9Hrl
3GriumO7eOxMilQNrvfy0x4muuMrP42tidIVMH5aPK4uWPAtdGu50AAvleq6
EqTuEGiOZY4ITZHuOZCx3cUDPZ9srZ4JHVEls52cWeWJonZLzd7ohV81bq2b
qjc0JbUFisyd2EH2p5v0+UDXLZvnlMsETSxUWX7KF8395eFPpcHQmhxmTyX4
QXXLhH9K4x60o35KsoqOBPOiHDeVCJw3tCnjavqjcDx6IKgmEpzuYMNqZQaO
6LMsIrWjwNCrsDs/yoBdadtj8+ZoxKaJG0SDAWjNcxWOORxA43FZC6+VCZ7N
XIdtRSwGQsuqWssCMaOlKS6piIOZl/VbFWoXDL0a4878cQiigcEH+od2w9pE
GDO97lswM7xXm3GCcHcP31A7LR5Dk+EezbuCESLTZ0quJIBUFX7QdWKhVz7c
VKdyGIfWF66MswiBmZN15ZqtSchUqejmLw1FjO0v196dPoLOiYSRC0phcM3d
W0jJkzFU+fZ3+rsw/LOi6HGNUSqMPpCslv49qPRU4ujGp6GY+9LLWLoXNpf/
jVl0KR1fXLd/Jt69D/4pepZN8gwE128js7v2oW7m5885ZUdxwXTD62GE45yP
qKvOORMl8Vd1V0rCwZ++rL9mJBPlvnozD6wi0NWzszv55DEsy5eBcS4Cb1wk
G/mmWdAp0Qi4or0fy5Xtt/T3ZWG940aDImo/jDlJaoOx2XhRNf5jkEIkMmZ8
bdrUcvBZ4rToUVwkZEObJ367lAO3B8NBI39HQiN09oDILxdai5j377DY8La7
k31+Lhef2usUHB+ykVAZHco5fRzJOwVzRz05GHct4i92yAN7myyEauSA6TC8
d3AoD9zOF/EGdlGQdLx4FXXkBBS0igMLqqKgdE5H22RtPiSTRi+TV3Fx9mrf
9zel+bhrr/DaUsiF1X7/v4y5J5HlybvWoRKNC5Gvmy+uKICf/Lrzw6PRMFFg
Bn9sKMDF4LLMO2+iMSb1YJpvL0SGncwxLyIGJeO7nudMF4Jd9OSH9qEYbBNZ
u1ScIvG0+9i9P30JLHbWaA4TkFBn24z2+xHolc44GwhJ2NGG5+7SCTCn2pyE
pSTSNzgadTIJRLowHApEJFZWy2NFoQSy7ydaJLWQoKf6LuEdIHBrvk2P/pxE
2zodT1oBgXx+zeml4yTGujtZ7wsJ+K09oSN9RWIpmxf7hiLQ783Q8pgmEVJ6
u3hEQEBeNLLMZZ7ErGrKyO3ShT5LNVUzLQo2Tx5l8n9a6Gt9lTOqQyEwtaA4
v57AKfp95Wp9Cqm6zuKsBgJ6iaWKXxpS6NnxXV/CZQLmHdbzmhYUOE1uhqxm
AhMBGmn3rCiQgVP2Aa0EJGMzH0kbClcmz3r4thHYtKT9vZojBUXT2YObbxCY
L69J6nKiYN5Znel0k0Cn3YnZ7E0U6GH0YtsuArm/Hkzc+jWFw3Ofas27Cfjs
Zrz715VCuUDcskZK4D/81TRo
         "]]},
       Annotation[#, "Charting`Private`Tag$3596#1"]& ]}, {}, {}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{0, 0},
   DisplayFunction->Identity,
   Frame->{{True, True}, {True, True}},
   FrameLabel->{{
      FormBox["\"b\"", TraditionalForm], None}, {
      FormBox["\"n\"", TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   GridLines->{Automatic, Automatic},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImagePadding->All,
   ImageSize->500,
   Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
   PlotLabel->FormBox[
     StyleBox[
     "\"cluster size n  versus  bordering b\"", FontFamily -> "Helvetica", 
      FontSize -> 12, FontWeight -> "Bold", StripOnInput -> False], 
     TraditionalForm],
   PlotRange->{{0, 1001.}, {0, 214.}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.02], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
 CellChangeTimes->{{3.733092106441432*^9, 3.733092114376828*^9}, {
   3.733095534079273*^9, 3.733095588293337*^9}, 3.733095691409423*^9, {
   3.733095841551009*^9, 3.733095867662918*^9}, 3.733095911439714*^9, 
   3.73309596097033*^9, 3.733096067376163*^9, 3.7330961608936*^9, 
   3.73309620250974*^9, 3.733096245069832*^9, {3.733097259191615*^9, 
   3.733097273378603*^9}, 3.7330987083844013`*^9, 3.733143503210012*^9, 
   3.733146779301908*^9},
 CellLabel->"Out[13]=",ExpressionUUID->"c9e48d74-3546-40ed-9446-b5c4202ab4cf"],

Cell[BoxData[
 FormBox[
  TemplateBox[{
   "General","munfl",
    "\"\\!\\(\\*FormBox[SuperscriptBox[\\\"0.23398677352908603`\\\", \
\\\"499.`\\\"], TraditionalForm]\\) is too small to represent as a normalized \
machine number; precision may be lost.\"",2,14,1,17856018701862124304,"Local"},
   "MessageTemplate"], TraditionalForm]], "Message", "MSG",
 CellChangeTimes->{3.73314350369208*^9, 3.733146779928314*^9},
 CellLabel->
  "During evaluation of \
In[7]:=",ExpressionUUID->"9c6428c8-3148-4b33-adee-c444c14ac186"],

Cell[BoxData[
 FormBox[
  TemplateBox[{
   "General","munfl",
    "\"\\!\\(\\*FormBox[SuperscriptBox[\\\"0.04806740972032362`\\\", \
\\\"499.`\\\"], TraditionalForm]\\) is too small to represent as a normalized \
machine number; precision may be lost.\"",2,14,2,17856018701862124304,"Local"},
   "MessageTemplate"], TraditionalForm]], "Message", "MSG",
 CellChangeTimes->{3.73314350369208*^9, 3.733146780638442*^9},
 CellLabel->
  "During evaluation of \
In[7]:=",ExpressionUUID->"14d17d28-d4df-4eb1-b528-44e9a4d8e4b2"],

Cell[BoxData[
 FormBox[
  StyleBox[
   TagBox[GridBox[{
      {"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Standard Error\"\>", "\<\"t\
\[Hyphen]Statistic\"\>", "\<\"P\[Hyphen]Value\"\>"},
      {"1", "42.34800601800603`", "0.7408754698894574`", 
       "57.159411721816866`", "4.8582345`*^-317"},
      {"x", "0.17999998799998798`", "0.0012803522451913995`", 
       "140.5863024616947`", "0.`"}
     },
     AutoDelete->False,
     GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
     GridBoxDividers->{
      "ColumnsIndexed" -> {2 -> GrayLevel[0.7]}, 
       "RowsIndexed" -> {2 -> GrayLevel[0.7]}},
     GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
     GridBoxSpacings->{
      "ColumnsIndexed" -> {2 -> 1}, "RowsIndexed" -> {2 -> 0.75}}],
    "Grid"], "DialogStyle",
   StripOnInput->False], TraditionalForm]], "Output",
 CellChangeTimes->{{3.733092106441432*^9, 3.733092114376828*^9}, {
   3.733095534079273*^9, 3.733095588293337*^9}, 3.733095691409423*^9, {
   3.733095841551009*^9, 3.733095867662918*^9}, 3.733095911439714*^9, 
   3.73309596097033*^9, 3.733096067376163*^9, 3.7330961608936*^9, 
   3.73309620250974*^9, 3.733096245069832*^9, {3.733097259191615*^9, 
   3.733097273378603*^9}, 3.7330987083844013`*^9, 3.733143503210012*^9, 
   3.73314678066078*^9},
 CellLabel->"Out[14]=",ExpressionUUID->"6fbad36a-c598-48e0-8271-fc579c6f6187"],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"0.17999998799998798`", " ", "x"}], "+", "42.34800601800603`"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{{3.733092106441432*^9, 3.733092114376828*^9}, {
   3.733095534079273*^9, 3.733095588293337*^9}, 3.733095691409423*^9, {
   3.733095841551009*^9, 3.733095867662918*^9}, 3.733095911439714*^9, 
   3.73309596097033*^9, 3.733096067376163*^9, 3.7330961608936*^9, 
   3.73309620250974*^9, 3.733096245069832*^9, {3.733097259191615*^9, 
   3.733097273378603*^9}, 3.7330987083844013`*^9, 3.733143503210012*^9, 
   3.73314678067282*^9},
 CellLabel->"Out[15]=",ExpressionUUID->"87052b97-38f4-40c5-87fc-778a4e6d04a7"],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{"0.7408754698894574`", ",", "0.0012803522451913995`"}], "}"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{{3.733092106441432*^9, 3.733092114376828*^9}, {
   3.733095534079273*^9, 3.733095588293337*^9}, 3.733095691409423*^9, {
   3.733095841551009*^9, 3.733095867662918*^9}, 3.733095911439714*^9, 
   3.73309596097033*^9, 3.733096067376163*^9, 3.7330961608936*^9, 
   3.73309620250974*^9, 3.733096245069832*^9, {3.733097259191615*^9, 
   3.733097273378603*^9}, 3.7330987083844013`*^9, 3.733143503210012*^9, 
   3.733146780685544*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"f3acf057-6316-4730-bef1-d7af6f6edf25"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Section 4: Develop a module that applies Fit[] to the log of sizeList \
\>", "Section",
 CellChangeTimes->{{3.731949496460588*^9, 3.7319495039435663`*^9}, {
  3.733097923686881*^9, 3.733097926942137*^9}, {3.733098354130898*^9, 
  3.733098357166833*^9}},ExpressionUUID->"0920a608-da80-466c-a68e-\
eeadb56acccc"],

Cell[TextData[StyleBox["In this section we will compute Log of sizeList from \
section 3 and a fitted model for logSizeList. We take out the first 50 data \
points because the power-law relationship may only hold from a large n \
values.",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.73309841640688*^9, 3.73309844180204*^9}, {
  3.733098629131445*^9, 3.733098631064431*^9}, {3.7331441435668573`*^9, 
  3.733144148131731*^9}, {3.733144293433251*^9, 3.7331443858658533`*^9}, {
  3.733144443887145*^9, 3.733144444246665*^9}, {3.7331595200941668`*^9, 
  3.733159524675722*^9}, {3.733159731298893*^9, 
  3.7331597793343077`*^9}},ExpressionUUID->"e9e64421-6314-477c-add1-\
b354663cf1ee"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"logSizeList", "=", 
   RowBox[{
    RowBox[{"Map", "[", 
     RowBox[{"Log", ",", "sizeList", ",", 
      RowBox[{"{", "2", "}"}]}], "]"}], "//", "N"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"logSizeList", "=", 
   RowBox[{"Take", "[", 
    RowBox[{"logSizeList", ",", 
     RowBox[{"-", "950"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot3", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{"logSizeList", ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Blue", "[", 
         RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"PointSize", "[", "0.01", "]"}]}], "}"}]}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<logn\>\"", ",", "\"\<logb\>\""}], "}"}]}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<log(n)  versus  log(b)\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "12"}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}]}], "]"}]}], ",", 
     RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"model2", "=", 
  RowBox[{"Fit", "[", 
   RowBox[{"logSizeList", ",", 
    RowBox[{"{", 
     RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot4", " ", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{"model2", ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "1", ",", "1000"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot3", ",", "plot4"}], "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.733098547137018*^9, 3.733098683030788*^9}, {
   3.7330987321939707`*^9, 3.733098867631537*^9}, {3.733098932305698*^9, 
   3.733098935061144*^9}, {3.7330990255743103`*^9, 3.733099031003772*^9}, {
   3.733099077065371*^9, 3.733099110421887*^9}, {3.733099726916074*^9, 
   3.73309972712746*^9}, {3.733099771008587*^9, 3.733099772761011*^9}, {
   3.7330999582277803`*^9, 3.733099981402972*^9}, {3.733100041523731*^9, 
   3.73310005070576*^9}, {3.733100380110415*^9, 3.733100389121875*^9}, {
   3.733100448618989*^9, 3.733100516430684*^9}, {3.733100670806221*^9, 
   3.733100692935993*^9}, {3.733100726034979*^9, 3.733100737605361*^9}, {
   3.733100877515451*^9, 3.7331008995329037`*^9}, {3.733100946910185*^9, 
   3.7331009508181067`*^9}, {3.733100982361478*^9, 3.733101027531411*^9}, {
   3.733101073685927*^9, 3.733101098265786*^9}, {3.733101204557225*^9, 
   3.733101258198453*^9}, {3.733101345184375*^9, 3.73310137098486*^9}, 
   3.733101454319676*^9, {3.7331015285877743`*^9, 3.7331015311799726`*^9}, {
   3.733101599765017*^9, 3.73310160074443*^9}, {3.733101679571824*^9, 
   3.73310168121236*^9}, {3.73310366339367*^9, 3.733103666599807*^9}, {
   3.733103973565669*^9, 3.7331040487644463`*^9}, {3.7331049445729628`*^9, 
   3.7331049448016977`*^9}, {3.733104993101305*^9, 3.733105020773637*^9}, {
   3.7331051159222307`*^9, 3.733105194038163*^9}, {3.733105348965473*^9, 
   3.733105350682754*^9}, {3.73314322570619*^9, 3.733143227027227*^9}, {
   3.733143315114744*^9, 3.733143340828122*^9}, {3.733143393703514*^9, 
   3.733143410703861*^9}, {3.733143452722068*^9, 3.7331434694525957`*^9}, {
   3.73314352174442*^9, 3.733143564948697*^9}, 3.733143605860244*^9, 
   3.733143640803145*^9, 3.733143678892907*^9, {3.733143729418879*^9, 
   3.7331438064128923`*^9}, {3.733146933674408*^9, 3.733146941517367*^9}, {
   3.733147390307851*^9, 3.733147456606257*^9}, 3.733158189244535*^9, {
   3.733158226618156*^9, 3.733158228307836*^9}},
 CellLabel->"In[90]:=",ExpressionUUID->"9670b60a-7a9c-422a-8f0f-6933a6b7ca80"],

Cell[BoxData[
 FormBox[
  RowBox[{
   RowBox[{"0.5794153502287143`", " ", "x"}], "+", "1.3377405307492307`"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{
  3.733098805918116*^9, {3.733098844148878*^9, 3.733098869397148*^9}, 
   3.733098953603833*^9, 3.733099729865653*^9, 3.733099774279573*^9, 
   3.7331000562684298`*^9, 3.7331003908160563`*^9, 3.7331004529560738`*^9, 
   3.73310050641639*^9, 3.733100698136389*^9, 3.733100739468399*^9, 
   3.733100953984741*^9, {3.73310100708687*^9, 3.7331010301580057`*^9}, {
   3.733101075952832*^9, 3.733101098964424*^9}, {3.733101228062871*^9, 
   3.733101261128573*^9}, 3.733101378374453*^9, 3.7331014581368113`*^9, 
   3.7331015320454473`*^9, 3.7331016016958103`*^9, 3.73310168234025*^9, {
   3.73310394981922*^9, 3.733103989761835*^9}, 3.733104945729154*^9, {
   3.7331050088554287`*^9, 3.73310502155128*^9}, 3.733105128640909*^9, {
   3.733105168756222*^9, 3.733105194876418*^9}, 3.733105351404951*^9, 
   3.733142657520768*^9, 3.733143227773417*^9, 3.733143414088873*^9, {
   3.733143456222357*^9, 3.733143470816703*^9}, {3.73314355313389*^9, 
   3.733143570238624*^9}, 3.73314360893747*^9, 3.733143654764827*^9, {
   3.733143738188299*^9, 3.733143807309558*^9}, {3.733144623315065*^9, 
   3.733144632371006*^9}, 3.733146781037992*^9, 3.733146944017826*^9, {
   3.733147395251775*^9, 3.733147459570654*^9}, 3.733158217471978*^9, 
   3.7331595064741697`*^9},
 CellLabel->"Out[93]=",ExpressionUUID->"b3f573d9-1e53-4a86-81de-6e9bb04ef31a"],

Cell[BoxData[
 FormBox[
  GraphicsBox[{{{}, {{}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.01], 
       AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJxNm3k8VU/cx5ElS3aXe+1ChRapiOSqiFYlInuRIlkiSyVLhbIkqYiIkBLZ
9+zRgpAlskRlabEkIfE8zx/PfH5/3df79f7Od+bMOWfOOTNzpY87H7ZjYmBg
KFjGwPB/v/rpmU/k4njojg1FL9qKOelPt8mcWl3DQ59ZUey/sJmLziVQ7eky
ykNnklKY53zKRf+1WGIawMZL2OCMwpvfVLBp7LOHQkq89PJ69Z+rh7jo+3Vf
1n3W5KU/1q6pe0ldQfdeUz9ndxB8Y6Pdx/vW4BVn+L8Fu/LSO1YEn5o7sIL+
7eOqK+sDwKIiQyORt8Dnnmpo5SShvMea21G3clB/6so93Zuq4dfyrZ6MbeGl
r3DbqNEWsIK+tak08eUAuJR2O754nJd+ICNcxy5/BV1a2z7EY5GXvlZz3ZGK
ryvorzpkp/9y8RGvdXTa/IAoHyn/db2Mytk1YPGnCx3mqmBztqRYUR0w/0hl
9dPDyLdt474gbmvwmgMrBXWc+Ej957KzQgx8+Oh83wUpExRu+i4/r3/KQXz0
R6+Sz3O6c9M5lz3zGYvio5/ZcbOc+yE33S4+lelCIjjg8/Gwzxl89IKJm+tp
jdz093Yz7KuL+egepQbb1Ga56RzJs5a6dXz0o8c2ZrrJ8NB1y0cd6C189CzR
jMLWvTz0Vq1bzMK98NauFn9fj/DRlxgOLVm7/+91cz1T2GKaj1xHvErKok1L
fOQ6emEuVyXByU+8+6m0FwcpYJeOX9U20uCYijRPIyV++v9fp0Nl//wUVeFn
H7H5ftbmJ9fl7N4vy/33gR3cZwsYj/KT+hXFKVR7G3iKUX5ptiNYuUR01aAH
P7muy+uD+mcv85Pr+ruAqd1MCHjw/pbNPVGIT1FXaUiLB0fZTUSZp/GT+8Bs
dwXH/HNwrUHlfb8SMKsRY+5EDfIr68YW72mEf/ioQSiig5/cR2G57Iov+vnJ
fWQeqHKsawR891MoY98kP7mvdDQ/X2ye5yf31Y2wFLlnywQIZ+4eP+zJJUDi
C11STdcJgQ07QkJbxQVI/lfeu7bZysPvU71cOLQO+fbmTu45pAo/esdN+5kW
eI0s2+c/u8HDa6a8lQ1Q3tSqWu+YCXjHOuvbbtYCZBzInqIl+JwSIOPAm4L8
dDcXATIO5Lqc+3XMS4CMA8cbZEpU/ATIOPDZYY36vyABMg4oTYRnFkYIkHHA
0YTt0Im7AuS+j5b7Z8qQIEDGBU+h+wwRqQJkHIj3a/fhyRQg48Cczbl//vkC
ZBzg5xQu+1oGbxmsO6pZCy9RfK4g6A2YTYL5SG0ruMnGZPLXB5S/9KCtjPIJ
3vz44JDSiAAZZ2Q9GaO2jMP/dKr9rTIjQMadY1M9W+T+gdcOi/pxMAsSjj3g
Pz/EIUjyGTKoNWXzwRusT5RxFxEk49aOKkMJRUlBMm75eCsOd8jBz8/uTPJU
grePZnJfoQJvvpzmd28r8mfPho+K0MEiVmWt4bqID6g/ajG/Dz7cxeOJmSHy
6+0J7c4xRfyFdu8VDNaIT/589fjOk/D5WidnL56Br1Tc8SPDDfmsOvMt27zA
/M/aT0/4gg+FKisyX/1P/WOdL7hvIL9YjwmdNxL5rzC/er38LuKDfGq9ZuMQ
zxvkatufBL/BIjCt/DF8fNE2s6hMsIaD2z3rPOTf1BrmKVcCNtGL4B2sQPyB
c7wX7tTBK0boNu14A15IZpYefofz/y6jPCSwA+Vzsq5JiHwUJM8lxXqR348+
gTPF22QUhtH+Ez1lZWnfBclzq5M1qFN8CtxnEh4U+keQPMekmXsHfi0IkudY
iYvSH0MmIfIc2/L9a28GmxB5jh0/sObxIhd4S1TtKX1+sJqD1sYwYSHyXEsr
7aS8FhMiz7WqlVdWM8iAT9v+9lq/CvHXGu9STZTAhnG3RX2UET9xqC0ieosQ
eS46M/699kRDiDwXU6LP8RXTwX+57bWqdMDuxdVSNXuEyHNTwUyt7sVBIfLc
vGSmtS3viBB5bkoNCUUnm4K3/a7qC7UED/nzr3Y5gfJ/XRMv7z+F/C0Zzb9k
nVB/qPC7WzOu8DHhDxyqz8MbfWAJCb4Az3Qv4a+eH/LvVGpvXnYV9e+KG+Av
DoFPd2lsORWO8s5h3lSBKPgvLHdmCu/C939L8zgah3xbaOzJk4mIfzM/cv1a
ClhSb/9WyhMh8h4R+5q/JDET/JTvp4h8LvK1rz5jk1oIX6P364F0mRB57wi8
UvPxbiV40e2PHHudEHkPefX+ylWPV2CrLZqsvY3g7+FNOVqtQuQ95XFPTXR8
B/h34dXimW7Eb5k6tnpvPzhS/PT32CFwrFQ+z9dhtIetamWc0nd4XyWf+2cn
4K/7WgpnTKO+xp0m1C+z4JD5d0+o/4TIe5P3TcF2fUYKeW9a1XfooQcLhbw3
VUjeF4lnh7d3/LqzcgX8LjlupX4+sKxry4c5IQp5rxp8+N6Il0Yh7113Ux8+
lpFA/AOPx90bZJDfz0NrWl0e/HYo+R9dgULaO751fnHHOjDfzQMM2hsRz754
m11jC9h9Pl1GWR3xk4WHDVdup5D+iPbTTObbQSH9mcjDRFvQAfuPSlcM6lNI
/8rkLI+t2w+f5M+V8+gQ8jmonuX2MwIXvzEtOGqK+PSEr88VLZCv8I0wy4I1
uO3Mx+IGW8SzTvB0RJ4C55v7WB89g/z+8e+sqC44vp76vO7Oc+hf+vGCzlue
4C3phy32XED/cLZNuy76gg0Oz4hmBSB+o3CCjfk1cFyznD7rdcSryz/vygiD
N2cwpBhEgl+U8C+fuI3zr2M3lBt6Dz7132NJ+TgwV7rqobIEcPl3wwMHk1G+
9EuOWH8qri+fy/NVjk8o5L09Nr1q5+9n8KWhN9MuZMMffzEwtZgHz+Iurni5
CKyUn2y0UAp+GVrteb4CfG9r6Z0f1cj306Kq0OYlhXwnOK7729/6Ciyd6C+o
3Yj4rmjrYxnvwA+47hYIvkd85G7hdT6dFPKd8Wv+Zd3Hbvjw85cDtvWh/LHm
FadjPsH7LVfx/fUZ5TM0Hr3cM0Ih3ymFlRy7E76Bs44wM0/+pJDvFjX6Sib6
FIV8t5he+qkT+huc8vFTW/ss4t+sOftcbIFCvmtiD871Wy+BT13oO5nEJEy+
c8a1zup/YgHHFf68Ic4uTOKNU+9vPMoFnx+bqh3OA1/JrVpezQ8vEeyQPi0k
TL6b3hYe4ZClgg1jloYMxITJd9O1SEf1C5LC5LupOjSCL1kGXCxm6tQgJ0y+
o15FZxz6thr+k6JdNacS8k1yG9WsWQ/+dU3HRGcj2I/nzxXLzeB/btL7PdSQ
z/bLubwQDfh4iegX97cLk+8+H/eVzk+1EW8wnfeqaBdY/mJnS81ulI9o473x
dg/KS9T2zbbuB19TLpToNPhP/ygvZ/hgCDaVUUvsMgZX1aQwdZiCq/c+Vmgx
R77lvhGir63g66qSuyqPgz1U5e3y7XC+DjQZ1KadAlubGs/ddcT5dXp0nO3a
WbBLcOykqyvie5ZES83c4UPVF0/u9AR7Zx+aXe0jTK5XirqUG9clsBvVqfPH
ZcQLCW9WbAwA34k46/LkKur7cY39ydVgcGjJeLflDcQr5TKzbglH/tMdcus4
I+G5vRSM+6Lg1b8N+WXdgVeJEs/2jUH+FdNxo3vj0H8CNdJrhRPAEzN2vgMP
wRmd4gNpj3A+TnPPHnFKg/+WF9+/4Qnyf7PNCZjKQP1bd7//33d6+IcycWIu
OfDpDvEiSvnwQgdOqn4tRH61ussXHpSAxRKdho+Uoz3nbmdeYq+Et9pbRi+r
xvXK4syp4lQHni3kPCrWgOvbY35V1qvXuB/dVV5qejSCQ6J3s0i+EybzHInR
y5fXt4LD/1jrObX/p3x3cgNfF/y/Wskbed3CZB7kfLbELaNeYTJPkk/92j/d
D/+FqfD8rUFhMk8ix95ksu4LeMnhRHDDsDCZN6kUfMpmMwaf9vR915/v4Kl7
nIuh48hfK+jpLTUFP75BfX/ONDhy6NyFHX+QP45Pi7VlDrwuJmfcYkGYzNO0
m33aOroIfqbX/dWNUYTM2xw9+Xh2fpkImacZ+njY1Z8VfGXp42FWdvBWY92E
EE4RMo+zN+26CSc3fA5Dkt8NXrBVlR+VQ0CEzAudShKWCxIC81QapjCJIL40
SSLuIg3ctNuS87cY4i0Vfk87SoLb91Yf+yQN9tn2cLuRLMr3rHRIqpcHS7+a
DFJbg3jd++yTaYrgrcWX+4XW4fgG5iSMAjagPM90puGPjYg/yjr90Xgz+vO1
ZdbUC1X4ZKaEO3Lq4M72S6+vb0O+yKRlN39uB0dJfRkx0Eb9IQttHdk7Rcg8
Gb3f14xPF/zSPdLLWU+EzJv5n8tb/3YPWDI5LGDVfhEyjxYwVOHmfxAcrvKN
sfuQCJlXGy7MUFY+An+v6sGyIGNwxnsT7x4TxGdnOt9eZwZeud/D2M8C8e9+
zVa+s0J7xFWC2iWPw+/8mR7tZIvy7jYjTCUnwROr6iVYTiOefUfD2AFHMFXU
8dRdJ7BFhFJMnzO4cfKar6wb6u/im6KedgdTtIfOZJwHj5o0ev30AnNWSmmt
vwAWcNlbdfYS+n9gvojh2WWcL8abrYuj/jifLExLZXJX4KveJG2zvgYOFJK6
HBOM/NtKOq62XP9P/++TPbw8DJ5NzPaLZgR82dZxulsk+iv+ktjJlCgRMi8q
c4Nu3BktQuZFFTqe8S+/B1+TlhenGgv/VSb5t12cCJk3ldSvk4h6AL9xyZZa
kQifyff662gSvPtW9WCBFBEyz7pBam5RI02EzLM22+rtO5EO35Ng7xryFL7L
/rZz5jP4UCMOvdYseJ3lyxams0XIPC2bQkEoJQ++7Jzp3y0FImTe1iBbSt+4
CNxQp+npXiJC5m11sgauRZaBOW+oej97AXaJOXuwoVKEzPvWaz/nGKxGvu0U
uYz5WrSnoPPXJv56cIS3durqV4i3k1Nm3P4G/u32z7qHG8H/Tvp42TUj3suP
765nCzglpy45uA1s5ZSdcK8dHHDmx420TrT/dOnN0/kfwMeLnm+t7hEh89TS
Ijb/GnvBUqO5+V39iF+6U3xi8JMImcdu/+G3/NsQ2FdcOGXqC8p/LAtSnxtG
ecHGj68WRxH/okziyLLviOd3PNLN+hOeQzHQgn1ChMyTX3yZ2ccxBb/W/4Ml
5zS4T5S5n2MG+Q52rrVmnwVv9TwyxDqP+N067qeXLYAZBYN/Lf4D52iHBMwt
obx+qbPwL0YqOZ71S8q535ZRyfVAlaozHmIBt9WLLetmQzyjumJxMzu4z/iT
dy0nldQ3Fb9Ot2gFlawDVGxglHzKA3Zo12GJ5wNbzw3PhQlQybrARG7Lv0tC
VLIuMNDwjc9JGH7xibiqGRX+se/+s3qi4INi1kWbxBH/h38rRUryP/FPSoI5
pMHDX+t4p2WoZJ1BNEQ966MslaxD/Ob/Y1srD25t6Nr4dDXi+fJKqZEKVLIO
4ZjiQzuvBN6lOr752Drwj+bJM5obwFPJppWSG6lkXcOwaGED4yZwxeesik+b
wa/adZ2qVFE/lee2euJW8K0fp1b7aoAvpUermmmCU01GT6lqoX73CeVyfm0w
u962zT92gPVp71te7kL9P9lqoxJ04cNj2i946oEfaLdcP7CHStZhXCx9q2T3
Ucm6i/K2ZKn5/VSyLtNNYUlvOoh4ObFTZkmHwF5OEWoehohvktbU3m2EfLG7
N54XOQpO/7qua8QEvMA4ZV90DOUTY7bKB5nDB2ZUCxhZgjXV9ZVlrKlkXejz
ruDLP22oZF3oTevWpeITVLIutGeCkn3FjkrWgS5ld0YfsAf/VVHJFD4NvpXY
P9/vQCXrSHNKty6mnaGSdaPL4ksqZ89SybpR0csP0ptc4OUi3uvOuSJf10Ri
cvk5eOEL42r+HuC5Nme2XZ7gv8XveFi9qWQdKpGl8nC9D+r7+2TsXdBFKlmX
2lG6cH23L9jy5EN/Vj+ww0JAfq0/eOf4nlUBgeBl7x92br9KJetaV2zX189f
A3f4xc7kByO+ckPISZfr8EopicIKofCJW69yDIWh/eke3Vr3I9B+fe6dRYcj
wZbD1u7sUWBOoXb3itvoD0qvVon7HfC+IrrOmnvo34+G5ym9MfC5t3zX3bwP
bg6fCdsRj/iYWJ/N0w/AyjtS5VIScb5r9glaGSXh+nlbbjnI/Aj+VIfE89wU
lL9H+dFgkwZ/ZbvNWp501P/iIH9f6RPE+4ff6bHPQPxYQKC8QCb8lnmnivIs
+N9sY6n22WiPNtuNXt5c+O/3h04U54ELc6+q2xSAzwcxHV9ehPbQ48a6M4vB
d2sbU46UgkOlNWrmynA+hlZWqTx4Aa8S9uOPdiU8u+t+3i9VOP9Lm50uBdWA
u+50aa2pw/UQl0wzef0SPrW6/KVDA66niuy9oRyvqWQddW3Q/vT0N1Syzlqe
s1dCrxHcvdj47UsTykdEmvEHvoO//cPnlmQr+Miv7DOlbcj/69PVh8bt8AdM
LihPdiDfohFF/EYX2mu0rdtWthvevFmfvbwH/HbLMIdRL/JxLKk4fO+DF6Fl
rwocAEcO9O2kDiL/XX7Niswh+NquPbE7v4Al6661d35Ffob3Z10dR+BZnPzO
Lo4i358Tao03v8E/n+MPl/kBz9Z7Izv3JzhLUUJt1wTiJc7tknk/CXbM9nI5
8Qv9t95408qpaSpZ57beOaHuN0Ml69zKTbJFK2bhRyv078fOUcm69wvjjs/y
f6lk3fuxrnZ0zgJYo4ThmeYi2NvwulLDEpWsiycpGQoeZqQRz/ah2baHiUby
d9u5S9gy00h75o7laX9noZH2FJUyNZ9jg89la305vxwcrnF7jT8HuKal5A8r
F9i2pkA5dAXyDQbPdvLy0Ej/NNImxm7zwsdMcTmK8KN95pwvTO4LoP3hg655
4kLwNV53fBMoNHK8TodjCqVEaGQfQKTXsNVDKo3sA/D3/uEpLQp/97Qww0Mx
8NjQu1lJCRrZJ7DEcfX4A0mUvxryQFNMmkb2DTBEBV2PkaGRfQO+n9P1KbI0
sm9Ap+7shVty8NfYtKS4V4GZ2APVQlYjX1VuUu0yBfhcX5a6S4rIZ3WWX/2P
Eo3sQwgV2iDrso5G9iF8jq4NHFkP/yBxs6m1MryM58ijzo3gp9+07Q9sAn+3
SU2s3Yzymc3eh9RVUb9HuaRPlhqYuvsvTVYd7CjmvOmeBsq/9u99yamJ/PUG
j1t9t9PIvojC2Q0mk1pgo8tdR09o08g+iiP+Sy3vd9DIvonZjaz1OrsQ/9x6
n2qBDuJTnJVXyu+GN7QUvBGth/KidYbOzHtoZJ/F+47979z20si+CnMb68yB
fWDFphH2AwfA/TK7hkoO0si+jNeMzeqrDsEHPPnCG3UY3mF+yG7REF7bSU79
tBF8qjhLeJsxOL2u10bTBPH7KtgLUk1pZJ8HL8vQdR4zsNfGmgFPc8TfvLKs
qN8C3rV7gXe3FXxdPtv3Z9bwezUu6gsep5F9IsYqOWt8TtDIPhEDDsawfluw
RkvFuV0nEb9uu07fY3sa2XfSljH8lus0eEzin5aLA8q/2Naj2uZII/tO0ow6
8zY7wftb6xbePYv8axYuas05g6edm/cfc0X8mozIwRI3sM8e4UWaO/h8eG6i
jweY7WTx6w/nka//9IWLal5o72Tv3oI73ogfzr7oNe2D9kb4W9UcughvGK1z
N/MSjeybSRm8PM1xmUb2yZznsuo56Ucj+2gST6jsq/aHPza8b7d4IHySzt+3
nleQL7snrLflKrxNgaG3YhB8/K37iVeCwU29uYa9IeCdBWO3Nt9A++O+hVmG
hcJvWRor+BwGjpA9GKcRAc7VZuO5dRPlH45b8YxE4vhzfFLiNaPglwZXlty6
DX75Ud5uOBr52t7NPdC4C+77OGwfcQ/xk20alYMx4Lfb6Y8338fx7zbXkgqO
Q//lSAeu7o6nkX1IqhpGVYoJYA2upbGLiTSyL6k861VK40P4vNl/k+LJyOfs
9rnZ6RH4l3WdVnkK6j9aNaLNlYb2r9xzp+PYY/h0R6HFx+nwJwyeFc08wfHs
633EuisDXjNnx8jNZ/85XwqNFr2Z4EWz6/ZrnqP8FLWSxSMb/GUqV60yB/Gb
Nz5l4MyD3zAwZmGUj/NlV11okFCA65/J3PjDSCGuf9mfHLPKxeC1XtRnPiVg
kaK62epSjB/NCoY9HOXgi1wsRodfIL7bROJkTAU49M0n7oFK8JmByIPy1Whf
T5vL6jM1aJ9AcmVMdi3ifbhqkmfq4MciinU16uE5nL8EXG5Ae6KfBZnUvEJ8
Bt+netY3YGnZjS36b1G+YUOpe2gjylfrluU0Nf3Hq3nc4H0H9oyXWzrUAk5o
5uKOakX54J+ni9rawLv6LiwKtIPHd/sPGHaAmWOLLKM6kS9+ubV3axd8ZVj+
Br5ujOepozNXDvbAMyXbnA/7CN7HpsLyphf5mBUfKy/vB8+N/vu7awD9ETns
Yuv/CZ4rQtOxfBC+SSCRd24IPs71/ZFNX3D+Rs9Q1Jy//uf8RT0sSR8GP1XP
/jg0gvJFgX4PxMfgfx/dw3D0GzzVUZ/p5nd4YY2iRw0/wFa8LcMM46hfcrHx
jdoE/JXBEQOXSeT7d9LMM23qP+ePwUqr7xfiH4wrpQv+Rj4J7cW8PTO4nwYb
OU/4/cH9tl8jPDd/Fuz/tzx1bA4cvb9tm+RfjA+dmqznDRcwvrClhh8O+ge+
u/ZRa8ki4j/p+878WAKf2n7khRSjKIkXErdQMGQCO71sVr+6TJTEN48Ojhcw
i5L2ZLE06o+wwLfY1+pR2eAvr136rr8cXL8zd4sPOzjHniL/hAN8Sd6y9AMn
8rX1VswsX4H2MJe5dKlyi5LxdltlivVJHviKkqiI27woT93vbFfNh/wWY2cG
x/nhjTsrWcUFwZUKCa36Qsh/qVFD9zxFlOw7fevz2iZJWJTsO71Y6SvXJCJK
9pnKLd2InKOKkn2oRx3kUmVFEd+x3PfEQTHkq89reOUtLkqeJ2Oiah+SJcCv
KxYjGyXB478PTM9IIR9l5BiDlAzy2dw+UqC3Ej640YniKgvm1uiSipFD+1xM
69or5eHNfE5uHlmFfF2137V41sD7bw75vVkB/vmDM8bmimhfzuP64wFKiJ9s
eU17vBb1rXVO9W1cB5/y9Hr41Hr43SFpe4SV4YcGlHI1NqK+o8f3NFipwJfd
Xn89cBPK3+/gnEvdDM9vyS30eguYc+vJvu+qYJv9O4x4tuL8CZm2eCuri5J9
wp5T2vsMNeCPBFe/dd8Gljzu/zdaE7wqMaOrYDvaU6DhZtupBdZnY0j4Qwfn
ZQYGC+9AeYNCBWnVnf/J5ydz2ngX4sV33LP20EH7xXiL2W7ron+2d2Zb5uzG
+Zi1zLJ7pwe/U79P7Kc+uI7f8QrnXuRb4g+IW70P9TtfPGinsx/9UWW+1Gdz
AH4yrpzd9yDYh7HgS4wB2uulIHgu/xB8R9SynHeHwa7RWUnfDME6J3bqsBqh
fMSKsRgpY/h5v6ZE9aPg6muix46YgHtq/9Y5maK82IuYz9eO4fisWyULEszA
z7cXqReZg+3Fw9zeWYDNTzZYjFiCGcejFpaswOEO1IPCNuhPWYcAo3XH0f9T
5rM8OifAFjcf+ZrZYnyR3PD8oasdxi8X3V3eQScRr1gWwhJvj/zfTdJ1c07B
a9a82VZ/GpztRx3pcUB+npkSvQlH+LP2XZbMTvBWszEKImdRP1fGqgxFZ7CF
XfrQdhfEU4d1Ww65wrO9X+Vh64b8Qlnnms+fQ7xXz6GBYHfw67jRtFgPxKsI
H5fPOI/jCzw4aFbuCX9MMVKvyQus4xr9vc8b+T41COmN+4DnS9TMli6gfVnZ
svI8l8Dhs0uPJXzBJWUTQ2svY/xnWSHeuc0PzwtK95OAvf7wgTMVX00D4Gco
ActOBcLrt1I+eFxB/scPk+0Dr6J9mQP7s25ew/FebVR/Hh8EZg6JcHgSjOu3
heLaVxACvmjBzFVzHfdjnPzZX003RMn/GKYl2m53h8KHdhyd+BKG+6OxUop1
Mhz5Ujcd7vgbAZ98hteGNRKeO+hyKu8t5OO4VpNKi0J83kGG47K34TUyTbvW
RoMZDX6xq95B+wrqJ2a07oqS/1mcD/VI0LsnSv5nMTj3aNmhGFHyP4udBkkr
TWPB8Xk3l6zvi5L/XbRtiLp7Kk6U/O+iMvftN+d4+MN81gvnH4DTmFxeXUpA
/Cp5ScMriaLkfxmJSxEx1x+CbWc/x91MAv8P9qeDnw==
        
        "]]}, {}}, {}, {}, {}, {}}, {{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJwVxX08lAcAB/AjIpXsKq8lVCy60JVCdj+9LKuOjVYWLuqQ8n6ct+Puea7y
Es8jN9uhmj5WXVdKL/O5ZJJKfazT5m3yUqtpjuUlCkVs/fH9fK0PxPiGajMY
DO7/Pv2o67fgxYxhTnGrnb7q8jRHFd4bvaD3LScnZ7ds6C4DihHt9DmN7zlf
2hfaSFna2Kpz+vG0Wgs1EcaeSh1dLDdv0YyUzEWFxjRj4roB3ticzR8oMcFX
o4aiBr2F8L+Q2dNfshIObxaOhfEtkduXSvWVrEPKmOG7J5Q9+hV7af2hrdA3
4VOX5a4oaBJbNnR6g2WZkJSeuxlKcUbwvp27oSzbNHzP0Qusc1cy9xd8B4uJ
4UQe1xvxD0IL3MU8KLKmhr/Q8kWUW+at5v4QVIuE5ROm36K3qWJQls2Hu5Jt
uNnGH3Wa04WcyTCYV5m/5I7vw5p6qd2NsxFIuJZULZ3Dw6QVy4FpEQn5thmd
GkEwqOpLgsnMaDj49HwkOkPQQZ/0a/s+Fmm7dD50hBzEUobj+RVP4jDldstx
cSMfq1q8BtuWC+C29OKhHv8wyL3U/t7xCWCdqhKkdYVDlcYWuSsT4TsvRDW6
LQKlTOPC3Ckh2LGbGVqNh+F6qCPMb2Mytj1QEUPukSi2tliTQaVgkeHh2YJH
UeCl1ZtEP0vFj6Obgmq4MWiW+VU1m4hgG9okTAyMhZFRkO/pI+ngNZRhj3Uc
pnd4lXUpMuCk4FywexoHH9szxi+eixGw+/VYYF48bq5V3flZLkFoqsujrq8F
KHIu+H2hCwF9s50uXL0EHDhg33XzLwJZzh9X6qsTEPT3ePLGLBJ809ppWzoR
FVtcRoQrpGhP4ia4bBVC07o2rFcthcM3yUfy9ZOQvL20khl5FCZlUTMdd5Lg
OmeL3ErvGALvCUTZgmQ8nypLCL56DJIgkzTSJQXdTxJZAbuOwyP2mHHk6xTc
adF9kTl2HO1T+ezSa6l4r5odF1OYCQPt8IH7EWkIHdhU4OqcBXkrN59rKcIv
fs5X1ndkQWd5002rlyLwin44/zopGwZ9VsIeeTr+NaoqarXJwcDn82+N8DJw
sZ1pebYhBw4bDF7xrcWoFXdl2IWfgFbNH+UP28S47c/veG6Ui4D9FqIqrgSG
nwUx/a7mQmvJjZjsixLIzE50KnbkoZXTvDdnFoEpv6SfnEbzILeqkziFEvg1
qoUjyKeg8DCe8K4lEL1PGe8soyDui259dZeA1XbJuaFCCv6FD6+L7hEgrVbP
jSyhoD+QHKWsJ+DVLG3jn6cQcarzpa6aQNNGdtSeagqrJkvV1U8JaGYVFLtq
KFxSOZTZjxIoGg5Xj/dTkPKPSu6+JbCj24NROUAhwKg7yH+MQHllX+jaUQrz
DlGmx98TiI/wdHaYphBp8ibvxQyBj43DD5csosESVibL55G4crv+Q4cxDV2b
BXscDUkEK06xisxodKvD2fULSNQRXrLFy2jkrTQbHGWSyFpXyjO0pzHYknbQ
x4yEq7Xw5OPVNOolLfjHnET//F33cxxpnFnNskxfQoLbO75Kbz0N76PP2pXL
SEw3qwMfbKBh67Sh0tOaREVtGS11ozHTScvabUiElKfUeXrQ+DNTExuzggSz
2OfdDIfGVban92xbEv8BjH86FA==
         "]]},
       Annotation[#, "Charting`Private`Tag$225710#1"]& ]}, {}, {}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{3.905032608257462, 3.3726288895002337`},
   DisplayFunction->Identity,
   Frame->{{True, True}, {True, True}},
   FrameLabel->{{
      FormBox["\"logb\"", TraditionalForm], None}, {
      FormBox["\"logn\"", TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   GridLines->{Automatic, Automatic},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImagePadding->All,
   ImageSize->500,
   Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
   PlotLabel->FormBox[
     StyleBox[
     "\"log(n)  versus  log(b)\"", FontFamily -> "Helvetica", FontSize -> 12, 
      FontWeight -> "Bold", StripOnInput -> False], TraditionalForm],
   PlotRange->{{3.9512437185814275`, 6.90875477931522}, {3.4657359027997265`, 
    5.327876168789581}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.05], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
 CellChangeTimes->{
  3.733098805918116*^9, {3.733098844148878*^9, 3.733098869397148*^9}, 
   3.733098953603833*^9, 3.733099729865653*^9, 3.733099774279573*^9, 
   3.7331000562684298`*^9, 3.7331003908160563`*^9, 3.7331004529560738`*^9, 
   3.73310050641639*^9, 3.733100698136389*^9, 3.733100739468399*^9, 
   3.733100953984741*^9, {3.73310100708687*^9, 3.7331010301580057`*^9}, {
   3.733101075952832*^9, 3.733101098964424*^9}, {3.733101228062871*^9, 
   3.733101261128573*^9}, 3.733101378374453*^9, 3.7331014581368113`*^9, 
   3.7331015320454473`*^9, 3.7331016016958103`*^9, 3.73310168234025*^9, {
   3.73310394981922*^9, 3.733103989761835*^9}, 3.733104945729154*^9, {
   3.7331050088554287`*^9, 3.73310502155128*^9}, 3.733105128640909*^9, {
   3.733105168756222*^9, 3.733105194876418*^9}, 3.733105351404951*^9, 
   3.733142657520768*^9, 3.733143227773417*^9, 3.733143414088873*^9, {
   3.733143456222357*^9, 3.733143470816703*^9}, {3.73314355313389*^9, 
   3.733143570238624*^9}, 3.73314360893747*^9, 3.733143654764827*^9, {
   3.733143738188299*^9, 3.733143807309558*^9}, {3.733144623315065*^9, 
   3.733144632371006*^9}, 3.733146781037992*^9, 3.733146944017826*^9, {
   3.733147395251775*^9, 3.733147459570654*^9}, 3.733158217471978*^9, 
   3.733159506529107*^9},
 CellLabel->"Out[95]=",ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztfWu4W8WRoLjve21sQkhmspMAmcxkdh7JspkMk00g5DUJszu7M9833+zu
L14GbZjEIptklH1ECwaZQMRLJERJUIYgHkomERkULD9kQBhhEPYFhLFsZBAg
QCCDbGQf27Itu7aru6tPS2rdq3t1dc/h+6zvu9JV65w+VdXVVdXVVdX/+eJv
f/WylRd/+/IVF5/5lW9e/I2vXr7iW2d++Ypvsqbhkzyekyz293dnevB/8Hjo
jb8O4Jv8spT/PwzNRgOsRpN/q2i/T/H/x9hbA7JRP3i9XvEXTEIN4HXt0jH+
P97QBMuyoAny1ShBNBSDkoVfLEgGvBArWPqtw/z/90IhFoBYvg6NchYC+BQf
Ps0P6VKdX1fLJ8HH2/38M5wpg1WIQSCaF4+y8hDwx4Fd7XmpS//R/GsQZ/dG
p2v8t3o+zp4RhVqjAH6vD7JVBLwGCQamL5bn/fv1/n2xmfvf+jSD3Q+FhkS/
WeT95vA5vgRQcynBcIgVRP8x2X+jh/4RfkaXcJaPE1RzUdZvnI2egD/XAj/2
L/CrYvflFPu/W/+nyv4Z/avTEPZL+vtCMF0VUFdycU53n88n6J8qsdFlY8LG
w+8TnBFNl17u6HpYdY3f9s7h9wPO/v+XEjyNqT3n62110fY5amuyeWRxannG
ZFuD3SemFnjOM9xruo760yaSZ7SP674g0DmLfyyBRjEOAcYc/MWY2scYJ4/z
s16AEE1xbwjYgMhHdL+nXkwzxqPJGoBs2fJ8lt8zBcV0hPfFJ20gBmVOGQti
Xi4D8IGQDPmG8CIGcaIg5iTefI7swtB9axf1IkR8Xuoili13dIFYhX2ElR8y
JWuuUMzaRXcoBP08H+Mfy6FZL8N0bhpy7C8TD4EvGIN8GRm/AbVanY1iE7Ih
LwSSJX5/KRVkpEuh3JDQ9NKHBQkGbSRX5X3UpqNq3lMfwOWFF3IS2emoD7xS
zlm1GliMc6rZMHj9SS6zhOwIQqmp9dH9OaqPYpxpjUgO5AXENlofTUgHEd8i
v6Sej7E+mDybExyz99ENjs/rwzMGVjHOBzmWSkMiEmCPDUM6Xx3iP49zVokx
aZery+Gt55jk0ynbtYth+TOqWvWyEBCOzLnqZ6U4oY70FGM/JH/Ox/yM3MSg
NYgyZpy2tLuNndPd0GyC1j272wvJcuMchZ3FdWOEywLsP6P4ZhYA5tJFFyha
RgKFWAPSbCKEMoy5ygmmJAueYW0Y4n5tGBiqfl9cG4bZ7meTLIxCJQypdEpM
7XCWi+NzFJQViKDsCschnRKKjyldGkZUWDYdkCsCrXTo/gSiQyUT4cIkzhgl
HmbTzxsBptQVADi7/ELqxRJxCDIAAvF8TwDMeneXZ8sx+Dj/GGFIlPiNkWiQ
yf68TT4xyEJ4AlcNON+Q/Oeabx2Ro1JKsgkRTDDRlYVcocLpca42YmgdBhMZ
yGWzUKhYnhGFrY+bM/zVFKZOnvje3O2w/A0tFJQYGfZbNlfQoBwTjMSsmVSx
AqV8HgrFElRr1vDMTz2nl7vND5aK+N/xD6ajSykuK/zxgodIhDPMF8kyYqDN
zQFQdoPpem7hMau3zCRdtVKBSo0U8CgUmNTzx3JMK9Sgwn6rWU0PISekZBQq
TGRUssKKZCCeO0O3dKtQHgHIlatQq7LfKjUkOVG1yGxbbyDJ5W8xyf73p/Dn
oV4ePOvdXZ4tOfcsSahiIoAzXn6T/NUoQ8wvlbg/itaIMqCM11vMGlJKn02f
REER1irKpYj8SxTq9n184ALytwBkyvaAmPpUVBV2iL3ICiSQCPb0qGk/ByFb
afT2yK73dX2kJCZn+3lazn8gJg3r9Bl8O0lr5VbSMklqfKp8eZaeaDvR1mfb
BL6fdILlTrQtWtsJljvR5gqWa7voZEPb+YvQ9iUjcPya/y4v37p1a0cXC9X2
BfW0GWizGICY6XACkBOAvFsBoSXnkaNN+JeNO+Gaf97MP9l3z9+6DUgE7IJV
a9Qf++4+IJGCOpDsu/uAdA0l3ycvsrStDdruaLI22sIcmQEZ1/Due/nHFIBV
hAh3S4h1eDQjPP+FGGvjLok6pCN+2lsIxLIgfavSk+gCYSFwWQbFuM/2stfR
/+znm0QNqwZ1qyl942I3Ev0SuH8SzddpF9AleKBfxHZQQ11sdk1btg+K/az5
6WE64oNAqsTRGMXu/is+hj0skUjwB+K38xe47XP4PqoBfxr/QD5vQCosNusD
yYLaeKhkeFM4noJUPMz/z1SbJtqP83fzXPkbJ8fFqpagWBHsX8lGwOvn/qoh
CS9utmPgQDybg2wmBxWraZQFhJ9JsDmKXyUdBNwuQbd+JRPmEQOW5vtEx7vP
C7HpMtRqVaiUK/jzTENoUjAOoojSutrhjiTZXEyGZLv4rWBht++C8TNpT9dN
IleblicAaQdkSr7r/BPfsAN+saEAq372GFwdZX/sU/LUl9wF97A+OZnownfR
Rtd9wW0Aa6JSAYxtrgH4nnXbYf+Bw+8GCnO7hAPn/d56+NdHitA4fFS4rTwn
z8rS1KaLzOWyX5M4/ZhjiPKBYfoLMf36DWnYkHsJjjaPLZPAmhTU8hl+WwxE
uvjHJvU5AU9ur8B3fpjhgH0r/BA89uxrcPz4cemMNJsVy2f4zbkRYrNFRwxf
x44dh8xTZbjipgc5gP/3x5vgmWJ1iYsHrSdBRpb+YTZF1jz2Ilx+/QYO5Oo7
NsOu8p5J/vP8Zt+IJGWPyuhP3UKSUUmSA4eOwK8e3AmXrl7LSXJzfAu8Wt3X
qyQVBvhI+/APYcOw4fo/dgv+4xL/vfsOwc8f2AYXXZ2CCxnwP/nNM/D2Owcn
5C2zqT7Cv21qD2HDsOF61+A/KfF/s2bBbb9+igN+SXAt3M1U6b4Dh8fkLZ3j
P8alQzuLX3HTxuF3x8AvkYi/XHkHvn93Tqni32SKcOjwUTGfJ9R8Ns1/RoAJ
juIEp4Wc2iYF7JrpfrLEenvpbbjq9izH+mshpZYn5S2zsftZhuc7gs9yic+W
whtKHX+TqeNs/jU4dvz4uLxlNvHlGnxOkfigBn5E18CRTfD081UxHXvWMu1M
rF+ucadQi0s6eH1VNNvyqfXs+V23UOxUSTHU6anNtk4PMp1eLO851cAB3BqV
tor+Qltl0sAd73cLrqdJXNuV9U1CWZ8ib9FnL+GqzV7+Qptz0jCzXYPr+yWu
e/c34OdrtsHF19iK+a29BylwxDSuJIyRTPjJvktcR8F348YOXj7NMZw5yr8r
MUUV/KPE07YKXstV8FL5s2kSM8yEkl3CkdINb0K0fRIzDc3I0Mnkp/aIyABW
Vq0D/28kvqiYb7jnSamY18F9Dz+PirnbIorpbo/J0iL3qD4pTNctBv49Mf4H
Jf6Fl5iKjtoqet0TJVTRJ0v82xdabGw9JoOL8Ncniuk61+B/usR/6w6m0m+z
Vfqjz7yKKl3YbZNdlZU+sUfncK1r8D9D4s9NgKfL8I2bhQnwf5gJ8NTON8nB
YFpoz7RAdw1+H5b4ocJe+3gJLv/+BvKCwPOv1KYkDiYnyZIZfnMNfh+R+KH2
+fVDz8Ol167jQN5475OopCdmHb/xjt+01YXnFANcjuD5hxJPVNCx1HNKQf/4
vqdRQY9JPE3uyCWyO/oNFbRJJi0zwOYIrn8kcUUVHZEqGvFFvOtW4/0SVxNf
uma8/lji8PIb9Ra9mhB69TSJg4kvXYPDn0ocdjDduEoaPFw3Ps5146kSBxPP
uQYHyqSc3vkm/K/bHhH67Ral3/7QLXB+XMKJeghh0/UQg/09M9DaNfP2LIkD
ArmuU9fM5LR2DgeOwick5AcbR/gMJS2CM7f8Zp2sQJPEWQzIe7LkPylxeKdN
Q6AE3b33AFly7dQ3LdGmeoRr4Bx1tsSpuucA13RtmoAcaLpmM622JgzPdwSf
T0l8XmFa4cZ7t3B8LmO8hpbLwcZRk0PQtHpyDT6flvjseLnGzX3SEGs3v4ga
wuQQJHyQ53DcXIUPpdI+pWmL/3nLg3xlwLSFWOGMtKxwrrhxo8mYGjeA4AhK
lMfKwG9RLP/7R4/gom+mwKzJHuEYOA7nSRy4YmEL839sVSyjkuIzGfKuGY/P
S1xQ1dyXKbarGnIgzGasjxqe7wg+X5T4oNpBHxqpnR/9+ilUO5QSbOIv14zJ
X/H3CdjN1Ax6PEnNoCe0Lurn2Dtskx0eT80Y45ey53v+oxN4fEXi8cqbde6l
JvWC3ut9DA/JUAKP8Q4nJvGXY/D/NX9vVSc413GngX0nuHHjsI2XFhluDvZ/
4g0juFPEpSmC8j+uW29DuW7RqNuTgfg3/H2UawLcsSNNcNE1KSliRlt2VheX
pG3A/hfZhoDgdmk7kPqccwzIv5Nt964rtG2LCGEtO5R87QhwtrFy5e1ZFAak
TP7fTx91Fr6/5++dWvvCq9dA6I6HUPB6SHuYlr2fdw7oJS0ARe/Pc3ejLm2Z
Mc9vwLWtSX2f5wTwIpZ4DO7fsJkv90i8YnCWEq+tKw1H4PwH/j7ZQuRb7nqI
78H5f5hRG7EXXW0vJ3Rbz0GgW9c6Swwjf66TwOkLyyVSmq6SC399bn1m0EB+
SF7UbMycJtedqqOoazv8SAMHXOT84YK1AQl/EPKywpmdHsf+T4YoPS6YKKi+
TOlxhB+xur2vOsYNCH3vTG5DD6EkHzHM1f8waNyFw1cUuCsXsOBfDqZzGQj5
fBBL56HetFProIrpXUEoIn1k8a1UmSR6KxH+G383OwYFU7Quub8eSg9LXYtW
+sobNnLRwC751OKQwK455g3GIJ1OQMDrg3A8DdWGXVDNysfAR9mHgGl5fpVd
yAPS/4FbLqPwRG4LXHPb/fCdWzfArXc/jN8Zalt5ih1+J8Yfl4N+c+whDiVe
cba8bmBpfJQoJSodhqHKBj/pDyDjeyjCl1eNjUwrTPPRwFzyKGktq0snCvW7
9ueP86Vle1SglFQLHUrWMeCUxSgKX0YghHPeAh1zqu3IRQKmw80dc128UZzQ
JcFUh9NMVyQfXyzUeQJnUFQaLdioY0ItL1FcaUCjkqXk2rmh3hktQJIOs0tQ
qTMyKPQvtIWDx2StDJwklPGNVd8C3gAv2UkSqoy142S2YDRbNmkzEx26RZhc
tnqt9GWPqYXOgONc552QRqVdEEbpVFDhYajF0KRHnL535+Nw+/15IzsvdKjy
vJGhDR2TEx6Z9K6122HljXZYF3ojyG5yHRLdPO/tKkUfiT9yFxLjRhGB39mn
kpa4jUXbCa5BgEId9VGYkACjYiOAP+ougFvj8FfdnsUwRi6MhVDeSNYnepAe
2voKblp7Fnonfd44/B5/79zovJwtFzDQE//H4PUnnqt4KFPf5F/4kKFvBxEy
S1daCmPoIoZpdPM3/J6hX0eQ+aBs65SuIlr41w/t5EY99wSp38bhsmvXci7E
a/iiSNuu/p1FQa5rxTBCqTM4ulXWThnG5X2G/gbuCO6OROs6i6luZXldvaYD
DzRQ5bprwfGYN3+dYRiMUwyEf6/bAEbCy/wTxT1IYSlyjekDrkHiTAPV28WQ
Lppu+9VTPE6g3SwhhegyxE7paoLQPolBgUDjyFG46d4tKMYWIwR83jiK6t9j
yiVMfyiGn31ht2vGgqqUm5gMVxaYwiBGy2UA61N7SOo6JO5FUqiGf7kVY11l
ep0LACevSielR+G7P97E1+GYuos7H3v3N9T6zzWAd7pTx/kIvMRY5Af/Ms2R
QtMpztTXgUNHVIi04wh8wkD5IbluIOGxq7yH59txo/b6DXD/pl0Y3r/Q0Tp9
49DJ9lP8O8lMfD334ls82/ICGa6zccvLuJpY6MiQeePy513HY6ojofH48eOQ
217h6TO0yHh82+tjhm4X0TolDDpHY1gfDQ+59Y42j/E1nVTPPG95e+ntEUPX
jhiqn+w6Iq3qdUTi0zh8FH6bfUEFjeC6m2mKIUPXjnAY4TPb+Ag2GsccRfhl
eodaykrlgX32FQQxbwT+oscBodoPb79zEP75t88qJXLdzx7EehD9xCD0DbuJ
+Lpfj+o2vLZ7H/zwV7b+uHc91x9fcQL2sw10N9mZBPtv05t56RW89qvXYxUk
rjb+yknYke7o7Ned/PoCR8A+oe7l6uLHpC42QvrJl1BdfNEJHOgwvtn4flK+
a/cqVfEtqSowqmHzs6/1VSKrT0QmDDm+6/R1jsm+oK1OVBoPT7/Ca07hrVdy
pfHW5wy3ODIuI13lUeu9utLAqBhSGjhvmO34WefwGWlZYqogzomWUTIEpJr0
3Jh8389UCV5GquTmX2yBN962znECy093HbVRU5yqx4TXhMQLNcwdD2xTGuan
//oMaphPG25ZNLx0DTMhsdOH7onclqUS/Nd372dKRhQHWsGUDHqArENH+olt
mDf4FE+jDwv5/zA0ga4j2He9ugdWMwvLVjJFVDJnG7oekL1LEOsEn5KM5A+v
73DzE+RYo+e7mmrBoFOmWmZKrO1nt2/e9i/NTn1ETubvp+kuNKVakQa609xU
wAbtSbZSUSsX5Lqb7t0yJGutUTcrmXTXIkgXas9w3swpSNG5W+W7Mc3nlvxs
LzIjcr/ttl8w9HTlvZJR60eiOJexHJMjuJKJow87bd5SoBj7VIjx1CTJ/tSm
e1gXFAmxeARoVKFYKOBBfwU8R65QKPGYMLzmA3PEU89/eZ9hM2WpgQn72Ubt
hlOjjAf++lB4+P2szecLgTxnsmechIk9pqYkTscVwTXEn2dIddAW1t7PhmoH
NsKwwQP/ypCORyESiUAsmYNa094bnA2N9gfhPPp9wzB8ZCEBV4d8WgVeS90f
jvETTCMBjLThpcllDZteERhRbVJFUCAHDgtKRr22jRa9oooarWTDRjW9MLlW
zzrV5cRA6NAeVMiPffbxQ2/nTQc5mO2Zs+17ybT3THRo93/h9dQHapTB0gGP
bgim1QH2ABUIe+dDhw+20GE+Gz1odqPeQX8++vWRkdBo61InbqYaSWgzYSW6
gdEN6RT1imC+cDgkTg0NpTGAcY50a20zbYITMwR+8ihfueMVH15IpEYkUs3a
drjnluvgqsAquOWOdbCbNdcKOchXGj0L6ZmQ0t0Uo1Ja06IJR7H0+t4ze+xy
DoM1KUMsvRCIxCGNQg/PZQmmoMBDDXlwtSoNOB/8TjcMGkWfoxX0y407mKIS
eGK6+htv7z+9x67nhGc9H2VCLQ/2q85P8i408YAJMakXAk8aRzSUdJdArX4Q
7lzzHLMAL+QpJVhqYM++Q/yufiJnuiLcrDAp4Y/xM4fx1ahm2UBHoFwvQcir
Idx1mWN6HPFgpzdzsmXj5fW39qvSd2jq37X2ObAOHlHzf0HWMyTlGqUEhcXi
klccGhSfphNtVPDOfMaVhIlpdwY1t4xO4Re/8Npe7vvnK9Pr1vMCd4ePNPup
rNlleNl7swH1eh3/2Dfxbtkn9PTFy4Qz4iVTTCUvd+5IYR7td5nwRZxRKWHu
fPPYsX6qTnbayaPil0YNbX8oFotQKpXYXew//l4Ga45rANNzf1+2UYlWXWl2
loYe4x7W6R1vwrdvFfkT//SDh7H8QT9xSp2oTwnMS3yJAH51kLYfgsGgOpC7
0IB+NCthjsw8k5+ffJdsgGHT06/y1Tpej4meTPv2E8bR3ZooMznWYoVVISLj
4vvHeUnLigkjAGlVf9Xtj6qYVApBU77bI6LSHaWHoeXFFPR7BoE/WuPeYEY7
nmu+VmhrG6394sI10eLDxja6jnDef/AwryFjK+snUVn3E1DSfcybzIL0STYP
BQXLh/u3IAlnfYYTzthG15H3ulY/xOsA4fIM/1Bps7Z+YlG6GJgiBaZSLEAu
m4NCqUZ24UKjiezdWhNnDCpMWVMRPBxeRJkp634iVubthJoJH108CRv5VD4s
6JjCBRK+Zltcvch0NEZykY5OCB3dT2RL98zTWQ5onBv+ox3+VrQxaQGEuOLL
XiF36uqdTFcHpK5G2xsLqDFR3k8kzDySVwdytiP5yU071bhe5oHtP3tMkmYJ
d7AgGTXpLvT5Tluf4yfTcp5+omv6SnDtcnakp9f8t5kI1RnpPywd2xu1DaRJ
Hh+M1+H8QrWP9Z1Wims4L217Ybfi7QUkUG/prwt2KuWfyDZMk70p9iBPCzhF
UgB3dOR1/Pfbf/UIXLpapA/4wxtgzcbHGUC9PWfB0mHncLRlr5zxZ7JNn0LL
JGeglKUtqSE5QExD8BUOGQRYXQo1ydY5RAnNLRe29+Mu54qzPhtIouqZoJRR
hZqGtCcxIRoJuLwlIyF4+4PY1muw0ZxzYns/EbNXKtBxToiZTAlWwhMFJM59
JheM2SUEFo48lU9EjrhzDbcheg1bmkMabO+HZvaKP+X06ps/7WE2qFdpjxdt
KN3eYDqZXY8++nGuTLD4yffufIL3hUsETDhiSwYJf/fUN6mIei2ysKA2lqDB
ycrRjEoT8UQns8kxrR1YoDYJxebtWo1/lsAtdz3M733+Fc32uEGUb0R9QjYN
u0/thF7B9Iukrwq1WVRi6BEVcm0odzgn4R+vXwv66+nn3+Qrf7IYbrznSZVI
Q0YZ9tZrfMKCIkKVYHTOntDkGy1syfTCEcG6XSu12U7mk+NIIEtcIc2TKelC
1laqarQYe9F2tL6Q/YseH7egsFM4F8J297oCd/5Oapux9GLT6mTDQtQRmCmk
zhR12CkiR7jSQ0bCGhpaLISuP+nSK3iUiFqR/XmPEA0AvfGW1RRJX2Sxq6KP
GW2AEfm+MrSB2wFYVfNr398wqU0TuRDt+fSpBcXry7KtMyppjGsa/VUs11DT
T2oCl6Z4rwUuFhR2MhNMRUp0a0MtEe3yjEqe6df1ekDlHAO/ZgId6fdPt6yT
BBesgospTbGTpapPokltEmn6tte6G/OO8ZppfYarHBRWjKWNLkFyj+FsQVww
TN3LZgMGnLAZQ5v3F1+zpn3C9xwTtKDc1d6GhWF1+5ZYf6k2lPH1Bb4Xq4Kd
NMu/17CZgSLRLQ3fVJbDZLj/QQ/PGBAS46ptpnMd8LdV4kRDhYAmx5xEoLVN
LFU63XMrVq+14+euVUsG9wBOm2U0j23RNQyXXGNn3vcaJTJwgGmjXOdt0+mP
vUaALBLAbSUPrl03qgkVulQIn8kOLjLtQy3t8fEDx4+iHpCDvNfqO8UjPG9U
A3pCu5Sm9bjUFHc+sI0zIKp16WzsuWD/wHE8Q7ZpueuK6XStSIXWcADJi0oh
pHrpoF7dyYuGmGk2YegRJgnoiHWGg8+vXvyiIUZTCe2YcclwaNFjeh0mcOAL
r4mltvNCr1TE4zLbv+MepEim6VON5Agl0iISQ5gh9V75w1WaA6JXL/0C2MQz
xRhokYye5WQpBlNwV+o5GWG+nm9skflPjjZxd2/PGohVrLepHVa5E2LP9ZkP
odUXWr0iM3DOImRWacsSMnsxy67V4S2m1ZVSgpNfYMXqzkqSKPRaNw1chGxn
vbRRvh9MhaHwUKz8rt0q/KHdUpXOKFrp6YPOHmfa/1g01LgDQCqgcYkgMR6e
0bjxyZd5EPAFMnJn16t7aAsYsWr3d+i+A6dQo7QLvhMhgo5UmRNEAVfOuHNB
DhjalcDrsXbFClnD7LrYE6jWSAWvkqex2lw7ok9Rx9AVq2RzRhjlWu/dd4jn
hmP2ISq1HwilRlIeMaesAj3ERSYCGnfoBo7XX0rY0bLAwi5Ieyxnj0HA1qEj
JBF1/wipuK9el1I62RHBIXZlJnkJIBVAwmbRA9kXMICEX8R1M9NdplWfaSds
4ECLTMVxlGQtmYeYSI2+ftRN5D/TNyu+bOhr4MCeI4HFOIzN215XcRgyc51v
sxFl9U0JR4AVjn2Rlo4C1Q6FzMKtv9w6Lp17+v4Evkwz+jwn4CcvOKaho4Sk
UEY8uBLLH+CWF8kN3QtO8xH9g9SX6byBgSMg6hiM8gxzlIMr7LqbatmFgNP8
o81j0wh80tD/gMzgL8k2XH/czkQgbgjRYgrNYrqONBj6gy/QEtOlQWbaORm4
1UsTDcW3PHhTknpMN9P5C8nMJqmpKsK/75naCwi7aYeBTnnAuUvZdrgWvFTV
xR7mVU9dATgSlEKTTBOTQnVQn5JkdARg0TbF22iyoe/k5vgWVYwEfUO4ac4Q
ILh1V4pp/2yR4G5to+K9KOUf3GJLebbEVgTX92tdAziZgRipnnyUi3cjh7gG
YDLJMcwco6l0ZxQq+7f2HmzRSZfZuTeqCKLjSAif5RQvlkXHtV0gK4RgYAwF
UuKLTQmT69q0c+kILlMSI8zkulXWM1wh6xli6B++GOBj2niQqnUNDkv5+zhP
zNKDvnFGWIzRkKHGtFlBCPyZWxAQm3ljGGXbkmWF0VOXrk6NSQW2IijWgLjg
QwMBl8CmyiGO4LBc4oBm/RPPVZRZ/y32edHVa/T4CYwgwYtNu8SOwE4ZPbhW
enDrKzwuaonmUqEX7aBpGximPagB2ZSmNsrFQgWAFXIuCa5Vu6i4AkTZyVjF
M2xg/zMM/Q3cuOy+syqUwr3ruSEsfcaty1VTxrAj/EKZpjIKsmWrBCnsWoB1
l4sKRLMZ3G0Aj3XEqf6QqaiLg6kpAzu7BnbK6ULA5HLVGCpqOhbCEYAJEB4W
99NHWwDWEp8VYo4DfAZ/F/vtv0jv4LqGsQJ5jXSnhSnp3hGYaWMMJTV6Yi5c
1ZqDQGzsGoDpGXfbrgm1WUArUvR8MpP9A4bbO2AmKd+wLMrmUduusx3eOCjc
kIEuk7usJp+j6TyQDsSEDmaGdL0IEZW57YVYlldfme+BjguNM+0bofUlt7+U
Uxg3V4gBTYntXXBeBsW4H7xURKg2DT6vD3M4Zju8sZ80jrluJnVuAYoxpmXl
aYZ7u+DLVq/NppaiXecFcJLlRs8HNfaaqdZr20dlG8Zy4sG2F2hRdOhGoutO
NdzbJSNtDCqZCGA6UjyVhnjYx/6PQAVAIdnfGY39jOm/pWdcu65jTDFZj64z
nYjRZUxHoVFKch7N5LKQzRXYqOKtNq59nMrYP66j3CVC+K64Bp1q9s7T9M43
1PWmsgTd+bhRAL83ALlyFWrVClQqNeRrlXHZ34mM/aBNXgW9YAGl8+LiHhf5
9JohC9wYWNadHGzQkyGfEt3eQMLOtpv3gYwLQQYKXLK17xiXXrTVjZV20Xdz
z7rto5pRTEtsUxSaI9aEfoKJTKv10F584Ceb4GfJvNjjvnoNMbe+MDHFZzmC
B3lSfVpIEC2uMWEai79izBNJJTkQ/C5TaJajSMiqPy0bC8hcrZtpLgJYX78u
l1YbOutRAFzJuIjCenB0Nj2WE1tZTgEr82k5mFR4FF13RF1HKEl7RTgPaa1n
W77ikI6VoQ0tlV7p3i86AbB+Bgwu7Gx9IFJP9AxdzC6l7F3MVmnfAP68EwjQ
rjPq8VaJMdoRN8CkhZqIuCVMfZznBOAiK7HTByNzl5QW1IQIcY2+oDXDPhBv
6HmyTfIJB2ZIk2vURgyvx82eo6BaRBILgM1JiiLBb6O2Zlmj7vuskaYDBpbO
99Dr8xAl91m2mjvXSeD0sHwCDqupYinU48ePOwMcRVTglJB7y2pY8cwqTC7H
Kz7jBHAi8kA407AQPwpMNjEo5JWqk2It/l2v7vm0EyCeL0GkNj7IN6RpX4wn
Ja62i4u+vnu/51NugBNfbMzHpd2OfKjZubyQxtvvHDzbOVBHW0DFlx3XOM7n
0avVuqrwiayAm/b7DxwWWs0tQDONJOfThFKbeEDf9XflhJvye6JKZ+Pw0U84
B/dYB9zcCIs+Ni7dNprnkb+w6KaMwOfn4WDRzaPNY2c5h0Jrm54p8PMHWoM5
9Dgf8ipgFU2M8uTbwuGHeJEpJpM/7hZ8aJEp89a0Nfcw5/v2kkDHjh3n9TFx
oxWvx4DK6Z1vfsxt+HC1E9qgmc6tGR/k7MKsAax3SVkDeFbRrvIeqspkCqn8
iBHXgRh3pjZyb8gwxJbFAc4YFAY6uMSxFqMIr3K5eq3a6dKjKk5XWDg8hJSx
oZ+abaqEQIjRKe76Fh5aHDTUpxvHy0HE9NQqQgz1jhw1Dz855mSDUPmQ2zBB
Tw/Nsg9qiw4pOeTJsGPdbP0pt+BDNTJRatCcInwwe6/9Oi10WS229IGacAti
U7JNX78s1ya/zJzw6JVtaN7oopTcJmPa/HItsuQ2QZ/kMgk2FUrV3efsd7Kg
cNDRDWQ6By7wk036Hprr0F6qoY3wXhXNKqfEG29Z6jraycBJ2y4vvx6yq87Q
Xhqa7q5DlooyX6nZK3ac49qOVPQrf7pJXfd+zctE1427C7HRjni7O377rGdY
Q5B2OfS8Rpqdv2Nwo425BUHaejSZVLzsNA3PJZrTxzXQUySB7lKjXBzMtcS6
6rqtqa1phCB1ExKo0trdly+9/g6/CIeDJhRuwrkUic60Zfx+5wPbqPBGctML
PC6XXr4b0govfXBGjHgtmjX/J/x9mbEAJW7t3nTvFnutFdkEW3eI3W8UA1+7
XpVb6Vh6usGQF/su43yNtf6Jl3hqDOIRZGusYnkP/3VUKaANswyJAwiIVDSx
G3BfpqhWGqG7c7iM0nOVXAc71Vfcyxjl7rXbeaUQzv+uEalU/xRTSmR9XY9e
YMF1JBV50mP8NC5ZTJj/IXkPHT7qGgn5GQknzrGrZB0ArKOCZzIfbR4bdguc
Irl2Ap5+vgrfuU34yDBJN/NUGY5J0e08lCIdepx7vrCA7DduFkdPYBmNLYU3
eGSjwyCKjGdhWqU2v0gZrG4Ajc5kP9g4wpM1XACSIzEI70JAOBx0IIwO0vmL
0NZTUo5TwJnpxYOclxkuX3qi7URbn21mB8wJljvRNrA2M8txk+wZrZW3Hej8
fxf/GOXnmuczKUgk4hCLJyBXrHrsox6rkE2lIF9tKI8dxiVPp1MwXbFUJZR+
+viNAKeoLi9DWMQ1C3zZTcFEXmU4FJMB1uqD6TqI3nD5QceV8BOzP7dgPd0n
QNsuO2xaNahWq1CtWXJzawmAVYJELAbxeJqfd9qo5vgxMV5/HIP2WUORJ5h4
I1kV/m0VxKGh8XxdlRxmnUON9V2rW/K6EX7QZq2Gz6zhKZPDCow61EVsk32E
W+ulyxYIYtH9OFj8uE8VTtWCgSTSM5IViglxaF6iiJ7tJhSSYUn9D4AvlICK
JH+jkoNwQIS3+0MRCPrUiSpQm47x9th0TYyU7DNZtFq+J/BiAjfMDyq0iVfL
8dOFMWkkV5NZBJNQzsbV0ZX80eEUVNl9VlFgFEwUJKUmoNJ2LTJLPFfhIFRy
sSHJU95AUBBQBOjTbUOK8dRd0ik8E5HeM08iWZJv2XV+n4dAC8Ry4vTKViIl
9CH7sIImXmCjnI9LXAMQTyUhIp/ui07bveBTgiEIqMSyAIcoHw8QuVuGCfvt
/G5Bgp9yH4Zyk4bto1BK2QfshFIlfl+zKk6c8vqCkEhnIBaUFEkU1bgF4pS5
VqFJD7FMHgq5hBxDzBeqytPBWd+xBESDembFTPfB+EBJRaQQgxuPRyXnYl5P
rYNUM4xegXfL04EkAxUEazIsnt8mwQ2lJE8UJIgIEj1CgDf76DUhw7OwxPFI
gjxYO6EGcb82Y5hMQQRqOZGs4o9mod6woFouwjQeP1esGkbQgnIxD9lMFqan
2WiH5SmciP6+gpJVHMxGXpt6M9zXJBgHRS8LknIIk0UhJ2l++jlerfTqOoRv
yctINOggBeHZ5ySpkiVFqmRg9iFMtMlNe0gDdL0kzzg0K2k5VBnISfKkmTwq
JAIenpiH51z7YvAGDVuiKOTlcjZ3w5qsFKI3nivblJMCsv07u09Kv2H7HjCN
2EKSh34LQpGkLU+6szlRJ88Mk47+D6dJXGXF/OVUSkoRExdnKdJvvNsGpEPm
EYtkhJqZjvlbRCaho086YuIIEyaNYlyhYBUkAwdTUGs2oZQWw+OLZGGfHD2v
PwKFqgXNZl1N3EgyB8VCjhlzSchMl5hpMdPwzXCfcQgXkl72pItkRF5zPh7U
Jl0rvVqGcKQFpGYlI1TmBJPW9tnWQgHZSsPr84nLxqlTKCWlaVJqEM94qKtA
wK/mgvi9BjG/QJRJxTE5bVBlRaR5wfWBYmkmuA4UIUR6JBBQkCHQTHDKafM+
hu82yIYl1NEYRKSK5HnQFgnNhDDwWr43u92nLKlBEIqPjCKUz+/TCGW1E6pl
4M7skGL1YgZCqgc/RFMFZeM3a3mml/3sCX6GXxyiYaaLQzEoN3B5kNaAByhn
YhAIRiGVjkOAXe/3a33VshxZfzxP065ZRdssAMFISpCVaYlMlHUfCEOOmZSN
ap5ZGOzRPh97PLNcZOp9o5yBIDMZI3G0QHwQSrKlj1WGZIT0lY/1mYASP4+1
zKwdfEZaTD22pImHtO/m++TQDYpUNK38EEswkvkRvwAkhNTWSdWyNBjWoZHh
FLi+qEOCDzW3FChqoizNwGy1qU5+67x9vO12W8Qzka4yBfW+ukDEtIlaiVlS
BDFG1NarVT5FfbFptS3Wee+Ydu9uSPHc3aB9gGZrN3LZ6znp/wNADXzI\
\>"],ExpressionUUID->"4c390618-0fde-4523-a57d-43dc3e557838"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Section 5 : Create a list of {log[c], d} and find their statistical parameters\
\>", "Section",
 CellChangeTimes->{{3.733144462325074*^9, 
  3.733144585047222*^9}},ExpressionUUID->"6c70c45b-21cc-46c9-b1bc-\
03386cc6c6e2"],

Cell[TextData[StyleBox["We generate sizeListList below as a list of 100 \
{log[c], d} vales, then fit a linear model to the values and calculate their \
mean and standard deviation values.",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.733158027007984*^9, 3.733158071269161*^9}, {
  3.7331594386238403`*^9, 3.733159477126528*^9}, {3.733159595070384*^9, 
  3.7331595955648823`*^9}},ExpressionUUID->"c97e94f9-b0c7-4336-aca0-\
f618bfdb2128"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"sizeListList", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"m", "=", "1"}], ",", 
   RowBox[{"m", "\[LessEqual]", "100"}], ",", 
   RowBox[{"m", "++"}], ",", 
   RowBox[{
    RowBox[{"oldCells", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"{", 
         RowBox[{"0", ",", "0"}], "}"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"1", ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"0", ",", "1"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"0", ",", 
           RowBox[{"-", "1"}]}], "}"}]}], "}"}]}], "}"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"sizeList", " ", "=", " ", 
     RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "\[LessEqual]", "1000"}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"newPic", "=", 
        RowBox[{"growCells", "[", "oldCells", "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"clusterArea", " ", "=", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"newPic", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"borderArea", " ", "=", " ", 
        RowBox[{"Length", "[", 
         RowBox[{"newPic", "[", 
          RowBox[{"[", "2", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"AppendTo", "[", 
        RowBox[{"sizeList", ",", 
         RowBox[{"{", 
          RowBox[{"clusterArea", ",", "borderArea"}], "}"}]}], "]"}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"oldCells", "=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"newPic", "[", 
           RowBox[{"[", "1", "]"}], "]"}], ",", 
          RowBox[{"newPic", "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";", " ", "\[IndentingNewLine]", 
    RowBox[{"logSizeList", "=", 
     RowBox[{
      RowBox[{"Map", "[", 
       RowBox[{"Log", ",", "sizeList", ",", 
        RowBox[{"{", "2", "}"}]}], "]"}], "//", "N"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"logSizeList", "=", 
     RowBox[{"Take", "[", 
      RowBox[{"logSizeList", ",", 
       RowBox[{"-", "950"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"model3", "=", 
     RowBox[{"Fit", "[", 
      RowBox[{"logSizeList", ",", 
       RowBox[{"{", 
        RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"d", " ", "=", 
     RowBox[{
      RowBox[{"model3", "[", 
       RowBox[{"[", "2", "]"}], "]"}], "/", "x"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"logc", "=", 
     RowBox[{"model3", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"sizeListList", ",", 
      RowBox[{"{", 
       RowBox[{"d", ",", "logc"}], "}"}]}], "]"}], ";"}]}], 
  "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]", 
 RowBox[{"sizeListList", "[", 
  RowBox[{"[", "1", "]"}], "]"}]}], "Input",
 CellChangeTimes->{{3.733145423904855*^9, 3.7331454419573507`*^9}, {
   3.73314547905374*^9, 3.733145513251049*^9}, {3.7331457205262423`*^9, 
   3.7331457917594624`*^9}, {3.733145823007183*^9, 3.733145885993135*^9}, {
   3.733145990294484*^9, 3.733145997400737*^9}, {3.733146073095085*^9, 
   3.7331460752201843`*^9}, 3.7331461066488113`*^9, {3.733146179815926*^9, 
   3.733146187502722*^9}, {3.733157395829384*^9, 3.73315741582018*^9}, {
   3.73315764395319*^9, 3.7331576461609907`*^9}, {3.733157730393697*^9, 
   3.733157859882216*^9}, 3.733157894922194*^9, 3.733158051263834*^9, {
   3.733158120908431*^9, 3.733158122508792*^9}},
 CellLabel->"In[23]:=",ExpressionUUID->"a32f40f2-4dc9-4745-9e47-686810c3c85c"],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{"0.5928797144668164`", ",", "1.3441558256070871`"}], "}"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{
  3.733157420533369*^9, {3.733157591879827*^9, 3.733157620559038*^9}, 
   3.733157658377184*^9, {3.733157849211298*^9, 3.733157874876786*^9}, 
   3.7331581536650667`*^9},
 CellLabel->"Out[25]=",ExpressionUUID->"b3abc91f-76b9-43dc-a5bd-3348a3e3604b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot4", "=", 
   RowBox[{"ListPlot", "[", 
    RowBox[{"sizeListList", ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Blue", "[", 
         RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", 
        RowBox[{"PointSize", "[", "0.01", "]"}]}], "}"}]}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<d\>\"", ",", "\"\<logc\>\""}], "}"}]}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotLabel", "\[Rule]", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<d  versus  log(c)\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Helvetica\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "12"}], ",", 
        RowBox[{"FontWeight", "\[Rule]", "\"\<Bold\>\""}]}], "]"}]}], ",", 
     RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "500"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"model4", "=", 
   RowBox[{"LinearModelFit", "[", 
    RowBox[{"sizeListList", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "x"}], "}"}], ",", "x"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot5", " ", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"model4", "[", "\"\<BestFit\>\"", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot4", ",", "plot5"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"model4", "[", "\"\<ParameterTable\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"model4", "[", "\"\<BestFit\>\"", "]"}], "\[IndentingNewLine]", 
 RowBox[{"model4", "[", "\"\<ParameterErrors\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Mean", "[", "sizeListList", "]"}], "\[IndentingNewLine]", 
 RowBox[{"StandardDeviation", "[", "sizeListList", "]"}]}], "Input",
 CellChangeTimes->{{3.733158107718308*^9, 3.733158316062747*^9}, {
  3.733159354139402*^9, 3.733159364375424*^9}, {3.733159408715228*^9, 
  3.7331594233176003`*^9}, {3.7331599130579157`*^9, 3.73315992407051*^9}},
 CellLabel->"In[96]:=",ExpressionUUID->"f8649f83-d5f8-4efc-b6f5-7ff69c5ca5dd"],

Cell[BoxData[
 FormBox[
  GraphicsBox[{{{}, {{}, 
      {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.01], 
       AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJwVkns81IkaxieXKApZOVMrKlmNo2KXUHpcilNCWAdRuyxyiXHUWndbTUWF
xakthckuoywVsbHaJepMrIjSqS0zv99vRJtlhrn8xqAz54/38/z3Ps/7Pt/1
UezAGA0GgxGtnv+r5G+axChN4ptDqoIbJ2SYLjesIsaFaHIw0mSmq9BRuCFg
oymFeHtq05bzs8jafGapvi4JP5pRMJQlx+B5n/B1bwkUFa/ID+YpEHpJ9uOB
e0J8KBYMr55TIVuE7OKPKYTmTJct95ZC0nzt1txTAgbSo71e62nkv/hvg6Y1
iXMHH0UkRSvwF+KVQTYi3Ld/NfxuSIxcsxpK44YQDxMeasfeV0HqebUwHyQw
nTZJ5ssRtqKuSDJHgSWw+XJaRwxbeeeo/QMh/BI2iI0HVbDNfNnruoSEtXvg
Onc2DXPdC917m4RgVchplygVIk2eeQ8sCpHKsRjR3zEPZDjphd4kUHmg6p3B
uAI6qprSXg0B4uND/K3TPmCj5r4fJGwCwUd3DP47U4nfi/lVQboiGOdqzzqH
ziDgU4stx30omNtETERzZuF1JbtyIIBC9MSLQg8zKbj7LfVsLUUQh36Xbeg1
g6DI2pXpTwhUryllENo0Qp1Y5+uDSEzO9NgneMpRt35z54tVBDQSd/ucTVah
OqcuRKedgLb/BV6KlRIH0vpPRvuJYFr1zVqTZAmajb+93nqCxOMmgXSJtxyC
vdE/H9tK4E+vsdLG20os0/CI7M4SIt0qeH1Y1gKS4mJ6EjkEgvzzTBh26nsM
42a5DhTELkUfCpKlqDjW6svKpzB5J6fxK/4sesrLQ9JaCPAKuIOv1tII9u1p
uztAoDbkg3V3nQId/+QFuBmr+Rio5DXtUeC09rUdlp0EmJtSI7T6adx3ypvq
cROhrrHKw+29BI+Xk329LRT0isSTnqkzMNDcasg4LsSiT9XrHEKFl9Zdb051
i8ApkwmZbtMoftqarCwlsffva7wsVsgQXfTkoiSDwk2qfE58ToauT0bICEcS
8XbFqzY3y/Hp0MUjXYcpjGdKTvt8PYshO8l21wgChb+cHuSk0YiTPLPVf0OA
pZFgMLaSRsAbs2F3fSG+rnEM1zNbxLpcXzZtSaHjsz0695hSBLc5hXCpMQzu
1HHt8nmH1oG7fXHv1X6BZuJLWjLw5Eltns9JTOfcea7rLMOR2HJOOkmgf0Hx
etCURkluaJvmJQIBOx81mNmr/83ksBO/JWEVnhFhUyODwO69RYW3CLtNLzoy
LWdQ7d+1llD3L7dWzK361xxO1Cfq3KglcL1KaBX7HyUCK6R8Q3cKTtf7z9B5
M3CXMV2SXAi4Rdr7Xrwyh9TOuqRqNV+8vw7LF/k0nDnbl1JZBMgMm2adu0rE
lqSk1Kv7NDlkodfYTqNxC+1h0UPi19Zw2f2fZFgo9/uthU3C12iabxckw+9j
Ji2j/QRuZsKopJaGF/OBVgxfiJOtucFTkfOY2rVxKnWSxJuJz3/ZIJJBZu+0
8nsGhbNJ7xfGvpch/eHw6OxjAttHA3zEHTTodJetVAqB+jzC/Nw7Gn+43Sn7
rlGIhC8iXp9NmEdrVBGbp+a9b1felphYGlxmStfWORJD1do+7iNSiD9mrVjT
KcQ93adWmdfmcbvdryacR8J8160SO3V+hx8rb75cSyLAJemn/SoFHvyx6Zp4
nxDMw6rDvgcXsHo67PbUUhKJ8b/5X7agIexjdTqYCCDLrX7q7bsI/u0FrU+M
SJy6craWy1PioJFTJLuCQPuysh8KvlJin8a2flcuif5dfs8mwuS44OFx3XNE
CL7pR6zjwfOIKn8lDhwm4bhOq3CbmRwzh14UP4khMbv4dvvmJTSiuM5liCIQ
zXZm8LRUuCrpfl6xmkTZP86Mm7kpUM1yyF2mKcKjo6/cd7ZJcKtuR1BQL4kM
t/kJ80kpPvN6e3L3ZQLj1uGOY8eUaHhW3FS5jULy/oi9BUelsI3Jbu/fQ2JP
Un1on64CzQcv6xb8OQrO5/zKyOpFVHa0NPjuJ3DVtUV+ZpzGKTSUfOQ8juXu
r98axwmgX7pgeWpEANBaYTHdi1jUO2FQek69jxtZEXVEhv8BXfc3QQ==
        
        "]]}, {}}, {}, {}, {}, {}}, {{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJwVzXk41HkcwPFxlStmTMnIsYoUjYiUrs+3ZXOmRaUiZkVqtjZabfTY2rGj
JeeWR5M7V4cn40yL+vYYW48aSUbKamaaBoPhNzPreEhm2z/ez+vPt03UueAY
TRKJtP9r/+sXI+t9OnZoz+D9+aOzqVSk7WkeJ6K5QdWgk70Jm4pcSzlFjTQv
OKcfPU1PoSL/7oGGYtpB0GK+yo65TEWh0XOPr9Oigb6xgif4hYo+bP7DM4uW
AFfuBNAbY6mobNWBq9dobFhXWUqK86Uilm6jfTotD84Uf3dnzJCKjhy9zVUo
K0HDPuNYbboJMlCzwr2tuPCkPD3MztgE2QVNyXgZTSAvOtygdZOCLr7yAM70
X2DlZHmucSUFrW6iCDU0MCTWFyWtKiSjhxJPfs+uDhC8XDXus5qMOpjFa47K
O8Grq5Xx3NEYzazfbfzg7HNwEbWl/bDDCGnO15DX2L0A23V1ySVuK5C7mvEw
pJYPsf6mhggZIj7b7dUzhx7IiqN5131vgKxzJ+bvTb2GRzzjQF6APpqYEwn7
nd+A2+eLih3H9FDleNsLo+w+eMjoizgUoouMLO8QUqkALHgvE/cELUeZp5tP
ltm+hV49d6dMxjI0i2qkp64MQNblG9zkMB2UsKWzyZT3Dm4NRPF7Y7WRqLYz
o9dmEGoi7Jq7IrXQZ6agcOOlf6BUwvbNDddEjnuXL51sG4Jhv8NJbUwN5F1u
E1u/WgjOH80vpEeT0EXLyLUfw0VQ0FJulihagk3yQwYD+8XAqfYtuaReBKqx
QrdzQQxlyrdtHbMLoI0LueuKPoKfjg9ZQJ6HH93a3uv5S0DzQumig+Yc1PhG
lT2ZlUCy0BlZf5mGmQlz8kj+J+CQq1uHzP+F5hBvrTwvKXSUaHSLjJVwxSe1
tGlECmYTrE1HDkyBzQe/dxWpw7A9M+Gs4/Vx2LGNnUfbMgJxw67C9N9HINyq
Z16jbwQUPa3tl1olMPMiSxLHGgXSsEX+JHwAjnZEQrCtDMoCQrL18/uBKmwu
1+yUgdj90zMPOz7o/Mo66BE/Buh0FcO0vgVI9yK4EpNxIHkuJnZPFeDzBdws
k45x+GY4oocyyMOPHdN0KT9NgNhjtICr/QZT/M5vSjGSQ5kw57zLn+/xEP11
X1K7HDjpBy4EisW4rgRspRGTwLmew4tvkOLjXHeK5dIkmN0NO3EsXoZ/S73N
Gro/BYk6YX0tLDl+ILMsNAoggH+qi7l3vQKnFacyhIEE2I4G9jxyUuDo4Em7
2iAC2DECV+dtCmzR3l63P5QAFCVatPJR4Gs5Yc8yowioC5vJXjytwKfcbyn1
kwjgBNo0tjxQ4LXslT7LqwmIdkv6TN+qxEvbk1cM3CXgacNSZNVuJX4/+elN
dQ0BFi7sTot9Spwb2nh8Xz0BAnpOlkGoEqsdgn9ObSfAx77KeixRiYd6c0q0
+ghwNu/1rGxX4kdX50709ROQyTl8d83fSpy3M3JjxTsC5KZDhje6ldi/yqnp
WyEBldTRfpZQiVsTu7tYYwSYrfgSG0lS4Xz61pwgOQEJGSn8t7oqHC8pOmhD
EPBaT9clkKLCGwLOiJ5OE5C2jLqwa60Ka2sIKnPnCJCxORFNDiosbt7JZCwQ
4KVlxXN0VeF2ZsVm5y9ff6wK+4qdKnzT2mBGrSZgUb0h09xLhf8DF9toMg==

         "]]},
       Annotation[#, "Charting`Private`Tag$226997#1"]& ]}, {}, {}}},
   AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
   Axes->{True, True},
   AxesLabel->{None, None},
   AxesOrigin->{0.47927344606913114`, 0},
   DisplayFunction->Identity,
   Frame->{{True, True}, {True, True}},
   FrameLabel->{{
      FormBox["\"logc\"", TraditionalForm], None}, {
      FormBox["\"d\"", TraditionalForm], None}},
   FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
   GridLines->{Automatic, Automatic},
   GridLinesStyle->Directive[
     GrayLevel[0.5, 0.4]],
   ImagePadding->All,
   ImageSize->500,
   Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         (Identity[#]& )[
          Part[#, 1]], 
         (Identity[#]& )[
          Part[#, 2]]}& )}},
   PlotLabel->FormBox[
     StyleBox[
     "\"d  versus  log(c)\"", FontFamily -> "Helvetica", FontSize -> 12, 
      FontWeight -> "Bold", StripOnInput -> False], TraditionalForm],
   PlotRange->{{0.4830648979126792, 0.7257178158997565}, {
     0, 1.9663470362925983`}},
   PlotRangeClipping->True,
   PlotRangePadding->{{
      Scaled[0.02], 
      Scaled[0.02]}, {
      Scaled[0.02], 
      Scaled[0.05]}},
   Ticks->{Automatic, Automatic}], TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.733159926232314*^9},
 CellLabel->"Out[99]=",ExpressionUUID->"a8ef70d7-84ab-4af9-897a-8723d89102f8"],

Cell[BoxData[
 FormBox[
  StyleBox[
   TagBox[GridBox[{
      {"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Standard Error\"\>", "\<\"t\
\[Hyphen]Statistic\"\>", "\<\"P\[Hyphen]Value\"\>"},
      {"1", "4.880822691670382`", "0.06328475395891828`", 
       "77.12477945065253`", "1.5346862834673565`*^-89"},
      {"x", 
       RowBox[{"-", "6.0224496947470145`"}], "0.1111168253028416`", 
       RowBox[{"-", "54.199259908058245`"}], "7.121598062934457`*^-75"}
     },
     AutoDelete->False,
     GridBoxAlignment->{"Columns" -> {{Left}}, "Rows" -> {{Automatic}}},
     GridBoxDividers->{
      "ColumnsIndexed" -> {2 -> GrayLevel[0.7]}, 
       "RowsIndexed" -> {2 -> GrayLevel[0.7]}},
     GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}},
     GridBoxSpacings->{
      "ColumnsIndexed" -> {2 -> 1}, "RowsIndexed" -> {2 -> 0.75}}],
    "Grid"], "DialogStyle",
   StripOnInput->False], TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.733159926247053*^9},
 CellLabel->
  "Out[100]=",ExpressionUUID->"f68c795d-ef46-4b88-9bd2-ac51d78958ed"],

Cell[BoxData[
 FormBox[
  RowBox[{"4.880822691670382`", "\[VeryThinSpace]", "-", 
   RowBox[{"6.0224496947470145`", " ", "x"}]}], TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.733159926257661*^9},
 CellLabel->
  "Out[101]=",ExpressionUUID->"304cf539-ab4c-45b3-9052-84068c9366d2"],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{"0.06328475395891828`", ",", "0.1111168253028416`"}], "}"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.7331599262634277`*^9},
 CellLabel->
  "Out[102]=",ExpressionUUID->"4ca753dd-0e81-42b0-982b-ee2d349bee95"],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{"0.5680667719460142`", ",", "1.459669134368186`"}], "}"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.7331599262675858`*^9},
 CellLabel->
  "Out[103]=",ExpressionUUID->"61cb03cf-75c8-4bb1-8595-d6903b776769"],

Cell[BoxData[
 FormBox[
  RowBox[{"{", 
   RowBox[{"0.04105361694266624`", ",", "0.2513336539127534`"}], "}"}], 
  TraditionalForm]], "Output",
 CellChangeTimes->{{3.733158207839028*^9, 3.733158240475546*^9}, {
   3.7331582841150713`*^9, 3.733158317085606*^9}, 3.733159366391346*^9, 
   3.73315942555824*^9, 3.733159926271749*^9},
 CellLabel->
  "Out[104]=",ExpressionUUID->"6e9a90af-5e94-4e1d-9ca2-2bc1ac9b52e9"]
}, Open  ]],

Cell[TextData[StyleBox["From the above statistical measurements, we could \
conclude that the best value for d is 0.568, and the best value for log(c) is \
1.460.",
 FontSize->18]], "Text",
 CellChangeTimes->{{3.733159800993657*^9, 3.733159809751629*^9}, {
  3.7331599384950953`*^9, 3.7331599585759478`*^9}, {3.733159996191429*^9, 
  3.73316004328229*^9}},ExpressionUUID->"2a1fa5e3-7911-4d9e-a197-\
696b5288c5f7"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Conclusion and Speculation", "Section",
 CellChangeTimes->{{3.7323866135820007`*^9, 
  3.732386622545092*^9}},ExpressionUUID->"bd860b5c-dce7-48c3-b5d4-\
d33c964b366d"],

Cell["\<\
Based on our initial model speculation of log(b) = log(c) + d*log(n), the \
relation between log(b) and log(n) is log(b) = 1.460 + 0.568*log(n), \
resulting from a sample size of 100 generations each on a sizeList of 1000.\
\>", "Text",
 CellChangeTimes->{{3.73238670161837*^9, 3.732386831081956*^9}, {
  3.733160073156723*^9, 3.733160241087892*^9}, {3.733160275436922*^9, 
  3.733160322631507*^9}},ExpressionUUID->"348fe2a1-b9d6-4c74-8d88-\
baa7778ea4c1"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{808, 813},
WindowMargins->{{Automatic, 0}, {Automatic, 0}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 368, 9, 210, "Title",ExpressionUUID->"4ff34a9b-c9e3-4e77-b2ad-1e53638bf927"],
Cell[951, 33, 295, 5, 64, "Text",ExpressionUUID->"60176971-891a-498f-b05e-1c7b25d8f0be"],
Cell[CellGroupData[{
Cell[1271, 42, 380, 6, 67, "Section",ExpressionUUID->"9463184e-e0ea-4946-b109-adcbb5c7ef73"],
Cell[1654, 50, 1006, 20, 116, "Text",ExpressionUUID->"4fad08c4-ae96-4325-9755-20c424ac9425"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2697, 75, 271, 4, 67, "Section",ExpressionUUID->"79597f95-25b4-4b2b-9164-b6fa06ab1942"],
Cell[2971, 81, 1297, 34, 94, "Input",ExpressionUUID->"06f8333f-1096-4432-b651-057b0df4e38c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4305, 120, 443, 6, 67, "Section",ExpressionUUID->"cb4ce077-be6b-4232-94f6-9320176a7592"],
Cell[4751, 128, 202, 4, 38, "Text",ExpressionUUID->"b3c68032-dd45-4fbf-8e7b-cd945415abf5"],
Cell[4956, 134, 2550, 64, 367, "Input",ExpressionUUID->"b7df351a-bfac-4d07-8ad5-7f5fafa20ca2"],
Cell[7509, 200, 335, 6, 38, "Text",ExpressionUUID->"15d372f7-6c98-4b21-9ceb-2509e7a9fbb5"],
Cell[CellGroupData[{
Cell[7869, 210, 5038, 116, 409, "Input",ExpressionUUID->"35ad2452-f27c-4c52-b62d-92a2a587d859"],
Cell[12910, 328, 7507, 143, 351, "Output",ExpressionUUID->"c9e48d74-3546-40ed-9446-b5c4202ab4cf"],
Cell[20420, 473, 519, 11, 26, "Message",ExpressionUUID->"9c6428c8-3148-4b33-adee-c444c14ac186"],
Cell[20942, 486, 519, 11, 26, "Message",ExpressionUUID->"14d17d28-d4df-4eb1-b528-44e9a4d8e4b2"],
Cell[21464, 499, 1391, 28, 65, "Output",ExpressionUUID->"6fbad36a-c598-48e0-8271-fc579c6f6187"],
Cell[22858, 529, 666, 12, 30, "Output",ExpressionUUID->"87052b97-38f4-40c5-87fc-778a4e6d04a7"],
Cell[23527, 543, 670, 12, 30, "Output",ExpressionUUID->"f3acf057-6316-4730-bef1-d7af6f6edf25"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[24246, 561, 320, 6, 105, "Section",ExpressionUUID->"0920a608-da80-466c-a68e-eeadb56acccc"],
Cell[24569, 569, 684, 11, 90, "Text",ExpressionUUID->"e9e64421-6314-477c-add1-b354663cf1ee"],
Cell[CellGroupData[{
Cell[25278, 584, 3939, 79, 220, "Input",ExpressionUUID->"9670b60a-7a9c-422a-8f0f-6933a6b7ca80"],
Cell[29220, 665, 1491, 24, 30, "Output",ExpressionUUID->"b3f573d9-1e53-4a86-81de-6e9bb04ef31a"],
Cell[30714, 691, 29998, 509, 353, 13898, 247, "CachedBoxData", "BoxData", "Output",ExpressionUUID->"4c390618-0fde-4523-a57d-43dc3e557838"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[60761, 1206, 231, 5, 105, "Section",ExpressionUUID->"6c70c45b-21cc-46c9-b1bc-03386cc6c6e2"],
Cell[60995, 1213, 440, 7, 64, "Text",ExpressionUUID->"c97e94f9-b0c7-4336-aca0-f618bfdb2128"],
Cell[CellGroupData[{
Cell[61460, 1224, 4028, 101, 388, "Input",ExpressionUUID->"a32f40f2-4dc9-4745-9e47-686810c3c85c"],
Cell[65491, 1327, 412, 9, 30, "Output",ExpressionUUID->"b3abc91f-76b9-43dc-a5bd-3348a3e3604b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[65940, 1341, 2259, 52, 262, "Input",ExpressionUUID->"f8649f83-d5f8-4efc-b6f5-7ff69c5ca5dd"],
Cell[68202, 1395, 5876, 116, 353, "Output",ExpressionUUID->"a8ef70d7-84ab-4af9-897a-8723d89102f8"],
Cell[74081, 1513, 1200, 26, 69, "Output",ExpressionUUID->"f68c795d-ef46-4b88-9bd2-ac51d78958ed"],
Cell[75284, 1541, 428, 8, 30, "Output",ExpressionUUID->"304cf539-ab4c-45b3-9052-84068c9366d2"],
Cell[75715, 1551, 415, 9, 30, "Output",ExpressionUUID->"4ca753dd-0e81-42b0-982b-ee2d349bee95"],
Cell[76133, 1562, 413, 9, 30, "Output",ExpressionUUID->"61cb03cf-75c8-4bb1-8595-d6903b776769"],
Cell[76549, 1573, 413, 9, 30, "Output",ExpressionUUID->"6e9a90af-5e94-4e1d-9ca2-2bc1ac9b52e9"]
}, Open  ]],
Cell[76977, 1585, 413, 7, 64, "Text",ExpressionUUID->"2a1fa5e3-7911-4d9e-a197-696b5288c5f7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[77427, 1597, 173, 3, 67, "Section",ExpressionUUID->"bd860b5c-dce7-48c3-b5d4-d33c964b366d"],
Cell[77603, 1602, 466, 8, 58, "Text",ExpressionUUID->"348fe2a1-b9d6-4c74-8d88-baa7778ea4c1"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)