(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 67654, 1408] NotebookOptionsPosition[ 66737, 1374] NotebookOutlinePosition[ 67146, 1392] CellTagsIndexPosition[ 67103, 1389] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[StyleBox["Monte Carlo Statistics for a 2D Ising Model with \ Periodic Boundary Conditions", FontSize->16]], "Subsection", CellChangeTimes->{ 3.417433398298102*^9, {3.417433790647161*^9, 3.417433800308227*^9}, { 3.417825918718843*^9, 3.417825919003312*^9}}, FontSize->10], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Set", " ", "the", " ", "Size", " ", "of", " ", "the", " ", "Lattice"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{"rows", "=", "16"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"number", " ", "of", " ", "rows", " ", "and", " ", "columns"}], " ", "-", " ", RowBox[{ "best", " ", "set", " ", "to", " ", "be", " ", "an", " ", "even", " ", "number"}]}], " ", "*)"}]}]], "Input", CellChangeTimes->{{3.417429727748459*^9, 3.417429763088395*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Generate", " ", "\"\\"", " ", "Initial", " ", "Configuration"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"maxj", "=", RowBox[{"rows", "/", "2"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{ "the", " ", "last", " ", "row", " ", "to", " ", "be", " ", "filled", " ", "with"}], " ", "+", RowBox[{ RowBox[{"1", "'"}], "s"}]}], " ", "*)"}], "\n", RowBox[{ RowBox[{"nplus", "=", RowBox[{"rows", "*", "maxj"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"running", " ", "count", " ", "of"}], " ", "+", RowBox[{ RowBox[{"1", "'"}], "s"}]}], " ", "*)"}], "\n", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"0", ",", "0"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"A", "[", RowBox[{"0", ",", RowBox[{"rows", "+", "1"}]}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"rows", "+", "1"}], ",", "0"}], "]"}], "=", "0"}], ";", RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"rows", "+", "1"}], ",", RowBox[{"rows", "+", "1"}]}], "]"}], "=", "0"}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"the", " ", "four", " ", "corners"}], " ", "*)"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}], "=", RowBox[{"+", "1"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "maxj"}], "}"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}], "=", RowBox[{"-", "1"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", RowBox[{"maxj", "+", "1"}], ",", "rows"}], "}"}]}], "]"}], ";"}]}]}]], "Input", CellChangeTimes->{{3.417429750830261*^9, 3.417429757611227*^9}, { 3.417429804516849*^9, 3.417429847194579*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Generate", " ", "\"\\"", " ", "Initial", " ", "Configuration"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"nplus", "=", RowBox[{ RowBox[{"rows", "^", "2"}], "/", "2"}]}], ";"}], " ", StyleBox[ RowBox[{"(*", " ", RowBox[{ "ok", " ", "if", " ", "rows", " ", "is", " ", "an", " ", "even", " ", "number"}], " ", "*)"}], FontColor->RGBColor[1, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"1", ",", "1"}], "]"}], "=", RowBox[{"+", "1"}]}], ";"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "1"}], "]"}], "=", RowBox[{"-", RowBox[{"A", "[", RowBox[{"i", ",", "1"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}]}], "]"}], "=", RowBox[{"-", RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}]}]}]], "Input"], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Impose", " ", "Periodic", " ", "Boundary", " ", "Conditions"}], ",", " ", RowBox[{"Calculate", " ", "Averages"}]}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "0"}], "]"}], "=", RowBox[{"A", "[", RowBox[{"i", ",", "rows"}], "]"}]}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", RowBox[{"rows", "+", "1"}]}], "]"}], "=", RowBox[{"A", "[", RowBox[{"i", ",", "1"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"0", ",", "j"}], "]"}], "=", RowBox[{"A", "[", RowBox[{"rows", ",", "j"}], "]"}]}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"rows", "+", "1"}], ",", "j"}], "]"}], "=", RowBox[{"A", "[", RowBox[{"1", ",", "j"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"means", "[", "0", "]"}], "=", "0."}], ";", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"means", "[", "0", "]"}], "=", RowBox[{ RowBox[{"means", "[", "0", "]"}], "+", RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";", " ", RowBox[{ RowBox[{"means", "[", "0", "]"}], "=", RowBox[{ RowBox[{"means", "[", "0", "]"}], "/", RowBox[{"rows", "^", "2"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"corr", "[", "0", "]"}], "=", "0."}], ";", " ", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"corr", "[", "0", "]"}], "=", RowBox[{ RowBox[{"corr", "[", "0", "]"}], "+", RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}]}], "]"}], "+", RowBox[{"A", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j"}], "]"}]}], ")"}]}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"corr", "[", "0", "]"}], "=", RowBox[{ RowBox[{ RowBox[{"corr", "[", "0", "]"}], "/", RowBox[{"(", RowBox[{"2", "*", RowBox[{"rows", "^", "2"}]}], ")"}]}], "-", RowBox[{ RowBox[{"means", "[", "0", "]"}], "^", "2"}]}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.417825879768565*^9, 3.417825912938848*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Generate", " ", "and", " ", "Display", " ", "Graphics", " ", "from", " ", "a", " ", "Configuration"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"ipl", "=", "0"}], ";", " ", RowBox[{"size", "=", "0.1"}], ";"}], "\n", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}], "==", "1"}], ",", "\n", "\t", RowBox[{"{", RowBox[{ RowBox[{"ipl", "=", RowBox[{"ipl", "+", "1"}]}], ",", RowBox[{ RowBox[{"draw", "[", "ipl", "]"}], "=", RowBox[{"Graphics", "[", RowBox[{"Rectangle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"size", "*", "j"}], ",", RowBox[{"size", "*", "i"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"size", "*", RowBox[{"(", RowBox[{"j", "+", "1"}], ")"}]}], ",", RowBox[{"size", "*", RowBox[{"(", RowBox[{"i", "+", "1"}], ")"}]}]}], "}"}]}], "]"}], "]"}]}]}], "}"}]}], "]"}], ",", "\n", RowBox[{"{", RowBox[{"i", ",", "0", ",", RowBox[{"rows", "+", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "0", ",", RowBox[{"rows", "+", "1"}]}], "}"}]}], "]"}], ";"}], "\n", RowBox[{"Show", "[", RowBox[{ RowBox[{"Array", "[", RowBox[{"draw", ",", "ipl"}], "]"}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"size", "*", RowBox[{"(", RowBox[{"rows", "+", "2"}], ")"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"size", "*", RowBox[{"(", RowBox[{"rows", "+", "2"}], ")"}]}]}], "}"}]}], "}"}]}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.417826021382011*^9, 3.417826069737455*^9}, { 3.419002823863179*^9, 3.419002824125122*^9}}], Cell[BoxData[ GraphicsBox[{ RectangleBox[{0.6000000000000001, 0}, {0.7000000000000001, 0.1}], RectangleBox[{0.8, 0}, {0.9, 0.1}], RectangleBox[{0.9, 0}, {1., 0.1}], RectangleBox[{1., 0}, {1.1, 0.1}], RectangleBox[{1.1, 0}, {1.2000000000000002`, 0.1}], RectangleBox[{1.4000000000000001`, 0}, {1.5, 0.1}], RectangleBox[{1.6, 0}, {1.7000000000000002`, 0.1}], RectangleBox[{0, 0.1}, {0.1, 0.2}], RectangleBox[{0.9, 0.1}, {1., 0.2}], RectangleBox[{1., 0.1}, {1.1, 0.2}], RectangleBox[{1.3, 0.1}, {1.4000000000000001`, 0.2}], RectangleBox[{1.5, 0.1}, {1.6, 0.2}], RectangleBox[{1.6, 0.1}, {1.7000000000000002`, 0.2}], RectangleBox[{0, 0.2}, {0.1, 0.30000000000000004`}], RectangleBox[{0.5, 0.2}, {0.6000000000000001, 0.30000000000000004`}], RectangleBox[{0.6000000000000001, 0.2}, \ {0.7000000000000001, 0.30000000000000004`}], RectangleBox[{1., 0.2}, {1.1, 0.30000000000000004`}], RectangleBox[{1.6, 0.2}, {1.7000000000000002`, 0.30000000000000004`}], RectangleBox[{0.5, 0.30000000000000004`}, {0.6000000000000001, 0.4}], RectangleBox[{0.2, 0.5}, {0.30000000000000004`, 0.6000000000000001}], RectangleBox[{0.30000000000000004`, 0.5}, {0.4, 0.6000000000000001}], RectangleBox[{0.7000000000000001, 0.5}, {0.8, 0.6000000000000001}], RectangleBox[{0.8, 0.5}, {0.9, 0.6000000000000001}], RectangleBox[{0.9, 0.5}, {1., 0.6000000000000001}], RectangleBox[{1., 0.5}, {1.1, 0.6000000000000001}], RectangleBox[{1.4000000000000001`, 0.5}, {1.5, 0.6000000000000001}], RectangleBox[{0.2, 0.6000000000000001}, \ {0.30000000000000004`, 0.7000000000000001}], RectangleBox[{0.9, 0.6000000000000001}, {1., 0.7000000000000001}], RectangleBox[{1.2000000000000002`, 0.6000000000000001}, \ {1.3, 0.7000000000000001}], RectangleBox[{1.3, 0.6000000000000001}, \ {1.4000000000000001`, 0.7000000000000001}], RectangleBox[{0.2, 0.7000000000000001}, {0.30000000000000004`, 0.8}], RectangleBox[{0.30000000000000004`, 0.7000000000000001}, {0.4, 0.8}], RectangleBox[{0.7000000000000001, 0.7000000000000001}, {0.8, 0.8}], RectangleBox[{0.9, 0.7000000000000001}, {1., 0.8}], RectangleBox[{1., 0.7000000000000001}, {1.1, 0.8}], RectangleBox[{1.1, 0.7000000000000001}, {1.2000000000000002`, 0.8}], RectangleBox[{1.2000000000000002`, 0.7000000000000001}, {1.3, 0.8}], RectangleBox[{0.1, 0.8}, {0.2, 0.9}], RectangleBox[{0.2, 0.8}, {0.30000000000000004`, 0.9}], RectangleBox[{1., 0.8}, {1.1, 0.9}], RectangleBox[{1.1, 0.8}, {1.2000000000000002`, 0.9}], RectangleBox[{1.2000000000000002`, 0.8}, {1.3, 0.9}], RectangleBox[{1.3, 0.8}, {1.4000000000000001`, 0.9}], RectangleBox[{1.4000000000000001`, 0.8}, {1.5, 0.9}], RectangleBox[{1.7000000000000002`, 0.8}, {1.8, 0.9}], RectangleBox[{0, 0.9}, {0.1, 1.}], RectangleBox[{0.1, 0.9}, {0.2, 1.}], RectangleBox[{0.2, 0.9}, {0.30000000000000004`, 1.}], RectangleBox[{0.4, 0.9}, {0.5, 1.}], RectangleBox[{0.7000000000000001, 0.9}, {0.8, 1.}], RectangleBox[{0.9, 0.9}, {1., 1.}], RectangleBox[{1., 0.9}, {1.1, 1.}], RectangleBox[{1.1, 0.9}, {1.2000000000000002`, 1.}], RectangleBox[{1.2000000000000002`, 0.9}, {1.3, 1.}], RectangleBox[{1.3, 0.9}, {1.4000000000000001`, 1.}], RectangleBox[{1.4000000000000001`, 0.9}, {1.5, 1.}], RectangleBox[{1.5, 0.9}, {1.6, 1.}], RectangleBox[{1.6, 0.9}, {1.7000000000000002`, 1.}], RectangleBox[{1.7000000000000002`, 0.9}, {1.8, 1.}], RectangleBox[{0, 1.}, {0.1, 1.1}], RectangleBox[{0.1, 1.}, {0.2, 1.1}], RectangleBox[{0.4, 1.}, {0.5, 1.1}], RectangleBox[{0.5, 1.}, {0.6000000000000001, 1.1}], RectangleBox[{0.7000000000000001, 1.}, {0.8, 1.1}], RectangleBox[{0.9, 1.}, {1., 1.1}], RectangleBox[{1., 1.}, {1.1, 1.1}], RectangleBox[{1.1, 1.}, {1.2000000000000002`, 1.1}], RectangleBox[{1.2000000000000002`, 1.}, {1.3, 1.1}], RectangleBox[{1.3, 1.}, {1.4000000000000001`, 1.1}], RectangleBox[{1.4000000000000001`, 1.}, {1.5, 1.1}], RectangleBox[{1.5, 1.}, {1.6, 1.1}], RectangleBox[{1.6, 1.}, {1.7000000000000002`, 1.1}], RectangleBox[{1.7000000000000002`, 1.}, {1.8, 1.1}], RectangleBox[{0, 1.1}, {0.1, 1.2000000000000002`}], RectangleBox[{0.5, 1.1}, {0.6000000000000001, 1.2000000000000002`}], RectangleBox[{0.6000000000000001, 1.1}, \ {0.7000000000000001, 1.2000000000000002`}], RectangleBox[{0.7000000000000001, 1.1}, {0.8, 1.2000000000000002`}], RectangleBox[{0.9, 1.1}, {1., 1.2000000000000002`}], RectangleBox[{1., 1.1}, {1.1, 1.2000000000000002`}], RectangleBox[{1.1, 1.1}, {1.2000000000000002`, 1.2000000000000002`}], RectangleBox[{1.2000000000000002`, 1.1}, {1.3, 1.2000000000000002`}], RectangleBox[{1.4000000000000001`, 1.1}, {1.5, 1.2000000000000002`}], RectangleBox[{1.5, 1.1}, {1.6, 1.2000000000000002`}], RectangleBox[{1.6, 1.1}, {1.7000000000000002`, 1.2000000000000002`}], RectangleBox[{0, 1.2000000000000002`}, {0.1, 1.3}], RectangleBox[{0.5, 1.2000000000000002`}, {0.6000000000000001, 1.3}], RectangleBox[{0.6000000000000001, 1.2000000000000002`}, \ {0.7000000000000001, 1.3}], RectangleBox[{0.8, 1.2000000000000002`}, {0.9, 1.3}], RectangleBox[{0.9, 1.2000000000000002`}, {1., 1.3}], RectangleBox[{1.2000000000000002`, 1.2000000000000002`}, {1.3, 1.3}], RectangleBox[{1.4000000000000001`, 1.2000000000000002`}, {1.5, 1.3}], RectangleBox[{1.5, 1.2000000000000002`}, {1.6, 1.3}], RectangleBox[{1.6, 1.2000000000000002`}, {1.7000000000000002`, 1.3}], RectangleBox[{0.6000000000000001, 1.3}, \ {0.7000000000000001, 1.4000000000000001`}], RectangleBox[{0.8, 1.3}, {0.9, 1.4000000000000001`}], RectangleBox[{0.9, 1.3}, {1., 1.4000000000000001`}], RectangleBox[{1.2000000000000002`, 1.3}, {1.3, 1.4000000000000001`}], RectangleBox[{1.3, 1.3}, {1.4000000000000001`, 1.4000000000000001`}], RectangleBox[{1.4000000000000001`, 1.3}, {1.5, 1.4000000000000001`}], RectangleBox[{0.1, 1.4000000000000001`}, {0.2, 1.5}], RectangleBox[{0.5, 1.4000000000000001`}, {0.6000000000000001, 1.5}], RectangleBox[{0.6000000000000001, 1.4000000000000001`}, \ {0.7000000000000001, 1.5}], RectangleBox[{0.7000000000000001, 1.4000000000000001`}, {0.8, 1.5}], RectangleBox[{0.8, 1.4000000000000001`}, {0.9, 1.5}], RectangleBox[{0.9, 1.4000000000000001`}, {1., 1.5}], RectangleBox[{1.7000000000000002`, 1.4000000000000001`}, {1.8, 1.5}], RectangleBox[{0.8, 1.5}, {0.9, 1.6}], RectangleBox[{0.9, 1.5}, {1., 1.6}], RectangleBox[{1.1, 1.5}, {1.2000000000000002`, 1.6}], RectangleBox[{1.4000000000000001`, 1.5}, {1.5, 1.6}], RectangleBox[{0, 1.6}, {0.1, 1.7000000000000002`}], RectangleBox[{0.6000000000000001, 1.6}, \ {0.7000000000000001, 1.7000000000000002`}], RectangleBox[{0.8, 1.6}, {0.9, 1.7000000000000002`}], RectangleBox[{0.9, 1.6}, {1., 1.7000000000000002`}], RectangleBox[{1., 1.6}, {1.1, 1.7000000000000002`}], RectangleBox[{1.1, 1.6}, {1.2000000000000002`, 1.7000000000000002`}], RectangleBox[{1.4000000000000001`, 1.6}, {1.5, 1.7000000000000002`}], RectangleBox[{1.6, 1.6}, {1.7000000000000002`, 1.7000000000000002`}], RectangleBox[{0.9, 1.7000000000000002`}, {1., 1.8}], RectangleBox[{1., 1.7000000000000002`}, {1.1, 1.8}], RectangleBox[{1.3, 1.7000000000000002`}, {1.4000000000000001`, 1.8}], RectangleBox[{1.5, 1.7000000000000002`}, {1.6, 1.8}], RectangleBox[{1.6, 1.7000000000000002`}, {1.7000000000000002`, 1.8}]}, AspectRatio->1, Axes->True, PlotRange->{{0, 1.8}, {0, 1.8}}]], "Output", CellChangeTimes->{3.41900283088627*^9, 3.4190028974517813`*^9, 3.419003140428793*^9, 3.4190032997140636`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Evolve", " ", "the", " ", "System", " ", "and", " ", "Calculate", " ", "Average", " ", "Spin"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{"Tinv", "=", "0.30"}], ";", " ", RowBox[{"ncount", "=", "10000"}], ";"}], "\n", RowBox[{"(", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"m", "=", RowBox[{"RandomInteger", "[", RowBox[{"{", RowBox[{"1", ",", "rows"}], "}"}], "]"}]}], ",", RowBox[{"n", "=", RowBox[{"RandomInteger", "[", RowBox[{"{", RowBox[{"1", ",", "rows"}], "}"}], "]"}]}], ",", " ", RowBox[{"(*", " ", RowBox[{"pick", " ", "element", " ", "to", " ", "flip"}], " ", "*)"}], "\n", RowBox[{"S", "=", RowBox[{"-", RowBox[{"A", "[", RowBox[{"n", ",", "m"}], "]"}]}]}], ",", RowBox[{"B", "=", RowBox[{"Tinv", " ", "S", " ", RowBox[{"(", RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"n", "-", "1"}], ",", "m"}], "]"}], "+", RowBox[{"A", "[", RowBox[{ RowBox[{"n", "+", "1"}], ",", "m"}], "]"}], "+", RowBox[{"A", "[", RowBox[{"n", ",", RowBox[{"m", "-", "1"}]}], "]"}], "+", RowBox[{"A", "[", RowBox[{"n", ",", RowBox[{"m", "+", "1"}]}], "]"}]}], ")"}], " ", "2"}]}], ",", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"B", " ", "is"}], " ", "-", RowBox[{"beta", " ", "Delta", " ", "E"}]}], " ", "*)"}], "\n", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"RandomReal", "[", "]"}], "<", SuperscriptBox["\[ExponentialE]", "B"]}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"A", "[", RowBox[{"n", ",", "m"}], "]"}], "=", "S"}], ",", RowBox[{"nplus", "=", RowBox[{"nplus", "+", "S"}]}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"n", "\[Equal]", "1"}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{ RowBox[{"rows", "+", "1"}], ",", "m"}], "]"}], "=", "S"}]}], "]"}], ",", "\n", RowBox[{"If", "[", RowBox[{ RowBox[{"n", "\[Equal]", "rows"}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{"0", ",", "m"}], "]"}], "=", "S"}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"m", "\[Equal]", "1"}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{"n", ",", RowBox[{"rows", "+", "1"}]}], "]"}], "=", "S"}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"m", "\[Equal]", "rows"}], ",", RowBox[{ RowBox[{"A", "[", RowBox[{"n", ",", "0"}], "]"}], "=", "S"}]}], "]"}]}], "}"}]}], "]"}], ",", " ", RowBox[{"(*", " ", "decide", " ", "*)"}], "\n", RowBox[{ RowBox[{"means", "[", "ict", "]"}], "=", "0"}], ",", "\n", RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"means", "[", "ict", "]"}], "=", RowBox[{ RowBox[{"means", "[", "ict", "]"}], "+", RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ",", RowBox[{"(*", " ", RowBox[{"updated", " ", "mean", " ", "spin"}], " ", "*)"}], "\n", RowBox[{ RowBox[{"means", "[", "ict", "]"}], "=", FractionBox[ RowBox[{"means", "[", "ict", "]"}], SuperscriptBox["rows", "2"]]}], ",", "Null"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"ict", ",", "1", ",", "ncount"}], "}"}]}], "]"}], ";"}], ")"}]}]}]], "Input", CellChangeTimes->{{3.417430161360611*^9, 3.4174301868307867`*^9}, { 3.417430231667152*^9, 3.41743023602277*^9}, {3.417430269220998*^9, 3.417430374391692*^9}, {3.417430609403017*^9, 3.417430675168221*^9}, { 3.417430765480627*^9, 3.417430799521707*^9}, {3.4174311006872272`*^9, 3.417431104570904*^9}, {3.417431254181736*^9, 3.417431263404743*^9}, { 3.4174314309521017`*^9, 3.417431431292149*^9}, {3.417431518121952*^9, 3.4174315182807302`*^9}, {3.417431775552869*^9, 3.417431775709869*^9}, { 3.417433614616404*^9, 3.4174336312152042`*^9}, {3.417440679087316*^9, 3.4174406796975594`*^9}, {3.417440992034153*^9, 3.417440994278307*^9}, { 3.417826165084106*^9, 3.417826198033606*^9}, {3.4178262824528713`*^9, 3.417826283152513*^9}, {3.419002083194409*^9, 3.4190021096948566`*^9}, { 3.4190022731183968`*^9, 3.419002274758841*^9}, {3.419003059891561*^9, 3.419003060272457*^9}, {3.419003190815522*^9, 3.419003197056642*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Display", " ", "Evolution", " ", "of", " ", "Average", " ", "Spin"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Array", "[", RowBox[{"means", ",", "ncount"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.417433472322674*^9, 3.4174335017760344`*^9}}], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointBox[CompressedData[" 1:eJxc3X+Y1lP+P/Bb0rYJsyQjyUgYSUaSJBlJkmQkGUlGkpFkJBlt296SjCQj yUgym2QkGUlGkpFkJIwkI8lIkiSztrWztm2/1/cz53Rd98M/ez32+Tqvc97n /eO+59fd8cNuH3BTo0QikXZyIvH//7fhv7rzw/++Hf6P7FQ3xWm4BU7HrXEG bovb4ZNwJm6POzS48pfgjjgLn4E74TNxZ3wW7oLPxl3xObgbPhd3x+fhHvh8 nJ3qJE5cQI4TPcmDkzG/MNVJnOhFjhMXkeNEb3KcuJgcJ/o0ODvmOHEJOU70 JceJS8lxoh85TlxGjhP9yXHicnKcyEl1Nk7iSpy4gvE4GZyI43FiQKqzcTK4 cncYjxNXpjobJ3ElTgxkPE7iSpy4ivE4iStxYhDjcRJX4sTVjMdJXBmcjONz U52Nk7gSJ65hPE7iSpwYzPjBrB9XUp+4lvE4iStxYgjjcRJX4sR1jMdJXIkT QxmPk7gSJ65nPE7iSpzIS3UGzsZ5OIlLcSWuxYkbmB9n4zycxKW4EtcGZ8f5 h6U6A2fjPJzEpbgS1+LEjcyPs3EeTuJSXIlrcWI48+Ps4EScHydxKa7EtcGV P4f5b0p1Bs7GeTiJS3ElrsWJEcwfnIzzk+fhJPWluBLX4sTNzI+zcR5O3sz8 uBLX4kQ+8+NsnJefOj6JS3ElrsWJW1KdgbNxHk7i0uDsOD+uxYmRDU7E+XE2 zsNJXIorcS1O3Mr8OBvn4SQuxZW4FidGMX9w5a4wf3AyOA8ncSmuxLX0T9yW 6gycfRvzB2fH+XFpcCLOj2txYnSqM4Irfwrz4zycxKW4MjgZXIsTt6c6Izg7 zh+ciPOTJ3Ep9ZW4NrhyZ5i/INVpwcngDJyFs3EOzsMFOImLcSkuZ72VuBrX Mr4OJ+5IdRrOwFnB2fH4cQ7OwwXBiXj8uBiX4nJcGVz5Yzj+uN7gWlyHE2NS x6eNSe2fgbPGcPw4B+fhApykfzEuxeW4ElfjWlzH/Ik7U52GM4L3n3+cjXNw Hi7ASeYrxqW4HFfialyL63BiLMePM3AWzsY5OA8XjOX4cTEuxeW4ElcH77// cV1wfB4m7kp1Gs7AWXel9ssmz8F5uAAncTEuxeXB8XlfiatxLa7DiXEcP87A WTgb54xL7Z+HC6hP4mJcistxZfD+5z+upb6OPHF3qtNwRvD++x9n4xychwtw Mnj//Y9LcTnjK3E1ruV46nCikOPHGTgLZ+McnIcLcBIXF6auvxSX40rGV+Na XMf4xD2pTgve//wjzyLPxjk4DxfgZPD+1z9cistxZfD+139ci+twYnyq04Lj +4MMnIWzcQ7OwwXMlyQvxqW4PDg7uDI4EY8f1+I6nPhzWM+OcPw4IzgZnBWc HZyNc3AeLsBJXIxLcTnrqcTVuJbjqcOJCRz/hNTxGTgLZ0/g+HEe9QU4iYtx KS7Hlbga17KeOpz4S6rTcAbOwtnBiXj8OC+48odw/MHJ4CQuxqWML8eVuBrX 4jrWl5jI8U9MzTNw1sTUftkTU9ebg/Ni/3j8OImLcSkuj+uJx4+rcS3j68gT f011Gs7AWcGV28Px4xychwtwEhfj0uBkcHlwdnAlrqa+FtfhRDLVTXEaTscZ OBNn4a7J1PVm4z7BibifOBfn4XxcgAtxEhfhYlyCS3EZLscVwZXfh/OHq3A1 rsG19N+B63A9Ttyb6qY4DafjDJwZHM9nFnlXnI374Byci/NwPi7AhcFx/5K4 CBfjElyKy3B5cDK4Alfiqrh/8fzjGlyLd+A6+tfjxKRUN8Vpk1L7pQcngjNw JvVZuCv9s3EfnINzcR7OxwW4kPUkcRH1xcH7zz8uxWW4HFfgSlyFq3ENrmW9 O8jrcD1O3Nfg/fc/TsPpOCM43u+ZOIv6rjgb98E5OJf+eTgfFzC+ECdxES7G JbgUlwXvP/+4AlfiKlyNa3Bt8P7zH48/nn/yesYnJqe6KU7D6TgDZ+Is3HVy 6nqycR/qc3Bu8P7nP3k+LqB/4eTU/UniIuqLcQn1pbgMlzO+AlfiKtZfjWtw LeN34DrWU48T96e6KU7D6TgDZ+Is3PX+1PVl3596PH1wDs4N3n/+cT4uYHwh 8ydxES7GJbgUl+FyXMF6KnEVrsY1uBbvYL46XI8TU1LdNHj/+SdPxxlTUufP xFm4K+OzcR+cg3NZXx7OxwW4ECdxES7GJbg0OL7/KSMvJ6/AldRX4WqOvwbX Ur8D1+F6nHgg1U0fSO2fhtOD9z//cSb1Wbgr82XjPjgH5+K84Lif+XE928L5 x4U4GdcXXISLcQkuxWW4PK43uAJXUl+Fq3ENrqXfDlyH63GiKNVNcRpOxxk4 E2fhrjgb98E5OBfnBSeC84Mrvwvnn7wQJ3ERLsYlzF8avP/84/K4nnj+Y794 /+Oq2D/e/7iG+lr678B11NfjxIOpborTHkztl44zcCbOol9XnP1g6vH1wTnU 5+I85svHBbgQJ3ER/YvJS3Bp8P7zj8txBa6kXxXzV5PXML4W76C+DtdTn5ja 4P33P06bmlqfTp4RHO+PTJxFfVfybPI+zJeDc3Ee4/PpX0BeiJP0KyIvxiX0 Lw2Oz5sy8nJcgStxFa7GNbgW72C9dRxfPU48lFrfFKc9lNo/HWdQn4mzHkqd ryvOpr4PzsG5OA/n078AF+Ik44vIi3FJ8P77H5fhclxBv0pchatZXw37X4t3 4Lo4f7g+62P/4MS0VDcNzg5OC04Ep+MMnImzgiu3hvOPs6nvg3NwLuPzcD71 BbiQ+iQuor4Yl1Bfistwedzf4ApciatwdTwfwTW4Nq4veAd5Ha7HiYdT3fTh 1PnTyNNxBs58OHU9WeRdybODK78N5x/nUJ9LvzzWmx/HBxfE8fH842TsF1yE i3EJLqVfGfOX4wrqK8mryKtxDfPXku/AdbgeJ6anuun01PWk4fTguN8ZwfvP P86K/eP5J8/Gfeifg3OpzyPPxwW4ECdxES5mvhKOpxSX4XJcgStxFa6O+x/u hxpci3cwvo711+PEI6n1TR9J7ZeG03EGzsRZwcngruTZuA/Owbk4j/Xnkxfg QuqTj6TuRxEupr6EvJS8jPnKcQWuZH+qcDWuwbV4B66L6w2ux4niVDcuTh3f lLw5TgtOBLfA6dS3pn8GeTuciTswPgt3pr4r7o6zcS/ch+PpF1xZG67POH/w QJyLh8T+wXl4eJwvOD+O/6bBo8gL8FhcSP8JrCeJJ+MiPA0Xsx8zyUvI5+BS 6ueTl7H/izg/5XgprsAr6FdJvpq8ivWsI69m/Aaux5o4PngzrsXb8I44X3i+ 7MJ11O/B9XgvTjya6sa4KW6O0x5NPd4Wj6buRzp5a/IM3I76TNwBZzG+M+6K u+Ns3Cs4nt8+5P1wDh6Ic/EQnIeH43yOdxQu4HyMxYXUT8BJ5puMi/A0XMx+ zSQvwXNwKZ5PvzK8iPpy8qW4Aq9gfCX5avIq9msdeTXeQL8avDmen/A8qcXb qN+Bd+E6vAfXs7695IkZqW6Mm85IXV/z4Pi8TsMtcDrjW5NnBMfXp3Y4k/oO 5FlxvcGdY314veyKu+Ns3Cv2C+6D+8X6LeH+j3nwwJh/He5/PCTWB+fh4Tg/ jt8c7n/yAvKxMQ8uJJ+Ak3gy44ti/lW4/3ExnolL8JzYP7g05pvC/U9ehhfh csYvxRV4Ba6k32ryKryO+mq8AdfgzbiW/dlGvoP5d+E66vfgeryX+RKPpbpx cDz/TXHzWB+chlvgdNwaZ+B2OJP5O+As3Bl35fi642zcC/fB/XAOHohz43qC h+A8PJzx+XgULuB4x+JC6ifgJJ7MeorIp+FiPBOXsJ45uJT6+biM9SyKebj+ y/HSWB/vf7wCV+LV9KvC66ivxhtwDd6Ma+m/jePdgXexP3V4D67He+mXmJnq xrgpbo7TZqauvwVOj/XBrWP+Zbj/cTucGccHd8BZsb4m3P8xD+5K3h1n068X 7oP74Rw8EOdyPEPYjzzy4Tif+lG4AI9l/wvxBM5fknwyLsLTcDGeyXpKOJ45 7E8p52M+568ML6K+HC+lvoJ8Ba6kfjV5FV6Hq/EGjq+GfDPz1ZJvI98R8y/C /Y/r8B5cH/sF7yVPPJ6aN348NW+Km+M03AKnx/7heFrjDOrbkWfiDjgLd8Zd 6d+d483GvWL9xnD/436xPjgHD4z1n4f7P+bBQ3Ae9cPpl08+ChdQPxYX4gk4 iSfjIuabhovxTFzC8c7BpdTPZ/4yvAiX46W4Aq/AlZzf1biK62Mdrub62oBr 8GZci7fRfwfexfrq8B5cz/HuxYlZqfWNcVPcHKfNSl1fC5we64Nbk2eQt8OZ uAPOwp1x1zhffP2PeXz9x71ifXz9J++Hc+g/kPlz8RDq8+g3HOeznlHkBXgs LmS+CeRJ8snkReTTcDGeyfgS1j8Hl1I/H5fhRcxXjpdSX4FX4ErGr+Z8VeF1 XL/VeAP3Rw3eHOuDa/E26neQ7yKvw3twPevbixNPpLoxboqb4zTcAqfj1k+k ri8Dt8OZuAPOwp2ZryvujrNxr9gvuA95P/IcPDDWh+spFw/BeYwfjvPxKFxA v7HkhXgCTuLJHG8R+TRczPwzyUvwHFzKfPNxGV6Ey/FSroeKmMfXf/JKvJp+ VeTr6FeNN+CaOD68v9mMa6nfhnfgXbgO72G99eR7yRMlqW5cklrfFDenPo28 BU6P9eF4W5Nn0K8deSbuQH0W7oy7Mr4768mOeXCvmG8I9z95P5yDB+Jc5huC 81jfcJyPR+ECPBYXsh8TYh7vf/LJuIj6abg41gfPjHl8/cdzcCnj5+MyvIjx 5Xgp9RWsdwWuxKs5/iq8jvpq5ttAXoM306+W87UN78C7GF9H/z24nvq9OPFk qhs/mTq+aczj6z9Owy1wOv1a4wzcDmfiDqw3C3eO9fH+J++Os3Evxvch74dz 8ECcy/qH4Dw8HOezn6NiHt//xzy+/uNCPAEnY7/wfnwyLsLTcDH9ZpKXkM8h L8XzcVkcH7wIl+OljK/AK3Al61uNq/A6XI034Bq8Gdeynm14B97F+Dq8B9fj vVw/idmpeePZqfM1xc1xWhwf3CLm8fv/uDXOYHw78syYB3cgz2J8Z9wVd8fZ HE8v3If96Rfz+PU/HohzGT8E5zHfcPJ8PAoXMN9YXIgn4GTsF9//kxfhabgY z6RfSczj+3/yUjwfl+FFuBwvxRV4Ba5k/avZ3yq8DlfjDbgGb8a1eBvewfp2 kdfhPdTX47048VTqfjR+KjVvipvH+uA08hb0S8etY318/cftcCbuEPvF1/+Y x9d/8q64O87GvXAf+vfDObE+vv/HuXhIHB/f/8c8fL0xPObB+XgU9QV4LC5k /ATyJJ6Mi/A0XIxncnwleA4uxfNxGf0X4XK8lPEVeAWuZPxqXIXXsZ/VeAP1 NXgzrmX8tph/Fu7/mAfvIq/De3A9/fcyf2JOqhsHx/1pipvjNNyCfum4NfUZ uB31mbhDrA/Owp1xV9w99gv7k4174T64X+wXnEM+EOfiITgPD8f5eBQuwGNZ XyGegJN4Mi6i/zTyYjwz1q8P9z+eE+uDS/F8XIYX0a+c+ZfiCryC46nEq3EV Xoerub42kNcw/2bWX4u34R0c/y7yOryH+nq8N9Z/Gu7/p1PdGDfFzXEabhGc CE7HrXEGboczcQechTvjrrg7zsa9cB/cD+dw/APJc2NeHe7/mAfnkQ8nz8ej cEH0J+H+j/443P/RH4X7nzyJJ+Mixk+LXhfu/+gPw/2PS/Acxpfi+biM8Yui 14b7n3wp4yvIV+BKvJrxVRz/OlyNN+Aa9nMzeS35Ns7nDvJduI5+e3A9x7MX J+amunF02I+m0WG/m5OnMb4FTseto8P6MnA76jOZrwPryYr+INz/0VXh/sfd cXb0++H+j14T7n/cL/q9cP/jgdTn0n8IzmM9w8nz8Sj6FzD/2OjV4f6Pfjfc /9Grwv2PJ1NfhKfhYjwTl+A5rK+UfD7rKcOLcHn0O+H+J6+g/wrmr2T/VuMq 9nsdeTX9NjBfDd7M+mpZ/7boynD/413U1+E91NfjvdErw/3/TPCKcP9HLw/3 P3lz8jTyFuTp0RXh/o9eFu5/8naMzyTvgLNwZ9yVft1xNuvvRd4H96N/Dh7I 8eXiITgvemm4/3E+9aNwAR7L+ELyCTiJJ3M8RXga+1FMPpN+JaxnDi6NXhLu f/Iy8kW4HC/FFXgF/SvxalzF+HW4Gm/ANXgzrsXb8I7o8nD/Ry8K9z/eg+vx XpwoTXUj3Bg3iQ7raRod1tsMN8eHMj4NH858Lchb0i89OpyvVuStydvgDOrb krcjP4n1ZOL21HegX0fyLNyJfp3Ju+CuuBvju+Me1Gfjnqy3F+4dHe7vPrgv 9f1wf5zD/APIB+JBOBcPZj1DosPzaSh5HuOHkQ/HI6jPxyPxKDya4y3g/IzB Y7kfxpEX0m88+QTGT8RJ6ifhydRPiS4Lrw/kU/E06qfjYjwDz8Sz6F+CZ7P+ OeRz6VdKPo/x89nfBZzPMrwQL+L6WYzLqV+Cl1K/jOu5Ai/HKxi/ElfiVYxf Tb4GV7HeteTr8MfUV+P1eAP7vxHX4E14M/22kNfirXgb47fjHXgn3oV3sx91 +Fe8B/9Gv3r8O/V78T6c+Fvq+EbR4fgb4ybR8fUfN4sO91NzfGh0uP/SyA8n bxE9P7z+43TcCreOLg2v/+QZuC1ux/iTcGZ0SXj9jy4Or//RyfD6j7Nwp+jw X+f9+fn/9z9dcFfqu+Hu1Pcgz8Y9qe+Fe1Pfh7wv7of74xw8YP/+NXgg+SCc y3oGkw+h31Cch4fh4fvPb4NH4Hw8cv/10eBReDT1BeRj8Fg8Dhfuv14bPH7/ /dPgCXji/vutwUnySfSbjKfgItYzlXwano6L8QzWMxPPwiUcz2w8Z//zqcFz cSmeh+fvf/41eMH+52WDy/BC6heRL8bl1C9h/qUczzKOvwIvZz9XcH5Wcv1V kq/Cq/Ea+lcx/1q8Dn/M+Gr6r8cbWO9GXMP9uon7fzPeQn0t/baSbyPfzvp2 4J0c3y7y3eR1+Ff2aw/5b/SrJ/+d8XvxPuoT81LdCDeODuObkDfFzahvHh2u 30NxWnS43g+PDvdDi+hwv7TE6dHLGtwKt6a+Df0zmL8t62vH8ZyEM6lvjztQ 35E8i/k7sb7OuAvH0xV3iw7H3x33wNm4J+6Fe0dXNLhP9PIG9yXvR94f50Sv aPCA6JUNHkg+iPG55IPJh+ChrC8PD8PD8Qj65eOReBTjR+OC6LC/Y/DY6HB+ x5EX4vF4AvNNxEnWOwlPjg77O4W8iHwq+TQ8nfpiPCM6nP+Z0ZUNnkVeQj47 +p0Gz4le1eC5uBTPw/Ppt4C8jHwhXsT6FpOX028JXoqX4Qq8nP4ryFdGv9vg SrwKr8ZrcBVeG726wevIP2Y91ax3PfkGxm+kfw35JvLN0e81eAt5Ld5Kv22s Zzvr3YF3Ur+LfDeuiw7Xx6/R4freg3+LDvdLfXS4v37He6nfR7/Es6luhBvj JtFhvU1xM9yc8YfiNHx4dFhvC9wyOhxfOm6FW+M20eH5l0HeFrdj/pPIM3F7 3IH5OkaH53MWeSfyzrhLdHz9x91w9+jw/qEHzqa+J+6Fe+M+rK8v7sfx9cc5 1A8gH4gH4Vz2ezDnawj5UMbnMf8wPJz6EeT55CPJR+HR1BewvjF4LPXjyAvJ x5NPYD8mcr8l8STqJ9NvCi6ifir9pvE8mI6Lo8PzcAb5TPJZuIT62XgO65lL Xorn4fnMtwCXRYfn/ULyReSLo8PrSzleEh1ej5biZbiC8cvxCuZfSV5Jv1XR 4fVyNfka8iq8Fq/DH+Nq+q8n34A3Rq9pcA35Jvptxluor8Vb8TbGbyffwXp2 4l14N66j36/ke/Bv0e83uJ78d/rtZf37cGJ+qhvhxrgJboqbRYf5m+NDcRrj DydvQd6SPD06HH8r3Bq3YXwGbkt9O3xSdNj/zOiq8PqPO+COOIt+ncg7k3fB XanvhrtHfxBe/3E29T3p34vj703eh/F9cT/cn/E5eAAeiAfhXPoPJh9CPpQ8 Dw/Dw/EInI9H4lHMP5q8gHwMHovH4ULO53g8IXpteP3HSeon4cnUT4n+MLz+ 46nUT8PT6V+MZ1A/E8+ivoT9mE0+B8/FpYyfh+dzvhZwP5ThhXgRXozLef4s wUt5/i3DFXg5XkG/lcxfiVfh1Rz/GvIq+q9l/nXkH+Nq+q3HG/BGXMP6NpFv Jt9CXst6trL+bXh7dHj/tCM6vL/aGR3e/+2KDu8nd+M6/Cveg3+jfz3+nfq9 eB/1iedS80bPpeaNcRPcFDfDzaPD/hwaHfYvLTq+/uMWuCVOx63o35r525Bn ML5tdLge2kXH1//o+PpP3p7xHcg7kmeRd6J/Z9wlOjyfuuJuuDvugbNxT+br Rd6bvA/r70vej7w/zsED8EA8iP3LJR+Mh7CeoeR5eBj9h+MR1OfjkXgUHk2/ AjyG63EsHsf1W4jHUz+BfCJOcn9M4v6djKfgIp4fU8mnkU+PDl+PF0fH7/9H h+8vzMSzosP3O0rwbDwHz6VfKfPPY33zyReQl5EvJF+EF7Mf5XgJXsp+LuN8 VeDlnN8VeCWu5HpZRb6afA15FV7Letax/o9xNfXr8Qb6b8Q1rG8T3oy3ML6W +bbibdRvJ9/B8ezEu/BuXMf5/pV8D/lvuJ7r63e8F+/DiQWpbhQd+jcmb4Kb 4maMb44PxWmMP5y8BXlLnE59K/LW0eH+bIMzqG9Lv3b4pOhwfjLJ25N3wB1x FuM7kXeODtdbF9yV+m7k3XEPnM34nrgX9b3J++C+HE8/9rc/+5+DB0SH5/NA 8kHkudHh+T44Ony/fQgeivPwMDwcj4gOP//IJx/JekaRj8YF1I8hH4vH4UI8 nn4T8EScZD8n4cl4CuOL8FQ8jfVNZ/+K8Qw8k/GzyEuiw8+zZkeH36eZg+dS X0q/ecw3n3wBLsML8SL6Lcbl1C9hfUvJl+EKvByvoN9KXMn+rIoOv5+0mnwN rqJ+LV4XHX//D1dHh98fW483RIffd9sYHX//D2/Cm/GW6Pj7f3gr9dvIt5Pv IN8ZHX7fcFd0+G83rqP+V/I9+Lf99W//3//U49/3r6fBe8n30S/xfKobRYf1 NMZNcFPcDDePDvtzKE6j/nDW0wK3xOm41f5+DcfbGrfZP3+DM6LD7++3xe3w SYzPJG+PO+COjM8i70TeGXfBXTm+buxHd/a3B/ufjXtS3wv3pr7P/vWH13/c Lzrcv/1xDh4QHe7/gXgQ9bnMNxgPwUNZfx7HN4z9G04+gvH59B+JR7He0bgA j8Fj6TcOF7Ke8XgCnsj4JJ7E/JPZ/ym4iPqp9JvG/NPJi8ln4Jns/yzOTwnX /2w8h/q55KV4Hp7P+AWsp4z1LiRfRL6YvJx8CfuzlHwZrqB+OfkKvJL6Ss7f Ks7varwGV+G19FuHP6a+Gq+nfgPeSH0N+SaObzPewn7Ucj624m14O9fDDq6X nXgX9btxHfW/ku9h/t9wPf4d78X7OP5EWXDYn0a4MW4SHcY3xc2i4+s/PhSn 4cOjw38tcMv99Q37lb5//ga3wq2pb4MzcFvGt8Mn4Uzcnn4dyDuSZ5F3wp2p 70LeFXfbf74a3B33oD4b98S9cG/ch/X1xf2o7896cqLD31MOwAOjw99jDsK5 1A+m/xA8lPo8PIz64XgEzscj6TeKfDT7UUA+hvFjycfhQjye8RPwRJzEk+g3 GU9h/UXkU/E0PJ35ivEM6mcy3yzyEvLZ5HPoP5e8lHwens/1tgCX4YWMX4QX 43LGLyFfynqX4Qrql+MV9F+JK/Gq6PD316ujw99nr4kOf+9dRb6WfF10+Pvx j8mryddHx7//jw5/r76RvIZ8E95M/RZci7eynm2sdzv5DryT+l14N/tVh3/F ezgfv5HX49/xXryP9SReSF1/o+iwH43Jm0THv/8nb8b45vhQ6tPod3h0/Pwf 3JL6dNwKt2Z8G/IM1tMWt2P9J+HM6HC9tccdosPnaXSMDp+flUXeCXfGXXBX 3A13Z74e0eHzzrLJe+JeuDfj+5D3Je+H++Oc6PB5dAPwwOjK8PqPc6kfjIfg ocyfh4fh4RzfCPJ8PBKPYv7RuACPYfxY8nEcfyEeT/0E+k0kT5JPwpPxFPaj iOttKp6GpzO+mP4z8EzWO4vjLcGzo8PnIc7Bc6kvxfPwfMYviI6fTxkdPt9x IV6EFzO+nHxJ9Orw+o+X4Yro8HmUy/EK6lfiSuZfRb4ar6G+Cq/l+Nbhj6mv pv961r8hOnz+5sbo8HmoNeSbyDeTb6F/Ld5K/Tby7XgH3sn4XXg366vDv+I9 +Ddcj3/He6PD58vuiw6fl5tYmJo3ig7jG0eH9TfBTaPD8TeLDuezeXQ434eS pzH+cPq3IG+J03Er3Jr52uAM6tsyfzv24yScSX173AF3xFnM34n1dcZdcFfc jX7dyXuQZ7OenhxfL/LeuA/9+uJ+1PfHOXgAHogH4VzmG0w+hHwoeR4eRv1w 9m8E13N+dHi+jYwOrxejyEczvoD+Y/BYPI71FbL+8XgC53Mi93uSfBKejKfQ v4h8Kvk01jud4ynGM6ifSb9ZzFfC8cyODs+3OdHh88Ln4tLo8Pni8/D86PB5 5AtwWXT4/PKFeBFeHB0+v70cL8FL8TL6VbCe5XgFXokrOd5VeDVew/gqvBav wx/jao5nPd5A/UZcgzcxfjP5FlxL/Va8jfrt7McO8p14F/124zrqf6X/Hvwb rmf87/Tfi/fhxIup4xvhxtFhvia4aXS4n5pFh/uveXS4Pw8lT2P84eQtyFvi dNwKt8ZtcAbH0xa3Yz9OwpnRYT/b4w7Ud6R/FuvpxPF3xl1wV9wNd8c9cDbu iXtx/nrjPrgv7hcdntf9cQ4ewPiBeBDryyUfTL8h5ENxHvXD8HDqR+B81jMS j+L8jsYF1I+h/1g8jvUV4vF4Ap6Ik3gS801mfVNwEZ7K+Gl4OvXF7McM8pl4 FvUl3E+zud/mcD/OjQ6vp6Xk8/B8vID+ZXghXoQX43K8hPmW4mXUV+DleAVe Sb9KvIr61eRryKvwWryO8/MxrsbrOb8b8EauhxryTXgz/beQ19JvK/k2vB3v wDvpt4t8N66j/lfyPaz/N1xP/e/028v9tw8nFqW6EW6Mm+CmuFl0fP2PDus7 FKfhw3EL+rXE6bgVbo3b0D8jOr7+43bR8fUfZ+L2uEN0uH86RofnTxbuRH1n 3AV3Zb5urL877oGzcU/ci/3qzX72Ie9L3o/ro390fP2Pjl//kw9k/CDyXDyY fkPIh+I86oeRD2f+ERxfPsc/Eo/Co9nfAvIxeCz14zjfheTj8QQ8ESfpN4nr azKegou4fqfiaXg6LqbfDPKZ0eHfQ5uFS6ifTT6HfC4uZf557Md8vACX4YV4 EV6My5l/CV6KlzG+Ai/n/K7gelqJK6lfhVfTfw2uwmsZv475PsbVeD331wa8 kfuxhvt1E95M/RZcS/+teBvejnfgnXgX3s3x1uFf8R728zdcj3/He/E+nHgp db5GuDFugpviZrg5PjQ6/vuB5IeTtyBvGR1f/3Gr6HB+W5O3wRnUt8Xt8EnR 8fv/uD3ugDtGx+//4064M+5Cv67k3aLj9/9xD+qzcU/cC/fGfXBf5uuH+1Of Qz6AfCAeRH0uHoyHMH4ozsPDGD+cfAR5Ph5J/Sjy0bgAj8Fj8ThcGB2+fz4e T8ATcRJPig4/n5hMPgUX4amsbxr5dFyMZ+CZeBbrK4kOP1+ZjefgudHh5zOl eF50/Pk/+QJchhfiRXgx6yknX0K+lONdRl7B+OWsfwVeGR1+n6IyOvy+xSq8 Gq/BVdHh9znWkq/DH+NqvJ71bWD9G3EN3sR+bGa/tuBa9ncr3kb9dryD+p14 F96N6/CvXO978G/U15P/jvfifYxPLE49nka4MW6Cm0aH/W6Gm1N/KHkaPjw6 nM8W5C3J03Gr6HD9tMZtqM8gb4vb4ZNwJm5P/w6svyP7kYU74c64S3Q4f12j 48//cXfqe9Avm/X1ZP29OL7e5H1wX9yP8f2j4+//kQ9g/EA8COfiwfQbQj6U 483Dw/Bw9msEzscj8Sg8mv4FeAwei8fhQjweT8ATcZL1TcKT8RRcxPU1letv WnT8+T95MZ5Bv5nMNwuXcDyz8Rw8F5fieXg+XsD8Zax3IV7E8S1mP8rJl+Cl eBmuwMvxCtazkvVX4lXUr8Zr6F/F8awlX0f+cXR4f1wdHb6eWR8dvt7agDfi GsZvwpup30Jey3q2st5t5NvxDup34l3U72b+Ovwr691D/huup/736Pj9f/J9 5ImXg+PX/7jxy6n1TcibkjeLDvM1jw7rPRSnUX84eYvosH8tydNxK9ya/m1w Bm7L+HbMfxLOjA7nuz3uQH1H8izcKTr+/h/uEh1f/8m7kXfHPajPxj2p7xUd ng+9yfvgvtT3I++Pc/AA1jOQfBD9c/Hg6PC8HoKHRof3A3l4WHR4/zCcfATO xyPxKPqNxgV4DB5Lv3EcTyH5eDyBfhNxEk9i/GTyKbiI+qnk0/B0XIxn4Jl4 Fi5h/tl4Dp7L+FLyeXg+9Quiw/vXsujw9e9CvAgvxuXR4e8jluCl0fHf/8UV 1C/HK5hvJeuvJF/F+NXMtyY6/P1TFflaxq/DH+Nq5l+PN+CNjK/Bm6jfzPFu Ia8l38r53oa34x14J95F/93kdfhXrr89+Dfq6+n/O/levI9+ifJUH1CeWt+I /MDo8PxpTH4QeRPyP9C/KfkfyZtFh+M7ODqcv+bkh5Afig/DafhP+HD6H4Fb UH9kdLgeW5Ifxfh0fDTH3wofQ31rfCz1bciPI8/Ax+O2nJ8TcDt8Iuf/JHwy zmT8Kbg9PpX1dcCncbwd8emcjyx8BvWdyM8k70x+Fue/S3R4Pp5N3pXx5+Bu zHcu7o7PY3wP5jsfZ+MLWG/P6PC8vzA6/H1sr+jw97QX4d744ujw97h9yC+h f1/mv5S8H74M98eXM18O+RV4APVX4oHUX4UHUX81zqX+GvLB+Fo8hPHX4aH4 epzH/t5APgzfiIeznpvwCOpvJs/Ht+CRXC+3Roe/hx8VHT6/4Lbo8PkHo/Ht 0eHzEwrwHdSPwXfisYy/Kzp8HsU4fDcuZPw9eDz1f44On68yAf8FT8R/jQ7/ JcnvZb5JrOc+PBnfj6fgB3AR8z2Ip+KHGD+N8/Ew+XT8CC6m/6N4BvWPMd9M /Dj1s/AT1Jdw/T5JPpvxT+E5+Gk8l37P4FL8N9YzDz9L/Xz8HPMvIH8el+EX 8EL8Iv0X4ZfwYvwy/crxK3gJ41/FS/FrXD/LuL9e5/6rIH+DfDl+E6/Ab9Fv JX4bVzL+HbwKv8v41Rzve+zHGvbzfa6nKvwBXos/5Pm/Dn9E/cf4E1zN+E/x euo/wxvw53gjx/8FeQ3+kvpN+Cu8mfFf4y0czze4Fn+Lt+LveD3eRv493o5/ wDtY7494J8f7E/ku+v2Md+NfWH9ddHi/8Xf8K/X/wHvwPxn/W3R4//MvXE/9 v/Hv9P8P3ov/y/h9+H848UqqD8CN8IHRYb7G+CDcBP8hOpyPptHhfP4xOpzv Zvhg3Dw6PG8OiQ7Po0PxYTgN/wkfjo+IDs+/FtHhvyP35+f/3/+0xEft79fg 9P3rb/DR+4+vwa3wMbg1PpZ+bfBxzJ+Bj8dt8Qn0a4dPxCfhk3Em6z8Ft6f+ VNbTgf09jbwj+emcryx8Bu6Ez6RfZ3wW7sL4s3FX6s9h/d3wuexHd3we9T3I z8fZ+ALck/NxIe6FL8K98cX070N+Ce67//nQ4EvJ++HLqO9Pfjnz5+Ar2L8B 5FfigfS/inwQ+dXkueTXsP7B+5+XDb6WfAjjryMfiq/HefS/gXwYvpH5huOb qB+Bb6Y+n/wWPBLfikfh2zie0fh2XLD/9anBd5CPwXfisfguPA7fjQuZ/x48 Hv8ZT6DfX9iPifivOMn5uBdP4nq9D0/mfrmffAr9HiAvYvyD5FMZ/xD5NMY/ jKdT/wh5Mc/fR/EMntePkc8kfxzPwk/sH9/w/qGE/En6zyZ/ivFz8NP7j6/B c/Ez1JfS/294Hut5Fs/Hz+EF+Hlchl/g/CzEL+JF+CX6LcYvc3zl+BW8BL/K /i3Fr7Hfy8hfxxXUv0G+nPnfxCuof4t+K/HbuBK/g1ftv38a/C5eTf17eA1+ H1fR7wO8Fn+I1+GP9j8vGvwx+SfMX83+fUq+nvGf4Q3Uf06/jZyvL3AN/hJv ot9X5Jvx19Rvwd/gWsZ/y/2ylfw7vA1/T//t+Ae8A/9Iv534J9a3C//M82A3 +S/kdbw+/J38V8b/A+/B/2S9v+F/UV+P/41/Z/x/2K+95P8l34f/hxNLUn1A dPivET6Q+sb4oOhwPzSJDvfLH3BT/EfGN8MH4+bMfwjrPRQftr++YX/T9vdr 8J/2r6fBh+MjcIvo8P7kSNyS+qNwOj4at8LHsN7WHM+xuA3HfxzOoP548rb4 BOrbsZ4TWe9J+GScyfhT6N8en0p9B/LTcEfqT8dZ+AzW14nzeWZ0eP/cOTq8 /z6LvAs+m35dme8c1tON/Fzy7vg83AOfj7PZrwtwT+ovxL2ov4i8N+u/GPdh Py7BffGluB/9LmP+/qzvcpzD9X0FHsDz5UqePwPxVdQPot/V5Lnk1+DB+Fo8 hH7XkQ/leK/HeezXDeznMPIb8XB8E/1H4JtZXz75LeQj8a0c7yh8G+djNL4d F/D6c8f+66vBY/bfvw2+c//93uCx+K7o8P3uceR340J8D+PH4z9Hh+/HT8B/ wROjw88b/kqeJL83Ovw+yKTo8PtA90WH3xeaTH4/46fQ/wHmL+L4HmQ/puKH 8DT8MJ6OH2G+YtbzKPkMxj+GZ+LHGT8LP4FLGP8k19tsrsen8Bz8NOPn0v8Z 8lLyv5HPo/+zeD5+jvEL6P88eRnjXyBfiF+kfhH5S8y3mP1+GZdT/wpeQv9X 8VLqX6P/Mvw69RX4DeqXc32+yf20Ar/F/bcSv019Jf3fwauofxevpv491r+G /H1cRb8PWO9a/CFehz/CH+NPmK+a9XzK+tfjz6jfQL/PyTcy/gvOdw3+Em/C X+HN+Gu8hfm/YX215N/irfT7Dm/jfvgeb6f+B/Id5D8y/07W+xPexf7/zPne Tf4L4+vw36n/lX7/4PVuD/4n/o3x/yKvJ/83/p31/AfvZf3/xfvw/3Di1dR+ B+BG+EDcGB8UHdbfhPwP0fHfD8J/pL4Z/Q7GzaPDfh4SHd6fHIoPw2nR8ff/ yQ8nPwK3wEcyviXrO4o8HR+NW+Fj6Nea/FjcBh/H+Ax8PPVt8Qkcb7vo8Pvn J0bHz0sgP5nxmfQ/BbfHpzK+Az4Nd2T86TiL+jNYbyd8JsfXGZ+FuzD+bObr is/B3fC5uDv9z8M98PmsLxtfEB3+3qInvhD3ig5/z3FRdPx7UXxxdPz7UXwJ 7osvpV8/8stwf3w5zokOf296BR6Ar8QD8VX0H4SvxrmMv4Z8ML4WD8HXsT9D 8fU4j/N3Ax7G+b4RD6f+JvqPwDdTn0+/W7geR+JbqR+Fb8Oj8e3MX0B+B/kY fCcei+/C4/DduBDfg8fjP7OfE/Bf8ETG/5U8yfVzL57E9Xkfnkz9/fSfgh+g voj8QdY7lfwh8mn4YTwdP4KL6f8onsF6HyOfSf44nsV+PREd/t6/BD9J/Wzy p8jn4KepnxsdPn/gGfJS/Dfq50WHz0N4Njp8XsJ8/BxeEB0+H+J5XBYdPk/i hejw+Y8Lo8Pni76IF1H/El6MX8bl9HuFfAnrexUvxa/hZfR7HVcw/xt4OX6T 8SvwW9SvxG9Hx89/x+9Eh8/HXYXfxavxe3gNfh9X4Q/w2ujweb4fkq9j/R+R f0z+Ca7GnzJ+PflneAP1n5NvZP+/wDWcvy/JN+GvqN/M9fY199cW8m/Ia/G3 eCvjvyPfhr/H2xn/A/kO/CPPj534J7wL/0y/3eS/8Dyrw3/n+fcr/gfPyz08 f/9J/hv+F67H/6bf7/g/vL7sxf/F+xj/P+ZLLE3ND8CN8IG4MT4oOvRvEh32 7w/RYb+bRofz80fcLDqcz4PJm5Mfgg+NDtffYTgN/4nxh+MjqG+Bj8Qt8VHR 4f5Nx0dT3wofg1uzvmPZnzbs93Gcjwx8PG7L+TyB890Onxgdrr+T8MnUZ5Kf Eh3eX7XHp0aH928d8GnR4f18R3w6zsJnRIevNzqRn0neGZ9FfRd8NvVd8Tm4 Gz4Xd6f/eexHD/bvfPY7G1/A+emJL8S98EVcL725ni7meuyDL8F98aVc3/24 /i/j/uiPL+d+y8FX4AH4SjwwOrw+X4UHUX8168nF1+DBjL+WfAi+Dg/F17Nf eeQ34GHMfyPHN5z8JvIR0eH90M04Pzq8f7qFfCS+lfpR5Lfh0fh2xhfgO6LD +80x5Hfisfgu5huH78aFjL8Hj8d/Zn0T8F/wRPxX+iXJ78WT8H2Mn4zvx1Pw A7gIP4insl8P4Wlcbw/j6fgRrtdirvdH8Qz8GPfPTJ5Hj/O8mkX+BHkJfhLP 5vn5FJ7D8/ZpPJfn8zPkpfhv9J9H/iyej59jvgX4eVyGX8ALeX16ES/CLzF+ MX6Z+nL8Cl7C6+ereCn1r5Evw69TX8H63sDL2d83yVfQ7y28Er/NeirJ36H/ Kvwu61mN38NrGP8+81XhD1jfWvIP6bcOf0T9x+Sf4Gr8Ketfjz+jfgP+HG/E X+Aa1vslx78pOryf+4p8M/6aflvIv6FfLf42Oryf3Iq/i46fjxcdfr70Pd6O f4gOP7/aER1+Hvcj+U78E96Ff6bfbvJfWE8d/jvH9yvH/w+8h/p/kv9G/i9c z/z/xr/j/3A8e8n/i/fh/zF/4rVUHxAd1t+I/EDcGB8UHeZrgv+Am+I/Rofj a0Z+MG6OD2E9h+LDqE9jvj9Fx5//4yOob4GPxC2Z7yjWk46Ppr4VPob+rfGx 1Leh/3E4Ax/P+W8bHZ4PJ+B20eH5cmJ0/Pofn4wz8Sn0a49Ppb4DPi06PF87 Rofn9+nkWfgM+nXCZ1Lfmf5nRcd/r4n8bNyV+nPIu+FzcXd8Hv16RIf3k+dH h/ej2dHh/ewF0eH9bk/yC3Ev6i/CvfHFuA++BPfFlzJ/P3wZ7s/xXk6eQ/8r 8AB8JR5Iv6vwIHw168nl/FxDPpjx1+IhrOc68qH4epzH+BvIh+Eb8XDG30Q+ At9MfX50+HrsFjwS34pH4dvw6Ojw9d/t0eHrxwJ8Bx7D+DvpPxbfhccx/m5c SP09eDz79WfyCfgv1E/Ef+V8JMnvpd8kfB+ejO+n3xTyB8iL8IOsbyp+iPpp 5A+TT8eP4GLGP8r9N4P78zE8Ez+OZ9HvCeYrYT1P4tns31NcP3Pw09TPxc/g Uub7G56Hn8Xz8XN4AX6e+cvwCxzPQvwi9YvIX8KL8cu4nPv9lejw/acl5K/i pfR7jXwZfp36CvwGx7ec/E28Ar/F+JXkb7OeSo73HfJV5O+Sr6b/e3gN9e+T V7HeD7he1uIP8Tr8EffXx/gTXI0/pd968s/wBvw59/tGngdf4Jro8H7tS7wJ f4U30+9r8i3k35DX4m/xVvxddPz3cfD31G9n/h/wDup/xDvxT3gX8/+Md1P/ C67Df2d9v+J/4D34n/g3rod/4Xr8b8b/Tv4frre9+L94H9fz/8gTy1J9wLLU +kb4wOhwvzbGB+EmjP8D8zWNjj//x83wwbg5PiQ6/ntD+DCcxvg/4cPxEay/ BT4St8RH4XR8NPvViv08BrfGx+I2+DicwXzH47b4BNyO9Z+IT8InMz6T9ZyC 2+NTo8PrRwd8Gu6IT8dZ0eH16QzcKTq8/p0ZHX4e1hmfhbvgs3HX6OoGn0Pe jfnPxd3xebgH/c7H2cx/AXlP+l1I3gtfRH1v8ouZr09wMvgS3DfWf9rgS8n7 kV8W8+D++HLqc8ivwAPwlXH8+gYPxFdRPwhfTX1uzIOvwYPxtXhI7PdZg6/D Q/H1OA/fgIfhG+P8wcPJb2J9Izjem8nz8S3UjyS/FY9i/tvwaHw7LsB30H8M +Z0c/1h8Fx6H7479NjS4MObB9+Dx+M94Av4L/Sfiv+Ikvpf1TuL47yOfTH4/ noIfwEX0exBPxQ/haaz/YTwdP8L4YvJH8Qz8GPs9Ez+OZ8Xxnzf4iZgHl8R8 Y4OfJJ/N+KfwHPw0noufoX8p/huex/qexfNjffBzMf+iwQtiHvw8LsMvMH4h +Yvki5j/JfLFjH8Zl+NX4viaBi+JefCreCl+DS+j3+u4gvo3yJeTvxnzLxu8 Ar+FV8bxwW+TV5K/Q76K/F28mvr3yNfg92P9pgZXkX9Avpb8Q7yO+o/wx7E+ +BNcjT/F62O/rxr8WcyDN+DP8cY4fnODvyCvof+XzL8Jf4U3M/5rvIX6b3At +/Ut+7sVf4e34e/ptx3/QP0O/CPeiX/Cu7j+fibfzfy/kNfhv+Nf6f8P7sc9 +J/U/4b/hesZ/2/y3/F/qN/L8+K/eB/+H+MTr6f6ANzo9dT5D8SN8UG4Cf3+ QN4U/5H6ZuQH4+b4EHwoPoz+afhP1B+Oj8At8JGxX3BL8qOYL538aNwKH4Nb 42NxG+Y7Dmfg43Hb2C/4BPJ25Cfik/DJODP2C6+/p+D2sT74VNwh1ofX/9Nw x1gffHrMw/ubLPIzyDvhM3Fnxp/F/F1Y79nkXcnPIe+Gz2W+7uTn0a8H+3k+ eTb5BeQ98YW4F/NfhHuz3otxH/bzkpgH98WXxvrw/rpfzIMvi3l4f98/5sGX xzx8fZeDr4j1wQNiHr/+Jx+Ir6LfIOa/Gufia1j/YHwtxzuE/bgOD8XX4zz6 3cB8w/CNrHc4+U14BP1vxvn4FtY3kvN/K9fLKK6n2/BorsfbcQG+g/FjmO9O 1jMW30X9OHw3/QvxPaxnPP4znoD/Qr+J+K+sJ4nvxZM4vvvIJ+P78RTmf4C8 CD9I/VT8EMc7jfxhPB0/govp9yiegR/DM+n3OJ5F/RM8P0t4/j5JPhs/Rf0c Xk+fJp+Ln6FfKfnfyOfhZ/F8xj9HvoD8edZbFvPgF8gX4hfpt4j8Jfotxi/H +vD+qJz8FfIl+NVYH7wUv4aXxfHh64/XcUWsD34j5uHrt+Xkb5KvwG/hlfht XInfwavi/MHvkq9mfe/hNfh9XIU/YH/Wsp8f4nX4I87Xx+Sf0L8af4rX48/w Bvw5823EX+Aa/CX9NuGv8Gb272vOzxbyb3At/pb+W1nfd3gb/p7x2/EPeAf+ kX47OZ8/4V34Z+7n3TwffuH5UYf/jn9l/D/ovwf/E/+G/8V66znef+Pf8X/w Xvr9F+/D/8OJilQfUJHavxE+MNaH89UYH0R9E/wH5msa87A/fyRvRn4wbh7r gw/Bh8b6cD4PI0+j35/w4fgIxrcgPxK3pP4o8nTyo1lvK/JjcGvqjyVvw3zH 4Qx8PG5LvxOYrx0+EZ+ET8aZ9D8Ft8en4g6s9zTckfrTcRbrOYO8E/mZuDM+ C3fBZ9O/Kz4Hd8Pncnzd8Xm4B/fT+dyP2fgC3JPxF9K/F74I92b9F+M++BLc F1+K++HLmL8/vpz6HM7PFeQDyK8kH4ivwoOY/2qci6/Bg/G1eAi+Dg/F1+M8 fAPrH0Z+Ix6Ob8Ij6Hczzse3sN8jYx58Kx6Fb8Oj8e24gPnuIB+D78RjGX8X Hkf93bgQ34PHx37x6/+Yx6//yf+CJ+K/0j+J76V+Evl95JNZz/14CvUP0K8I P8j+TcUPcb1M43p7GE/HjzC+GD+KZzD/Y3gm638cz8JP4BL250n2bzb5U+Rz 8NN4Ln4m9os//8d/o34e8z+L51P/HP0W4OdxGeNfwAvxi4xfhF+ifjF+mfpy /Apegl/FS2P/4NfIlzH/67iC+jfwcurfxCvwW5yfleRv40r8Dl6F38Wrme89 8jX4fVyFP8Br8Yd4Hf6I9XyMP8HV+FO8Hn+GNzD/5+Qbyb8gryH/Em+i/ivy zfhr6rfwvPkG1/I8+5bn31aen9+Rb8Pf0287/gHvYD0/4p34J45vF/6Z/diN f+H+quP+/Tv+lfp/4D34n/g3+v0r5uHnG/Xk/yb/nfw/eC/z/Rfvw/9jfOKN VB/wRmp9I3xgrA/72Zj8INyE/n/ATan/I25G/cHkzVnPIeSHMv4w8jTyP5Ef jo/ALZj/SNwSH8X4dHw0boWPwa3xsbgNx3dczMP1lYGPx23j+OATcDvqT8Qn UX8yeSY+Bbdn/Km4A/WnkXfEp+MsfAbuhM/EnfFZrKcLPht3Zfw5uBs+F3fH 59G/Bz4fZ+ML6NeT/ELcC1+Ee+OLcR98Ce6LL8X9uL4vw/3x5TiH++UKPID6 K/FAfBUehK/GuRzPNeSDWc+1PE+G4OvwUMZfj/OY7wbyYfhG6ofjm/AIfDPO p/8teCS+FY/Ct9F/NPt7O3kBvgOPwXfisfgu5huH7+b+KsT3UD8e/xlPYP6/ sB8T8V+pT9LvXjwJ34cn4/vxFOZ7gLyI/EHyqezPQ3gafpjx0/EjuJj5H8Uz 2L/H8Ez8OJ5FvydwCX4Sz2a9T+E5+GnGz8XP4FLG/w3PY3+fxfNjffh9n+fw Avx8HB9chl/AC+P48PtUL+JF+CW8GL+My/EreAnreZV8KX4NL4vjg1+Pefh9 sgr8BvXL8Zt4BePfinnwSvK3cSV+h/6ryN/Fq6l/D6/B7+Mq/AFeiz/E6/BH +GPW+wn7U40/jfXh7wvXk3+GN+DPmW8j/gLXsN4v8Sb8FdfXZq7Pr/EW/A3j a/G3eCvzf4e3Uf893o5/wDvwj3gn6/+JfBfr+Rnvxr/gOvx3/Cvn6x94D/X/ xL9R/y+ul3r8b+p/x//Be/F/8T7W8z/2L7E8NT9geWreKOZh/w/EjfFBjG+C /4Cb4j/iZvhg5mse8+BD8KHUH4bTqP8T+eH4CNbTAh9JfUv6H0WezvijcSvO zzG4NT42jg/nvw35ceQZ5Mczf1t8AvXt8In4JMafjDPZj1PI29PvVNbfgfw0 3BGfTv8sfAbuhM+kX2fys3AX6s9m/V3xObE+uBv5ufTrjs9j/h7s7/k4m/oL cE/6X8h6erHei8h744vp1wdfwvx98aXU98OX4f6Mv5w8B1/BegeQX4kH4qvw IHw1zmW+a/BgfC0egq/DQ/H1OI/13MB+DcM34uH4JjyC6+1mnE/9LaxnJOu9 FY/Ct+HR+Hau1wJ8Bx7D+DtZ31h8F8czDt/N8Rfy+nAP+XjyP+MJsT7+/T+e GOvj3//jJL6X8ZOY7z48mfXej6fgBxhfxHwPkk9l/EN4Gvv7MJ6OH8HF+FE8 g/keY30zWf/jeBZ+ApfgJ/Fs/BSew/l7Gs/Fz8Tx8ev/mMef/+N51D9LPj/m 8ef/eEGsD34el+EX4vj49//4ReoXMd9L5IvJX8bl1L+Cl1D/Kl7KfrxGvox+ r5NXkL9Bvhy/iVfgt1jPSs7f21wPlfgdvAq/S7/V+D3mX4Pfx1X4A45nLfmH 5OvwR/hjxn+Cq/GnHM96/Bn7sYH78XOeDxt5fnzB86WG/EvyTfgrvJnXn6/x FvwNrsXf8nq3FX+Ht+Hvef3cHvPw/ZQfyHfgH+m3E//Eenfhn9mf3eS/0K+O /O/kv+J/4D2M/yf+jfX8i/Ndz/Xyb/w7/g/X217y/+J9+H/Mn3gzdX0H4EZv ph7PgTEPx98YHxTrg5uQ/wE3pf8fmb8Z+cHkzfEh1B/KfIeRp+E/4cMZfwRu gY/ELel3FE7HRzO+Ffkx5K3JjyVvw/k4Dmfg43Fb+p1A3i7m4f4+MebBJ+GT o8PnwWXiU6hvT34q7hAdPm/uNNwxOnye3+nR8fP/8Bm4Ez4zOnweYWfys5iv C/nZjO+Kz8Hd8LnR4fMTu0eHz1c8Lzp8XmUPfH50+PzLbHxBdPg81J7R4fNU L8S98EW4d3T49xIuJu9Dfgnuiy/F/fBl9O9PfjnOwVcwfgD7cSUeiK9i/CD6 Xx0d/j2JXHwNHsz4a+k/hPmvIx+Kr6dfHr6B+YfhG6kfTv+b8AjWdzN5Pv1u wSPxrYwfRX4b+Wjmv53rvQDfQf0YfCcey3x3kY+j/93kheT3kI/Hf2a+Cfgv 7MdE/FfOZ5L8XvJJ0eHfX7kPT6b+fjyF+geiw78fUxQd/r2ZB/FU6h8inxYd /v2ah/F0/Ajji+n/KPkM/BieiR/Hs/ATuAQ/iWfjpzieOfhp6ueSP4NLqf8b nkf9s3g+fo7xC8ifx2XR4d8regEvjA7/vtGL0RUNXkT+Evni6GUNfpm8nPyV 6KUNXhK9pMGv4qX4NcYvw6/jCvwG61lO/iZegd9iPSvx29RXkr8TXd7gVdGL Gvwu+Wr8Hl5D//dxFf6A9a1lPz7E66j/CH9M/0/Iq/Gn1K/Hn1G/AX9O/Uby L3AN/pLj24S/on4z+dfkW1jPN7iW8/Ut3oq/w9u4Pr4n3858P+Ad+Ee8k+P5 iXwX+c/ku8l/Ia8j/zv7+Sv+B96D/8nz5Tf8L1zP8+vf5L/j/zDfXtb/X7wP /w8nVqT6ANwIH4gb44NwE/wH3DQ6nI8/kjfDB+Pm+BD6HYoPw2nR4Xr9Ez4c H4Fb0O9I3BIfhdPx0fRvRX4Mx9s6OlwPx+I20eH6OY48g37Hk7dl/AnR4Xpt h0/EJzH+ZPpn4lNwe9Z3KnkHfBrzdWQ9p0eH9xdZ+AzcCZ8ZHd6/dMZnUd+F /Gzyrvgc1tsNn4u7M/483AOfz3qyyS/APfGFuBe+iPX1xhdzfvqQX0K/vvhS 1t8vOrzfvCw6vJ/vjy+PDl+P5JBfgQdQfyX5QOa/Cg/CV7P+XHwNxzuY/bmW fAjjryMfyvjrcR71N+Bh9L8RD6f+JvqPwDdTn0+/W/BIfCsehW+j/2jmvx0X cD3egcfgOxk/lvwuPI7n0924kOfbPeTj6fdn8gn4L3gi4/+KkxzPvXgS+3kf nozvx1M4Pw/gIvwgnkq/h8in4YfxdO6/R7h/i8kfJZ+BH+N5MBM/Tv0s+j9B XoKfxLMZ/xSeg5/Gc+n3DHkp+/U38nmMfxbPx8+xHwvIn8dl+AW8kH4v4kX4 JcYvxi/jcsa/Eh2+f7UkOnx/7FW8NDp8/+018mXkr0eH7+dV4DeoX47fxCsY /xZeid9mfCXrfYfjX0X+Lvlq/B5eg9+nXxX+AK9l/Iecr3Wcz4/IP8af4Gr8 KV6PP2M9G/DneCPjv2C9NfhLvAl/hTdzv37N/byF5+c3PN9reT34lnwr478j 34a/p992/APegX/EO/FPeBfz/4x3s/5fcB3+O/4V/4P93UP+T+b/Df8L1+N/ c3y/4//w+r4X/xfvY/z/mC/xVur6D3grNW+ED4wO/RpHh/kOIm9C/gfypvT/ I27G+g4mb06/Q/ChzH8YTosO77/+hA/HR+AW+EjcEh+F01nP0eSt8DG4dXT4 /sCxuA31x5Fn4ONxW8afgNtRf2J0+P7FSfhknBkdvn94Cm6PT8UdossafBru SP3p5FnR8xt8Bu4UXdrgM3Hn6JIGnxVd3OAu+GzcFZ9Dv27Mdy7r647P4/h6 4PNxNuMvIO+JL8S98EXsd2/O98VcL324ni4h70t+Ke6HL2N8f/LLcQ733xV4 AL4SD8RX4UH46ujwfMrF10SH59tgfG10/PofX4eHMv568jx8A/XDWN+NeDj1 N+ER+GbG55PfgkfiW/EofBsezXy3cz4KuB7uwGPwnXgsvov+48jv5vosxPfg 8fjPeAL+C/fbRPxXnGT8vXgS9ffhyfh+PAU/gIvwg8w/FT+Ep+GH8XT8CPMV 40ejw/NrRnR4vj1GPhM/jmfR/wlcgp/Es/FTHM8crq+n8Vzqn8Gl+G/MNw8/ i+fj5/AC/DwuY79eYL8X4hejw+vVIvwSXoxfxuX4FbwEv8p6luLX8DL8Osdb wX68wflYjt/k/K6IDs+bt8hX4rdxJePfwauiw/PzXfLV5O/hNfh9XBUdntcf kK/FH1K/jvwj1vcx/gRXM/5TvJ76z/AG9vNz8o30+wLX4C85vk3kX5Fvxl/j LYz/Btey3m/xVo7vO67Hbfh7vJ3xP9B/B/mP5DvJfyLfRf4z8+/Gv3C/1eG/ 41+5X/9Bvgf/k/l+I/8X/erxv/Hv0eH58Z/o8Dzai/9L/T78P5xYmTrfAdFh vY3wgdQ3xgdR3wT/ATdl/B9xM3ww45vjQ6Ljz/+jw/VxWHT8+j863A9/ig73 z+H4COpb4COpb0l+FE6n/ujo8H68FT4Gt8bH4jb4OJzB/MeTt40O3+86AbfD J+KT8MnR4fuVmeSn4Pb4VNbXgfWfhjvi09n/LHwG10cnfCbuzPizmK8LPpv1 d+X4zsHd8LnsX3fy88h74PNxNr4gOnx/uSe+EPfCF+He+GLcJzp8f/wS3Bdf yvh+rP8y8v74cvrl4Cuiw/fvB5BfiQdSf1V0+PnFIPKrGZ9Lfg35YNZ/LfkQ fB39/l9Xdx22VbG2Yfy2AxQDUTFAbAUDCxDFRMUAsTsAsUBRQEpYSNiB3S1g t2JgN4piN3YH2O23D5+Z+Y7123/sfZz7nOuee2bNs55635f96ecA/IH4g+CD M6fvi3rhe1OvD/4QfF/4UMYfBh/O/EfARzK+H/X7sz9HsZ9HwwPgY+BjOW8D 8YPwg/HH4Yfw+B2KH0Z+OH4E/nj8SOqPgit4NPenE+Ax3M/Gcv8bB4+HT4RP gk+m/in4U+HT4NO5P58Bn8n4s5hvAnw248/Bn4s/D38++3cB+30h/iLyF1P/ Evyl8GWMvxy+Ar6S/FXw1fR3Df1fy3m7Dp4IT4Inc36v5/F6A/5G8jfRz83w LYy/lXq34W+H76DenfBd7M/d+Hvgexk/BX8ffD/8AP09yHqm4h+CH4YfgR+l 3mNcj8fxT8BPMv4p+Gn4GfhZ+Dl4Gvf35+EXeD6YDr8IvwTPgF+GX2G+V+HX 4NfhN+A34bfgt5n/Hfy78Hvs1/vwTPgD8h/iP+L6fsx5+QT+lPP6Gef9c/wX 8JfwV9T/Gv4G/pb8d/D3jJ/F+mbDPzD+R+r9BP8M/wL/Sr3f4N/hP+jnT67P X/Df8D/wv9SLh+s8BzwnPFfm1N/cmdP65smcrve8+Png+Rm/QOb8/T++Cb4p fiF4YcY3gxeBF82cnh8Xwy+Ob45fgvlb4JeEl2L80viW8DKMXxa/HP0tD7eC W5Nfgfpt4BW5fivhV6beKvCq8Grw6vS3Bn5N5msLt2P8WvDajF8HXhduD69H vfXhDRi/IX4j1tcB7pg5vV7thN8Y3zlz+nxrE3jTzOnzsC74zfCb47dgvi3h rehva3xX/Db4beHt4G7w9tTbAd6R8Tuxnu5wD9a/c+b0eWLPzOnzx13wu8K7 wbuT3wPek/F70c/e+H3I7wvvB++fOX3ee0Dm9HnxgfiDyB/M/L3orze+D3wI 9frCh8KHwYfDR1D/SHw/uD/jj4KPZvwA+Bj4WHggPIj9HMx+H4cfQn4oPAwe Tn4EfDzzjcSPwlf40cx3Avs1hus/Fj+O/Hj4RPgk+GTqnYI/FX8a/nT4DPhM 8mfBE1jf2dwvzoHPhc+Dz6feBfCFzH8R/mL8JfhL8ZfhL6efK+Ar4au4f17N /fUa+Fr4OvIT4UmMn8zz4/U8f96QOb3euBF/E/mb8bfAt8K3wbfDd1D/Tvxd 8N3wPfC91JsC3wffz349gH8Qngo/BD9MvUe4Po/iH8M/jn8C/yT+KeZ/mv14 Bv8s/Bzjp3EenodfyJzen0yHX4Rfypze38zInH/+P3N6f/xK5vR+/NXM6fPo 1+DXGf8G/Cb8FvO9jX8H/y78Hvw+65nJej+AP2R/PsJ/TL1P4E8Z/xn+c/wX 8JfwV/DX1PsG/pb1fwd/z/hZ8Gzm+wH+kf35Cf8z/hf4V87rb5zv3/F/kP8T /xf+b/r5h/X9i49HalwVTuPnwBdO882ZudFPBcdcmRvrq+bCz12vV8ExD36e ej8xL/3BMR8ejvnx8zP/Ang4FqzvVwVHE/JN8E3xcCyEh2Nh9mdh9rcZ16cZ 61sEv0j9esWimRv346pw/nwED8fi9fFV4fR80bzOVXPqL0Eejhb4FtRfEr8k /S2FX4r80vRXOP9+Rp2rluSXIQ/HsswPx3J4OJbHw9GK+VvhW5NvTX8r4FfA t8G3wa/I/hROr09Xytx4/VqthF8ZXzi9nl8lc+P9S7UKflU8HKvhC6f3R6tn bvz+RQXHGng41qzXq+Boy/xt6b8dvh1+LTwca9fHV3Csg1+H/LrUX5d8e/Lt ya9Hfj3y63N91+f6b8D5gWNDfOF0/jbi/G9EvgO+A74j9TtSvxPnvxN+Yzwc nZm/M34T8ptQf1PycHSh/y7s72Z4ODanPhxb4Leg/y3pH46t8HBszfMLHF3J d6W/begPjm3x27L+7fBwdGP+bvjtycOxA34H5t+R6wPHTvidqN+d+t3J9yBf OH8eWudqZ/I9ycOxCx6OXam/K343+odjdzwce+D3oP6e9AfHXng49qb/vfH7 kIdjX/KF0+c5+3F/LZzur/tzf4XjAPJwHEgejoPIH0R/B3P/hqMXvhfz96Z+ b/J92B84DqE+HH3J98UfiofjMK5f4XS+Dud8wXEEvnC6vx2JPxLfD98P3x/f H38U9yc4jiYPxwD8gPp+xDHsDxzHsr9wDMQP5PoM4nzAMRhfOJ3f4zKn119w DOH8D8EPJT+U+sPIDyM/PHN6fTgcP4L6I+rj43j88fiR1IdjVObG7xNXhdPv G1eZq4av8KPJj6b+Ccx/Ausbg4djLPXH4seRh2M8+fH0fyLrK5z+cxIejpPr 4ys4TiEPx6nkT8WfVu+3Og1/OvnCVZf//ueMkv+Pq8IXNvyZZb8avvC1DX9W nSs4JpCfQP2z61zBcQ79nUP/59a5guM81g/H+fjzyV9A/QvIX8j1g+MiPBwX 4y+uX8+4BA/HpfR/Kf4yzgccl+MvJ38F9eG4Eg/HVfQPx9Xk4bgGfw3X51qu T+HGeYnrOF9wTOR8TaT+JOafRP+TWd9k9vd69rdw+vsMN3D/geNG7l+F8++v 8vwBx808v8FxC3k4biUPx23159sKjtt5fobjDvJw3ImH4y78XdS/O3N6fXE3 /h76g+Ne6sMxhf2Zwv7dx/7CcT95OB4gD8eDeDimUn8q/iE8HA/jH2b9j7A/ cDyKh+Mx6j+GfxwPxxN4OJ5kf+B4ivxT+Kfrj7cKjmd4fMLxLB6O53j8wzEN P4388/QHxwusH47p+OnkX6R/OF7CwzGD/mYw/8vs/8tc/1c4P3C8in+V6/8a 9eF4nTwcb+ALp/vDm9xf4HgLD8fbmdP7IzjeIQ/Hu3g43uP+9h7+ffLvM/9M +ptJ/gN84fT+8MM6V4Xz34fDf8T8H9MfHJ/g4fiU9X+K/4z8Z8z/Of5z/Bes H44v8V+yP1+x/q/If838X+O/IQ/Ht8wPx3f479if7/FwzMLPYv7Z9Deb6/8D /gfq/8j6f8T/xPw/Mf/P1P+Z/C94OH6lPhy/kf+N9f1e56pw+n70jzpXcPyZ OX2/+if+L3zh/H1+5vT9Lxz/1MdX/+D/JQ//73/rvnD+/h8Px5z18RUcc2Vu /PxFBcfc5OGYh/nnwc+Lh2M+6sMxPx6OBei/cPp5nAXrXMHRJHPj54GqJvim 1G/K/Avh4ViY+ReuzxfNMjf+/msFxyL4RcgvSv1F8Yvh4VgcD0dzfHPWtwTr h6MFvgX1l2T/4ViK+Zei/tKcDzhaMn9L/DLkl6H+sng4luN8w7E8fnnyrajf ivlb41vjV2B9K1C/Db4N+78ifkXyKzE/HCuTh2MV/CrkV2V/VsWvxvyrsb+r 41fHr4GHY008HG3pvy2+Hfl2rH8t8nCszfleG78O9eFYl/pwtCffnv1ZL3N6 floPvz7XZ/36+NiAPBwb4jek/kbU34h8B3wH1t8R35F8J3wn/Mb0tzH9dWb/ OuM3of4m1N8UD0cXPBybMT8cm+M35/pvwfkonB7/W3L+4NgKD8fW1IejK/11 ZX3bsP9wbIuHYzv8dtTvxv7CsT15OHYgD8eO5Aunn+/bCb8T+e7Uh6MHvgf9 7Ux9OHrie1J/F+rvQn5X8nDsRh6O3bn+cOyBh2NPPBx74eHYGw/HPvh96H9f /L74/fCF89934fFROP+9Fx5fhdPrnwPxcBxUH18VTq8fD65zVTj9e6a98HD0 xvemvz7M34f+DqH/Q9ifvvi++EPxcBzG/h7G9T2c8wnHEZzvwvnfr6lzBUe/ zOn9ZT98f3zh/O/d1Lk6qj4+jmb+o/EDyMNxTOb0/hmOY/FwDMTDMQg/iPkH 4wun9/fHZU6fP8AxJHP6fGkIfii+cPp8cVjm9PkkHMPr46vh9XoxgvpwHE8e jpH4kcw/iv5Gka/wcIzGw3ECHo4x+DHMP5b+x7L+cfhx1B+fOX3+DceJzH8i /iTyhdPn5yfXuYLjFPJwnIo/lfxpmdPn74XT9wenZ07fP8JxRn18dUa9XpxJ /3CchYdjAn4C6zub9Z3N9TmH/T+H/Ln4c/HnUR+O8/FwXICH40Lmh+Mizl/h /PuZ+Ivxl3D/geNSHv+X4i/Dw3E5/vJ6P3EF/V1B/krycFxF/3BcjYfjmvr9 uCqcf38UD8d1eDgm4idSfxJ+En4yHo7reX6B4waeHwun588beX6G4ybycNzM /HDcQv9w3Mr1uRV/G3k4bmf+2/F3kIfjTjwcd9Ff4fz9P+ezcDrf9+DvqdeL e/H3kp/C438K/j7ycNzP4+t+/APk4XgQD8dUPBwP4R9i/Q+zv3A8gofjUa4f HI/hC6fz+TjnG44nOF9P4J8k/yT1n2J+OJ7Gw/EM638G/yz5Z9mf5/BwTMNP I/881w+OFzhfL+Cn0z8cL+LheAkPxwz6h+Nlrh8cr+Bf4fq+iofjNTwcr+Nf p7836B+ON/FwvIWH423mfxv/Dnk43sW/y/reo/575N/HwzETP5P9/wD/AfU/ pD84PsLD8TEejk84f3B8yuMDjs/whdP9+3Pu74XT678v6lzB8WXm9PoXjq/w X5H/mvpf1/uJb/Df4L+lfzi+w3/H/N/THxyzmH8W9WdTfzb1fyAPx4/4H5n/ p8zp/U/h9P7n58zp/Q8cv9THV3D8iofjN+oXTj8f+jv+d/J/0D8cfzI/HH/h 4fib+f/G/4P/h/7/zZz+X/h/T+B1D8ccmdPPx86Bn7P4Lv/995z4uahfOP18 7tzkC6ef/50nc/r5XzjmxRdO+fnw8+Hnr3NVOPW3AB6OBalfuPHz8tEkc/r5 ejiaZp7c8E3xC9XrVQvhF6Y+HM3whVN/i1AfjkVZHxyL4Qun/VkcD0dz5odj CfJwtMDDsWS9n6pwur5L1bkqnP6zdJ2rwinfEt8Svwx+Gfyy9AfHcvX1VHAs T355fCvycLTOnPYfjhUyp/MDRxvyhVP9FemvcFr/SpnT/wvHyng4VqnXq+BY lfyq+NXIw7E6vnD6zxrUL5zuv2tmTvdnONriC6f7e7s6V+3q42MtPBxrZ07P T3Csg4djXeqvi2+fOb2+KJx//7/OVeH8+/91ruDYAF84vb7bsM4VHBtlTq8P 4eiAh6Nj5vz7//hO5Aun+TfGw9GZ/jvjN2F+ODal/qb4LpnT6284NsucXv8X Tu8PNs+cPr/aHL9FnSs4tiS/JX4r5t8KvzX5wun7r66Z0/drXfHb1LmCY1vy cGxHHo5u9NcNvz314dgBD8eOmdP3i3DshIeje+b0/XF3fA/yPerjY2fycPTM nL4fLZy+P92lzhUcu+ILp+9/d6tzVTh9f7w7Ho498IVTfs/Mje+rq8LPNPxe da7g2Dvzsw0Pxz7kC6fvx/dl/sKpv/3qXMGxP+uD4wB84bT/B9a5guMgPBwH 4wun79d71bmCoze+cDpffepcwXFI5nR+C6fHT198X/yheDgOq4+v4DgcXzh/ /0//cByJL5z2p1/mtL9w9K+Pr/rjj8IfhT8aD8cA+huAP4b8MfXxcSy+cJp/ IOuDYxC+cDrfgzOnx8dg/HHk4RiCL5zWN5T+4RiGL5zWP7zOVeF0PkZkTudn BP74OldwjCRfOP97P3WuRtXHR4UvnL6fGl3nCo4TyBdO9cfg4RiLL5zqj8MX TvXH48fjT8TDcRLrg+Nk8nCcUp+vKpyuz6l4OE6rj6/gOD1zOj+n488gD8eZ eDjOwhdO53cCHo6z6a9wenyckzk9vgqnx+e5da7gOA9fOD3/nY+H4wJ84ZS/ EA/HRfRfOK3vYjwcl9THV3BciofjMjwclzM/HFfg4bgSD8dVzF84nY+rM6fz UTid72syp8fPNfhr8XBcVx9fFU7nbyIejkn4wqn+ZOaH4/rM6f4Axw14OG6k fuHU3030B8fNmdP+wnFL5vT4g+NWfOH8718xPxy34+G4g/XBcSe+cMrfRf3C qb+78XDcQ3047iUPxxQ8HPfh4bgfXzj/e12Z0/l4AP8g/RdO9adSf2p9fDxE /cL5733XuSqc//43vnA6v49Sv3Ca/zH6f6w+Ph6nfuH898PxhVP+SeaH4yl8 4fz3xqlfOM3/TJ2rwvn7/8zp8xc4nsPDMa1eryqc5n+e+Qun/l9gfYXT/k/P nM7HdPyLeDhewsMxAw/Hy/jC6fH7Sp0rOF7NnO5Pr+Jfo37htL7X61zB8QYe jjczp/2F4y184fzvs+MLp/rvMD8c77K+wmn979W5guP9zGn/3sfPJD+zPj4+ yJye3+H4MHN6ffAh/iPycHxMvnB6ffVJnSs4PsUXTq/fPqtz9Vl9fHyO/xz/ BR6OLzOnz2cKp89/vsqcPv+B4+v6+AqOb/BwfIuH4zvmL/xcw3+PL5zys1jf LPxs5p9dHx8/ZE77Vzjt7491rgqn15c/ZU7n5yf8z3Wu4PgFXzjN/yv+V/xv 9AfH73g4/sic1v8H/k/8n/i/6A+Ov/GF0+PrnzpXcPyLL5we//F44vz7/3WO OTKn+9Mc+Dkz59//z5zmn6vO1Vz18TE3+bmpPw8ejnnx85KfDz8f/c1Pf3As gIdjQTwcTTKn69cE3xTfFL8Q/cOxMB6OZlzfZvhF6B+ORfFwLIZfrL6eWDxz evzA0RxfON1/l6hztQT5FuRbMP+SeDiWwsOxNPPD0ZJ8S/wy5Aun++eyeDiW o/5yrG95PBytOF+tuH6tOT+t8StwfeFog29DfysyPxwrkYdjZfIrs75V2D84 VmV/4FgND8fq1F8dvwZ5ONbEr8n627K+tuTbkW9Hfi3ya5Ffm/za5NfBr4Nf l+tTON0/2nN+2nP91qO/9ehvffz65Ddgfjg2ZP4N8Rvh4eiA70C+I/11xHfC d2J9G7O/cHTmfMGxCXk4NmX+TfFdyMOxGfMXTvenzfFwbIGHY0vm3xK/Ffmt 6G9r1gdHV+p3pf421N+G+tvit2X92+Hh6IbvRn77zOn1Mxw74Aun9w871rmC YyfycHTHd6e/Hvge+J3xcPTE96S/Xei/cHp/smudKzh2w8OxO77wtIbfAw/H nvg96X8v+odjb/Jw7EMejn3ZPzj2q1+Paj/8/uT3p/4BeDgOpD84DmJ9cByM P5h8L+r3Yv7e9NebfB/q98EfgoejL/PDcSj5Q/GHkYfjcDwcR+DhOJL5j8T3 Iw9HfzwcR1EfjqPxhdPjZ0CdKziOIQ/HseThGIgv/HzDD6pzNYj8YPKDyR+X +YWGLzy94YdkfrHhC7/U8EMTVw2uhuKH4Yfhh+PhGFEfX43AH08+c8xo+JH4 kcw/ivqjqF/h4RjN/ozGn0AejjH4MfV6MZb6cIyrX68KjvH48eRPxBdunIc4 6f/Py/8BDuOAlQ== "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->{-1, 1}, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.419002969554741*^9, 3.419003166086637*^9, 3.419003308034099*^9}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztPFtsXMd119xdcinKkiyrUai3LUs0VYmSIlMS9aLWlGXJskRZiay0pl8Q SltBk2wdP0QTTV0KqcHWqbJpEsnxg0hRFrH1kQpFwQAGXZiFAaIoYYAfFozw QynAH34w7bbFAuXH6ZyZc+aeO3fu3qvUQOwgV+Du6tyZM+c1Z845M/eefPr5 Z//o608/f+H80xsefO7p8rMXzn9rw9FvPqdAuduC4LZngyCY2RDk1W9QP+UH XhfwI/q7SX/lYGFhQYMXagtBgWHVeZivabDpc072aYEbY8OIG4avTUJVd/4F 9Kn/j/yiBtMjZegfHoHhvl4YvVHz9M+pjyoMd/TCZA2CBg1rgvnZWajNXlet yjCLsOo4dARD+NvgeNTgCBZBbWYMyv1lGBgYoL8ylIfHDCkwD5VSL0xUIac7 NML8xBB0DEzom1CbglLQAZOmrUWq2tbmYVbRoP4UTfhZgNm5edMsRBow1qlK CTqGp+j+nObHDpqD7oOvA6sgB6yHk3E5uHirk8OK2nGXWhev+lZk4v/wOmXw NuqvPCpEyXMK+oPesG+jktuIFu+cptiI+iZY1B99NJsFdQ1ujI8osvoQdZ7A WuuVUaj09ymtV9mQCJ2DkqSwlBrVlLVVjREGsuPhQ6Zjjj4HX/45SZss4rjB s1L0+d7lf9Gt83awE/FGeD3Q86YajBs9JBvlbSO68kkEPCi7FdxuSIwgxEMt I2sgBMSx/bv27rRpdCy5r9Yei+jE8ZEohUdFxxwZetj8kd6/TW7uG0ex88df +wdL3/2Hf5xGXt6Ol72lGAVHUFLQLY6k9G0krlCKEa7uT+nYRB0ffsiRXiml Y9Hy1fOptVSwp554186aAtGGkzNC2+FkPAXJVE/pjeSOul+BOpW6jVbxGy8l jiDUXB0hMgKcfDgcfuOlaM6GYJGhxRCsZmaU4O6UjvkkER3K2PGJvnduraPi FaeOnKqK71S9tNhBM7VU9oAWKUdJ6bbYyjF7SzVKqfvH1uZw4uF1QPTVXXUP NiqclqRlg7MrZSgFE44j4KmKpkbT23TYk4zndvxsFh0jKtud0tEhoEn4CnLr pkNnMp4l+NkifAV1Ni12pXSk9YMJaLFaqtPRrJAFLaLekz+xnhAXlmIWBCsI duzo23ZktvpjR9+KinBnOh6cYhHPlEbA71kC3tLNlcisS4stp9vT8OQ1ARdf HNPf6i+kQRKvu7biL/99B7duyRbJvssStS2lo+OvWbpCXabD1mQ8q/CzSB1Z x7bjlpSOCnbk/jdiCwYaJ3nEzHjYzuQUjZlJu8ETxrAq1lRx8oJYwFQyAdWq DrEpcIyOt4ZgfnP6fTOA6ai8R3UKyqV+GBkdgt7yNdB5CSUcb37w9xpx39Aw DJQrcKMGBc94a/GTp3xsmbjXjGc6LtEMjfUH0D9m2k0MdMAA/dZJivo3T9kR zF6DUt81ZJ55j41LcnUN5MIzP/OSoTyDP7/5KQx2lGCcMpLpSgcEQ5M2YIYF lSp0DKtEJVSDJGUdfrLDF/YCTzz+jjQPbybUCLP//s8q5C/ZvAlzqtLQJIf/ sDAP15Wcgr5R1I8G63h4jRxWiMGu8luk8FHl8yqTUmzOMZslZJMDZEyYrpVL lFn1wcTcQuBT+Dqi6zSt0qdNoBsujbFhazDaizo34xr93+TMaPpqPwxcn4aq ypF99myGa4p4J7yQSWI2MgcbmR1t2mUYvTasTbsKHGere7Vplb8FhlUl1mqC Yu8m2KNnR/VY6jsyVjhP0WyrRKmBLdRqsEB4ckm4OS4Tobhp0WYG8DmUu40B xCfcZWosJBV8PyMe6Yn+mhSMeFjK586O5j3o/yoFfaOHzO8S+gvP/kzGX2zu KAoGq25/mTKCWP+kB5ADslXxUqK+g0ue0ZgIxiFXi1fTyRD8eB1RgYakNJr/ LoVD5iMGLrQjAv/gL9IpQQzYJaSkoE0YGUJh4R+mvzk79CU5GmeeRnbxmPA7 ggJNwGb8pdfYnMTD4S62/jaRJpi3QsKL3BfFjnlNPsYwJJChZKbDwdUN9EWk ZkP6K9RD6DkyLsFV6hHOA17ZX0kZVDgG7LzEjidz3khM6ky0ZsIiQtDgW+mD ypU9HJRdiS/B5ewT/6dFxTmQLIqo0Z9LHryNYHKyNQtlCYuOkIMXDW/JkUbA msOcj8xcwdgu6xCkl/QWogCvehTQomwoyIc9o4FyPvO4bEKCc1ODowCAXYfi PCd0wDgG9WdjLNHFSwlYyom7NBApUnbNHh587RqysKTIkc4QucBLcZKzM9uH V8GE1wpeluIhBOGsa/RgabY6Yb7RkzK+b6aRnXfdGktAuvsoeuGUL+rPQqTy xp04MHK1hCbKqzwvhhiFNJGJcd4r819fxP+NZNY6CSbTapY2B5WW/Hyon7y+ K4gVtYiotfK6wAs8/ykbCHT7sqBOE7dbalauES8lmQbTJLwk/i1D8BLPIvH1 ZIHsJpiM9lkgcq1/MSSyIHHr68RxKrppCpbirzrmtYfUfy5KPZqPkbsUK3dq EDbJHdSipAfk6So9tMbxpGNMGFndnkzYXjk4yjbCfE7OM4lxmUdNbJZ4eUli k8ESXBpJznqiTOQlAosVC/84JZZrSYNe/Yw9Ys2MfR+G3D5BP255rUNWl5SU tDSTBre409TKC0fHFZk9Ik7vLcSMnDGqXYMcQ3rsPNkP4sQqKO4fGFf1tuGI xSCV5CrO2vkTHsPKwLyeCtvsL68TiM6gNtkYL3bgPN3kQiZXG0SMZtLgGUjD NklJySEfI95QQr8Wb7yYyy0ajpZ8HgMvDgkS2eLpKl3dRqJURjxqeXicVO2W lFGjizNzsZgGlDtHm6ibMDpLMVWTIjC5PxNysdkv4Mc8xsb6wf8zokVZWNji IDNOTN3gGgL7UfUdrEfYIrJBsSmBvRZJih1zPyGUyhFzo2gqvEckFZERSovg RjPTRdpzwljNQTsNyMsquiPNFc9uUXe2ndZZm+T5IXNDHfFyOM0OE2u3zBxX kNS3Zc6hLWRuSYpuyAFJEpV9Gb6WkQLEXLTqdOKkNTGWRKBhK/5xy41KlYv6 XV5bAbl4SddEHjlIY3gfDekGbIj+XmkOMhtstvLkIU8LLriayjMbHVxRikGM ZcRwL92IzVgOGTmmU98BTxqxu2d9Kg+JW4WhcKhIwAuUiHclU74NyjrulAXn lGS0tbc55ITrRl4u4DZREEZtxbta6Jt8UGC8XaPMTVksiICVIVfUzVaHTaE8 ok2MkIt6JCYDNcExF27W6L0NiPkk3hByUwHFpGFYC5AlyxGDndGkcvYyGczU nZfoA9rsGFwPkxu/qkeLlKfcm1gWWmXe9UrtUcFpbrj4IMJbGw/ItTByqEGW xPTCxoukXM71csEyIRlFrFI2riOn/XJkzkeBXUujowi8wdGTDD2KVmBc40G7 IGs2ATDntiIatI5flk4aQr2ESxA2DgVa9FDVYGXuGGuLlKD0tqwbX6guYVqg d2SUoDyBsckaRBA9FFAMBWVMT6yoK0OxN7pVTI6SNEU7NKr4TDJnFiDmrULx ObWAyLa0zHR04+V4Y3kK946X6ubQoil2IIgnglrYjAlziiqN/k4aR9bgbCr2 FUOgPZpB5k9isLNLynu7tQ1243Ij21flZ9cl/TFXTbkdLZq6Ox8RkHGw3KGP iPQOuoHukm8sS5bwPoLhmcYPJmfDSaJg//bhx/acY4NoxzDev/OVOGRcYILP xa4W1d8V2L/3MhwyaZgLCwPiK3YLlZUnq7a+fCQn9RFzAI5FsM9RzMF8tZbI b1bYyqhFFMkQOEtFco4eeTOs39MpOA/rAmQpH3z5n+CrZ9+OwGQ7tha5shTJ uvASq5FMVaOGVYRDB65o2KH9V+xMv//wVfjg/Y9RE+wE5P52Q0ZYLtkczVmJ FQBzE9DL23v6rwwzqq88Ptu+6weaqvs298DGbc/wkVqNyDlmu9ZqxJkyJpNY Dh9N/iOsWHoODqAU9g3CmqUHoe3eU3ZHXcJWrH4GDop2Yuf9Ds9MZR2hEDm0 CCsOicVe7caVvc9PjUBv3zCMVBRHVyf1hqTDnY9jH4wk3ujz6uj6dFYUrNK0 sGs9tPcFWNmyB9raz8OaO3tgZ9frdYWxbdvpGIzbTaphRNrXRNS4QRQ2UZ68 kcwRL+noTyivSeVluxTFM5lW6Nx6RpvO3e0X4UP0RuW/gcc6+6IHJfDyHaDI BiPr9lryfWTJ1ekR6Bsa03qrzozD2Ac/BXmIuyU4BV17vqZgPdCFuvjwXVgS nCHYI+5h7zXSkuXJPGMveBCoBrMzM3DjhvqbGoU7giOwc/d3YfuqAHac/Tst bXNm45d6H79167AW3Y61q2Hd1ldhW6uB4UVnOyjvKMQ28DG91bSER/eaLBkz 0Bd0wLjeUq/CcFCCiY8959d9Z9oTzrnz+WN0nmwNaiXmsCru8FbArrvaoHWL 4XD3PW2wbPOfWQ59XOMS5cLCdsMseLm4sgt3yvAFO8Udf5Cj5UAUpfSFSQrn h3iF+c5bDWTS3QcrcHfLDuh+GNPbpHM4PvgvM8JmIyatLbqTHNHN6wNQHp3R jd1nBpYFq6Fj6yloWfuC4XP/YASGgZk4rR9w9uPbipLLuDXqhWnVFw8yDcFQ uU9PtB09fwr9m/Y5R37GYKijBNu7jB7u27gagvXfgA0tbRoWtpskow63W/Bb jc0qlrXmZmoq3ZD6zba+cPMaBB0VesCiBlc7Anjl6kDsqYr3fh5/0sIHQwmx S/HVaCSMDX3ftuPQsv4FONg1CCuVdDbsuGS5RlZcSbjS+dKhNwTsYoGM1Kkr 6ZglHzNvzC6YoMhWlizjqN8m8F1FoMtwz9Id8PaHOMKcVptzeisBPpYRNlnH O5vNqOUwPtALV6fmY4+rtN2+E9qU72i5a5CofU3PPAkTj7XwPoDPnmW6aCLh lTA3XoHKBNqLWg1LFahpzTvnx779Ew3buc/Ib3/HKVi27ik1r9osjM6ZRYqp i4gWGfWKNCK22Rpa8XV9aGxyrgZzU1fV7xKMvhZ/iufN4WywCWHF8QAzCmMr xjh2yxdMyNdy1wsq+BzU/D71vMnF//DYozEZSLk89NV3QhhbcbhlizWiRmu9 TpmryaMvf2oRN2KtSp8O0Yjr6DYVVteId7HmwnNqUJuuQNB7DX7lPCfVGbTD 9m1nlDf+k4iHljDpf9iiY7mTcZt44nkBrg8MwMTMBPQpE573Hk18T8MOHDNu oGdbu3IuZTWTlBPaaxC/++I+nKta0tp03fIcAJ9vxnXghj7YGvQOwLBaC65P z3mfCHt/PBsMOf0i/r/oWApWcZvJKA92nldiXw2n/uCCCmtfgqMP/wj61ndG OXW4um9jm4adWvUleP39WQG76NuwLUij9FWB5Q6W2B5mK6yEVug/IZoAfy8j LJMVKv8xNGJ8wM2xsl4R52ZGwH2krmvvUzrY3YeU7zWBb2fnkxYmHr0zVij2 DFT6vzADQyVlAaWriiPfsdRPNOzLP5zWoeoPv4xB3CVovxO/XzMZRDve/9fA V332pcSagFay+5vXkTV9kthYkofFTzLCFIva+tgk0DNz0URsYK6A3Zt3QMvG PzcaUb52eGouxvlGtURJLk34esm2Qy4N7FW5icJ7wwVBQ2hpEFvWfIWniEVi feLA3kHYGXSRL/MeHk6Af5IRdjPVIo2Lmqj02SJCZcJ4cedpTA3DIkLrxvNw T+seaN9VMROWYNSOg9bYrrjZdlhuTKGjYkwjfnY5BjuJ0ZHOkXugvf0MrG1/ Dg6Fx82ihzjN2HKvypeMfJFImRmrQEklq6Pj03poH8sZYK00+GkRx/GKoKJK TZVdu/dehDVfOK64uwp9yj/rGoRHCpJjrAp0I5uiHcIOoDkZ7qS0tWSK4YxR SdeRSNIV1lXYGp0Hcs3l001GnSXA6lgjHwRbqFWhRs/x8lPd4tlefWFF+OC+ y3DgQFjWxP8feUAvS/aRGllh9nkxrkPJU+08poDZhL6n9CPY31UhJ9AU2c+X 7axj0l5gjTXNGszenLVyJK7Ye8gHmBs8MMOUKX3jFgoVd7Semy2LVxTRsqS4 OMIIXkYgP0BG9O/Qk+l2S6zx8MTBaj9zxVk5dhKnu/RuMVdUpTRzHlhDRlid +udhgjEJOLwmzURAJiHjOjpNfO8eg9ZNeKDb2aXl3UjlYvRMOoNtcqR854hJ 5ABZkXTlnFEUDydHD8cVrdDpAViJOGpNztqzIVlIPQSTJy9YPk9GCdNjaPms DwfK291ROnDC24rykP4jUiQoWXa9+C3zLRYJFVH54U/eZJT702I7IOq2w83w nNPNOvqN6fIQj/PYJwtl/M4u9YTYRlwrScKLimvedf4rlkw+VCFLETmPFHkb lE3W7MQDOfA0nqRFsbiXS4JRSlxT10d2CiQweQ5a0aMV/CgNI/eImQ95ukeP KvfceOgY4cuI2CDQpeamZHZu/5RhOUuK76EnUwKISkHdD0u1TvslSSOImi9v VuLRI3Y48v0VmpgmD6JPm3UfzMcAM8p1WXVPZrA+bqNbsrkk8foQ1Dkj+Gkz G26uJ+o+51FQw61ov2jHEFjsS0GC6Akr3yOGv1mtF6xjYtrR8SQK03eItsHT +LOveGcfQj+w2h0rpzR7ume2AruiBv5TDZ9lK0AvXt8KsnJRR+0HkwZgqclj PEy9XKc1IWs9WIoE84X/rF9zqiCu880efIcSRRd9tlNFhDqe8GWmd3pw1BGO fkuNXivk2xtkxIYzlWvn58SBIX5WyDn0w55bPuB7kZo65/b5LQlyhijwPSnU egWeZ2FzdMBycw65dBDY4VK+VIpTkIoHtimZtgdCkzKrtHxs3dQvQ3yRReA/ U5DqaImPLcvnK1RsmM+CQARZMpSMlZaKKcgeJJisiPoed9bIfKP8SuDWqBlj rLwXuDVW8L/04j+SqT1GMFEsMbh5KYm94uW/0pHFCA3rgnV6H5cCYSeNE6uY paMshrMOWRYya6mjuePCDJh8fbzkf5L7mLeM5eWpaRQ9J2dyQ7YOB4zFOavg O+X638lYThFMxoNaBLXkPr0Ek/mvdGTdlAVKFBrDaSLZ8RH/mzzUIwSTGXCG 5jKZqdNcf/b+Ft0Pb3+WqPq83//cSZXDp7R38vhiQInHZLm3/OodX9Au0ZoV 99d+w47vDSVx9PjpPWzhe1NOZpT+re+E19/Qk0+M9KRL3i2/4iabYG/pNTbZ TODWX1GTJtHQROu/jsaXnUo8v+mp9tvpQD4X938n1d9J7bNy/3MnVfa/aa/9 TltyDJ5bexV4Npy3/nrw7GtO/DXeCS8H//+grPdm8FhUIHnO+CpwfZ013Zus sm7ljd9+DFle7K170kvT+WUdcqtb3te/6W3zwW3/B6SNq/U=\ \>"]] }, Open ]], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ "Calculate", " ", "and", " ", "Display", " ", "Correlations", " ", "of", " ", "Spins", " ", "One", " ", "Site", " ", "Apart"}], " ", "*)"}], FontSize->14, FontColor->RGBColor[0, 0, 1]], "\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"corr", "[", "ict", "]"}], "=", "0."}], ";"}], "\n", "\t", RowBox[{ RowBox[{"Do", "[", RowBox[{ RowBox[{ RowBox[{"corr", "[", "ict", "]"}], "=", RowBox[{ RowBox[{"corr", "[", "ict", "]"}], "+", RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", "j"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"A", "[", RowBox[{"i", ",", RowBox[{"j", "+", "1"}]}], "]"}], "+", RowBox[{"A", "[", RowBox[{ RowBox[{"i", "+", "1"}], ",", "j"}], "]"}]}], ")"}]}]}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "rows"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "rows"}], "}"}]}], "]"}], ";"}], "\n", "\t", RowBox[{ RowBox[{ RowBox[{"corr", "[", "ict", "]"}], "=", RowBox[{ RowBox[{ RowBox[{"corr", "[", "ict", "]"}], "/", RowBox[{"(", RowBox[{"2", "*", RowBox[{"rows", "^", "2"}]}], ")"}]}], "-", RowBox[{ RowBox[{"means", "[", "ict", "]"}], "^", "2"}]}]}], ";"}]}]}]], "Input",\ CellChangeTimes->{{3.417433509273056*^9, 3.4174335215139837`*^9}, { 3.417433703731361*^9, 3.417433710367592*^9}}], Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"Array", "[", RowBox[{"corr", ",", "ncount"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}]}], "]"}]], "Input"] }, Open ]] }, WindowSize->{796, 556}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{1, Automatic}, ShowSelection->True, FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 292, 6, 36, "Subsection"], Cell[885, 31, 610, 18, 44, "Input"], Cell[1498, 51, 2267, 72, 100, "Input"], Cell[3768, 125, 1504, 51, 86, "Input"], Cell[5275, 178, 3310, 102, 100, "Input"], Cell[CellGroupData[{ Cell[8610, 284, 2406, 70, 100, "Input"], Cell[11019, 356, 7709, 131, 373, "Output"] }, Open ]], Cell[18743, 490, 5442, 133, 189, "Input"], Cell[CellGroupData[{ Cell[24210, 627, 549, 17, 44, "Input"], Cell[24762, 646, 40088, 663, 235, 32450, 536, "CachedBoxData", "BoxData", \ "Output"] }, Open ]], Cell[64865, 1312, 1592, 49, 72, "Input"], Cell[66460, 1363, 261, 8, 27, "Input"] }, Open ]] } ] *) (* End of internal cache information *)