(* 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[ 64024, 1353] NotebookOptionsPosition[ 62202, 1286] NotebookOutlinePosition[ 62665, 1304] CellTagsIndexPosition[ 62622, 1301] WindowFrame->Normal ContainsDynamic->False*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Here is some data (time, voltage) from a Cavendish Balance", "Section", CellChangeTimes->{{3.4432982782783012`*^9, 3.443298311589781*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"data", " ", "=", " ", RowBox[{"{", " ", RowBox[{ RowBox[{"{", RowBox[{"20", ",", " ", "1.236"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"40", ",", " ", "1.181"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"60", ",", " ", "1.111"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"80", ",", "1.053"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"100", ",", "0.981"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"120", ",", " ", "0.889"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"140", ",", " ", "0.800"}], "}"}], ",", RowBox[{"{", RowBox[{"160", ",", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"180", ",", "0.654"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"200", ",", " ", "0.580"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"220", ",", "0.527"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"240", ",", " ", "0.491"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"260", ",", " ", "0.473"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"280", ",", "0.436"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"300", ",", " ", "0.436"}], "}"}], ",", RowBox[{"{", RowBox[{"320", ",", "0.455"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"340", ",", " ", "0.473"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"360", ",", " ", "0.509"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"380", ",", " ", "0.563"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"400", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"440", ",", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"460", ",", " ", "0.790"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"480", ",", " ", "0.854"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"500", ",", " ", "0.905"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"520", ",", " ", "0.944"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"540", ",", " ", "0.999"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"560", ",", " ", "1.035"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"580", ",", " ", "1.053"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"600", ",", " ", "1.072"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"620", ",", " ", "1.072"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"640", ",", " ", "1.072"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"660", ",", " ", "1.053"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"680", ",", " ", "1.016"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"700", ",", " ", "0.999"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"720", ",", " ", "0.963"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"740", ",", " ", "0.926"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"760", ",", " ", "0.870"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"780", ",", " ", "0.832"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"800", ",", " ", "0.781"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"820", ",", " ", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"840", ",", " ", "0.690"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"860", ",", " ", "0.654"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"880", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"900", ",", " ", "0.580"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"920", ",", " ", "0.580"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"940", ",", " ", "0.563"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"960", ",", " ", "0.563"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"980", ",", " ", "0.563"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1000", ",", " ", "0.580"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1020", ",", " ", "0.599"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1040", ",", " ", "0.617"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1060", ",", " ", "0.653"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1080", ",", " ", "0.672"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1100", ",", " ", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1120", ",", "0.745"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1140", ",", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1160", ",", " ", "0.817"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1180", ",", " ", "0.854"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1200", ",", " ", "0.875"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1220", ",", " ", "0.908"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1240", ",", " ", "0.926"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1260", ",", " ", "0.926"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1280", ",", " ", "0.944"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1300", ",", " ", "0.944"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1320", ",", " ", "0.926"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1340", ",", " ", "0.926"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1360", ",", " ", "0.890"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1380", ",", " ", "0.871"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1400", ",", " ", "0.853"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1420", ",", " ", "0.818"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1440", ",", " ", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1460", ",", " ", "0.772"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1480", ",", " ", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1500", ",", " ", "0.709"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1520", ",", " ", "0.690"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1540", ",", " ", "0.672"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1560", ",", " ", "0.654"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1580", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1600", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1620", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1640", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1660", ",", " ", "0.636"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1680", ",", " ", "0.654"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1700", ",", " ", "0.672"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1720", ",", " ", "0.672"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1740", ",", " ", "0.690"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1760", ",", " ", "0.708"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1780", ",", " ", "0.726"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1800", ",", " ", "0.745"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1820", ",", " ", "0.781"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1840", ",", " ", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1860", ",", " ", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1880", ",", " ", "0.817"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1900", ",", " ", "0.835"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1920", ",", " ", "0.853"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1940", ",", " ", "0.853"}], "}"}], " ", ",", RowBox[{"{", RowBox[{"1960", ",", " ", "0.854"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"1980", ",", " ", "0.853"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2000", ",", " ", "0.835"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2020", ",", " ", "0.817"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2040", ",", " ", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2060", ",", " ", "0.799"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2080", ",", " ", "0.781"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"2100", ",", " ", "0.763"}], "}"}]}], "}"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Length", "[", "data", "]"}]], "Input"], Cell[BoxData["104"], "Output", CellChangeTimes->{3.441842277492064*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"gr1", " ", "=", " ", RowBox[{"ListPlot", "[", RowBox[{"data", ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0.", ",", "2."}], "}"}]}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"PointSize", "[", "0.011", "]"}], "}"}]}]}], "]"}]}]], "Input"], Cell[BoxData[ GraphicsBox[ {Hue[0.67, 0.6, 0.6], PointSize[0.011], PointBox[CompressedData[" 1:eJxdlHtIU2EYxg8SEhEyIkIqbIaUIKJBEhHhGyUVhtiFqBBZYmEWaeKl0mxq Zl5qm3NeNnVHossfERIREiEjSqIgDAIJMc2pcxd1Zhf6Q+js23P+eL8DYzy8 73ee7/d87/lOlBQkXlAUJUb7FUX+abcinsyVkSdLfS9S6RD0z1nTgeY9xygH ehn1U9Dhb5uaUmLW01nopcJzW5+t3UX50KGyQeM712M6b45qf0/kcVIh6r59 YgFdRH1mwVZUNphIl1CftuQl/35rocuoT02kbXYPV9FV1CevtB501MVSKUX1 96yBj56GDXQN/eM+0UBlHq7LsX4iYeh4sHEbVUjrK3W/o6/zMzJvURXWT6+J t5e/KaYbWD8jcHdStZnz1KDff/j5xps/jlAt+oMn/37WIqPb8AtpcHnJf8gM vZiRWaNFQOYp5CnaU6gO7w8jj3ojz78efss7/Ne/rv6iBhPXdxRJq3x9I/zD 2mlnjzylRsn/LvwXR8VAUBP8F0RcJdQE/+DYuv2lD8foHvwD7uHKDy/3UjP8 51dXXmUNfKJmlefVAv85LS0tMmqBvz4PrWaedxv89flo83B938TP64EiaZVr C/H1Fvh7tXC0ESAr6l6Bc4asVuwH34cV/bN1sUkTaVvIls75bNi/D/nYRqPa H8Ezvqd28ASSxIBTeynyxLy0gy8U5SC7AfnHndYIKshu4udhH+S6A/z6fHXk SlqV+sNcO8AfEsdTTQ7wBycFIDnAH8wWFwR1gj+QKjZInWbO2wn+eS1MQ0EO dRl5Xl3gn4tcH4+81OXh89Ft4Hl3g1+fl27w6/PSo0g6V9KqpMNcO4m/32nl /k7p/F3pfL8u8M+J6zOOXKOct1fnx3z0gl//Xno9PL8+g6TBr89PH/gDkest YYj6FX4+/bmSBr8+b/1hXncTf5/byv3cU3w/arykifOoRdDFtf/Gt3/5Dwx1 GOM= "]]}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->{0., 2.}, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.441842278070369*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Fit to a damped cosine wave", "Section", CellChangeTimes->{{3.443293787709642*^9, 3.4432938069532213`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ff", " ", "=", " ", RowBox[{"FindFit", "[", RowBox[{"data", ",", RowBox[{ RowBox[{"amp", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Gamma]"}], " ", "t"}]], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "per"], " ", "t"}], "+", "\[Phi]"}], "]"}]}], " ", "+", "const"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"amp", ",", "0.5"}], "}"}], ",", RowBox[{"{", RowBox[{"const", ",", "0.75"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Gamma]", ",", RowBox[{"1.", "/", "1320."}]}], "}"}], ",", RowBox[{"{", RowBox[{"per", ",", " ", "660"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Phi]", ",", "0.01"}], "}"}]}], "}"}], ",", "t"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"amp", "\[Rule]", "0.4670321740100197`"}], ",", RowBox[{"const", "\[Rule]", "0.7774220086179631`"}], ",", RowBox[{"\[Gamma]", "\[Rule]", "0.0008217998126268101`"}], ",", RowBox[{"per", "\[Rule]", "657.7365277741015`"}], ",", RowBox[{"\[Phi]", "\[Rule]", "0.2112458980109271`"}]}], "}"}]], "Output", CellChangeTimes->{3.4418422790535173`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"f", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"amp", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Gamma]"}], " ", "t"}]], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "per"], " ", "t"}], "+", "\[Phi]"}], "]"}]}], " ", "+", "const"}], "/.", "ff"}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"gr2", " ", "=", " ", RowBox[{"Plot", "[", RowBox[{ RowBox[{"f", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0.", ",", "2100."}], "}"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0.", ",", " ", "2.0"}], "}"}]}]}], "]"}]}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVmHk8lN8Xx5U1ISSlVEILSZKtzTkJbULFtywhFZGQfckuhLFvYxnSQlSS Eoo2RdoIkd0wlhnDhCgVv/v7y+v9eu5znns+53POneuEi+2GC1xcXIu5uLhT //8X+HxUNPd+6IKpmtdBS+qVYDA7MtuNXQFTz3RvfNphA9PZ5vfEg7Rh6umT oBMX9ZEnW3cn/dAimHp47sayDbooka0126B0Dqbypni6bg2hfPZOs5p1KTAV yPxV8UoDNzJm1H0hD6b8rH3C9WNw8/ZKMTWbezDl+Yam+u80KrzRri/KfwtT ThZrZJfKozLziGb6phmYOl2fH/WZjmrZKtEVK8l65XtpW6Seo/p+sV+phaow pSC3Z++JCVRnTNm5a5H9ykWn9gYfRc3t5TrKZrUwtWrrY+extbjrzZ7ft7K6 YYqrMlPhxkXUZh68mCgjApNfJCjeubGom71dTl1CCybPHk9M/f4EdX/GiwnM +sCkmV06R80Y9QzH57s6KmDyuFP63e/iqL/oQUf4DbJ+/9HJBx51eNBeObFF mbCMhWJnnDUeUVXmcj+kBT86Wx6lUbXRqF6pu/TqLvhxmHuNdaEVmvqtMJCc XAsc7nXi75+ZofVIrJKFAS8wX9t+CVsrgDYCVWBlPADM+4Z5A2+r0Wat8sW4 P7XAzNhT4qd4BG30Jas4xyOA6bz6hM6JGrRJH7Z+9G8JMKUmDJT7/uLZPVF3 tUyWwaib/CszwXG0Df2w98DiNTCyaR2tfc4IL4gZ2562UgPGnbnNFhJ/8MLm YDMdoxZgJK12db/+Ay/se2istN8DGAFHbMVbHPGCo6g2l/xjYJzsWGvzoB4v vG6UKmCR9xfV/OlUeIl2HkaNsz5qMGjz6VWlTTnatxlBWpIaDGxz+lY2/hLt 2UEawWEtMLB6VGyqWAEvcpdsc/T0gAEB/4vP16TgRZVla7RPPwb6wIC8SsUX vBj1eWZovRrQs970+byQR4fdhg/U7xFe9uNbyKNodMwxlG55pw79SxWF/tv8 Eh0fcfMJHqyAvr9vBX4s3oeO755OwPu90Dfu3SqcxIeOnA2vij7qQt/Xjcfb S9zxks5Pu+CvJtCXeXPJspkxvDSS/VCpxw36VPYzYhPc8NK/41Rb6xnok7OX qaxwRydxvtD0fj/ok7w9qDTwC532uJouHgyF3rm+RtZ9a3SKO/C7nZkIvRk+ 28oO7MLLO1kHwmceQk/3NbwXEI2XD+UpVRF9eupWbrac1sTLZ0xXTPyugJ6H vYvyKUJ4OaJmxPzvC+i5JqZeLvINL39PouxY/Bl6VLj/yXxKRuegXd+6RdjQ naCraNSfh87J4y/FE1yhO+CgnW/2BnQuvHn3kNg0dDuGMiavG6HzV2H/Mom/ 0K13U/nAV390ke+Xub5GCLrmY8HbqQldGqIc1TcrQpe/1NG60nJ06fu9QW7C Gboc7W6+2boDXWYd2sWflkKX2a9VEiIX0VV8Y/4veWno2jVct3WNMboCv2jr 43zo/Ptvq6lCALoeAadjq2qhc9z5FH1hP7qaetfVXmVAZ5+887ETUeh6aSSw TG8LdNZal7RlWqNrxvuxhNYH0JkgURkZkYSutxYfEtjVCJ1hJc4hzGfoWrL7 ZlDOD+j0jr+zY2kgur4tNne2U4NOG5kOrytq6DoZW390pgo61SK/ihtdR9d/ b+XemHdBp6LxWqbaerwi8C9od80/6JSxMi62dsIr6501FCIROoUNKE4mrnjF wOgW/8p30MG8svmOHg9eORW5EOg/Ah10Bw9dyRK8YvvS/GffEujoKGijfPuG V3x3iA0WGkDHhxlmhMRtvBLucNlSyBk63vJnWxbH45X4G/XNrvHQUWPlsL3t K165Ix70WusrdJT2GG0PPo1XSo907sqeho57/Of5PxzFK8/DNEoXrYCOO14x al9+4ZWvU+zc96ehI3vTArf5FrzSo3R45TY/6Ei/iF3NCXhl9PytuMQs6Eia CPKPnMYrC60WAWY90BH9peXW9i3otlS4fLqGCzoit3z/81UM3ST1xZ1kN0BH WF3l42e96LYh8PJAhA50BN9z0FQzQDel8vfmzPPQEfCtez3rA7ppTsh/NYyA Dv+DS9DsMLod2Bx8+FEBdPj87ctdrY5uhtadryTrocPrl6mPSRm6mWVoavkx ocNzl5m29Ap0O9+YVNIrBB3uL7oy2mLQzVVgfPOBbdDhFvPFvlgE3fz3H6YV GELHlRyZl3Q/dIvwvS251BU6XKdqJq2a0C3xERfFJYlwfHoo72d0y2Za8jQ/ JuyarrwzAt0KZZ9e1fxGOO3uuPJRdCuzEJ/KmiXxFlVNZt9HtxprITmx7YTv 9+v5p6Lb+3O8JyLsyfezzyS4BKJbi918yJ9csr9PUjP7hNCt13G21LWN7B8W 6CO+6MZ05vQNLSP5Tf8SUmSj2083pqjFQZI/i3ll9A26c3kNQGMQ0Ue+xqqk CN2X+nW56D0l+uWaX9x8F90lA1tpVRNEX4tbLplC6L4h9PNnlc1Ef7OQ/Esy 6K4UUffvtjV0RGR888uio7tm9Kttq9NJ/aQiXhh4oLtOXJVl/BfooHR6XkPC x5LKYnn5oSOx59r2C4rofi7zNnPCGzoyCyh6v3jQ3ZmWu/p8CXTkXj7ITVdF d9/8jMPfh6HjlvutIyL70D2+KLqw9hR0lOyvjRLoQ/fMB2Htu4g/nyw5cPrq JLrffhTA/6AOOp6Jdp+7rYPuz6pc7Kia0FHX2+zoHoru72oupom4QMfnO4x2 q2J0b3p99m0Y8UPr/WC76n/oPtxgIu+8EjoGfcrc+nLRY3mHFl1nFjp5F5dt jaxDj3U9qmIV20m/vbq/TUsUPRToW/dvs4dOyZrR4+rn0AOYa/NWtkGngrtP fgUdPS7NLT4z9hQ6j7/k4Q0l8bz+/aGcnYBOSyf5vY0H0CN00fTzts3QaW+6 8mhDG3pkLBle8yodOgO6Hpn0vkKPWqkP39O8ofN+Nao17UKPL2vfLllaAp1V R4LuvwpFj44NNVrBw9BZJz734doa9OAoPEy/dAo6+00+ttTxoKf0rlRT1ISu ValdnHJu9PQws2ocnYWuxM7CX78b0TMwt2nTZhnoymOmHUp9h55RQwf8zx+B rocKo3kGM+iZ475lYw8Nur5ozcApVfR8F/vDr1kfupcl/ZGwOI5eq16Gyz9P g24q//wLg+/oJcc34/v7JXQX80R17NVEr23HLn7RYEJ3tbLU4xu30Gt/xzHf Um3opvM7as3No5fjtOTn2wzo2fbxUmHYBvSq3nzXJ56cJ58ujWAPeb/Oec2n j1bQQ1+760WdH3o1PY6TWxIFPTOzX4/vikGvof1un8K6oHfD5q70n9fQe5nF XlmfMOj1drJ64euK3mfjPn88+xX6VhQoBpz5hT68P6dl1K5AX5vxokTpHvQR W33shMQc9DGW5Y2pjaGPtPadsJ/h0DdJjyqnMtBHNfL0cHkG9Iu88j23+Sv6 WK+qLtGqgf6DVh+DfHTQp2JPpPY+Qeh/Ni4rmD+Jvo4hqy318oFeusyy/58N +n0Rwoz/mmHwpf/LQwZxGHDx5bdQzxIYaXzSt/uMIwa4TeknbUcY6YvycN04 jwEBG8tvjHyFEc5/H46sqseAxOi0VxYzMLpsSHF081IMqDI5xaWjDaOGH7Id 4CgGLh1uCxT5DKMfEhjDB5swsESk3bdgDJgv+X76N1/BwMr9h66fDQTmF8Ex N+tMDHzjXkFdIwrMHj65XzLLMbCtPaMiTg2Yf/o/5FUoYhBX/ulZD39gaSrq vtM7jkEn1L977hcEVjE9tLWqAoN+Wna4dijAGEVIN8M6EoO5NDIeZFFgLKv/ MPthMwbz1tTKN/DBWFEBDizTwGDRL2vENv2CsXreS3pLD2Hwpon60b5OYC/m S5buccXgkzvks0xuANvtz48vf+IwuLisc37vNhg/+Hyh7jB5XrqX38PxDoyb ON78kB6LweW1qkzqOhg/y8Nw8NbH4Fct17/NiMK473zdmqYTGNz2U7OkZBrG i3mqjnRzYQiPRvJZuSqYWAIjjlWOGGJZeaR2qT5MPHkS1eP9GUPO0h99WlsB Ey9W7HK2u4wh9kKrv6kowkT92U2fXXQx5IoVc8RUGCY6ql1dnXdgSARPjBCt BSbmr30zzH2HISVGH0xVbYGzQe5e605uDHnsu9NKtwU4G3e4b314FkMq87Ps T+kDR0Gt9XCtK4a8+XnJ96oCcHZIGy8+yoch7VlLc95ygIOXLZIXf8fQRUNH h8wCgXNGfIlxy2sM5V9WNnGJPLfxOv+j9xaGCu1a8yuQfO/cl58vFCcwVDJm TPCmHnAcbW+srgvC0K07YpXZQsDx5o2tMa7AUBXzac0FEs9PfL6/ygRD1cMt 94uReAErApJ892IotCmd1GgGTujvuMvRtRhqcvWjd3AmcOL2XVu2dBWGmt1R C04m8RIm+79v3YyhVl+yr98m8ZJzFF7kKmGog6xTVsNZ4GQ0Oh1deIShAfVC LyW2ACe/1JJafQ1DQ394vN9E4t0WLi3jKGBo1Oqur7uWAqfAqu1eRS6GJl2+ N3hmHDjFvTHabzwxtGD5MYGCJ8B5bH7kpDoTQ+8VFwYv/Qac8siEZ7YSGPpI l3vW5RdwKgqyvZcZYGi1R9WQ1l7gPH8xl27wAEPfCEueybYCTk3FoXKox9D6 21dauYKA8/LmyRwXIQxt/rblTf1r4NQez2Af6sXQ787he5QGgfNOovHYsQYM 7eXrfZTIC5y6hprz6nswlKmZdsPsEHA+iFSdVtXGMC4qT+DoZ+A0zq45b1mF YXw7rGeOET2/Gv1wW6aEYUvfP7v8SAw4zVm+fhK2GCY552bpZwKcb4J5vFYz GCad9Km5xws4bcqgkrcSw2S3KhzVSQdOu27iHfE4DNtm0bdbsAM4HYdPOk1J Ypjq9J5S57/A6dSspTlwMEwrNn1LM7lfda0a6+DXxzBt+ck8TW3CY59cB95h 2IFqw5VZNsDpLrPvaCLrD5vejVsIBU6PU3W54TkMMxzn5Tt3Ezi9qxrrXcj+ TSJsAupqCT/N/9v6E8PM1z37uXUIOH2HNS5QQzHM+ulKpwQBwp8oe3Z9wbAL Rm4D08S//bp3u2K8MMxx+LPF6SOEH8R9vv0Pw1yDFb4+dwIOfan286ZNGOa1 6toRmTjC5s99rrRjmP/DvlfhJYRzBH/N5GFYyOG9WiONhL9us3jzEsMi+tMf GkwSnts0uK4Dw2J9pzaXSgBnQGLuz2pVDEsSN8xdoU5Y9v72oCAMSy8qkvT9 j7DcfvUqGwzLOcBL6fYhLFk1/NoCw/I7z/LuJ/6j/1unxxOAYYXuz6/efkb4 m4ubQTyGPRBaOb2ki3D+o8r7DAwru+V+6fI84bNj0ZdeYVjl3s/0rzKEJTb4 9J/GsJpWRXON/STfyv+0HZwwrPbytaZM0o/9x+Mr9fwxrIG3/9B8ONGr6+tY eBqGfaHtfWl7h7CF/AL5vRXWqpGh+a6O6P0xBr4uxbA+e6NN8YKkPlFF27jP YBhjoShnSonUr9lZQYrox8rgW3HqGGExywkJot9MfTXP+gTiB7tmbno+hv09 u8o/7BHxS6Dda7GjGL7ot/vkMJk/369rb7mtjeHCior9D1cSv/l8/c+0EMPF C03VpUm9v9kOe428wPBVm4KirrsCpxUbbUPFMVxOtmW7bQPxc7tW4YZdGL4l bz70izdwmhLfdgjnY7jy+i3f9soT/4NORLI1hu9aExCwgsyDT16bbTfGYrg2 taAxdCvplwWaY74Khuuu/Co33g6chgDZPAsOhhtJbGqo30H6y+TKrc0PMPyC UOMqP9J/Lz0Fz6/KxnDH6LlLw4mkf3MbDCOcMdx1ifwLE+Lv5zXt5QIJGO7P 62u3jfRLRQf3j8VKGJ4wL1vWR/z38Fq86wwLw59PeBrpPQROTuXTYfFSDH99 OTf/0RngZG0rLypbieF1rPcz64neGWk5x6JcMfzryLrsX+fIPDtoEGFP9Bvp q2MWSwInauNZPbY+XpP4ujpC9CpwXHpUuzrD8NpqY93vAcT/Trd7FT+O4bX1 n12URsl8cjinoLzvDl5TbHjT/GY7cGxr6tw1JPEavrm8wZMOHJM/0Tb6bXjN 6fGL552HgaOZZu/5YS9eu9Ky6eZWMv/UuKKZDkvwmtcU5frVFOCo2JprLlXH ayGqFv+t7SPnh7DCWt9wvJZaMjNhQ+bHaq63BkWP8FpNsfIGxi2YmG5R2bN+ E0aI3siO4HDBxCXVX4tW3sWIFa+4L+MdmLCVWZy7hIURq/scTyYawMTpBecn mTIYsXG9loxqOkzop7ZonDHHiF3ZzZWeyjCxQa4ie3QOI2wzBFlzljDe3Gm+ iVKJEY/jfI/zVMC4it85/1JujKiSbr5dIALjGxd/0OSnYMSL4m1zh8/BuFQQ apqHYkRDHf1mHHm+2ITxTWItRvQvOjqz8jywm723o0cDRgp7rMtSFAO2e0WQ xa5YjLQze8MwcoSx+1orrh9DjFq1UcQ3cw2whO4Eekll4PXzgvUW26pgeLaD O/QEHa873D1RylsBw8yxPp5HuXjd5WA3X3c5DPfw2XP88/G6X/iPR5QyGK41 XGqUlIfXExbWLGHfg+HEnUlWn6vxevXU5Yr7uTCs5FXvc/oGRq/sWi65LRyG HL4MCmXQMHqtX44TbygMnfkrmpBngtFyUptfdQfD0ImdKrQfzzFa+b+9lykB MLT7dcPhlN8Yrdt4oZbtBUNCjs0iP0ox+kpthfv9i8Ao/b1RYPUERn+4f6Zp 2zFg8HV9Sv1WjNFNO55yC3+Dwb9NSv7vb2L0t3JRtXEbGJxszmg8b4fR/TWv M0o8YLB38Z+KkHKMnm3cbLMjEwYrRCil/g8xRm6SM6E+BINObEtPj3sYo+B9 ZIOkKwye+1RIGxbBGOU/N0/M/IZB82r32QJRjNnNY1peIQSDhz50d7BsMMZ4 RVXgblUY3HjnSNUVsv6qVriwdiAM0M/p7LitiTEhz7thvQAMdOjJ8MwKY0zk fg3XhUQY+Kr5fJtoG8YkHR5peXkLBl4b+Ugu34kxhRbHsnUaYODmkfGPFT8x piVw5daDK2HgYnDRXY8QjOlY7Hpmcx4MnPUKNr6khDG9Ee/jBRRhwDxwnfp0 AMYw4/0n32vDgMFr9rNvXRjLlddXdfQ8DKjKO12/jxir9KbosNFDGOA6UdEj YoOxO0ZRREIV6HONZjuZ+zBWU6T1a/sToP+0UpQK+Ymx+80XWZ6tAjqzKFT5 vjLGmnJOu7jXAv0b/+L/VLkx1mLFuJqWPtCb+NdvlTbGWJs9Yb/+vgf6x5Wl drNFGOsYURJ67TPQXzuaHp2extigtXyp6e1AL/l2qc/oM8aG62SbW5oDvdiC T0X9I8Zev6i6XqYb6Hc4RpNV2RibXHam8C4d6LRjS3P+KGFs4eHHz56NAZ3y 5fOjgz4Y2+Rl29/HBXSXBUNbdW+M/Zb1686tUKA7JcauciP76XwVd8mBF+gO ah/fjwVj7JBQ5fSkINBtH9gtk/uLsXM3xQR4JID+n/Mpde8JjF2ov/OxPg3o J68LC/QIIIVnfG8iRQroxo97jZx2I2XZrotrJNcB/SjQK3MTkSL35YXyli1A 3y9Hkx2qQMqWnyaTY0VA16670ls5iRTl1cynpduAvveqa+mkLVK07CT3E7/Q taSmykSTkLIv5h7v/GOgayz2Db/2CCk6pfsb3mgCfeffbQ+TBJFi8PeyicFe oKvIyYYnFSLlhCy3lGgN0JWNz3kbhiDl1MGM7tb9QFeKazY/qIaUs0m19lak HgqgXXIBkWJXYa4kWw/0zaW7VkxuR8ql7gnOMNn8JjUb01t9SPFSWOPnagx0 OecdHpZHkOJffet+BNFPdgnjUbMoUkJObOvLJvuTqayt/LALKRFDT5c/GgL6 +oCP3rl0pMT6o/57os8607k8LRpSkpa99+01APpaHYP7G8yRkn7rxL2fgUCX 1nnRnU/W52h19i19CPQ1p45rdPEjJf/TheUbSH1XB/0L67RDSuHZcX0tUg+p 569MK7uQcn/G29fwIGFhKvcToldZDNe9835AX+UWXHysCCmV66/3+t0D+kqW h4pVGlJqHi8XT+gh7OOyD4h+tYey9QpECUtd1um+jpT33Rt9qnWALvnZwaeH 7PfLlQfFzZ6E0626n91BSiufZs9oIWGPQxp0XaR0ZL4SW+ggfEFGtMUZKX3b j+hJChF26HtxWAIpjDfNPkrahEPDzxxzQArrtGWxjivhEp5PdvVI4Ywxek7f JDxhUflEFSkzIS5izt/Ifg6EpkheRMrfFb90rwkQLnZLO5yHcYuKQryzdpP8 Nsu/Xt6KcfwgWFTqRLg8sf98McYJNSd319GIHmaPIvxnME784lrR7iain2iQ 0LIvGLfy750D0zyEOyZS3XUxbm3idm9BDaJ3xV+R6nyM21K1v1sjk9SnpM5Z gTxXNvwgavCR1K/eU6pCBON2DpgcsJ0n9Z3ORgtnjAMhu7vx54gfIlV3tDzC ON28ia7bpD/k2HEid0ox7oi677JnxF8b7fdFd1VhnKlVjOeIEtC3pLiozYhj nPmUROG8FdAVD3a9a2vDOJvInC6JBOJXEY+XXoYYd+nhQx38CfTtLb8cdfgw LnRxqwiV+FtLXzzOzxfjItOt9pdwgL47MeO93hDGUbYOe7yTJf00wV4pwI9x GSa/OyYjSf8NfhpKqMG4BwXrCo4SPxs84OcL+oNxj/cUdJwlfjaKS51614xx lY07hL3JfDsRcv3lhmqMq/19wP3WaqCfptnnlVpgXOfRi/B3AOgXHtmfqSJ6 9IUo886KAP3iqFf18nsYxyif+vCD1MdJRcEgRhjjOLJBpxgkH7cxPW0JE4zn +Z3m8pn0a8ha9vETBRivXFCbm0v0yrNQ8XHmwfid3dEXqK+AfitCI2lLP8bv Wm68NZnMs8KXPB+NhzH+QGDn06gDQC+1rpFzycH406aTX1xJ/m+O/c1q8cP4 UB6ZBZ0jQB9tW1D/LITxkbuG3u4l/h0fOeIYSOLFutyL0cgD+jSfzd8xVYxP 69RaqTgL9PlL9ypcuzG++JHR9uXkvFne3mueNoTxrTYBVoPzMLBf+b/yE18x QcLKqffAdxi4rct6ZHoHEyTHNFX01GCguAF3a1ZjgpTfohD9BBh4ZJ69Y2EH JqxLTZM7rA8DL7LXiXpwYcKWjy/sDcvIeZenvHVBFRP27RGbMKPAoHiWzYVH pphgt7J8wWU/DF5vEmoY3I0JT5sWrc+8C4yETL5+HmVMdJDlVtzsDiNLFYT9 t2dh4qWnDP83+TAikZ///t8sJjofq/tk/RVG1q25vbDsHCa6+8a6ZqjCyI6l 6S7D8ZgY0LTiqQAHRv7rS4zevQMTk4K3HGBehpH8I2Za5jGY+Kzn2Jn79jCq 6S/Mc1ACk0Qy0xN3WgDzP0dvY0MHTBKzn6L+sgWm9dsD8+a6mCShbpRf7QhM h1UDfL3imCTVyP/4oB8w/YtU/Yf1MGkjr3ebRQYw89xM71HkMGmP6+m1YS3A HAnjjd0gjEkX9VfdbTEA1pWFwybcDZh0abnHI6oJsPz0HHbfc8aky32Nz6ws gBUa4qi64RQmuftHfRxxBFZSR5Wo+k1MCnz4a+JvFLDKRHXL045jUuqa7xry tcD64Snl+GgTJr36QX3tuQ/GLtAwcekqTKot7d6lHQ1jTuaiSTflManOTeYh fzuMuS9peZaUhkkfJwtyMt1gLFj3j9PdB5jUNvXE58UdGMvcsmq3JnnO/tm0 XVAUxj6kvExulcYkTrnE7eYzMNYU1T3lcxKTprxOS2cXwVib6xJhzXBM+jXT K7hdF8YGVpocn8zBZO7ZMYaJL4z9rrfQb76MyVJzAlk0OrDl7ybvXuOOydJV BuIXtwNbof7t5EI9Jq/3i7++4yqwlTsO3Tgnj8nyfyS9aiWBrdUiIRMjgMnb /8ofHz0KbIMNhQ8LnmOy3jwK7HwMbFdpf6/qJEw+VBMW+G8RsD0sZVipaZh8 NPDdz3eGwPahhKZ9asBk44VjA6dHgR381sLbeQ6TLbksagLXATshKnlNlAMm uy329qiPBPb9UxvVX6tiskeymTp3I7Afeq1fcV4Jk73ld/+E1cAuizKVCtfG 5Kt6855PHgC7MuICbzUTkyOuX/PKbwF27UhB+9EFTL6+2l6jh8Svu7HCJvoA JscUH5qVsgd2g8ENI6l0TE74JOQdPwfsLwH8EybSmJwlmuwTsAHY33VemqVU YXLODU+tCkdgdyY3Ggh8xOQ81f9+TZF8u9tGpVwkMPm2yWpfx0PApu9SYr9o x+SS9Bt+p5yBzbRuwJhwTH60JWx3UgWwxw6NjZ5LwOTHlefnPhN9xjcuFT0Y h8mVnVv8dVOA/aNcPbVaBJNr15Vc3fEM2L/2Kn3UKMTkdyWJey/zAPt39fJP 0/6Y/B7d/hYSPf+o/FgxS/T7bKsRsI7U89+POL1l45jcOLVqn7kSsBe0Dc+J RGByc9jvf6meMM4VyKu+YiUmt9+uDhQSIPeXb6ffvr6GyZ0audoHT8A49/hM tVQRksMseD40C8Z55ijfa6cxmT6qG/hbGcb5mMWmHo8xmeG3CdR8YJz/y072 3TBMHlnKv+DyGsYF7lStvWqPyWylD0FDpjAuqPx2hRMNkyeq7+MGGuHeYys+ 82DypGE8l+UIjC8N/b5u71VMnu51fUnmxbiQ1KXVSj8xedb1RHDzVcI3eSWe XsHkuUVq+0Xewbjw+rta14lf/yZJLjqyjHC82Zd2Uq8F2dlX104TnpTSPT6H KYvKvoe8zIdxkYPM9aIJmMKj+2z/nzHCCV8ErqdiCl9rzmJNdcIfG+qDhjFl iV3ga7dAwn/69Ka1MGXprE3o/XoYX7ZOPHKFM6aIROnojIoT3nm+VKUCU8Sk 5Lnlyf1x2Z4OgRvnMGV5Ee8b6zuENT05L89giuTu4bDMCcKb9rCavTBFiltn bkGUMO+yeoPnmLLmA83tHMlXpKXFxDcTU9Yl/x6tMyGcGhf+jMTfYGl6dqvX //d/wLvVF1Pk5R+2x2eQ/EZ/eX0uw5RNbCGjqUrCAY+m2wcxReGJfd2pDsKL feJV/2GKUsAb7Wd/iX4+RsJCVZiyXX9ducw6onfvXmuHa5iiKuJHroOENfVv LGZhitq31lsjNqQ+wR6SitmYopm7Q9ogFMaXVH1QKp7HlN32sckPb5H6Dp3w v38SU/ZtH1kq8ZYwt+ibhAuYovOS9rub+IlP1FCI/QlT9KLmruxXhHFerk5L 9WJMOXTcdOT2UeIv+nOoKMeUY3ShtstxML7IL31daCimGBddNGx6SPypknfv 9AdMOen25p16E7Dn26WndI9hihm33+N5CeL3P3ExIs8wxeLDNyVbDWDPeac8 mBDBFKuUHTffnSL9wti+I1ccU87LjyTFZQL7Z6ruve9NmOKi/5/rehnSb4d/ BH0SxRQ3kdLhMB3Sjyep15+S+nm0CVsNnyP9aiSTNV2HKX72tcdK7gB7aM30 UkYIply7rroVST91Cze937MGU6KOU/Jvkf7rKEiq2P4NU2KkRlcLuAC7fWcE s8YJUxKK8gQaHwG7WVb93yvAlKyPwkNntYBdb/FDgoGYQktxOPPWHNhvoyj/ iTVgyg3L2lYFf2C/Lqgy8a3GlDts/zc/aoD9vHr2h0orpjwSGb0Rqkvm39Fz WiuXYEr98beWN48DO+Wnz5IDsZjSwHyx1tcV2InxGVNRhzHlU1hFj1EcsOPW VY90FmDK1yf3rP80ADtS+gAXXxqmdEslnzMh+fukRkU+c8OUH/1nHXlVgW0W U9vaS/w17Weh2GkMbFP/zizJaUyZXW7CKiXz7PjZagfpTZjyV+/Q5TNFwD68 bCXfOTtM5bu73bVcFtiaP6j2a90xdY3rP297cWAv7+9wGHHC1HVLZrX2qQB7 mfqidIk7mCqTz/m9nOgnGNi0VjYAUze2DPq9iAY218zNU7ftMVVF60PgqsUw xvYpefUyCFP1uKjh9ZMwVtscbPqPfM85QT1J8SuMXTZ2COPuwNQr0pKveDVg zL5MvF6fH1PdC2YneqkwZityY628Dqb6vKg6lmIDY/89LH60cimmhk6AwL8J GIN81UD+PZiabngoqEkIxkQPzRc8PoupL4UsLnvrA+v+yhs3+SYxTSwq5NCb KmC227/qKTiOaa+td79PyoYRNuPk9qM7Ma12audv/hYYGXrHWHFxFabVRW7b EiAEI71FpqmMKEz7WLI+0v4qjDR5r7x68SGmfVvg0dtrCSNPFLnODPZiGov2 6eWQNIwE8BjPWIxiukSPTcUuEl9k4VfSf56YvtLNbPgBiS+gLfjJ0wTTV/Od XClP4i8KXqb8sR/T1yvreYlcheFp0W1h+wDTFQIVdw5YwvB3W2+5yS5M1147 VUKRhuFbXa9njacx3f5MxJ2+HBiGRIU/TT2Y7rjlNyXwHgxrDQ14Ol7FdKdp Jw/p5zCsesDlpWIFprvFnNh/uguGN4rpPyX36/SAyrWdTSTe0kVf1zRkY3qK xBPRV9kw1K63xT84BdNffhjwy8+EIZ+3yStjvmD6m/T/zmIRDLmviPcr1sX0 d7bvD/ZUwdBlZ3PqQyamf/hdsmJ1BwzZ7k4VjNTE9LbNfiXJq2HIAKRFQ4cx nR0mNhhGhSFZnlXv7fswY/VuNDyfAYyWljhVmhtmSH87HuhxCxiNN7nTlyhg xno32wfhpcD46OczvOkxZsgXXRO51QCMWt2wSa3vmLF99fvPg3+A8VhvrqJj FjMO/DlpdMEKGGm8yzb/rsQM/fTzQZ6OwEiy37AQ0okZh3d6llzzBkbc1+I3 LwUww/BSxrLbCcCIfKd96IE0Zpzu7PrCeAUMPxkx/kPVmHH5uZ2xnTwwbG7e 3lV4AzNcT3sHe6kA40ze4oFkecxwm456GLEPGOZ31BzobzHDR6lY9M5/wDBp 3FQkqIEZodkTjUORwDic+cHvvSpmXNu1iGsmBRj6Qwa6wl8xI7JVfAdfHjB0 97x1W6KGGRRhtcRNFcAAnlK9ID/MSA/wPW7HBIa6SMeUkCxmZEpFh3jNAEO1 KOZt/B/MyH6SVRrJDQwV4xPtSPZ3g10jViANjK1VKiWOBphxz2px07AhMOT4 g/o7WZjxYE5i0aw5MGSmtZTVxzCjNG2jKr89MNax1g72v8KM8s8HkzYHA2P1 rL/CBaLPK+3YE/ZE/+XRVWeST2HGm0/hnKs8wBBralYoo2DGO8uAuASSr6hs rQD3dsz44OfcUDEHDCH6TnrlXsxoLT9+QFAHGHxtrlwhAZjRrne4fx3Jn/fy 5bmFq5jR0aITpDoEDB7hY/+tFcGM3smdzyyigbHIsWyx2FbMYCpL7rz3FQb/ Zm+hBV3AjLFqkaZX8jD4JwSfPu3BjAkDfpdWLxicc91neWQ5Zkw7/CqeXw2D vy4NCoWT/c78+nFE3AkGZ30o1Y3HMeN35OjIphoYnIlfo/eB6D9/u1PeyBYG p7/3SbUpIXWRWvPrc48Ji0gePhGEVO43H218+GBwykBZ8rwYUvlO1M7HnobB yZRN5YarkCrQX52dVwSDPxgL+juTkLrUtXz347+E95d/kBlHqvDCg/b3hjDI KTBQleAgdfnaPIkfkzA4kbC0fz4bqSvuZZTx6RJevpPCdQ+pq/YkHl+dBoPj eTt71MORurrh+oTyCGHNpUcG65G61iyUcmA3DLLbKlRCPZG6fsR/66lYwqFa LZVGSJX1dn9/qYewRswFz3ikbuRzsg9WgcGx6fu/Jp4jdXPqBd6UUMLPaZSn r5GqKH/mZmEL4TgLwZf6SFUq+29/9SbCDv3gKI/U7TqGvU0+hI2UBxK2IXVH 08EARgNh7f3pyeuQqmaDa+akCWtKDl6oRarGhFaViDPhXff2UM8idVfgjtOy Lwnrzc867kXqXmGFGU1xwpZC5h89kKqdLZty9Dxh/6ZEn39I3b91jap1OeHb R0vmupB6oGp5o4cA4XbvsueJSNU/LOQcZU7yW3EibEU1Ug+18wpl3yNs0c4+ PYPUo3b/ih7OE77HX9TfidRjP2cOvTUm+vF2npEzQKpx+PjQ93zCDsZxT1hI Pbl8OJw9TfjbuTI7Es80v09usT6pxzGxLfqk3qd3tL+SzCD8yXRneRxSzV82 WSsyST1PKfFeDUTqGaP3/2AvYWZi+af9SLXueZ11Mo7UPzL0iDU/Us//LWu7 qkr80j1Vtug+Uu3fbFKqtCR+KkjcqkjiOV7PCJmJIH67apL78yBSXSUDlFw7 iD8PajwL24VUt66JkAc8xL/7DpeNiSPV8+bZbyxl4m/tSwaeT5Dqr6IfcoHU c+7sq7IjCUgNmKn4lk/0+BM22jfyH1KDq7du7f3//+9L+S+HqiE14siyb+aK MLigJCC1mKxPPPd9q2ET6b9b1m4GkkhNUTwaHEP6l2+tu2ytG1LTONWt78n8 47+hu9aBgdTsgPzgAz7AEKxfeFbXjNSC9EutWjJkHmRwd083IrXoTI+i5xEy L+y2jGVpIvW+vHHQIw9grND+MHXICamPStUVleqBIbXMcVP9KFKrG+aDNrgC Y8M2ibrhZ0h9mejaYpVJ5pWh5j9ZUp83pwYUsmqBsdE36bvTVqTWD9a1SK4C hsJ4g5lFK1Kb5xMVBF8AY+cJDZsniFSmqtzXaRFgHJy2aA6zRurYr7QtO7SA ceSxfabpIFInXggEONsC41igoJ+QKFJ/HmVvHikHxsltoyP/5jCT60L51W5y fljNt09F38ZMSerhTe/uA8PHwmNMcQ1mSqmYvfVsAYa/V+T8ASHMXFPneF6e 6BVINeVWb8JMmZ8xN8L0gRH+c/qFCAsztx7/LA195PwRVDYe0MJMFDwpVikB jPupLv9dKsdMB3/Lv8V+wBhlKVIeFmGmk/jlLHNyHrBHeoV5ezHT+W7g7iV1 wPgxNt5aFo2Z7m253hfJ+3P8nGfXbDEzcGf/pPx9GBKMUj6/bRIzk8bshnP6 YEjpDrvQXBszn59x/ZqgD0OefML14ucx88XWP1su6sKQnzRX3MYUzHw1dy0Q dWAoSDNRUuc5ZtalZypwAIauBx3a9mAWM782vQkx1oIhmpH6IaoCZo7qS6qI boWhevq6z7PcmCWp/IySIArDG8zpQ7E/McuVi+9Q/HcY7nzaYvm5CrNlLL0h rwVGjTf8PbBzHrNlxyMNJiRg1NRGVOpiGmZvDE4332cKo+ZZaRN/1mO24q2n Ht+/wegFwY9nzOYwW21stki8DUZ964NybKMx+3CAj0R4B4zml3t4DAVgthvN Z/RiH4xOCu4q2dSD2Z4qUTNPZWB0pjdTJJqN2d6vqTy8Z2H0z8Nf0o99MPvq UOW6G3Rg8h4MeYaumB2p/PtkxwAwV+nkrOpWxOysF741BkPA3GdfuvSYPmbX 9vklq7CBGb5pmXprNWbXFfGcKuQDZtTrG9b1s5jd4B4nvX4dMClmi/JtyPca +W7cFjYEZprbvLxmG2Z3KtWXjz4EZqF5wQvvAcz+4SvZke8GzIYJwYmIdsye 1snLlYoG5hc52LCzHLNnhRTPJeQDs9k4SbvgLWb/zd03FtgMzM40pxhBXczh f3t+3kINmKxP3283jGHOOrEy2eUzwBIc2Ch73hpzZDr3DsWIAEvkVtkQSGKO 3K13xYs3AUvc6n3Y7iOYo6DZsZNzClhSz+0orErMUbNapPexAlhbpP22xqVh zpEiI4dr/sA6cJHFKyeKOcfcv2/7mwIs/RKBErHTmGO8z3bS/R6wDrO3RTws xBzTRi//s13AMj7uqnvRFHNsZnIo2trAsqxcpW7ZgDneOqyHs1zAcl/QWq5f ijm+3ctr+lYAyyt5otq2E3Oueu/5+F4RWL4bLD6zP2JOSHHMSNZJYAVugp1m IpgTK6EkA3eAFdlWWVhmiDl5DKf4KANgUVun/jz2xZybQSk5V84CK2ua62Fe MubckXpebO4FLJpgGK+RK+YUGwrVKeUCK1+mM1PSHHOelN+bb5wE1t3JCOXo UcxpiGQ7S1GBVbH5uAy3F+Z8kl1xdfEDYFV+sPzX8xdzvjzfd531GljP7Msa HRYwp4VDuV09Bqya67sv1S7DnF6zbd22CKxaF58EOzXM6Z8yZR01Bdbb1sKs 1YOYM0gJ+K3mCKw6tRXtki8xZ/TV5xV8ycB6/z0/ksyLnGlFZ4OiIWB9rjvn vdgCc2Zq08yT/wDry4TFgdQKzPltVXMxQBRYTcsioxuFMWc+WSTcaBewmrX8 ihPXIW3RNo0kLUNgtezHQ86hSOOus8rbYAusVh3NoQwzpAn8efBsKhZYbUrP rg++QdrS1G/vu24Aq11SfzJJFGki2+fb3pYT/sUtl05H2vLzhlPpvcDquCUc Yp2CtLUfJ7buVQVWN+3S73fqSJMxmy9QIPF69NV3HLFDmuywiJzkMOGhwP1h HUjbsmjb6h/ZwOoTTRzhe4s0xbi9aT3En32ZZ+5+/4u0bWsMxD8aA6tfujjL wRBpquqXBAt4gUUXbFJ6twppaq/9Isj9jkX3Krp+MRZpmkbRi0OqCHcuiT8n jbRdXdRAFwlgDWhOHDHORtpeh8I5S2fCFMu2ryJI056p8DpcT7jzpH3kMqTt D6ub1JQD1uCGlj18RC9d0TZn+UDCNn23W54gTT9niCnWTjjDT9h6L9IOK85c WCD5DtbfGOF8QdrRCt7+MQrhH2avirqQZqi34sz3EWAxxG+qXs1AmvFX+e91 OoSVgkZWH0HaSWs108ckf4b2j7xpZ6SZjh1oukHyZxye26ZcjLTTviePxZP8 GceynNPuIM2Cz/b91eL/P287nVmItDPJbnqOfP9/v4SDJF8bmZBXp2z+H3+z gdQQ0mzvJ+7Trfr/9yFFOwFpF3bnVaoSPQY5fxZrGyHNvu6h+nrn/+/fon8V qaejyYtSoff/z+88Zace0pz6vyjNyRM+u9LjmDbSXJx7Cof/r4es9wqvOqRd +TMu30L0GOiKWMRL1ntE/ct7vZNwvH57VAHSvCWFpUviCO9+bviS1M/3pnRG NtGD3jO6Of470q6qKElEEz3o/o2aU++RFli9J8Gb6EEX87jwcwppYW3mUSeO k/pv5uL/NYq0iPOOPEDmSd8ddvaQCdKifvgGbyP5962/lbhTDmlxSzN8+Em+ PX/8k6XmkZaQUTD1k+TbY5V/aGsJ0pI3PnWhuxB/VuUWuxF/Z8A3+2qSX9cZ pdAaFaRlfmQMFJH8Om++kVN5gbQcs5/W6d+Jvwe0GzLvIS3fQ+KUG8nnu8H0 3rXKSLtXeOLg5l+kfz7UZsVtQVqJ+tk3K06Q/vo+MbmB+KX09RVYfJ/0X7+M WPplpJV3JWh0k35r6sitHo9H2ivRz5uSPgDrk196XZ0/0t60njIrnAfWR80W iZ1PkfaO2hf7QgVYDWzu7TzTSPsoOznJIvO1Hk4ojugjrU195QtdK2C9Drc9 y/8Sad9/502aJwDr1aKBB78eIK2rRnGjK5lHL7yjC6YEkUY/uDcmZzOwnh89 n1ZO3mebnT31kwOsJ2lL3uawkcZZy4wRkgVWWauv1ultSJuku9fImgDrkbC+ dWoq0n5dipA3rATWA7uVht7BmMsdUMQpCAXWnWq7+IYSzJW6MX3djOidtiRx 7Ye/mCttF/DcRR9YyV99rbXDMHf9Vr6Jaz7ASkzRXl2mjrnyj6VMS7uBRVkc pL1bAHO3v9OWFSTzOyy6uU2mDnP1mJFVz8n8cjmR/WQhHXPddq4eW0/y1eH6 suFfM+Z61OtcLCD9BluPC2z/jbneZxwY258Da6/xbTuKNOZejazoQSlgacQM 91ANMDei679GmxZgKb59d1RBHnOzriU9zj9CzjtD2i+vHMytbRcM2KgBzO7p C+e492He0poA/hWRwHQofC/8UxrzhOU/5Y0LAvNCYLiN+EnME4lZt/tdHDDP Gs+/phVjntipF87eKcA8PXahU6UU8yR/LGol9wfmwXk3L6fNmCcrF3Er5zkw N2ZpK/nwY55WVILOxkkY7ZXccffsfcy7cOJmgIYljBoKmRzXmSZ85nT7+xMw emzUSTXjGeGLImqWhwm/C1/1PZJwgA8rVINwcNm0nwbhQkOzz8tg1ODncmaD LOF/vzUuvIbRo90tI2q2mGd313AidTOMHi6xaLhLI/yY5+iWdYRjPbyGyf7s airvPCO/tw47UBQXrhFukT/Tz0VYrqZ4XSrh+d8flcjvqUMZssIFNphnb3Kr sDYWRg8GM10vVxO2Nuc5FUrYivu22SbCjsusmb6E90gPpI0SDvJbKWoHo/o/ ja7aWxMuMrpmCYQdnmwY+0P4CU/fBMlHX//zzb/qhF9U7Q1VIiw3oiIVSLhV fqpwFYzq9ayx61pCuLfDaI8I4WfqYg/VCDPjiz7zEs4wfMfZS/inHp/NX8Ke FyPKUggvzJ2dInrrnQg9Lvb8f6iLUMw= "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{0, 0}, PlotRange->{{0., 2100.}, {0., 2.}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{3.441842280584812*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{"gr1", ",", "gr2"}], "}"}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{ {Hue[0.67, 0.6, 0.6], PointSize[0.011], PointBox[CompressedData[" 1:eJxdlHtIU2EYxg8SEhEyIkIqbIaUIKJBEhHhGyUVhtiFqBBZYmEWaeKl0mxq Zl5qm3NeNnVHossfERIREiEjSqIgDAIJMc2pcxd1Zhf6Q+js23P+eL8DYzy8 73ee7/d87/lOlBQkXlAUJUb7FUX+abcinsyVkSdLfS9S6RD0z1nTgeY9xygH ehn1U9Dhb5uaUmLW01nopcJzW5+t3UX50KGyQeM712M6b45qf0/kcVIh6r59 YgFdRH1mwVZUNphIl1CftuQl/35rocuoT02kbXYPV9FV1CevtB501MVSKUX1 96yBj56GDXQN/eM+0UBlHq7LsX4iYeh4sHEbVUjrK3W/o6/zMzJvURXWT6+J t5e/KaYbWD8jcHdStZnz1KDff/j5xps/jlAt+oMn/37WIqPb8AtpcHnJf8gM vZiRWaNFQOYp5CnaU6gO7w8jj3ojz78efss7/Ne/rv6iBhPXdxRJq3x9I/zD 2mlnjzylRsn/LvwXR8VAUBP8F0RcJdQE/+DYuv2lD8foHvwD7uHKDy/3UjP8 51dXXmUNfKJmlefVAv85LS0tMmqBvz4PrWaedxv89flo83B938TP64EiaZVr C/H1Fvh7tXC0ESAr6l6Bc4asVuwH34cV/bN1sUkTaVvIls75bNi/D/nYRqPa H8Ezvqd28ASSxIBTeynyxLy0gy8U5SC7AfnHndYIKshu4udhH+S6A/z6fHXk SlqV+sNcO8AfEsdTTQ7wBycFIDnAH8wWFwR1gj+QKjZInWbO2wn+eS1MQ0EO dRl5Xl3gn4tcH4+81OXh89Ft4Hl3g1+fl27w6/PSo0g6V9KqpMNcO4m/32nl /k7p/F3pfL8u8M+J6zOOXKOct1fnx3z0gl//Xno9PL8+g6TBr89PH/gDkest YYj6FX4+/bmSBr8+b/1hXncTf5/byv3cU3w/arykifOoRdDFtf/Gt3/5Dwx1 GOM= "]]}, {{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVmHk8lN8Xx5U1ISSlVEILSZKtzTkJbULFtywhFZGQfckuhLFvYxnSQlSS Eoo2RdoIkd0wlhnDhCgVv/v7y+v9eu5znns+53POneuEi+2GC1xcXIu5uLhT //8X+HxUNPd+6IKpmtdBS+qVYDA7MtuNXQFTz3RvfNphA9PZ5vfEg7Rh6umT oBMX9ZEnW3cn/dAimHp47sayDbooka0126B0Dqbypni6bg2hfPZOs5p1KTAV yPxV8UoDNzJm1H0hD6b8rH3C9WNw8/ZKMTWbezDl+Yam+u80KrzRri/KfwtT ThZrZJfKozLziGb6phmYOl2fH/WZjmrZKtEVK8l65XtpW6Seo/p+sV+phaow pSC3Z++JCVRnTNm5a5H9ykWn9gYfRc3t5TrKZrUwtWrrY+extbjrzZ7ft7K6 YYqrMlPhxkXUZh68mCgjApNfJCjeubGom71dTl1CCybPHk9M/f4EdX/GiwnM +sCkmV06R80Y9QzH57s6KmDyuFP63e/iqL/oQUf4DbJ+/9HJBx51eNBeObFF mbCMhWJnnDUeUVXmcj+kBT86Wx6lUbXRqF6pu/TqLvhxmHuNdaEVmvqtMJCc XAsc7nXi75+ZofVIrJKFAS8wX9t+CVsrgDYCVWBlPADM+4Z5A2+r0Wat8sW4 P7XAzNhT4qd4BG30Jas4xyOA6bz6hM6JGrRJH7Z+9G8JMKUmDJT7/uLZPVF3 tUyWwaib/CszwXG0Df2w98DiNTCyaR2tfc4IL4gZ2562UgPGnbnNFhJ/8MLm YDMdoxZgJK12db/+Ay/se2istN8DGAFHbMVbHPGCo6g2l/xjYJzsWGvzoB4v vG6UKmCR9xfV/OlUeIl2HkaNsz5qMGjz6VWlTTnatxlBWpIaDGxz+lY2/hLt 2UEawWEtMLB6VGyqWAEvcpdsc/T0gAEB/4vP16TgRZVla7RPPwb6wIC8SsUX vBj1eWZovRrQs970+byQR4fdhg/U7xFe9uNbyKNodMwxlG55pw79SxWF/tv8 Eh0fcfMJHqyAvr9vBX4s3oeO755OwPu90Dfu3SqcxIeOnA2vij7qQt/Xjcfb S9zxks5Pu+CvJtCXeXPJspkxvDSS/VCpxw36VPYzYhPc8NK/41Rb6xnok7OX qaxwRydxvtD0fj/ok7w9qDTwC532uJouHgyF3rm+RtZ9a3SKO/C7nZkIvRk+ 28oO7MLLO1kHwmceQk/3NbwXEI2XD+UpVRF9eupWbrac1sTLZ0xXTPyugJ6H vYvyKUJ4OaJmxPzvC+i5JqZeLvINL39PouxY/Bl6VLj/yXxKRuegXd+6RdjQ naCraNSfh87J4y/FE1yhO+CgnW/2BnQuvHn3kNg0dDuGMiavG6HzV2H/Mom/ 0K13U/nAV390ke+Xub5GCLrmY8HbqQldGqIc1TcrQpe/1NG60nJ06fu9QW7C Gboc7W6+2boDXWYd2sWflkKX2a9VEiIX0VV8Y/4veWno2jVct3WNMboCv2jr 43zo/Ptvq6lCALoeAadjq2qhc9z5FH1hP7qaetfVXmVAZ5+887ETUeh6aSSw TG8LdNZal7RlWqNrxvuxhNYH0JkgURkZkYSutxYfEtjVCJ1hJc4hzGfoWrL7 ZlDOD+j0jr+zY2kgur4tNne2U4NOG5kOrytq6DoZW390pgo61SK/ihtdR9d/ b+XemHdBp6LxWqbaerwi8C9od80/6JSxMi62dsIr6501FCIROoUNKE4mrnjF wOgW/8p30MG8svmOHg9eORW5EOg/Ah10Bw9dyRK8YvvS/GffEujoKGijfPuG V3x3iA0WGkDHhxlmhMRtvBLucNlSyBk63vJnWxbH45X4G/XNrvHQUWPlsL3t K165Ix70WusrdJT2GG0PPo1XSo907sqeho57/Of5PxzFK8/DNEoXrYCOO14x al9+4ZWvU+zc96ehI3vTArf5FrzSo3R45TY/6Ei/iF3NCXhl9PytuMQs6Eia CPKPnMYrC60WAWY90BH9peXW9i3otlS4fLqGCzoit3z/81UM3ST1xZ1kN0BH WF3l42e96LYh8PJAhA50BN9z0FQzQDel8vfmzPPQEfCtez3rA7ppTsh/NYyA Dv+DS9DsMLod2Bx8+FEBdPj87ctdrY5uhtadryTrocPrl6mPSRm6mWVoavkx ocNzl5m29Ap0O9+YVNIrBB3uL7oy2mLQzVVgfPOBbdDhFvPFvlgE3fz3H6YV GELHlRyZl3Q/dIvwvS251BU6XKdqJq2a0C3xERfFJYlwfHoo72d0y2Za8jQ/ JuyarrwzAt0KZZ9e1fxGOO3uuPJRdCuzEJ/KmiXxFlVNZt9HtxprITmx7YTv 9+v5p6Lb+3O8JyLsyfezzyS4BKJbi918yJ9csr9PUjP7hNCt13G21LWN7B8W 6CO+6MZ05vQNLSP5Tf8SUmSj2083pqjFQZI/i3ll9A26c3kNQGMQ0Ue+xqqk CN2X+nW56D0l+uWaX9x8F90lA1tpVRNEX4tbLplC6L4h9PNnlc1Ef7OQ/Esy 6K4UUffvtjV0RGR888uio7tm9Kttq9NJ/aQiXhh4oLtOXJVl/BfooHR6XkPC x5LKYnn5oSOx59r2C4rofi7zNnPCGzoyCyh6v3jQ3ZmWu/p8CXTkXj7ITVdF d9/8jMPfh6HjlvutIyL70D2+KLqw9hR0lOyvjRLoQ/fMB2Htu4g/nyw5cPrq JLrffhTA/6AOOp6Jdp+7rYPuz6pc7Kia0FHX2+zoHoru72oupom4QMfnO4x2 q2J0b3p99m0Y8UPr/WC76n/oPtxgIu+8EjoGfcrc+nLRY3mHFl1nFjp5F5dt jaxDj3U9qmIV20m/vbq/TUsUPRToW/dvs4dOyZrR4+rn0AOYa/NWtkGngrtP fgUdPS7NLT4z9hQ6j7/k4Q0l8bz+/aGcnYBOSyf5vY0H0CN00fTzts3QaW+6 8mhDG3pkLBle8yodOgO6Hpn0vkKPWqkP39O8ofN+Nao17UKPL2vfLllaAp1V R4LuvwpFj44NNVrBw9BZJz734doa9OAoPEy/dAo6+00+ttTxoKf0rlRT1ISu ValdnHJu9PQws2ocnYWuxM7CX78b0TMwt2nTZhnoymOmHUp9h55RQwf8zx+B rocKo3kGM+iZ475lYw8Nur5ozcApVfR8F/vDr1kfupcl/ZGwOI5eq16Gyz9P g24q//wLg+/oJcc34/v7JXQX80R17NVEr23HLn7RYEJ3tbLU4xu30Gt/xzHf Um3opvM7as3No5fjtOTn2wzo2fbxUmHYBvSq3nzXJ56cJ58ujWAPeb/Oec2n j1bQQ1+760WdH3o1PY6TWxIFPTOzX4/vikGvof1un8K6oHfD5q70n9fQe5nF XlmfMOj1drJ64euK3mfjPn88+xX6VhQoBpz5hT68P6dl1K5AX5vxokTpHvQR W33shMQc9DGW5Y2pjaGPtPadsJ/h0DdJjyqnMtBHNfL0cHkG9Iu88j23+Sv6 WK+qLtGqgf6DVh+DfHTQp2JPpPY+Qeh/Ni4rmD+Jvo4hqy318oFeusyy/58N +n0Rwoz/mmHwpf/LQwZxGHDx5bdQzxIYaXzSt/uMIwa4TeknbUcY6YvycN04 jwEBG8tvjHyFEc5/H46sqseAxOi0VxYzMLpsSHF081IMqDI5xaWjDaOGH7Id 4CgGLh1uCxT5DKMfEhjDB5swsESk3bdgDJgv+X76N1/BwMr9h66fDQTmF8Ex N+tMDHzjXkFdIwrMHj65XzLLMbCtPaMiTg2Yf/o/5FUoYhBX/ulZD39gaSrq vtM7jkEn1L977hcEVjE9tLWqAoN+Wna4dijAGEVIN8M6EoO5NDIeZFFgLKv/ MPthMwbz1tTKN/DBWFEBDizTwGDRL2vENv2CsXreS3pLD2Hwpon60b5OYC/m S5buccXgkzvks0xuANvtz48vf+IwuLisc37vNhg/+Hyh7jB5XrqX38PxDoyb ON78kB6LweW1qkzqOhg/y8Nw8NbH4Fct17/NiMK473zdmqYTGNz2U7OkZBrG i3mqjnRzYQiPRvJZuSqYWAIjjlWOGGJZeaR2qT5MPHkS1eP9GUPO0h99WlsB Ey9W7HK2u4wh9kKrv6kowkT92U2fXXQx5IoVc8RUGCY6ql1dnXdgSARPjBCt BSbmr30zzH2HISVGH0xVbYGzQe5e605uDHnsu9NKtwU4G3e4b314FkMq87Ps T+kDR0Gt9XCtK4a8+XnJ96oCcHZIGy8+yoch7VlLc95ygIOXLZIXf8fQRUNH h8wCgXNGfIlxy2sM5V9WNnGJPLfxOv+j9xaGCu1a8yuQfO/cl58vFCcwVDJm TPCmHnAcbW+srgvC0K07YpXZQsDx5o2tMa7AUBXzac0FEs9PfL6/ygRD1cMt 94uReAErApJ892IotCmd1GgGTujvuMvRtRhqcvWjd3AmcOL2XVu2dBWGmt1R C04m8RIm+79v3YyhVl+yr98m8ZJzFF7kKmGog6xTVsNZ4GQ0Oh1deIShAfVC LyW2ACe/1JJafQ1DQ394vN9E4t0WLi3jKGBo1Oqur7uWAqfAqu1eRS6GJl2+ N3hmHDjFvTHabzwxtGD5MYGCJ8B5bH7kpDoTQ+8VFwYv/Qac8siEZ7YSGPpI l3vW5RdwKgqyvZcZYGi1R9WQ1l7gPH8xl27wAEPfCEueybYCTk3FoXKox9D6 21dauYKA8/LmyRwXIQxt/rblTf1r4NQez2Af6sXQ787he5QGgfNOovHYsQYM 7eXrfZTIC5y6hprz6nswlKmZdsPsEHA+iFSdVtXGMC4qT+DoZ+A0zq45b1mF YXw7rGeOET2/Gv1wW6aEYUvfP7v8SAw4zVm+fhK2GCY552bpZwKcb4J5vFYz GCad9Km5xws4bcqgkrcSw2S3KhzVSQdOu27iHfE4DNtm0bdbsAM4HYdPOk1J Ypjq9J5S57/A6dSspTlwMEwrNn1LM7lfda0a6+DXxzBt+ck8TW3CY59cB95h 2IFqw5VZNsDpLrPvaCLrD5vejVsIBU6PU3W54TkMMxzn5Tt3Ezi9qxrrXcj+ TSJsAupqCT/N/9v6E8PM1z37uXUIOH2HNS5QQzHM+ulKpwQBwp8oe3Z9wbAL Rm4D08S//bp3u2K8MMxx+LPF6SOEH8R9vv0Pw1yDFb4+dwIOfan286ZNGOa1 6toRmTjC5s99rrRjmP/DvlfhJYRzBH/N5GFYyOG9WiONhL9us3jzEsMi+tMf GkwSnts0uK4Dw2J9pzaXSgBnQGLuz2pVDEsSN8xdoU5Y9v72oCAMSy8qkvT9 j7DcfvUqGwzLOcBL6fYhLFk1/NoCw/I7z/LuJ/6j/1unxxOAYYXuz6/efkb4 m4ubQTyGPRBaOb2ki3D+o8r7DAwru+V+6fI84bNj0ZdeYVjl3s/0rzKEJTb4 9J/GsJpWRXON/STfyv+0HZwwrPbytaZM0o/9x+Mr9fwxrIG3/9B8ONGr6+tY eBqGfaHtfWl7h7CF/AL5vRXWqpGh+a6O6P0xBr4uxbA+e6NN8YKkPlFF27jP YBhjoShnSonUr9lZQYrox8rgW3HqGGExywkJot9MfTXP+gTiB7tmbno+hv09 u8o/7BHxS6Dda7GjGL7ot/vkMJk/369rb7mtjeHCior9D1cSv/l8/c+0EMPF C03VpUm9v9kOe428wPBVm4KirrsCpxUbbUPFMVxOtmW7bQPxc7tW4YZdGL4l bz70izdwmhLfdgjnY7jy+i3f9soT/4NORLI1hu9aExCwgsyDT16bbTfGYrg2 taAxdCvplwWaY74Khuuu/Co33g6chgDZPAsOhhtJbGqo30H6y+TKrc0PMPyC UOMqP9J/Lz0Fz6/KxnDH6LlLw4mkf3MbDCOcMdx1ifwLE+Lv5zXt5QIJGO7P 62u3jfRLRQf3j8VKGJ4wL1vWR/z38Fq86wwLw59PeBrpPQROTuXTYfFSDH99 OTf/0RngZG0rLypbieF1rPcz64neGWk5x6JcMfzryLrsX+fIPDtoEGFP9Bvp q2MWSwInauNZPbY+XpP4ujpC9CpwXHpUuzrD8NpqY93vAcT/Trd7FT+O4bX1 n12URsl8cjinoLzvDl5TbHjT/GY7cGxr6tw1JPEavrm8wZMOHJM/0Tb6bXjN 6fGL552HgaOZZu/5YS9eu9Ky6eZWMv/UuKKZDkvwmtcU5frVFOCo2JprLlXH ayGqFv+t7SPnh7DCWt9wvJZaMjNhQ+bHaq63BkWP8FpNsfIGxi2YmG5R2bN+ E0aI3siO4HDBxCXVX4tW3sWIFa+4L+MdmLCVWZy7hIURq/scTyYawMTpBecn mTIYsXG9loxqOkzop7ZonDHHiF3ZzZWeyjCxQa4ie3QOI2wzBFlzljDe3Gm+ iVKJEY/jfI/zVMC4it85/1JujKiSbr5dIALjGxd/0OSnYMSL4m1zh8/BuFQQ apqHYkRDHf1mHHm+2ITxTWItRvQvOjqz8jywm723o0cDRgp7rMtSFAO2e0WQ xa5YjLQze8MwcoSx+1orrh9DjFq1UcQ3cw2whO4Eekll4PXzgvUW26pgeLaD O/QEHa873D1RylsBw8yxPp5HuXjd5WA3X3c5DPfw2XP88/G6X/iPR5QyGK41 XGqUlIfXExbWLGHfg+HEnUlWn6vxevXU5Yr7uTCs5FXvc/oGRq/sWi65LRyG HL4MCmXQMHqtX44TbygMnfkrmpBngtFyUptfdQfD0ImdKrQfzzFa+b+9lykB MLT7dcPhlN8Yrdt4oZbtBUNCjs0iP0ox+kpthfv9i8Ao/b1RYPUERn+4f6Zp 2zFg8HV9Sv1WjNFNO55yC3+Dwb9NSv7vb2L0t3JRtXEbGJxszmg8b4fR/TWv M0o8YLB38Z+KkHKMnm3cbLMjEwYrRCil/g8xRm6SM6E+BINObEtPj3sYo+B9 ZIOkKwye+1RIGxbBGOU/N0/M/IZB82r32QJRjNnNY1peIQSDhz50d7BsMMZ4 RVXgblUY3HjnSNUVsv6qVriwdiAM0M/p7LitiTEhz7thvQAMdOjJ8MwKY0zk fg3XhUQY+Kr5fJtoG8YkHR5peXkLBl4b+Ugu34kxhRbHsnUaYODmkfGPFT8x piVw5daDK2HgYnDRXY8QjOlY7Hpmcx4MnPUKNr6khDG9Ee/jBRRhwDxwnfp0 AMYw4/0n32vDgMFr9rNvXRjLlddXdfQ8DKjKO12/jxir9KbosNFDGOA6UdEj YoOxO0ZRREIV6HONZjuZ+zBWU6T1a/sToP+0UpQK+Ymx+80XWZ6tAjqzKFT5 vjLGmnJOu7jXAv0b/+L/VLkx1mLFuJqWPtCb+NdvlTbGWJs9Yb/+vgf6x5Wl drNFGOsYURJ67TPQXzuaHp2extigtXyp6e1AL/l2qc/oM8aG62SbW5oDvdiC T0X9I8Zev6i6XqYb6Hc4RpNV2RibXHam8C4d6LRjS3P+KGFs4eHHz56NAZ3y 5fOjgz4Y2+Rl29/HBXSXBUNbdW+M/Zb1686tUKA7JcauciP76XwVd8mBF+gO ah/fjwVj7JBQ5fSkINBtH9gtk/uLsXM3xQR4JID+n/Mpde8JjF2ov/OxPg3o J68LC/QIIIVnfG8iRQroxo97jZx2I2XZrotrJNcB/SjQK3MTkSL35YXyli1A 3y9Hkx2qQMqWnyaTY0VA16670ls5iRTl1cynpduAvveqa+mkLVK07CT3E7/Q taSmykSTkLIv5h7v/GOgayz2Db/2CCk6pfsb3mgCfeffbQ+TBJFi8PeyicFe oKvIyYYnFSLlhCy3lGgN0JWNz3kbhiDl1MGM7tb9QFeKazY/qIaUs0m19lak HgqgXXIBkWJXYa4kWw/0zaW7VkxuR8ql7gnOMNn8JjUb01t9SPFSWOPnagx0 OecdHpZHkOJffet+BNFPdgnjUbMoUkJObOvLJvuTqayt/LALKRFDT5c/GgL6 +oCP3rl0pMT6o/57os8607k8LRpSkpa99+01APpaHYP7G8yRkn7rxL2fgUCX 1nnRnU/W52h19i19CPQ1p45rdPEjJf/TheUbSH1XB/0L67RDSuHZcX0tUg+p 569MK7uQcn/G29fwIGFhKvcToldZDNe9835AX+UWXHysCCmV66/3+t0D+kqW h4pVGlJqHi8XT+gh7OOyD4h+tYey9QpECUtd1um+jpT33Rt9qnWALvnZwaeH 7PfLlQfFzZ6E0626n91BSiufZs9oIWGPQxp0XaR0ZL4SW+ggfEFGtMUZKX3b j+hJChF26HtxWAIpjDfNPkrahEPDzxxzQArrtGWxjivhEp5PdvVI4Ywxek7f JDxhUflEFSkzIS5izt/Ifg6EpkheRMrfFb90rwkQLnZLO5yHcYuKQryzdpP8 Nsu/Xt6KcfwgWFTqRLg8sf98McYJNSd319GIHmaPIvxnME784lrR7iain2iQ 0LIvGLfy750D0zyEOyZS3XUxbm3idm9BDaJ3xV+R6nyM21K1v1sjk9SnpM5Z gTxXNvwgavCR1K/eU6pCBON2DpgcsJ0n9Z3ORgtnjAMhu7vx54gfIlV3tDzC ON28ia7bpD/k2HEid0ox7oi677JnxF8b7fdFd1VhnKlVjOeIEtC3pLiozYhj nPmUROG8FdAVD3a9a2vDOJvInC6JBOJXEY+XXoYYd+nhQx38CfTtLb8cdfgw LnRxqwiV+FtLXzzOzxfjItOt9pdwgL47MeO93hDGUbYOe7yTJf00wV4pwI9x GSa/OyYjSf8NfhpKqMG4BwXrCo4SPxs84OcL+oNxj/cUdJwlfjaKS51614xx lY07hL3JfDsRcv3lhmqMq/19wP3WaqCfptnnlVpgXOfRi/B3AOgXHtmfqSJ6 9IUo886KAP3iqFf18nsYxyif+vCD1MdJRcEgRhjjOLJBpxgkH7cxPW0JE4zn +Z3m8pn0a8ha9vETBRivXFCbm0v0yrNQ8XHmwfid3dEXqK+AfitCI2lLP8bv Wm68NZnMs8KXPB+NhzH+QGDn06gDQC+1rpFzycH406aTX1xJ/m+O/c1q8cP4 UB6ZBZ0jQB9tW1D/LITxkbuG3u4l/h0fOeIYSOLFutyL0cgD+jSfzd8xVYxP 69RaqTgL9PlL9ypcuzG++JHR9uXkvFne3mueNoTxrTYBVoPzMLBf+b/yE18x QcLKqffAdxi4rct6ZHoHEyTHNFX01GCguAF3a1ZjgpTfohD9BBh4ZJ69Y2EH JqxLTZM7rA8DL7LXiXpwYcKWjy/sDcvIeZenvHVBFRP27RGbMKPAoHiWzYVH pphgt7J8wWU/DF5vEmoY3I0JT5sWrc+8C4yETL5+HmVMdJDlVtzsDiNLFYT9 t2dh4qWnDP83+TAikZ///t8sJjofq/tk/RVG1q25vbDsHCa6+8a6ZqjCyI6l 6S7D8ZgY0LTiqQAHRv7rS4zevQMTk4K3HGBehpH8I2Za5jGY+Kzn2Jn79jCq 6S/Mc1ACk0Qy0xN3WgDzP0dvY0MHTBKzn6L+sgWm9dsD8+a6mCShbpRf7QhM h1UDfL3imCTVyP/4oB8w/YtU/Yf1MGkjr3ebRQYw89xM71HkMGmP6+m1YS3A HAnjjd0gjEkX9VfdbTEA1pWFwybcDZh0abnHI6oJsPz0HHbfc8aky32Nz6ws gBUa4qi64RQmuftHfRxxBFZSR5Wo+k1MCnz4a+JvFLDKRHXL045jUuqa7xry tcD64Snl+GgTJr36QX3tuQ/GLtAwcekqTKot7d6lHQ1jTuaiSTflManOTeYh fzuMuS9peZaUhkkfJwtyMt1gLFj3j9PdB5jUNvXE58UdGMvcsmq3JnnO/tm0 XVAUxj6kvExulcYkTrnE7eYzMNYU1T3lcxKTprxOS2cXwVib6xJhzXBM+jXT K7hdF8YGVpocn8zBZO7ZMYaJL4z9rrfQb76MyVJzAlk0OrDl7ybvXuOOydJV BuIXtwNbof7t5EI9Jq/3i7++4yqwlTsO3Tgnj8nyfyS9aiWBrdUiIRMjgMnb /8ofHz0KbIMNhQ8LnmOy3jwK7HwMbFdpf6/qJEw+VBMW+G8RsD0sZVipaZh8 NPDdz3eGwPahhKZ9asBk44VjA6dHgR381sLbeQ6TLbksagLXATshKnlNlAMm uy329qiPBPb9UxvVX6tiskeymTp3I7Afeq1fcV4Jk73ld/+E1cAuizKVCtfG 5Kt6855PHgC7MuICbzUTkyOuX/PKbwF27UhB+9EFTL6+2l6jh8Svu7HCJvoA JscUH5qVsgd2g8ENI6l0TE74JOQdPwfsLwH8EybSmJwlmuwTsAHY33VemqVU YXLODU+tCkdgdyY3Ggh8xOQ81f9+TZF8u9tGpVwkMPm2yWpfx0PApu9SYr9o x+SS9Bt+p5yBzbRuwJhwTH60JWx3UgWwxw6NjZ5LwOTHlefnPhN9xjcuFT0Y h8mVnVv8dVOA/aNcPbVaBJNr15Vc3fEM2L/2Kn3UKMTkdyWJey/zAPt39fJP 0/6Y/B7d/hYSPf+o/FgxS/T7bKsRsI7U89+POL1l45jcOLVqn7kSsBe0Dc+J RGByc9jvf6meMM4VyKu+YiUmt9+uDhQSIPeXb6ffvr6GyZ0audoHT8A49/hM tVQRksMseD40C8Z55ijfa6cxmT6qG/hbGcb5mMWmHo8xmeG3CdR8YJz/y072 3TBMHlnKv+DyGsYF7lStvWqPyWylD0FDpjAuqPx2hRMNkyeq7+MGGuHeYys+ 82DypGE8l+UIjC8N/b5u71VMnu51fUnmxbiQ1KXVSj8xedb1RHDzVcI3eSWe XsHkuUVq+0Xewbjw+rta14lf/yZJLjqyjHC82Zd2Uq8F2dlX104TnpTSPT6H KYvKvoe8zIdxkYPM9aIJmMKj+2z/nzHCCV8ErqdiCl9rzmJNdcIfG+qDhjFl iV3ga7dAwn/69Ka1MGXprE3o/XoYX7ZOPHKFM6aIROnojIoT3nm+VKUCU8Sk 5Lnlyf1x2Z4OgRvnMGV5Ee8b6zuENT05L89giuTu4bDMCcKb9rCavTBFiltn bkGUMO+yeoPnmLLmA83tHMlXpKXFxDcTU9Yl/x6tMyGcGhf+jMTfYGl6dqvX //d/wLvVF1Pk5R+2x2eQ/EZ/eX0uw5RNbCGjqUrCAY+m2wcxReGJfd2pDsKL feJV/2GKUsAb7Wd/iX4+RsJCVZiyXX9ducw6onfvXmuHa5iiKuJHroOENfVv LGZhitq31lsjNqQ+wR6SitmYopm7Q9ogFMaXVH1QKp7HlN32sckPb5H6Dp3w v38SU/ZtH1kq8ZYwt+ibhAuYovOS9rub+IlP1FCI/QlT9KLmruxXhHFerk5L 9WJMOXTcdOT2UeIv+nOoKMeUY3ShtstxML7IL31daCimGBddNGx6SPypknfv 9AdMOen25p16E7Dn26WndI9hihm33+N5CeL3P3ExIs8wxeLDNyVbDWDPeac8 mBDBFKuUHTffnSL9wti+I1ccU87LjyTFZQL7Z6ruve9NmOKi/5/rehnSb4d/ BH0SxRQ3kdLhMB3Sjyep15+S+nm0CVsNnyP9aiSTNV2HKX72tcdK7gB7aM30 UkYIply7rroVST91Cze937MGU6KOU/Jvkf7rKEiq2P4NU2KkRlcLuAC7fWcE s8YJUxKK8gQaHwG7WVb93yvAlKyPwkNntYBdb/FDgoGYQktxOPPWHNhvoyj/ iTVgyg3L2lYFf2C/Lqgy8a3GlDts/zc/aoD9vHr2h0orpjwSGb0Rqkvm39Fz WiuXYEr98beWN48DO+Wnz5IDsZjSwHyx1tcV2InxGVNRhzHlU1hFj1EcsOPW VY90FmDK1yf3rP80ADtS+gAXXxqmdEslnzMh+fukRkU+c8OUH/1nHXlVgW0W U9vaS/w17Weh2GkMbFP/zizJaUyZXW7CKiXz7PjZagfpTZjyV+/Q5TNFwD68 bCXfOTtM5bu73bVcFtiaP6j2a90xdY3rP297cWAv7+9wGHHC1HVLZrX2qQB7 mfqidIk7mCqTz/m9nOgnGNi0VjYAUze2DPq9iAY218zNU7ftMVVF60PgqsUw xvYpefUyCFP1uKjh9ZMwVtscbPqPfM85QT1J8SuMXTZ2COPuwNQr0pKveDVg zL5MvF6fH1PdC2YneqkwZityY628Dqb6vKg6lmIDY/89LH60cimmhk6AwL8J GIN81UD+PZiabngoqEkIxkQPzRc8PoupL4UsLnvrA+v+yhs3+SYxTSwq5NCb KmC227/qKTiOaa+td79PyoYRNuPk9qM7Ma12audv/hYYGXrHWHFxFabVRW7b EiAEI71FpqmMKEz7WLI+0v4qjDR5r7x68SGmfVvg0dtrCSNPFLnODPZiGov2 6eWQNIwE8BjPWIxiukSPTcUuEl9k4VfSf56YvtLNbPgBiS+gLfjJ0wTTV/Od XClP4i8KXqb8sR/T1yvreYlcheFp0W1h+wDTFQIVdw5YwvB3W2+5yS5M1147 VUKRhuFbXa9njacx3f5MxJ2+HBiGRIU/TT2Y7rjlNyXwHgxrDQ14Ol7FdKdp Jw/p5zCsesDlpWIFprvFnNh/uguGN4rpPyX36/SAyrWdTSTe0kVf1zRkY3qK xBPRV9kw1K63xT84BdNffhjwy8+EIZ+3yStjvmD6m/T/zmIRDLmviPcr1sX0 d7bvD/ZUwdBlZ3PqQyamf/hdsmJ1BwzZ7k4VjNTE9LbNfiXJq2HIAKRFQ4cx nR0mNhhGhSFZnlXv7fswY/VuNDyfAYyWljhVmhtmSH87HuhxCxiNN7nTlyhg xno32wfhpcD46OczvOkxZsgXXRO51QCMWt2wSa3vmLF99fvPg3+A8VhvrqJj FjMO/DlpdMEKGGm8yzb/rsQM/fTzQZ6OwEiy37AQ0okZh3d6llzzBkbc1+I3 LwUww/BSxrLbCcCIfKd96IE0Zpzu7PrCeAUMPxkx/kPVmHH5uZ2xnTwwbG7e 3lV4AzNcT3sHe6kA40ze4oFkecxwm456GLEPGOZ31BzobzHDR6lY9M5/wDBp 3FQkqIEZodkTjUORwDic+cHvvSpmXNu1iGsmBRj6Qwa6wl8xI7JVfAdfHjB0 97x1W6KGGRRhtcRNFcAAnlK9ID/MSA/wPW7HBIa6SMeUkCxmZEpFh3jNAEO1 KOZt/B/MyH6SVRrJDQwV4xPtSPZ3g10jViANjK1VKiWOBphxz2px07AhMOT4 g/o7WZjxYE5i0aw5MGSmtZTVxzCjNG2jKr89MNax1g72v8KM8s8HkzYHA2P1 rL/CBaLPK+3YE/ZE/+XRVWeST2HGm0/hnKs8wBBralYoo2DGO8uAuASSr6hs rQD3dsz44OfcUDEHDCH6TnrlXsxoLT9+QFAHGHxtrlwhAZjRrne4fx3Jn/fy 5bmFq5jR0aITpDoEDB7hY/+tFcGM3smdzyyigbHIsWyx2FbMYCpL7rz3FQb/ Zm+hBV3AjLFqkaZX8jD4JwSfPu3BjAkDfpdWLxicc91neWQ5Zkw7/CqeXw2D vy4NCoWT/c78+nFE3AkGZ30o1Y3HMeN35OjIphoYnIlfo/eB6D9/u1PeyBYG p7/3SbUpIXWRWvPrc48Ji0gePhGEVO43H218+GBwykBZ8rwYUvlO1M7HnobB yZRN5YarkCrQX52dVwSDPxgL+juTkLrUtXz347+E95d/kBlHqvDCg/b3hjDI KTBQleAgdfnaPIkfkzA4kbC0fz4bqSvuZZTx6RJevpPCdQ+pq/YkHl+dBoPj eTt71MORurrh+oTyCGHNpUcG65G61iyUcmA3DLLbKlRCPZG6fsR/66lYwqFa LZVGSJX1dn9/qYewRswFz3ikbuRzsg9WgcGx6fu/Jp4jdXPqBd6UUMLPaZSn r5GqKH/mZmEL4TgLwZf6SFUq+29/9SbCDv3gKI/U7TqGvU0+hI2UBxK2IXVH 08EARgNh7f3pyeuQqmaDa+akCWtKDl6oRarGhFaViDPhXff2UM8idVfgjtOy Lwnrzc867kXqXmGFGU1xwpZC5h89kKqdLZty9Dxh/6ZEn39I3b91jap1OeHb R0vmupB6oGp5o4cA4XbvsueJSNU/LOQcZU7yW3EibEU1Ug+18wpl3yNs0c4+ PYPUo3b/ih7OE77HX9TfidRjP2cOvTUm+vF2npEzQKpx+PjQ93zCDsZxT1hI Pbl8OJw9TfjbuTI7Es80v09usT6pxzGxLfqk3qd3tL+SzCD8yXRneRxSzV82 WSsyST1PKfFeDUTqGaP3/2AvYWZi+af9SLXueZ11Mo7UPzL0iDU/Us//LWu7 qkr80j1Vtug+Uu3fbFKqtCR+KkjcqkjiOV7PCJmJIH67apL78yBSXSUDlFw7 iD8PajwL24VUt66JkAc8xL/7DpeNiSPV8+bZbyxl4m/tSwaeT5Dqr6IfcoHU c+7sq7IjCUgNmKn4lk/0+BM22jfyH1KDq7du7f3//+9L+S+HqiE14siyb+aK MLigJCC1mKxPPPd9q2ET6b9b1m4GkkhNUTwaHEP6l2+tu2ytG1LTONWt78n8 47+hu9aBgdTsgPzgAz7AEKxfeFbXjNSC9EutWjJkHmRwd083IrXoTI+i5xEy L+y2jGVpIvW+vHHQIw9grND+MHXICamPStUVleqBIbXMcVP9KFKrG+aDNrgC Y8M2ibrhZ0h9mejaYpVJ5pWh5j9ZUp83pwYUsmqBsdE36bvTVqTWD9a1SK4C hsJ4g5lFK1Kb5xMVBF8AY+cJDZsniFSmqtzXaRFgHJy2aA6zRurYr7QtO7SA ceSxfabpIFInXggEONsC41igoJ+QKFJ/HmVvHikHxsltoyP/5jCT60L51W5y fljNt09F38ZMSerhTe/uA8PHwmNMcQ1mSqmYvfVsAYa/V+T8ASHMXFPneF6e 6BVINeVWb8JMmZ8xN8L0gRH+c/qFCAsztx7/LA195PwRVDYe0MJMFDwpVikB jPupLv9dKsdMB3/Lv8V+wBhlKVIeFmGmk/jlLHNyHrBHeoV5ezHT+W7g7iV1 wPgxNt5aFo2Z7m253hfJ+3P8nGfXbDEzcGf/pPx9GBKMUj6/bRIzk8bshnP6 YEjpDrvQXBszn59x/ZqgD0OefML14ucx88XWP1su6sKQnzRX3MYUzHw1dy0Q dWAoSDNRUuc5ZtalZypwAIauBx3a9mAWM782vQkx1oIhmpH6IaoCZo7qS6qI boWhevq6z7PcmCWp/IySIArDG8zpQ7E/McuVi+9Q/HcY7nzaYvm5CrNlLL0h rwVGjTf8PbBzHrNlxyMNJiRg1NRGVOpiGmZvDE4332cKo+ZZaRN/1mO24q2n Ht+/wegFwY9nzOYwW21stki8DUZ964NybKMx+3CAj0R4B4zml3t4DAVgthvN Z/RiH4xOCu4q2dSD2Z4qUTNPZWB0pjdTJJqN2d6vqTy8Z2H0z8Nf0o99MPvq UOW6G3Rg8h4MeYaumB2p/PtkxwAwV+nkrOpWxOysF741BkPA3GdfuvSYPmbX 9vklq7CBGb5pmXprNWbXFfGcKuQDZtTrG9b1s5jd4B4nvX4dMClmi/JtyPca +W7cFjYEZprbvLxmG2Z3KtWXjz4EZqF5wQvvAcz+4SvZke8GzIYJwYmIdsye 1snLlYoG5hc52LCzHLNnhRTPJeQDs9k4SbvgLWb/zd03FtgMzM40pxhBXczh f3t+3kINmKxP3283jGHOOrEy2eUzwBIc2Ch73hpzZDr3DsWIAEvkVtkQSGKO 3K13xYs3AUvc6n3Y7iOYo6DZsZNzClhSz+0orErMUbNapPexAlhbpP22xqVh zpEiI4dr/sA6cJHFKyeKOcfcv2/7mwIs/RKBErHTmGO8z3bS/R6wDrO3RTws xBzTRi//s13AMj7uqnvRFHNsZnIo2trAsqxcpW7ZgDneOqyHs1zAcl/QWq5f ijm+3ctr+lYAyyt5otq2E3Oueu/5+F4RWL4bLD6zP2JOSHHMSNZJYAVugp1m IpgTK6EkA3eAFdlWWVhmiDl5DKf4KANgUVun/jz2xZybQSk5V84CK2ua62Fe MubckXpebO4FLJpgGK+RK+YUGwrVKeUCK1+mM1PSHHOelN+bb5wE1t3JCOXo UcxpiGQ7S1GBVbH5uAy3F+Z8kl1xdfEDYFV+sPzX8xdzvjzfd531GljP7Msa HRYwp4VDuV09Bqya67sv1S7DnF6zbd22CKxaF58EOzXM6Z8yZR01Bdbb1sKs 1YOYM0gJ+K3mCKw6tRXtki8xZ/TV5xV8ycB6/z0/ksyLnGlFZ4OiIWB9rjvn vdgCc2Zq08yT/wDry4TFgdQKzPltVXMxQBRYTcsioxuFMWc+WSTcaBewmrX8 ihPXIW3RNo0kLUNgtezHQ86hSOOus8rbYAusVh3NoQwzpAn8efBsKhZYbUrP rg++QdrS1G/vu24Aq11SfzJJFGki2+fb3pYT/sUtl05H2vLzhlPpvcDquCUc Yp2CtLUfJ7buVQVWN+3S73fqSJMxmy9QIPF69NV3HLFDmuywiJzkMOGhwP1h HUjbsmjb6h/ZwOoTTRzhe4s0xbi9aT3En32ZZ+5+/4u0bWsMxD8aA6tfujjL wRBpquqXBAt4gUUXbFJ6twppaq/9Isj9jkX3Krp+MRZpmkbRi0OqCHcuiT8n jbRdXdRAFwlgDWhOHDHORtpeh8I5S2fCFMu2ryJI056p8DpcT7jzpH3kMqTt D6ub1JQD1uCGlj18RC9d0TZn+UDCNn23W54gTT9niCnWTjjDT9h6L9IOK85c WCD5DtbfGOF8QdrRCt7+MQrhH2avirqQZqi34sz3EWAxxG+qXs1AmvFX+e91 OoSVgkZWH0HaSWs108ckf4b2j7xpZ6SZjh1oukHyZxye26ZcjLTTviePxZP8 GceynNPuIM2Cz/b91eL/P287nVmItDPJbnqOfP9/v4SDJF8bmZBXp2z+H3+z gdQQ0mzvJ+7Trfr/9yFFOwFpF3bnVaoSPQY5fxZrGyHNvu6h+nrn/+/fon8V qaejyYtSoff/z+88Zace0pz6vyjNyRM+u9LjmDbSXJx7Cof/r4es9wqvOqRd +TMu30L0GOiKWMRL1ntE/ct7vZNwvH57VAHSvCWFpUviCO9+bviS1M/3pnRG NtGD3jO6Of470q6qKElEEz3o/o2aU++RFli9J8Gb6EEX87jwcwppYW3mUSeO k/pv5uL/NYq0iPOOPEDmSd8ddvaQCdKifvgGbyP5962/lbhTDmlxSzN8+Em+ PX/8k6XmkZaQUTD1k+TbY5V/aGsJ0pI3PnWhuxB/VuUWuxF/Z8A3+2qSX9cZ pdAaFaRlfmQMFJH8Om++kVN5gbQcs5/W6d+Jvwe0GzLvIS3fQ+KUG8nnu8H0 3rXKSLtXeOLg5l+kfz7UZsVtQVqJ+tk3K06Q/vo+MbmB+KX09RVYfJ/0X7+M WPplpJV3JWh0k35r6sitHo9H2ivRz5uSPgDrk196XZ0/0t60njIrnAfWR80W iZ1PkfaO2hf7QgVYDWzu7TzTSPsoOznJIvO1Hk4ojugjrU195QtdK2C9Drc9 y/8Sad9/502aJwDr1aKBB78eIK2rRnGjK5lHL7yjC6YEkUY/uDcmZzOwnh89 n1ZO3mebnT31kwOsJ2lL3uawkcZZy4wRkgVWWauv1ultSJuku9fImgDrkbC+ dWoq0n5dipA3rATWA7uVht7BmMsdUMQpCAXWnWq7+IYSzJW6MX3djOidtiRx 7Ye/mCttF/DcRR9YyV99rbXDMHf9Vr6Jaz7ASkzRXl2mjrnyj6VMS7uBRVkc pL1bAHO3v9OWFSTzOyy6uU2mDnP1mJFVz8n8cjmR/WQhHXPddq4eW0/y1eH6 suFfM+Z61OtcLCD9BluPC2z/jbneZxwY258Da6/xbTuKNOZejazoQSlgacQM 91ANMDei679GmxZgKb59d1RBHnOzriU9zj9CzjtD2i+vHMytbRcM2KgBzO7p C+e492He0poA/hWRwHQofC/8UxrzhOU/5Y0LAvNCYLiN+EnME4lZt/tdHDDP Gs+/phVjntipF87eKcA8PXahU6UU8yR/LGol9wfmwXk3L6fNmCcrF3Er5zkw N2ZpK/nwY55WVILOxkkY7ZXccffsfcy7cOJmgIYljBoKmRzXmSZ85nT7+xMw emzUSTXjGeGLImqWhwm/C1/1PZJwgA8rVINwcNm0nwbhQkOzz8tg1ODncmaD LOF/vzUuvIbRo90tI2q2mGd313AidTOMHi6xaLhLI/yY5+iWdYRjPbyGyf7s airvPCO/tw47UBQXrhFukT/Tz0VYrqZ4XSrh+d8flcjvqUMZssIFNphnb3Kr sDYWRg8GM10vVxO2Nuc5FUrYivu22SbCjsusmb6E90gPpI0SDvJbKWoHo/o/ ja7aWxMuMrpmCYQdnmwY+0P4CU/fBMlHX//zzb/qhF9U7Q1VIiw3oiIVSLhV fqpwFYzq9ayx61pCuLfDaI8I4WfqYg/VCDPjiz7zEs4wfMfZS/inHp/NX8Ke FyPKUggvzJ2dInrrnQg9Lvb8f6iLUMw= "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, PlotRange->{0., 2.}, PlotRangeClipping->True]], "Output", CellChangeTimes->{3.44184228083947*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "Calculate the ", Cell[BoxData[ SuperscriptBox["\[Chi]", "2"]]] }], "Section"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"chisqmin", "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "104"], RowBox[{ FractionBox["1", SuperscriptBox[ RowBox[{"(", "0.018", ")"}], "2"]], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"amp", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Gamma]"}], " ", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}]], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "per"], " ", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}], "+", "\[Phi]"}], "]"}]}], " ", "+", "const"}], "/.", "ff"}], ")"}]}], ")"}], "2"]}]}]}]], "Input"], Cell[BoxData["189.09190118053925`"], "Output", CellChangeTimes->{3.441842281111742*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"chisqPer", "[", "perf_", "]"}], "=", RowBox[{ UnderoverscriptBox["\[Sum]", RowBox[{"i", "=", "1"}], "104"], RowBox[{ FractionBox["1", SuperscriptBox[ RowBox[{"(", "0.018", ")"}], "2"]], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"amp", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "\[Gamma]"}], " ", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}]], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "perf"], " ", RowBox[{"data", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}]}], "+", "\[Phi]"}], "]"}]}], " ", "+", "const"}], "/.", "ff"}], ")"}]}], ")"}], "2"]}]}]}], ";"}]], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["Find the uncertainty in the period", "Section", CellChangeTimes->{{3.443293870429606*^9, 3.443293876811593*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"chisqPer", "[", "perf", "]"}], ",", RowBox[{"{", RowBox[{"perf", ",", "645", ",", "670"}], "}"}], ",", RowBox[{"PlotRange", "->", RowBox[{"{", RowBox[{"170", ",", "350"}], "}"}]}]}], "]"}]], "Input"], Cell[BoxData[ GraphicsBox[{{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV1wc0lf8bAHArFJVRSsOKZJWMkkpPy0gJIamsxvsa4apLZWS8yB6REcqW GaKi/EpJtn9WigiReW9mkvo/nOM4n3OP+37Hs15DR2vxKywsLGz4u7T8FyiF NdaZdzWB3sey/AOU4k4fjgdcQFcbdylhrQVKRbq1T+ID3HBg9c63ooBSW717 9W9FuF4425y7RwaoA1pCP6/IgvNYhWSYbD5QoFZL92gH54P9gcUJLUAdu8NR KCsItOB7woHeh4DSLtLdL/ULnIae/OFdeA7UyQfmRHMYOBn/nMz0rQNKv+v3 Cd4ycPTLjrRUvQ2UUYt1w3wEOPhdlo2lvQDKtLgjVYkfrj04P+eY8AmoCyOt CiySYP8ksclEeAtQlp8STd1+gN2wN2v99l9AXdk8uS9YA+zEFpbKm6qAstHV tR89ALa0okuP1xkAdU2JDnNjYNObze7WEAoUzUG6LlYebIyU5SbOyANFt8qc +S4K5EDF7Z7kl2jnmnWXM4H8thBqyWxF+w4rD8oD2afWa1VVgM6Udh5QA7K7 4nCVNq6HPprF/KYPZHvTP11aAlAu17PHe7yA/KB09lsC/r+rX85gez+Q+T8M xRcigLr1uKDtfTqQHnysz2Nw/7fK74iH2wDplu7WJdSDrj/tYLoTyNvqO6yq JNDjTK6RZ0C62myMmpYF6rbi7v28dUDS2o0+qAShy0pSDBlAXh3euzGzCyi3 t8+v9aoDacC8yfc3Ct12tzzrH5D6ucqQlor+fo7L6R2QpwkxgdRzQLlzLqSw nQLyZH8Y9Rnjw/2EepuUJZBaw69mTKzRLa/U7fyA1JCzCr9xASiPnirO+RYg 5WIjJIzXoSdo/kdPAikbnJ/mtwa9JL4irAZIGa82l+wnQHlu8eKQqgRS2uHg X/cQ9HkN1jN5QG4z9uKU/IP+9GKxMADIzcYVvsvxeqe1iEloAMm7ZTwicw96 wMqxuBxIHiVZ+mFb9DTf5F9VIFdp07i7ACgvAceJGLwf7psbd4Raog0URt9u ApJjoLG8dxTd9Pi7yBwQS+P/K+LlA8q7NqWrIx8IZrUwx6GH6M9pbr9cgGCM OkcX8KLHMkQ2AZqfzb/fHyif1dnWF9uAmLCqdrmE8ehjUDA6+BuIUf6OOhFd dFfF4rQ2EAN54n6TOUD5jnSIrB0EorNWRiuNif7d+Xp3AXr9EeNuE6Aoni7r MzeB6LCOmDadRSt0Z8byANHO0vSmDJ9P0QZ2iSgB8VG3fnrNHHqBCfJeQDRu ft/rvwiU/6pV1jpbgHibFhHAh/vxFx7bM5GIFk2b4ilEyzTwRG0Foirxnecr vG9/ndCnX0SAeJMgP5xyGh3Iz3VNDIj/UhIsvgYAFcC9IS9SEoiK/gGFlXhe ARt+eaumox3fyzy5iJbuMv4sBUT5UulPC6wPAVoJ/7ZtR2+paNi2Ch0gYlC6 A4jnhOUGDU+g7nJKznXJA1G6i/PIgDhaaEW9Rx4QTxtKLp89j5YaeiShgLaz X+jHeL57PFvHbicQJXkr7OgzaD+5B0uKQBTvq2sRCQYqkEPpsLgKEIUBqZ+7 76HXCa6vfores+4JK8ZT4LbpUVtVIAq+R6VyOaOPlt4r2YPWKn7fi/kV6Ks2 dFwNiPyNvrdbbwAVxHYo2OYgEDnjVm8D8fyC+FtCTf6g8zr415JocavwYy/R DrrlgPEbBL73xNSBeDyzn4cF8znIsyahE7//Ma/MemfM96Aw08R3s0Bkt2UV Nt5HJ40kF+P6spOl3f1j0K94UsPwPrNVZONIzK+gxdM5x3cBkWWvoaOA9SZ4 VV+e0iT6QPX31+PoTbQCMYyXrDV6wSYYD8H77hX9wfPMfGrd3ZqOdu18UYLn nckbyy3fgA4gKx4NA5Hxffv5fn70/YVXYZnoyrJwhd/oss1vbPF+M2503Wl+ h56xqBUXByJ9YscvoQygQtiZ9Wv60I2Vz8yxnocIejf8eYQuNJZ/q41WSmv5 hPGTTqekRM3RR1Q+VvegzbYklWJ9CTGobivB+Es/VPos4wjaafhT+Cb06pHy +ji0183P7l1ApM1R8clYr0MiVnbbxqH7xFaswPwOKZTr0xRCl5nav9mLrnzZ r9yOTp01L1peT9OpQfFodETUB9tp9ITj8JIg2qFp0nY9eollZPQj2uLaIx0e oEJ5I8c+RaENeEslsN6HbpGYeG+APponXYXxGSpXwnjKh96rOzWD+R26/9jP lGa03Dg/tQN9on06PAwtHuJ5YPlzs6uzHqfQGxWUf+Shbebn7XjRa5sUzpxE 37q7YFqP5na8ormEDhRe1AxCs/F1hC6g43KWVHSASP1T7Dalj87e/0+CG/3L SP9gJ/pZI+vaGvTMnJ76NnSNOduSP3oqzulJOLqDwTGmif6pXmy7/LwhL86u FWhmzxpF7F+h85fkkgSWP4/nc1MCKoxbU98K8zt12sdBFftNmLAMXQrzLXUu Xqy0HS3LGz+C8Zz6u0fcEOM/TJ1RmY/1K41Fz8ZSDq37cYBmjOYcn1LF/hJ2 oYxb1Rq95nmptCv6Wrz8ghN6Q1Fu0CDa093glefyeXU2PcX+HxZu4eIdjFaQ 2zCP9TLs0ZEHx5fve1+u12OsR2FFUq9XYnymaZ1j24rzRVgV92Aj5kuaiXJM xgf0x/GVkW/QV/er3HdAzxSfEe5GB/7v1m+Ml/AVMa49o+gkS0H+12ihm4kp 8+jizYmf16LVNIZk8TzSeoWJKie0jsSqSVH0/PkHFJ5/uNmKXcWYT+l8dRnq lWi3+lvqWmitXaLV2F/CQwqS/hqhL0v6qmE9Ck+KrKrC/af7aKb0Yv6E/2fK e8IdXcU3x4b1LYJt8KFZMebX6Vca0u5owZp3Iq/RLhv+e4v5FyGZ86O/EZ2U XMMmgNakKdn+QDOVdO/9RQcuVbthvc5Mc6CZuKAT+kYPyaJb5S9rY75H5L5b y471K4tDKGsP7jeiMcg02BBt48S2nN+R/EJjD+5ifTmiqh+E8Rsp/pvP4j7a zZjLwgut1KO6LQ1d+uiT6wG0UZpnbiXWL4WXBcG43sg4Rf6KGax38hwO7Bhf UaI6e7stgMizbDorcBm9iyd8G9aPvKzIF5zdaI3GIRtcfx7DTEgC62+UuWH0 PMZDvvdC7hCeb1TS+SkB7FcFJdfzlPD77m12yNNOBaJIS0+5aDdabjd7+FV0 snmZLvbve+rTZh2436K5692ntqLNbq68VIL1/3Hunwg1dLz3Vc932C+2GM4c wvkmekO02NPvQDxTcmEO4rwXLW3iupiDjm8WomM+RO/d2HzEAfsRm1zAdVxP tEmSRzP2/+ftP1WLVdAxmZ/HOLF/BRbKsPwEKkbwRYw41tNKcWac0SRQ99f2 rgy1A+L9euuFRzg/3heb9s3B+H6fna0bZY9W4lqqiQCi5qBIOC/Oa/eNFKdY MX4/2CQ12fWj43w+0/8CUdcyZOP6DahY8e255liPmj/FD07jfccpOZzaVQ7E 5wZTH2YT+qhPzckGIL74FZy7ivNBnFHsYZteILpB5dDZerTLa9V0DiB6nrds NcX+G1chILJRD4i+51rnSYyX+KOlDJZ+IIaYCrleX4BKMFqK+LgSiJmjWSOz WJ8TrlSpqZkBMSvypRCwXye4+PcmYz+eXVR+zMmBjuPbaYffN/9MX5QP8y/h i0QDB6538aC6qR3OYw+stLn24jzFdqdJ+D4+L/FalEcizlv8vW+jGp4Clewv bXtVD0iVlXbkSBs6pv7YaRMgVYVrFjVxf8kZjiJq5kDukRM/UbYd/a7s4yoH INX061aV4X4esmvuKwgF8kDW9I0k7J8Pva5wzzUAedwvZdNFXN+j26lZATpA nv0Lc4pYH1Mdtw5mHwPyzsCFyB94X6legt+IBCC9T2/eug3jKTVy5dftTCB9 XrIeJkzRxbOdGYlA+sUfU0vHeEida6xLnQUy2Cp9M6cmUGmeHk+ScN6Plal2 VfsBVHpoj/s9biALYy9UmmN/zMxLWueB6xkY0lNPwHzKSfnxnZUTbAzaHxtn Y38s8AWu1mGwaeRtas3CeaP4hWifuhnYXpoa/t6C9aIss+Gc3h2wY6nlzInr AKp8hDVgbg/YvdwdeHoA58PKAfazuslg72a/Seoy5ktVLd9NZj5c0+RM1J/F ea6GPVNfRg8clBLnznXj+1PdHPeHtSngKLq1nPn1DFDNB6ZnA/3BaaPyiYGA XqBaHUfpZwTAKXMw4yW7IVCdZvkxdfpAU442aivkBupL/IJduxjQ2owL5w05 geqN5O8PTQFnZ4HwkCc4Dw7cUM2NqYXrgh1JYVKYf8O5/bK7feF6vaaq1IUK oMYenNe5zYAbHu0P5R9i/jF6L4qcnQC6lmVtc43XsnuzG9LRtyxtvn5F91kr HzZDR85LKOL8wPhGaMu+R7+VnPDuQw/QaEuPgK4t5T3Zi/HPGParTjuDHj2c 6I7vGwxGgf3kC6CfuN6Zthr7D4OpOH/JCR3izGHxZ9kl3p+2ozN4fcdwnmf8 fH7/TRS643BP4Wb0dNXraHug66rlLw5hf2DMd65TF0UvBQ/w4rzA+HUuuaAD 6CeFdkj/jF12t8y2EPTOd8F8OK8xFr5p8P5CWyzlzuP7L2NxjOj9iK6imdJw vmP8uTZldBf9ZU341u5lM91rD6GncwfZNdBLs5HFeUA/JTn0wxXji/H35hZp a/RBKo5Lctm/MxM3ok0krjQvP++fx26+JrTDa4Nm02X/q/DzQweYX+DO/j8N BYiG "]]}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesOrigin->{645., 200.}, PlotRange->{{645, 670}, {170, 350}}, PlotRangeClipping->True, PlotRangePadding->{ Scaled[0.02], Automatic}]], "Output", CellChangeTimes->{3.441842282244458*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"fper", "=", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"chisqPer", "[", "perf", "]"}], "\[Equal]", RowBox[{"chisqmin", "+", "1."}]}], ",", RowBox[{"{", RowBox[{"perf", ",", RowBox[{"per", "/.", "ff"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"errorPer", "=", RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"perf", "/.", "fper"}], ")"}], "-", "per"}], "/.", "ff"}], "]"}]}]}], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"perf", "\[Rule]", "656.7654201729828`"}], "}"}]], "Output", CellChangeTimes->{3.44184228264686*^9}], Cell[BoxData["0.9711076011187743`"], "Output", CellChangeTimes->{3.441842282685212*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", RowBox[{"per", "/.", "ff"}], " ", ",", "\"\< +/- \>\"", ",", " ", "errorPer"}], "]"}], ";"}]], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"period = \"\>", "\[InvisibleSpace]", "657.7365277741015`", "\[InvisibleSpace]", "\<\" +/- \"\>", "\[InvisibleSpace]", "0.9711076011187743`"}], SequenceForm["period = ", 657.7365277741015, " +/- ", 0.9711076011187743], Editable->False]], "Print", CellChangeTimes->{3.4418422827817583`*^9}] }, Open ]] }, Open ]] }, WindowSize->{754, 692}, WindowMargins->{{53, Automatic}, {Automatic, 11}}, PrintingCopies->2, PrintingPageRange->{1, Automatic}, FrontEndVersion->"6.0 for Mac OS X PowerPC (32-bit) (May 21, 2008)", StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", CharacterEncoding -> "UTF-8"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[590, 23, 145, 1, 89, "Section"], Cell[738, 26, 8319, 212, 522, "Input"], Cell[CellGroupData[{ Cell[9082, 242, 62, 1, 46, "Input"], Cell[9147, 245, 72, 1, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9256, 251, 336, 9, 66, "Input"], Cell[9595, 262, 1267, 25, 324, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[10911, 293, 114, 1, 89, "Section"], Cell[CellGroupData[{ Cell[11050, 298, 912, 28, 123, "Input"], Cell[11965, 328, 409, 8, 66, "Output"] }, Open ]], Cell[12389, 339, 435, 14, 71, "Input"], Cell[CellGroupData[{ Cell[12849, 357, 327, 9, 46, "Input"], Cell[13179, 368, 19128, 319, 324, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[32344, 692, 108, 3, 46, "Input"], Cell[32455, 697, 20106, 335, 324, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[52610, 1038, 98, 4, 93, "Section"], Cell[CellGroupData[{ Cell[52733, 1046, 1105, 34, 194, "Input"], Cell[53841, 1082, 88, 1, 46, "Output"] }, Open ]], Cell[53944, 1086, 1190, 36, 194, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[55171, 1127, 119, 1, 89, "Section"], Cell[CellGroupData[{ Cell[55315, 1132, 283, 8, 46, "Input"], Cell[55601, 1142, 5199, 91, 306, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[60837, 1238, 513, 17, 66, "Input"], Cell[61353, 1257, 140, 3, 46, "Output"], Cell[61496, 1262, 88, 1, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61621, 1268, 200, 5, 46, "Input"], Cell[61824, 1275, 350, 7, 25, "Print"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)