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

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

(* CreatedBy='Mathematica 10.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    125369,       2409]
NotebookOptionsPosition[    121901,       2344]
NotebookOutlinePosition[    122258,       2360]
CellTagsIndexPosition[    122215,       2357]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Torus
JMC
natural
created: 11/2018
last edited: 11/08/2018\
\>", "Chapter",
 CellChangeTimes->{{3.741392324820384*^9, 3.7413923493322353`*^9}, {
  3.741443109839137*^9, 3.741443126207752*^9}, {3.741446680536162*^9, 
  3.741446681327345*^9}, {3.741453500214444*^9, 3.7414535107803507`*^9}, {
  3.7419567182029247`*^9, 3.741956721160833*^9}, {3.741968004781842*^9, 
  3.741968007327126*^9}, {3.742301752606923*^9, 3.742301770299365*^9}, {
  3.742397232261779*^9, 3.742397232792536*^9}, {3.742911384949459*^9, 
  3.742911385875022*^9}, {3.742913298217881*^9, 3.7429132984790287`*^9}, {
  3.7456841396848087`*^9, 3.7456841496020823`*^9}, {3.745685390088977*^9, 
  3.745685390327573*^9}, {3.750623575326748*^9, 3.750623580776507*^9}, {
  3.750623635507696*^9, 3.750623647658661*^9}, {3.7506887519796352`*^9, 
  3.750688752351657*^9}},ExpressionUUID->"ad850729-9755-4498-84de-\
cb71834c15f5"],

Cell[BoxData[
 RowBox[{"Off", "[", 
  RowBox[{"NIntegrate", "::", "slwcon"}], "]"}]], "Input",
 CellChangeTimes->{3.451259087372103*^9},
 CellLabel->"In[1]:=",ExpressionUUID->"653d9c4a-289d-4333-86f4-de334ba4329b"],

Cell[BoxData[
 RowBox[{"Off", "[", 
  RowBox[{"NIntegrate", "::", "ncvb"}], "]"}]], "Input",
 CellChangeTimes->{3.451259087375103*^9},
 CellLabel->"In[2]:=",ExpressionUUID->"a84267fe-8547-4235-8779-ea70f56be419"],

Cell[BoxData[
 RowBox[{"Off", "[", 
  RowBox[{"NIntegrate", "::", "ploss"}], "]"}]], "Input",
 CellChangeTimes->{3.451259087379103*^9},
 CellLabel->"In[3]:=",ExpressionUUID->"26f78377-a4c3-42d7-8ce1-ac8760bcab6b"],

Cell[BoxData[
 RowBox[{"Off", "[", 
  RowBox[{"General", "::", "Stop"}], "]"}]], "Input",
 CellChangeTimes->{3.451259087382103*^9},
 CellLabel->"In[4]:=",ExpressionUUID->"8e7fb54c-ee22-4c34-9bbf-77f7e8b4ef02"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Alpha]", "=", "0.6"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Theta]s", "=", "0.01"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", "=", "3"}], ";"}]}], "Input",
 CellChangeTimes->{{3.7506314814737186`*^9, 3.750631493625293*^9}, {
  3.750688642154299*^9, 
  3.750688642515098*^9}},ExpressionUUID->"3b14f5c5-ea64-4841-a7ab-\
3bf3afa6d1d5"],

Cell[CellGroupData[{

Cell["build BC\[CloseCurlyQuote]s into basis functions", "Section",
 CellChangeTimes->{{3.7300592267233686`*^9, 3.730059254107985*^9}, {
  3.730478686574498*^9, 3.730478690740705*^9}, {3.7450002228898153`*^9, 
  3.745000226683486*^9}},ExpressionUUID->"fdbf016a-a481-49e4-a4a1-\
7fef444f2077"],

Cell[CellGroupData[{

Cell["null space", "Subsubsection",
 CellChangeTimes->{{3.737730674156197*^9, 3.737730688827662*^9}, {
  3.743520912861556*^9, 3.743521014943858*^9}, {3.74352104660297*^9, 
  3.7435210770410748`*^9}, {3.74352284608762*^9, 3.743522849789524*^9}, {
  3.750624046427766*^9, 
  3.7506240472131853`*^9}},ExpressionUUID->"d797f85f-1f44-4aa6-bb33-\
319d3a419081"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[CapitalPsi]1", "[", 
   RowBox[{"\[Theta]s_", ",", "n_"}], "]"}], " ", ":=", 
  RowBox[{
   RowBox[{"\[CapitalPsi]1", "[", 
    RowBox[{"\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"j", " ", "\[Theta]s"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Sin", "[", 
        RowBox[{"j", " ", "\[Theta]s"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalPsi]2", "[", 
   RowBox[{"\[Theta]s_", ",", "n_"}], "]"}], " ", ":=", 
  RowBox[{
   RowBox[{"\[CapitalPsi]2", "[", 
    RowBox[{"\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{"j", " ", "0"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Sin", "[", 
        RowBox[{"j", " ", "0"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
    "]"}]}]}]}], "Input",
 CellChangeTimes->CompressedData["
1:eJwdzl1Ik2EYxvEXlaIFkbIk7ANcbZbsIEIJJaMnmgVKM3FQSoiRlEtMPNDV
zJXEEnJ9iisnC7WDRR5UhntFbGSga5JiITmlMYvIeizcYh/CHD3XfXDxO7rg
n33+SnltiiRJGjF4a7Tok/IdZzsn/fNQ7Vn5Cu/tCC/BibLyCNTHC6NwWXW2
MFOoTQZJt6upKktYXOYgFbLNXeflLJhiGIOJ44MxKE8dikPXnC5d6+Ps0ooh
Aza2ntgO6/NWSXfOsWRoirO7NxrIZ/bIubAwTe6ogVGL9Q5UXU2QPboxY0KY
XtNqhjMLMSvUVnZ3wgbP0v51YapxWy4cChgKYOk3GxnUKe2wTdFIvrTadyXx
b/aTkTO+PudHzkZ0+f0wWHkzDpevjZJNpvWS1RnOTq5l62F8+s8F2PHrQT20
/d1dERJKW/eSpk2WAHz/4zZ5XflUFRa+1SyQrjdHLIo5zuxb2tth1DxghWs/
8zthZFA/AbvUT8gNaeOL0KjxksPOgtTNwosWRnrbhg/DwAEP6ZCyHkFt0UEy
I2Z+DuX73eS+LyWzcPrVafL1P1vxi3nOHvsc5AfnojvHLzq/V8nw4Z7xU54Q
Z0dn1Zfh742fW2BvtWyC/wE9F11l
  "],
 CellLabel->"In[50]:=",ExpressionUUID->"6ad48be1-82f4-4c5e-848e-3955567c1543"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CapitalPsi]", "[", 
   RowBox[{"\[Theta]s_", ",", "n_"}], "]"}], " ", ":=", 
  RowBox[{
   RowBox[{"\[CapitalPsi]", "[", 
    RowBox[{"\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[CapitalPsi]1", "[", 
      RowBox[{"\[Theta]s", ",", "n"}], "]"}], ",", 
     RowBox[{"\[CapitalPsi]2", "[", 
      RowBox[{"\[Theta]s", ",", "n"}], "]"}]}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.7355739710418577`*^9, 3.735573992623913*^9}, {
  3.735574101539873*^9, 3.735574110372566*^9}, {3.737730504214737*^9, 
  3.737730512867444*^9}, {3.750624101943598*^9, 3.7506241337449923`*^9}, {
  3.750624199673057*^9, 3.7506242089691877`*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"5c3f117f-edd1-4683-ace3-bed538a9f32b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"null", "[", 
   RowBox[{"\[Theta]s_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"null", "[", 
    RowBox[{"\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{"NullSpace", "[", 
    RowBox[{"\[CapitalPsi]", "[", 
     RowBox[{"\[Theta]s", ",", "n"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{
  3.446467244532*^9, {3.730478840077458*^9, 3.730478847237978*^9}, {
   3.730483273972766*^9, 3.730483286252482*^9}, 3.7306507223573847`*^9, {
   3.7355741816543007`*^9, 3.7355741901825943`*^9}, {3.737730515384025*^9, 
   3.737730518682877*^9}, {3.750624144770982*^9, 3.750624159393696*^9}, {
   3.750624211273332*^9, 3.750624219713533*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"698cc67f-7248-4dd9-a3ac-3bd95fdd14dd"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Recombine basis functions", "Subsection",
 CellChangeTimes->{{3.737730674156197*^9, 
  3.737730688827662*^9}},ExpressionUUID->"1c31843d-c3a2-404d-b4d8-\
27c8c57f9dbe"],

Cell[BoxData[
 RowBox[{
  RowBox[{"y", "[", 
   RowBox[{"\[Theta]_", ",", "\[Theta]s_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"y", "[", 
    RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{
    RowBox[{"UnitStep", "[", 
     RowBox[{"\[Theta]", "-", "\[Theta]s"}], "]"}], 
    RowBox[{
     RowBox[{"null", "[", 
      RowBox[{"\[Theta]s", ",", "n"}], "]"}], ".", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"Cos", "[", 
          RowBox[{"j", " ", "\[Theta]"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "0", ",", "n"}], "}"}]}], "]"}], ",", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"Sin", "[", 
          RowBox[{"j", " ", "\[Theta]"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
      "]"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.735576507017578*^9, 3.735576541692112*^9}, {
   3.7355769686003017`*^9, 3.735576974295094*^9}, {3.7355770646772203`*^9, 
   3.7355770909145517`*^9}, {3.737480077384235*^9, 3.737480092071457*^9}, {
   3.737480331145009*^9, 3.737480334054442*^9}, 3.737730524521524*^9, {
   3.750624290326643*^9, 3.750624377502618*^9}, {3.750624994936388*^9, 
   3.7506249962549877`*^9}, {3.750625042949799*^9, 3.750625045441626*^9}, {
   3.7506252292416763`*^9, 3.750625256190674*^9}, {3.750626516176787*^9, 
   3.750626518805224*^9}, {3.750631458077841*^9, 3.7506314588680964`*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"9af92159-c78a-49e0-9fe6-46bac46cfd52"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"y", "[", 
     RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
    RowBox[{"[", "1", "]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "0", ",", 
     RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7506261175806637`*^9, 3.750626137056567*^9}},
 CellLabel->"In[55]:=",ExpressionUUID->"b2b8f0ba-19f7-451b-b670-3aed0e09fedd"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJxTTMoPSmViYGBQBmIQzbyf69xj00Y7BihQfbnPWm5yvD2M/2uDfsrNSflw
/lORjB/bJtXD+bfcWg4tndQP51+Nn6Q+d9J8OD8oR/LWjEnr4fyJn2I2T5u0
H863PurF5VR8Ac4vqPJL3FxwA87P3ZHax6T8AM6/8sKTc3nSEzi/S/8bw7bI
F3B+ofaJqQeWvkG4Pzl73d28D3B+0t6wgAlpn+D85I9XDhzo+ALn9xhsz1/W
9g3OZzqw3byl6gecL6jlfe7+jF9wvsd8p8UfJv2B88Xfn3Gr3/4Pzp8xcX/i
y/0IvoTJpprgE0j8yumb1G8i+FLMSXIXfiH4cuLfvyrY/Yfz1ewVlx46guAv
eyh8QOcsknwL2+1pVxF89ZOvBbKfIfhTNKP3cbIwOMDDs+rKaTsOBF/kjO/N
Yh4E322yVH6CAIIPAE4xjLQ=
        "]], LineBox[CompressedData["
1:eJwVl3c81d8fx+17r5E7cC8NK/K1RTRw3k2klJSVUbJSSRKVnVlJdoloGEXK
jqQT2UmFRMoOtyGZN5Hf5/fXfTwfn/N5n9d5vcc9H1mn0wdceLi4uKa5ubj+
/zvU5GORzOKC3VYO16gZHkiCpi1qvJoLljcypYOTA1BKiYJysQIXyPyrq1Xn
j0FdfJl6Y/9xwZfFYRddn3gkYcnauUqdC8Tedjp3jiWhGxzBI5G6XDCUmF3F
+yEN3TSYTLY14gKG5ZXJc8NZ6FbjM27e41xQr8+fXRBahCYHCi/8OskFUk9T
Y5wti5HRQs7vXi8uKDR5/zVUpQTNqCYOlp7ngrJBg9afXaXILOFEjWs0F8j/
drbzX1+BeOxXhbbkckFXLvvb1N8X6MRUIFfCKBfob7BoU/NoQPqrty7LOHMD
/+nD3aYO7UhlgDoS6MMDwlwOpUnmg6jgAd29zJsX1ESWOHReNgros+9QieQD
4Y7d01l+k2jGyOC2yxU+2DtyysslchKdLFrllhnLB4svPymuSZ5Eh8N7F+g3
+OCPwsk/fsWTaLOyjdxCLh+EL+80ff1jEv05d/BMUzMfYLWPDoOOv5GviKmo
iwg/POc9on0ATaEJX+WeDBo/iNr16nH2TCHXAcr9bnF+EN/I3ptkO4UsS5p0
96zhB89Pdicen5tCG2x2OWhr8INm991k9/wpNJO1tYDbnB9GNRuoo2LT6Iz+
RtOMJH7I+3IyM6F/GnF5JgVfv8kPNyRLdMjfp1Fs5mRJSDo/yPzsczg7N43y
eB6uOpbFD1FuzYbKIjNoqElyYl0pPyRM9DSxNs+gg4f+xhV18ENuZSTXRPwM
0vV80VVHF4DBk1fT5DbNorpMKaFyCQGoOGCdv377LDrw3hflSgmA7d/sRt29
s8hTR+PBZTkBODKlNi91dBblLmReMFsvANWvJSbMo2cRKyp0Vbe5AJiovjP0
6pxFC5k7nL7HCYCp7plamWNzyHDeIEErWQAunKkdOnByDoWa6dX6pRLreTqs
zp2bQ+Sl/+T47glA+g4u4atRc0jcVnRoVYkACOaPrufkzSFNsU9HzT4IwF3t
ENe2iTl09kRHfFKPAPwMTu1enptD5bWtNZ++CEBZTFSXzPIcMjiDZd1GBcBG
nblLS3Qe7W7LGgyeFwDVhZCHV9TnkUv06aNFkiToDS71EPeYRw/63ePnV5Pg
VP6o4uYz8+iHrlONgRwJVgdAxd7z88j760HZFmUSqIg6amyNnEch2zYPDm0h
uHPHLuu78+jWIt9RMQcSWAQoUA0751Gfxb84WycSSLzP1jf9NI9k8+df3nEl
wUZPzbRdA8R+Nt9k1E6TIFDYvYn0Yx6Vlb8d2BlCgitnRPqKeTjondetI373
SMB+s0DNUuUg5TjqX58cEnBCF2kGWhwU9iQy2TuPBDfvitEbNnCQ7oR3y6li
EhRu+jWTa8hBt0+Y6jjXkqBbuDE2bx8HzV2paTvaQII6q9D3tRYctD9P77hj
Cwk2zHf0t1hxEB9bPsO2nQTp7a925TlykIfrIsl8iAQHzoqlZHtyUF2E9z2z
URLctzx6bYs3B63OHtff840EDUECsjXnCL3DnWeMpkjwIJeLlBbAQXpHC3oN
eMjA+UPfcPkyB8WFyPtuESDDr0aSb1oMB33LTKVuEiTDSy53pdvXOSijL2Kn
Dp0Mafu3vvFMJvTYOTxRliOD0iIjdvcdDnLw7zRRUiRD/KPsdx/vcdDTW7tH
FJTJUNqzyuFgNqG3R1dSdj0Zru14bbQyj9DLeVSyRpcMG5T10lwecdAalrzZ
qs3E85+/12U+5qD3lqJhzG1kcPEK+9xfTPjrG7FafBcZVlKlF0dLOSg8+e9T
+m4yaDsk7vhSzkEbO8d+iBwgQ9tWK7jxjNBr/tKSz4kMliMhF6CGgzheur+5
Xcmg2H7j3MtaDjoQ9+jq8nEyyLjzBWrUcRD/25svF86QodDowp1PDcR5JlYc
5pwjw0XupTKJJg6qEImYnb1AhhLzN81bmzno5J4zyr9DyUAZFh9wfs1BDSfG
6iYiyBATf/uzXSsHSV+1d/xxmQx2L+Jeb3/DQefzOv6wr5FhMVEhT7KNg9qb
TZLG4skg+SXrXB/Bqmys/jWZDHmCHioJbzkokqzbPJRKBq7Nw0067zhoYN2j
YwO3ydC+Od60geBNRnL/vtwlg0AyqWjXew7Sud8S3ZhN6FF/8buEYPVlb0bR
QzKUbderorZzkNLhlRm3CsjQUekgZE+w3NNXSuFFZDjdoTVwk+BVjJMlp8rI
kK+H9eoIljgtZmhVSYaQOUVqP8HU18+boJoMrtduubEJFlrnYqFcQ4Zv3+30
hwkWCBPpY9QT+TlfG9NKMFd/mftSE6GfoWWVTfDCZofp0VYytKySuXWK4JkU
gaB378ggdoF1UIHgianH5GedxPrE0rBWQv+4mVXi/W4iX0fvrDtG8FDe8upr
n8mQcDzVkE2c/4vAgwe+A2QIq/dsdCC422m/9pERMugpcNe8Ivxrf8GpNhkn
g2/HNmUWwa1Sd421f5BBrVKPx47wu8HXpGPVJBm6j8waXSfy87L9t73ADBk+
/Lm3UETk75n6rfFf82TIPXxBso7I75PRb0u1y2QI0P58s5zIf962xOhHvBRQ
r17PSSHqIytjCyOFRAG2/q4a10YOumkZo3RclAKvVvyN7STqK6FYp+QAgwKP
VWYSvF9xUMyKLwb6TArs6bFZXibqM7RBzUJUmgIbn36hsl9wUIBc1xeOHAV8
ng2W7ajmIN+gIPchRQpc56x+db2K6Cfdt4Fl6hQ45tYmMfOUg5wT/MiZ6ykw
/LTcXqT8//UrnRitSwGrSCwiQfSLRY7Xg8OGFOBv7subfcJBZjyS2ju3UUB/
aP5mawEHGTvUVKvvosCaBt2phHwO0pegd3CbUcD92K/xnzkcpBBZspRrT4GV
gWX/PqdzkMzQ4ej4owTfZttp3OIgKUM+hr8LBUbCc9TO3OCgFXMHlcxOUWD6
NFZ5Fc9B8y5zB2YCKCDt6VdzJ4KDmnduegDpFPho8L1+uwfRH41B6sF3KMCn
Klyo70rUq0ldaXUWBSbMz4OcE9EPe/bVbC6gQKPB4SMlNhy0w8Llk84LCiQy
3s/9NCLydyRe+L8Bwm+PxPuNMhzkNdSV4DZCAbrfroSYlUQ/Oq+SzBmnwAdV
sWOGEhwU5JarIP+bApk/9lW7CxHn9aw2XM0jCIeFunq6Z+bRU3+2F22tIKSa
5Gio1M0j1yX1+X1KgiALvDdfVs8jiWCfwFhVQXgx3aK7/ek88rm0fFlogyCg
EC5Nat48Wn9Z/D7/LkEQ7GPUNMfOo8cpWz/8cRME7leKajUH51F2UeqmoXxB
CMi4enmmfQ55/D2s3P1EENxoe/o7Woj/051rVraVEM9nkx7er51Dz7vvLVZW
CcLALblKseI51LGcj+NfC0IF7VGyQdwc4jZ7sXPrd+J9JwalyGQO2X8bMr+j
LARa+1IDrItmkZxO9rYUdSFgvvs3dCh3Fo0HumnHrBeC0KWA5h23Z5EP7YfY
+c1CsFnpuv48cV+4unH2477dQsBvKt/4xHEWVUaSHJaPC8FNUTXe84KzSFxe
1cPhoRDYHxOYe3RoBsUlXLV8WCAEL+tGdH4bzyBBnu/bZoqE4Bz8mlLQn0Fc
gw9XXnkmBEUnAhSPys2gHxnr3pS9FoI108H8MhPTqE5KXlNkQghG5JZKWi9N
Ix+61FzVemGIfHl2p9ndKdTOTb4k+VwYzOSNauPKJ9FEQarD/pfC4EO/L/Ai
ZxJRbFW2RNUJw5rqa44DKZMIis1mZlqF4acSTwofcb987JTs+u6LMBhPvyG/
0Z1EV1/JmUYtE/GFxxJ63X+hHREGYrPbRMDhmvOJvv0/URn5bM67ZhGi3vt/
qf8aR6riVA+n+hUQUCXxh7R3AP2KnRPT7hIFo4FXMCzxGslHRWmc9qAC5b+y
+06J+di7SyfI6SQV+B6KzG+JL8AvFYZaD3lS4agY+8GW2CfYoU7fQ9+bCmoc
6puO2iJ8i3sqi3KRCh8Fs0R4ysowPcBOKusyEU+3/2FJWDXmOaPF3/OACt+3
z0f4OTXgIdveT9vGqBBLlf1y8EQ7Fjawyj3CpkL46ub753Pb8QbpjrNB36mw
OMLM9Bxux9HDLcLPflFB6kSC6zvbDqxx8pmhFocK7+qpj12MOnFAYOo9aQoN
bMO46/OZXVj8jpXHojIN/jGyrQ8n9mDDSx26Umo0uL5pQcq9uQe7O+/j3ahB
A0/DbMOtyz34uZJR+lltGqRN3uOz9viEnYp0275toUEq5xOTYdiLH9eKr+/Z
QwOvkt/v0ns/Y+OvHX/KPGlQFpRy7u/PfhyoEuv/xIsGz1Sf+1NFBnDxGeOl
B940uLg/57Wg/ABe+e85V7ovDc5EOIyWmQ3gCYkc0qUgGtScOlg4kjWAk4zP
i5vF0sDu8KDJrT2DeCB/ldZoAQ2+P0mflbw8hMWnuor7n9BgKc+miZw+hHdv
jNfpKaLB2Y7n3V8eD+GSOv6NrWU0MNZajNvcOYQjv0wYFlfTgFZlHVC9ehir
itbsCXpDA/fCzsSh/GF83tvFXeInDX63ahn7l47gaNMJt82/aCBX+u7d2foR
fGOtn5vDbxqIhAbI7+0aweVd0a45szQgOz6k3ZsfwdObHznr/qNBiINcZOfG
r9iTd/qIpSgd0jrfC1mXfcWBn/2PXKTRgWNhI3a47iuOKeM7ksGgQ/4DmxzD
jq84303CcZRJhxN5H+1yJ7/i8deb7H1l6MBzXLP733+j+FhSiE2KFh0oLkdx
Xcoo9j5FsXmmTYc7DTyJb++P4ku7Eqz7NtChZ/7N52eFo/gO576V4mY6+ISp
dK5vGcV9do2HyrfRQfFepfODv6PYWmHFgS4LOqhlJWvnW49h938p5guH6CA9
Ovq9+tgY9vsobb7Gmg6T9SHaRZ5jOOWy1n5XOzpkKL0N2BQ+hjt+HjSbc6ZD
kZOI1fuCMWxWnrZbwpcOTUdjCmX/EvGmnZbFz9Ph9sKBxGzecRymqVwqfpEO
yy7BRiuEx3FFXsVq8SA6DC5n6MauGsdymR8mGZF02LiXq9JYfxzrf07PZkTT
Qbt91ZPp7ePYStLZlnGFDo7ruDPDTMdxTOLUK3osHWbsHlw6YDuO56JEb9BS
iP0LZpYe+Y5jan2XKe0mHYTT511CAsexCk8GF+0WHaJmOPMbw8fxkQBVD2oG
Hf7Q7Utt4sdxi5eJgWgOHViCima2eeM4wyZsWLicDp57SS4OHwm9KbtvClcQ
/je3vIr6PI47Omh7hZ/RIUW191XK4Dgmm90pF3pBh+KXzc0u38ex97bqy4IN
dFjnlTTQtEToDQ43FGwi1pNz9WV42Dj3uek0pYUOZb+bQ+wF2Piz7ic7Shsd
/mo1jl5ewca7VOY1yF10mP0brZu/mo2PuL8YIXUTfmkY9q+XZWP/7IhU0ic6
HJmeqr6zlo0LpcV4SX10kPc/vktDhY1b7HqfCgzQgf/bPN9+dTb+mnrvpMAQ
HVYIS8tbabGxlJhWF/8oHfrOcz1k6bGxjjnnCv84Hd63H1Ps2MTG+2Ix4v9G
h9jXReo++mwcTjZ7yDdBB1upq8IeW9k4Y6e4A98kHXQqb3JebGfjykuf6XxT
dHh8MDl9YScbTyx6BPDO0WGIFLVZcTcb21ib/vdnng5vnTTfCe9h47oSla6J
P3Q49bxwZ+9e4rIsKhw28pfQ57kzJ34fG6d6/ND4tESHxBU7ljTM2ZivofXz
22U6rLWVsiw/wMaesgWX67kZMCBp9nztQTbe0X1quJCfAeua7PorLdn4ibZZ
XA6JAZIFz0oGrIjzXlc3SKcwILbGv+6XNRtHfFvxLV6IAQkZ6avHbNh4cuev
lCgRBvw6pPW6wZaND999uz1QlAFUG5ee2MNs3LD4ZNKbxoCvjv4HkR0ba1rH
3XZnMODGvjfG3QSnlXjtdhBnAJJtrbWxZ2MBUfN5CyYDrJuW214R7OWhlWUi
yYC9kqSLEg5s3FtPM0crGdDw5MY7c4J3yU4t6awm9u8M7D9HcFFAe56yNAMy
XdJLLhG8qrvYSkaWAX4VP2zOExylncgnIc+A92FB/ZYET8WeLRJSYEDeKz/z
NQTbf7Nw4F7HgDOZYtVtxP5NO3WE5pUY8GxfroY7wevvilX8UGZAL+/N8m+E
/tuLM85DqgwIbHR2tCSYbP2B1q3OALL7Xb184vzeJWUv3mgyYMlmfg+b8OfL
ipQTr9YzYHuGQgGVYGMPX1alDgP25XK5yhJ+ltRb1j/WZYCRoXfESsLvNbJ6
3lkbGbDsdkFymchPdABT+tZmBkyYKGu2HmLj6Y/zr6/rM0BI9PHHS0Q+HbS7
z0cYMiC6ddsaeQtCf2yFgj8wINF1C73A/P/1dyHYdQcRf52QfpAZG1Pu2qja
7WLAf8q/YuuJ+vJZ3NRjbsyAniv7m+eI+jMpWVhvsIcBJj3hGkwjNp6RCRhl
HGRATOu4k7khGzsG2CVSLBlQscEhgb2F6J+P+rBsxYDjSyayJ4l+yIxdSv12
mPBXee6Wpg4b714M2ltzjAHr9X4+qVBi4zIrx4VyFwbMbOuJyVZgY5kSlPvI
jQEuXSK7AuXYeO44F/fNEwxQblp97edKNr77MbTU8ywDRAt8NDtE2Hi0XrM/
9hwDyjqTRxsobKxa2k954scAj2NOVvf42bg8zsDxlz8Dbh49VE0l5kmryR+K
VzgDHEyNJ958G8f0jQ904iIZ4BPwopV7lJiXilaOhdGE36wF6hpiPg3xlJVO
xjDAs5MiRyXmGafKy/FMMgMEu33bPteMY8M86SvxNxiQunx0m+VzYl7fbCst
SmWA1EV2a0X5OF5xTlVw6jbRP349wZvyx7G8+nipdw6h55lOt0riODa74yjo
85Tg7i2e4w7jOCl2xYakSgYYP9rH989qHPcEVDuWVjGAW9Jf7t/+cexss7Js
BhP9NShv+WLbOL5A/+h4rokBYnoiG7cqjOOscLMy3x4GyHnMC7qOjuEFd4Mj
F/4yoPuUvn+3zRjeXaUM55YYUJc2XChqPobTRFgyZ5YZkB9SPaxrPIb1i6f6
3XnFwNokK+mQ7hgO+pvraC0kBv211zy7aGOYN5buqLdSDDYUrvHn1I9i4ZIx
+9nNYqAU7X3SQn4Ur16KP3zmohjspn00Unw6go2H9lbeCBADwyLRZe68EezT
SGG+CBKDQy3Fe9+nj+CW+NB2wTAxeNDwWM700gj2VTxjknVVDGKOyzY83jOC
3+3br/cxXQy8vqq0mvUP40v3VjAMsBiIR/7JrlgcwqNGV1vIfOKwtOHk8WNr
BvErxVa05qo4eHYcs30e/Al7kHP4ZCkSoHq4ypvLrgO/5om+tl1IAvht/OOp
hh1Ydem4hKuIBLRptVSrynTgX7/VlPJpElD4Y77n8kg79uktM90gJQGu+IxF
76l2HPi4PsFERQKoWvYHFy+9x3EHv8p475WAT7R9tSe83uKnd9YavIonnrMk
Tb7kN+G1yzv8axIlIGiQQUn3b8Jx9i6VOFkCIuTaTYJMm7C7VM6G56kSIOgs
zvX4eyOWSlqnXnZXAlp0gmlYrREHRClL5xZJgGFk5ok3FfV4m6cm99X3EmCw
/ahn5fda/Lh1P7rcIQHRFFOfTS9r8UqVM4FRHyRgp3W4wFBSLZ4ZK/oT1iMB
GZcyHAoMa3H2Ue2pgEEJOD02XKqZVIPJh3SHPX9LwKRC6GuS0UvctkW//gCN
CbbmffAktQpfqTj3dxeDCbN25LT3ZlXYaMMTrS3iTGi5rXTWhK8Kv1SXy5CX
ZELQv3K3zV7PcIks2W9GhgkL0T1tinsq8U1Sp1KKJhOetR2wshB7io91nIjp
2ceELkNfc42xYiyzP7v2jTkT3BpzhW1zi/Hn1j5OjQUTSAIjOf1uxdiywdwl
z4oJvhafeRfZRdi4Us/Q35EJ56XixAKnCrF6Ju/kqtNMKMSeey/89xhzPNIs
jsQy4Y7fFPe0QS6+fvE+KTiOSXwv+r243ZqDFa/kV91OYIJMKm+8ol0OPvSw
Sr43hQlUD/3tHYHZuHisd/pQJhN28fGszXh9H590XpVkWsgEpffBQjLxdzCv
z1pjj2ImNCzoudao3cG3wlQXo0uZcGNxBzsyJRM339N3bqhgQs+CnkUnLQMr
DtrpbK1hApxlbb6yOg0P2N/u1G1nwlp8TvXa5WTsdyo7+lAnE76t/d7wry8J
rwgs0PfpYoLDoOxdHd0kbJBenVX0ifBDdcPFn+wEfOvTFx/VYSZoF9oJpLjG
4YNWa8TlZpigOBoweOt4FP7uqtgMc0yI2XnqutjeSHzJVz3QkcMEVlXZMFMn
AhclG35NX2TCSX51zSrhMLyi06GMyc+CHcdKQobeBeKm/ZmHRCRYULFGpbhd
4RT2qPv6gsRiAbVS+qFr+3EsslFViUeKBSGpJdRtV13xAenKhbnVLHinmPax
RfYI7v3ZnjmgwAIIjdOUpCIcdFSS0ruOBVxPatdvkViDZT44en/4j2CbSIr0
Lg3k/PzHjhY1Il6s5dHDbmZo4orAt5INLCAXvAyWPeCE4v7tsXisx4Lu4OnJ
okFnpO2d+PzBJkJfR0kdzcUNnbeRuX7bgOD3W4Z2GZ1EPEqbdaJ2skBYJNW0
77k3ykoLuR1qxILxN3/K/3WcRUaijQIBJixYvKJTtuKrD4qZs+jx2ssCYy+7
eu5ZXyRefyrI9hALSs3sRy1aL6KKjaVjB61YEPfW7yy1xB/ZPlrYv8+G0BMx
9bE+MQBlJEbJ77BnwaoHfJoy24MQkNpiDB1ZUPfMrayRFoyGL4rNbTzKgnSx
8g2nPgUjJae7TWouLAhSsq7tfRSCWj6MaSm5scBH+3aUBjkUnTJRT5M7zoL1
sxoCEUdCUbFm1SnmKRZ81r5J2ki+hPSXO+nc51gwkiSX3K4dhvq8Vwb89WXB
lVi7GGO/MBQ6evTr7HkW1C5WaTQ8DUMNbyaefgtggeJTalavRjhy37pB9msQ
C3RU6kXPuYUjwTL/K/0hLCjfy7SXTA9H+9LJ9p3hLEj6vvAjYjEcTYnua2iL
ZEFwyJsv+5QiUFJYskZzNAseD1cKK5lHIN353puvrrCgcGzAUtQvAnV7yPG8
iGHBWo+6NIFbEehin/uJilgWRO/dUb3iWQRafeBJZ3EcCxpOp2UqfYxAL+tn
DQoSWDB03+6/Q78jkNMm/dzcJBaYPTq9PZkciY7/fDP7OYWIN8l4KCEdibzu
Ou6gpxL52K7WM6sVic4f+p1glMYCN5sxWZ2tkSiYEjYYcJvwm/9PZ/ueSBRZ
LaZZnMmCeEqmUcfBSBR7Jido7C4L0gZT2nVtIlGywsY3q7JYIBYv2f6P4PSe
5pUHclhw3K8rW9UqEt2/dtgj6gELcqUnCmv2R6K8rT8rnuexoE3ipw3eFYmK
Z4NIU49YMFrbILhuUySqfEi1XPeEyIfBO82FdZHopf29LLsiFpirV+lvZkSi
JprOdHwJC67mxd0bX4hAb+vrtzaWscBG/FGPQH8E6rpgFbf4lAWyvOV2KTgC
fVFj92k9Y8GEdObDm+kRaGTwoprbcxasZKZtFfWNQD+ShQPSX7CAGbLy4V/T
CDRtktHy/iWxf0uLm+OaCLSwpCFJesUCNc0osv7PcMRTXOOmX88ivp/tJRMr
whFV6itfbjPh53rrjvId4YjZRkza14TeewceRPKHozWXyPdobSzgG6Ds/Fgb
htS+qaCAdqKe2Cf3cRP1p5NRfa2okwU5iencfV+J+jyw7/NoFwvufY7k25Z8
Ce1+5n3BvJcF7f+ac3LYocjck68p8gsL9FPL7uXFhiJruRSJ5/0suBTl90Bf
MxS5XaksURwh+qeKJ592PASF23JN/P1B9As7a+uRY4Ho6ooEfa1fLBD809b1
0DcAJdTKX3X9zYLfP7n8VcL80V1lI6X3syxILrey+Rp1Ab1YuOaU84+YD6sE
g1PhHKp/LF3UyyUJtbnBvpQ1PqjVqWiZyisJLjTDz20z3qi3pTPdnyQJwYca
ScdjvRDn1sqP+6mSYCJcqHp1xBUtmxUoRNIloVM6cu2NDc5IgBf5VIlJQoxb
pc3VC0eR2AknmqKkJHC/Kum5UWqFVkrPOB5eKQn39KP09k+YI7mOiMdxqyWh
LyAkV3OVMVKOYi7VS0tC7KbdL7n2b0RaWx6a/pWVhPZiE77Tg1aG/wP7zZ8D

        "]]},
      Annotation[#, 
       "Charting`Private`Tag$405382#1"]& ], {}}, {{}, {}, {}}}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->
   NCache[{{0, 2 Pi}, {-0.7927869327569858, 1.551492123319656}}, {{
     0, 6.283185307179586}, {-0.7927869327569858, 1.551492123319656}}],
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.750626133930895*^9, 3.750626137445798*^9}, 
   3.750626538832128*^9, 3.750631341445098*^9, 3.750631466348077*^9, 
   3.7506314992443857`*^9, 3.750685151168363*^9},
 CellLabel->"Out[55]=",ExpressionUUID->"17986e9e-d1ea-4bba-8002-8d5443cda6fe"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["orthogonalize y for \[Psi]", "Subsection",
 CellChangeTimes->{{3.737730674156197*^9, 3.737730688827662*^9}, {
  3.750624507579715*^9, 
  3.750624517649906*^9}},ExpressionUUID->"ea08ebac-e4fa-4fb3-a8a1-\
1284274fd018"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Psi]", "[", 
   RowBox[{"\[Theta]_", ",", "\[Theta]s_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "=", 
   RowBox[{"Orthogonalize", "[", 
    RowBox[{
     RowBox[{"y", "[", 
      RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], ",", 
     RowBox[{
      RowBox[{"NIntegrate", "[", 
       RowBox[{
        RowBox[{"#1", "*", "#2"}], ",", 
        RowBox[{"{", 
         RowBox[{"\[Theta]", ",", "\[Theta]s", ",", 
          RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}], "&"}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7286640954500437`*^9, 3.728664100741357*^9}, {
   3.728664198703786*^9, 3.728664212748205*^9}, {3.728665340925867*^9, 
   3.728665345580411*^9}, {3.728665377799589*^9, 3.728665443643546*^9}, {
   3.72866547537993*^9, 3.728665482858376*^9}, {3.728665741465156*^9, 
   3.72866574250451*^9}, {3.7287642189319687`*^9, 3.72876421929893*^9}, {
   3.7288336574145184`*^9, 3.728833660747314*^9}, {3.728833734451893*^9, 
   3.728833745132497*^9}, {3.728836481264612*^9, 3.728836483391611*^9}, {
   3.730482499125739*^9, 3.730482585729575*^9}, {3.730482710143909*^9, 
   3.730482710212982*^9}, {3.730482929331092*^9, 3.73048292940926*^9}, 
   3.7304831089964542`*^9, {3.7304856248871813`*^9, 3.730485642150317*^9}, {
   3.73048608666818*^9, 3.7304860883391*^9}, {3.73065064795049*^9, 
   3.730650655888479*^9}, 3.730651057303834*^9, 3.732895154358745*^9, {
   3.735410050793686*^9, 3.735410061376081*^9}, {3.750624551397264*^9, 
   3.7506245985584*^9}, {3.7506248095074778`*^9, 3.750624812808668*^9}, {
   3.750624856109438*^9, 3.750624860113104*^9}, {3.7506253888206043`*^9, 
   3.750625389402116*^9}, {3.75062542075815*^9, 3.75062544409221*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"06f9ecf8-085d-4573-ae5f-34752dd38caa"],

Cell[BoxData[
 RowBox[{
  RowBox[{"plot\[Psi]", "[", 
   RowBox[{"\[Theta]s_", ",", "j_", ",", "n_"}], "]"}], ":=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"\[Psi]", "[", 
       RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
      RowBox[{"[", "j", "]"}], "]"}], "/.", 
     RowBox[{"\[Theta]", "\[Rule]", "\[Theta]\[Theta]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"\[Theta]\[Theta]", ",", "0", ",", 
      RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4437384361685085`*^9, 3.4437384404555087`*^9}, 
   3.446467244604*^9, {3.750625526433076*^9, 3.7506256022422037`*^9}, {
   3.750626894562504*^9, 3.750626894959093*^9}, {3.7506273912215233`*^9, 
   3.7506273918905573`*^9}},
 CellLabel->"In[57]:=",ExpressionUUID->"46f7fef6-fb20-4fad-b080-d1185973ae5f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"plot\[Psi]", "[", 
    RowBox[{"0.01", ",", "j", ",", "n"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.714747109333274*^9, 3.7147471500383444`*^9}, {
  3.7150790509470406`*^9, 3.7150790510650406`*^9}, {3.7506255844740868`*^9, 
  3.750625592680398*^9}, {3.7506272490318737`*^9, 3.7506272518619967`*^9}, {
  3.750627395356153*^9, 3.7506274074506273`*^9}, {3.75068864967819*^9, 
  3.750688669084094*^9}, {3.750688739682876*^9, 3.75068873995862*^9}},
 CellLabel->"In[87]:=",ExpressionUUID->"85373ec1-21a2-4953-9ad6-303bee858a0f"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   GraphicsBox[{{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJw9mHk4lF/4/+1mbDFjGRJSlihKJIVzUlFISpZPliQhS1lCIVpkKcmWSFmy
prJmy3IQZatsY80+D4pkX8b2fb7X9ft954+Z63U9z3POud/v+77Pc2a3za2L
15kYGBiS8a///WVGHD/GVB5oMvy/j9TvquNi0VfA/2d6nqJtb9St/2OM32G1
OCrg/zhMuIE7p/zR/3GGrEUSQ04weDGFzS7yDldVq84eNEp+Cow2dxHiZhar
+rQDa9OjnoN6yyB7pUgmRJfjkvOrjQRXnj3YsXqBCwUrdbmnl0eDCZeLZu2/
yIj/WHL5j08vgP1/0awStiIo5YQjy+rHl6Btb6bP7SEJpHBW+dzuzHiQlOwa
Gre8F5Ubbr3QTU4A+zQiZPy996EzZg2DHvFvAKvSTkMZGwVEvRIl8yYqCWgE
WzPo9R9C+YcVjo7VJ4PZv79E7T4rI2/V6F/etSnA6Wfv1vUfqogFWEqllKeC
yeYEMcoLTdSkVduoUpoGrCVdlMb3nECR2jI3mz6lA+PRz9WfhE4iMYO54oWP
meChoNiBfAMdRLtgYh6SnQXS+x6QhyPPovfG5duime+AQUWq5+N/euio5eMz
2snvQdii8ojIPUN00Vm4Ly4qF/C9M4h+I22CuI4CFKuTB0j+6mcEC0zRbJkb
30h9HmikOV1q0f0PFVd2FXnW5gOPRtcHc6kW6GR98lZieSHYm2HWEk2/iqR1
Ogx/q30Cs8eaj55dskEcjayph0s/gXUBonXi8jXU2uKo0/CpCIwbHLKT5bZD
lp3KEXMfS8Bb+yUGFOqIoIn96PH9peBGRt6CZ5kT2tsTrxyUXQrU9/37+P2f
M/rTv9UjklkGdM5bbl5xvYW8xxp2n0ouB/OyIj9Drd1R5LxFYWwUAi2jBj43
DO+gZwKv1Ku5q0ER573xIoW7qD3/sbjKxWowkbAjYZLXB5n/scDE+qvBXOYD
GcYRX3TTnNNtfroG2FZfrmfIC0AxGg4h8bx1AGCH7NNYAlFvr5ETj3EdkN+O
nWSIDkRiXsDgUXwdUHk3s+PCnscoM0eQ30WyHpzZzefYcjYIfRavT4LKXwFQ
qp5TzAlBQ4ySJROmDeB4afGrnM5naG8i9yuL1w1gwC1yVudhOHI4tubXNtwA
6vR3H6Yeeo7m3Vu1ym80At3ih7XX9SIQC83/R7hvE9Dx+/KZXBGJ5L72YypJ
LUDqCg+zQUsMOl6vy6Hl0Qq0OzRGD2QlIN7OjJ+q0a1gIMK33ELwNaKNMr04
UNgKim/7qF1//BqFM5SLCy+0AuWmZU5muzdo9Li8yqx7Gxj6xWLidigJhRRw
Wie6twNy7se1kDspqCuppYju1glUtIbKnoemo/c5sr6zkZ1AqbKYUe1HOrpf
GQjH8zvB+clftw3IGWhfv3pz21wn6FMbqqhMzEB+gh+Hs9yo4M9jq4kPnzPR
nmfPOE3dusBwg8GrerZs5OpjcLXQtQec91xri/mXg3RvqmcRInqAa6RMKgnm
oj02cv8sc3vA8eGtHSsRuajrLNs9wkwPMOq+X7MmnIc0hKviLJ16gUtxcKwF
ykMcpQda2e36wO7Vxd5P5AKUvsQFLMx/gRcHJG3cRj8hl9Lr4Ux7hoF6LGf+
R4FyvL4vSDAqDwOKZ++laYNyNNurXrB9ahgYPf8xfi2kHN1eJ3dt2A0D3x0L
9zbXy5GPZu2ulexh4IitsPdgFSioTuzjlNII8B833iNYV4WSW7ubO7VGAVOi
tyuTcQ3qnDxLzLShAQcbyvqLJ/XoieIyQ/F/k6BEzuxCre0P5Cbf8KI6fRpc
MxqW1w3vRPRrTjkDN2eBxxAvs4RIH7KpNDGMsJsHizFRoZ2Hh5BSMYUkfGse
XDv5Gzt1dggx5fZ1pHjPg+YiFa6PVkPobbKVaWHIPKB+2KVuGzqERh/ZWXVl
z4PXhru/Ow4OIRtdL2fRf/NA501ov9XJYWTT/SIky3sB3PjElNA6MYyuzXVW
V4csgjGuKqmaXaPI2ifCtDRiEYRk2VmbS48iC2b9mdy4RXAtNi9xRmEUmQh8
2ZmUtQj+oyzPb4NRdEYt3+tewyLIMDoS03J1FCk8eLb/GGEJHEvyDGdKG0Wr
vNpxhcFLwOxMp0+I5BhajGdUfP98CVzIz3ggJzeG5iQr69++XAKC4UnLXw6N
oT/KyguRmUugIVB87RccQwNmkgau35YA9+YSc4bVGKpN3mI+wL4MDjswWITG
jaGwgyW3MoKWQWPsPUfASkNNSt3jnuHLwCkiczyWk4YIKiuWp2OXQc6Hrlwa
Hw0Fqqnq09KXAWGzYsNGjIZ8tUpkJeuXgd+/i7q/j9CQ46WSkTfMK+D8td1C
ZDsayjLpNnPhXAEu2SxqTU40NG628lOdvALKT31Z9XKjoauWqlX9kiug2pou
ne9HQ2Z2Ja+EtVaAMyeRPzOKhnTulBi9uL8C0lryV55X4PP7dDfZhqyAduVG
xtc1NFTrt3JCOWIFfNMfZHzzlYY0H6ge7EhaARczZcrvtdLQkSclXHxoBQiP
yyyXjtHQ3tcl9WGbK2BKI6eJm4Ahm8RudQvWVRAy7X7PlgtDyckrhfLcq6DC
bzAolxdDoumqb5tEV8HvOlFvSWEMkXNK/Anqq4Av+EGNzz4MMVWXqAb6rAJ9
k+CaaR0MZXOTMjYfrgLv6R1c7XoYumjuzO/9dBUsaYlNZZ3HUOqyxNyN16vg
3vrEUXlTDJ068CT7fNUqqLt6nnr/OoamfWiUxq+rQMG+4vGcA4ZiGjSDtX6u
As+Ju65GzhjCri1cUxleBQRspeGPO4aCXlns2sm0BvJzCSZ2ARhSmCx+GsOx
Bu7/jsq68hBDXSp8dG7yGmh3PzZ09jGGZNrquxj2rgHZ2yw/Z59gqIFdMWL8
9BoY5FF9GfwCQ67GoVtXDNZAh+jvSbY4DFFSx5x7TdaAmnlVi9crDN3QjDvb
Yr8Ghrd375FNwhCHJyNzQegaeJ09sfozE0MFteZu8lFrwC7qQmP7Owxd5i0e
Snu1BjQtjQbr32Po/XvHipfv18CyblGeXy6GLq3VyfF9WgNvHjw3h/kY2tAW
j39SsQbOLS9JLRZgSH+0w/Pe9zXwMVx+Sb4YQ4uKCrRl6hr458A1mVuCodf3
Qi66Dq4Bpt1vuvaUYegvRUPx2r81oF5b/6i/HEOxdi/f/FpZA0cVOQ6JV2JI
89McpwkDHUQXmpZcrMJQ+PmMyTN8dBB0HxwIrcbQkTcMprXCdCBzdC9XWA2G
Bv9crj8uSQd6gplv/WtxvY8WHS6So4M016RJ6y8YUgza8VbhMB187lylKtVh
qKfjBm/WcToQlSi0WcT5/u46/92n6OB63GhQej2GZG+J/X2lTweHmfOOaX/F
UGvFHXN+YzqYqzTz6sL5DkdH4zNLOkihKiiZfMOQhNmBo+x2dGAnFm9Xj3ND
enDG/Zt0QNAjc+1twP1aGOGne9GBOYvCTnecKSfUH3kE0IE/w/OwXJyrw2Pn
poPpgEFkzuwXzg6/Zq/YRdDBtN3eu6s488np/RiKo4PC8zzjzI0YGmdP4mtI
poPQ8xfvbuPXP2Pzl/Ky6IB3f5TXNM7Pv2jHxeXRgafE8Z4mnG1TXvXfL6UD
akdT3Cuc1QJmxG5U08G2w88Cc5x5LLVsLjTQwcHwCQkenMeOxaarteLzhdUO
5OPxlFL+TO7uoQNnMwmaNs7PljX2cwzTgVlSm8J3XA+bzshb8xN0MNpTUnoK
5yMFWEHfPzo4Zlzv9xHXkzNCbbl2hQ6e6296suM87PJM7f02Heye8Ui6iOv/
SW/EL5p9HZjdMlgMx/2yYg9lthVaB9VWwtz9uL+HsV/a+uLr4N2vd9/+4P4T
vhx8oiyzDmS9B19NITyf/Xt4WVXXgZO80t1qPH+CLfZf+qu5Djo6fX1jKjBk
cez+S6r2Oujpqgr8D8831mUZsQyTdWARqp1QXIqhvg7fq+FW6+ClgFayAZ6v
Ofk/07zs1sHf72kJ3UUYMnPxltfxWgdcTq3XK/B8/0D7enTyxTqg6W3+asXr
5UGtiF/rm3XQoLmwbxavJ5Pkm6g0fR1QSh6dZ8jCEIOFoHZo0ToIlcAEF1Ix
ZNRx3UiOug5ia9BxyQQ83/LLYkkD62AnJjH/Da/njXDuPjptHVwT7Le3isVQ
hm6RdfPiOjjZQAmyjMQQvYblpjP/BojrZDyZFoShn0lm+ZdEN0CwxHzWz0cY
Srv3YVF97wZ4HZRYNnUfQwZql3y5lTfAlPLJ9iUfDKXkpYbkGG2A0vEnGRdu
YuhMklbqbNQGaO++dUb6Iq5fQ+vAw4QNEL/NAYgGGPo6Z0URSNsA7q7NfH1n
MXT6lO8ztaIN8ADFaZ8+gSGt34XeD7s3AM+BKz9PKGJI/bCUPr/oJuh+L/+t
nR1DvhaFQRl7N4FIo79YNBOer49P1Bw9sAl+fmk7e2KThtR6LI9Yam6C8RIV
Mft5fD+4FyuRYb0JdmScae3pp6GD39iXVNM3gZG7wxnj9zR0a/aFYtPHTRD7
T8pgIY2GcoT3OloUb4Ji7YTjAYk0dMAZDt3/tgkEln26rCNpSI7vbmPj700w
ND91Xs8L31/M/7wxV9gCbUVBP5A6DQn9a9EOKNkCDcERTbBiDMVFoqu/0RZY
K01+31Q4hijKBX5GDVvgxvZ0zMn3ON99WSDTuwWQjkos06sxJMJsI9ZK3wLc
XIYLRV5jSExoZUlCcxuM1LXdXNk/hqTB7vTaum0Q8Y6UJ/BsFGWMkKv3f98G
KsrXNZ49HEXSgWz9sdRtIK9mVrzoPYpkGqd4nca3AYvwq/Bom1EUs8+8isjC
AH3Ca7r4VEdRbKGUXIEUA5TRGuURHhhBr759ZmS+wQBLU8p643eNoNnhvLv/
nBngAW3nqDW+EaRDz5jrd2WAercZ3fTZRtDi/uiRT3cY4F1D+K15ZhgZRDnV
2IUwQInLerOH0DBishR90JTJAM/dm1BJtxxGTvP3GKLGGeD3zhiOCckhpL7r
xLaELSPkOcPb/TKnH52x6HS3cWCEBqS2FNnn/ehSgv14qjMjXHQLNH57qx85
CYd/l/ZkhA6MRRwXDvajVwL9CfuDGCGrduXnpLw+tMrtefRoFiP0kI1Q3fOx
F33aynI7P80IgwKdijpiu5H8MC/t3m0m+K/x8dVJlQ70MYvkUOTODNveuFyt
lKpHfoOWHfJBLFDrebpCqkcUWtTReHP9CQsU9lUd6FSKQM75ovZJ4SxQjzni
9WPsKTIP7KeTXrJAEBagUv37ITom958kPZMFProWXSpwHqI1z0tuDY0sUF1t
mk3sWCTw4tbbcZ2bFTJPPNiju/UOuKkf1UuMYYWN1T7zg/6VgOFmTMDzOFZ4
kBLHU7dYCcKTZgvvv2aFO7UusJk4VoFspnei19JYYSu/uZr/JQRGG4RnZD6x
QhbVb832UdXgkvF6RH4HK6zvbrDyvVwLjtys6qojscH0gnvDJd/qQV2SCGex
IBv0f1r7enqjHlxs8wKZImyQNxv+mTn0FdxUVswKlWSDom8r5c0TvoJMetJd
AyU2eDMvZ63J6RugBD8Q7bnABm+JHiH9YGoE9KRTNlMRbPDrgZjvKewtQHNF
I+rQCzZYcv1U9axaC3hgoFrrHc8GQw0zdqw7teDvkfskWd6yQarok5fwZwsQ
uLxjVLSQDVYfn5MOjvoODvL3XTWgssEIJnHHWs6fwMOpIzKmlw0K8Np+SlP7
CYprW2r6Btig0JlUqob9T6Dhhnbbj7PBuzlTodE1P4Huj7SRgBU2+FwjKd1G
rRVcD7l1NV+YHWJh10iBY60ga8ghcmUXOxScFTghutoKpo/Y1GhIskMmw+7v
z7jagDt2aXeTHDu80tL1eEO5DdzXOjYyepwdTntEibAGtoFXGyxX+a3YoTRD
+nSWSDtodX1l7f2WHfaqHhG+s68DyEXwrt/OYIdme4orNlU7wKPcoBfu2ezQ
OyGs1Eq7AxyZcW9yKWCHA8J7XApsOsAbJz1l21p2qLbQjRLjO4Cj3Qb7hVF2
6KEQftmXoRPUPXZ/azCO30+Fmx+5OsGu9El1/T/s8MVOtsHPlE7QOtbppjPP
DjlV6cLXDnYC1asf+zWYCFDj2MC3KctOwGJhlSsnSYBpH78wt+HnUCvfzrOy
0gRI/V7rmVzeCUpe6dKk5Agwospf42x9J3DsPSK8W4kA+3V4E6R7OkGbyY5H
QloE+HJcSHnHRieQ83q8S0CbAN2k+VczWagg8MV6CUmXAOWC2vOEuangaOfE
NPdFAvzjPtcRvIsKEi9Um7DYEKB33dvCZ+pUsOp6ZI7RjgD/ySj+mz5JBRcj
PjzdvkGAgcRHzXv1qID1Z1w13Y0Ahzztodh/VOCs7yY394AAJ/znd6p6UMFX
p4m6mccEaMy60/LyXSoQf2p5ZToUj+/PoVXjACpobzwbMxFJgH1h2aJdT6hA
TUdyayCFAA+RRs6tJVGBcmpTyLd0XB9+8sWvaVSgsO1Ozn9HgD70HGaXd1Qg
WfJFNjCfAJ93FZjqFFCBKNm50KWIAJtPDsz4FFOB4C1+TdMyAtz3Ylk05DMV
cMrg+28NAT74F79PsZYK2B5xD5LrCfDmC7ethnoqYBgqcthsIMDdDh0m6o1U
sBjL5t/aSoB5R+bu5/6kgpn5HMLnTgL8ORykm9tOBZMGptGpPQQYK3E5OIRK
BQNsWVlewwR4XzNc/1sfFfTYGB62phFgSYrxj/0DeLxVq5VnJwmwMEem2GmI
ClpEUs4cniZA5c7DjIEjuF5eZztEZwlQVbn/g+cYFVS3z1myLeLzZeYVQ4wK
Piu8mvy3QoBdtXdEsXEq+PREy6N3nQDlu/f9spmkgtzxP5u12wT4Y2Vqoew3
FWRrRYd8YCbC0HWKzZ8/VJCWeJwcy06E1wLUpRancL/pY28COInw8yRZvXua
CuJMwmRv7CDCnfbh6dF/qSCqQLnwIpkIvyy1WO+boYIwngENdSEipFpIusTh
HOT4uEFqJxHqKpGbhnF+8PWA0Q5xIgzqM/Zm+kcFfpJdA6uS+Pgh7bc38ete
/v4Oo9JE+F0ztbYNZ9c+6YVmOSJU/E6+GoCz45Gf94oUiFCuPcWQDWfbKG9C
khIR8pZtRNrh67GaEY8OOUKELGWs0sn4es10G3a5HyNCn46f3IV4PEYZrlnm
mkTo3p6km4THa8AkfPi0FhGaaU8P2OJ6nLGqqVTQJsKPjywbmXG9tD7fOEPR
JUKd5O88frie6oKkDkYDIrRwPZP3nUYFR9w/W05dIEKlR5x5q6NUcPCHzWSn
MRE+EQ/i3RymAqmgws1MSyL0fGa6FfmLCiRGzUMir+LPWz6OlMT9F9FkIfte
J8Jd2ZTIsG4q4Fm+JGvgQoRv68DYdBsVEC5uFqi6EeF4rKrG2A8qYMpJ19jt
SYRXxyxIhc1UsHJ9+eKiHxHyEIXBUB0VzNckDgzcJ0L6oeSk4zVUML1Lx+Fb
IBEG2Kk8c6ukgmHqy3uvwojQ/u6TPY5FVNB4Wi0Lvsbnu5HbEZ5KBXe++SsE
JBPhnnsWxucSqUD2bN2nyjQiLM+/3TIYh/upf77m2EcizHJ1r/B5RgWnjK73
KVfhep2+9OblbbweOrOtPWqJcOGyQu2MCxWkmsyO538lQkzHVVrYHl//Zd8F
hZ9EmKB2sHvCDM9f60iufcNE6Grhz3PzGO73aFeUPY0IYUbCVWslvB/Yigpn
TBJh1MdtA2k5KvC3z5TaM0eEE72s8geEcT9uVmruYuKAzMXBprJLeP/z/e3K
t5cDyqnJ9wyndQK7TYWV87IccIqt4gT5VScQDLh9L3w/ByzIPHtA5HknuP1w
O5RThQPWTWp/T7/TCZRCBVJZtTmg+KgbufhsJ8iJPUFds+eAvudqUprGOoAV
Jdj8qDMH1JDNCDnX1QG4X7WMeLlywPffXUUyGjqA8xvTmYU7HFDifnZkzwd8
/0hzYZ8J5YDKM6Tdpz06QHp+vNroew6ocNxij+FqO3BcN5fryeWAXLbD3vq/
28HB02I7fxRywMrGv9JSfe2gouftRlk5B3xZ95DmVN4OOrbfo8hmDjgS/U7w
gl87YDSoOn1iigP27tr7KWylDVj+Gb2QLMcJ9yuu9Nt2tgJJ5XStWAVOqHsx
/b5gTSuYvGd/OEyJE57PijfK/NgKbvNN8985xgmNbLj+2AW1gqdHl7rP63LC
ug7xbOsjraAsiN1q+wYnvEo4SHwv8RMI7NnvaPWOE8b/S7nKcbkFREQ9NXn3
EeeWdVgu0wI4mKa0FvM54fOEayoyS82AYeTdziefOeE3+97IQ8+bwXSizPei
Zk6YOXo6kVLdBOpE9hzknuGEyrwiNzwFGsFtkshyuRIX1F/5t5wZUQ/aGQkP
hSu4YOsa15nlknIw8zHeyrCaCzaxVm25HS0HxMvyx4PruKD7lzmR6NLPABYY
LC62cMECpgibhpIykGPzwq51gAsu8u+xJBeUgKdfJPWCt7ngrMt64/fYQnDq
sQb/khY3vMk/dXpTLwMUETwyWhu54b+pG2qmv+PRfgFeR5t6Hpis1U21c65D
/8KX+Q937YD32j6VNXt0oD3BwYq3HHlhWomB183L/ci9S9nfxpkX3l2rPgq9
+1G11GiL8U1eaJQt/W4ruh9Z1ak7qrvzQidx870nvuPvu4zzaUQfXihUZhvQ
q/kLkfwsRNJCeaF68EVOKDaAmNwOsfZm8cLLzcJuB74PogvVg0Yt2bxwrPK+
w9WJQZS8I+wt+sALVWY32vwZh5BmzjjIyOOF751qX3uqDCGfqQQfj1JeeCVc
IdDw9RCav842y93AC7WvMHZ0KA2j0cv9fVoTvJCocTdo98Iw4tIwzbT+zQuj
sQP2ZdvDSEW8w8N/ihcuKH1z0uIaQSFjTVyf//HCGXaZPBmpEaTo/Fnz0Cov
VIO/ujVMRpDfvfi34kQ+mNq8Nyy8aAQJJJs6bsjxQdXCFH0vh1Gk+bDjiMgB
PujslS25y2MUOdieZz6qyAdPvnvUXHhvFFXI6rz2OMwHtV/Eu+dFjSKb/CM/
/hzng8J+M9Lj5aMop1ZAqVefD96dXWVT5xxDZ7COtaKbfLBuzkjd6s0Yuicf
7pvrygc5BDNZrdLHUIHbmc0sdz74faNyweDjGNq5VcHw2osPDr3+tZ+lcgzN
CGawP/Tng6Wn5t+U94+hmDN3BAzC+SDvnBPFXJCGGsOVXupE8EEPTy9Dp100
tNU5TTkRxQeP3dRxc95LQ/ZXr4oqx/JB9hjy1RNKNHTMR3ePSCIfXJopy1I5
R0PD70UPjX/kg/uTeF+evk9DAvNdBUO5fJB27brwtSAa0j0aqdybj6/PUsnS
LYyGCutYj7YU8cEaS0VJ6zgaChqY0Syo5IPxNcyeQbk0tH9Hjb7/dz74wEV2
51IfDV019v3h/ZMPjqg7c4kN01BsgoqhWxsfVPY8+UEVo6FtmWwjWyofjDL+
tKT+j4baYfRl3QE+uO+2asZz/Lx9x/26g+BfPsjYWuBhLI2hEL0Z+2P/+ODj
evLKqByGXu71trea44NmMt6+Vvj5vbgrxC5jiQ9KB4s+FT+KoYVjH2yPbPHB
65mp441nMMTEr2J7mYEE6XOiUd/0MUT6W3nNn4kEt09vnyk0xNDBpJ82X9lI
0JFh//AlMwzdZF6wNtlBgqeFpaar7TB075evtQ8fCfLv3y+n6IihsCIW60Qy
CVbOpAQ+dcHQe3vBK+NCJMjWrB/PcxtDk81qll4SJFjHeytC8z6GVtNqLV5J
kmDkdRYGsUcYYvfXs6jaS4KxFs9Spx5jSPqgpTnbPhL8K7wn5sJTDF2Luf9f
7CESvPzh8f0zLzDk7kL87/NhEnx6/MBi/EsMPdSOMhtUIcEvh7PxAzSGkldT
TaWPkeDS2Kzu3kQM5bXtN9VVJ0GFJitxxWQMoewik5uaJDht8MdE9i2GBi2+
GRdrkeCaZT7jeDqG/qoYGvedIkFymrF/diaGNnh6L21pkyB3z6Sr5TsM7az5
Y6StR4JSg6XN4R8wJPfKw8jxHAmKeElocudg6JjHxsXw8yR4xKrlmE8uhsyk
eC52GZFgkJ/BpHgBhhy2Yi/QjUmQeZonzawQQ97d4hfEzEgweWFiO+AThoLy
sgy1LpOgVV8oY2wRhmJDDxnaWZDglmNFzqtiDKXbfD7/xIoE9w9acYeXYOjT
8ZPnc6xJMLHOfo9rKYbq+FsM2m1IcGqWtgzLMNTx95LBsi0Jel9df8bwGUOj
XwfOidiT4Pq+iZkcnOeT7M5p3iBBCct1cf1yDDHe/adv40SCHUspMj048168
ox/kQoLR1eZcRhUYkpBn1M++RYKBp6o7K3BWZHmi98ONBG98uBcs8L///w6Q
9OY9SPBZrpm8Jc4GxQm6gl4keNNn/7donB0WbLYF7pBgnnWnYwnOjw7KfRLw
wePpSoAtOCe6zDkI+OF+P9I3aMe5NLt0l4A/CeoOVsc34tw+EdDOf58E3wnT
ZAtw/rtXJ5j/IQmWkQsYw3Am2PCo8weSoNJpVTlTnCWTqLPkIFzvvIgsMs7q
v16nk0NIUPDH1MMafP2mwraXyU9IcDwgCFnj7G4iv4Mchvu5nmM7h8cfFj3/
hRSO52PplzseOGe0lt0hRZBguu1hpglcvxruBwdIUSTY2ejHro/zL90zo3wx
JJhdq/D8La7/cvCOl3yxJMgpPZn0B/eHt75Ljy+OBAuKTOFenOWZEhn4XpGg
Cpe693ncT2u//Y68iSQ4q/Phiw/uv0/ZghhvMj7eOvOED54fL5Y/d+x4S4IN
M/YVznj+NLme1diRgftDeRkvlY8h2kfeeZ4sElzlzM6dxvNv+093Bk82CcpE
m8Wk4/mpfN2OlyeHBFMj0fj0e1zv/x6NcRWToMbgw4B2PP9LY3XjuEpJEO6g
9sqm4fnTwXeO6zMJmt0t1nTD64VgkFzMWUWC3e5qOj14fblrVYZyfCXBS4JT
igx4PYYFBGpyNJDgnPn654kYDGVW6C0Qm0jwe4BKRE0Urs+RPgviD/z+CcZj
+uEY0pZfUSR04eM5bHhUBeLxO1TR2HvwfnOCflTtIYZ80x/Hs/eRYMgX9acZ
AXg9i/Mzsw+S4PUUMbqZD4ZE+A91sY7j/ctUcnnXLQzNbDj6MS+TYHDxrqMH
TTD0n5nevrUVvJ9J1Y5MXsTro1C+a2YNr8fvb2VfnMdQvOO0Yt8m3h9Y0/Mb
8P54qsdlLI+VDM1fT79/fQxDCYWuulYCZKj9ck6TQRRDbDsurBgJkeGvT3Bz
WQhDro6H0s4Kk+Hdgf6pETIez+75TeVdZMjlu/d+DCdeX+Ee+ZxSZCjc4hDl
QqehM45elDJlMhQjPJ5coeL7Sb1Jfc4RMvxQt3L0SBsNie1WdU87SobRUrUm
N1poaKF7pfm5Ohn63xiWzKmlocTTdwPsTpFhxb4SkcIcGlqU8BsnXyLDndxD
zwsDaSil+8Gnmx5kqNjsaNcjS0Pj9QeHwj3JsGus/AqUxPenT0PEXG8yvNRr
WvF6Jw0VR2hc+edLhlKVU91KPDTUcnaN6BpIhs9BhVTX/BhaLXe94vaCDOMs
tZNZP48hg+QrHLdLyDCMeLzIBY4huoOG9d11MvR5evDv24OjSLdcDnpu4vO7
MA4ekh1FCdwUCbdtMgwvffm3THwUqRfMDzkw88MrCZuvinhGkf965hUzTn54
uvZEmMX0CGIOJ11R3ckPRcoelnhkjCCuwgnLpWP8MIWgUigmOIJ2bUaau/nw
QzOLoUNaLUNoXOdpE4FFAAqqe4wZnetDX6RbgNhTAciZqdH6vKUTORIyWHYT
BWGhx1fPF00/UUnyXo0vkYJwPsPplKjKV7R3+5RvTbQgnGM6xzZH/ooiLK+X
oReCsPuuWGvffD1yEMlQqYgXhKTUksq/efVIJEZGoShFEP5NrW8gKdQjv2A5
8cx8QQgu7/Cy31+HtG4eZHzahvOnpuBUtVr047h6/UU+IfhF5+nl50qV6Emp
57o2WQgeZFxgyN2oQDoquYeOCwhBo7/6pza/VqBqBcnEPcJC8GRHyKSURQUq
3E3wXpQQgnFlAmHvQspRHHunbOxBIZjUa0gN+VOGrnU4hfWeF4JOi/QL5PZi
JGGYXvv9ghCUUzzwWuptMfrVMrhaYyQEx+JW7Rzci5HJ1wvXs02FYK7e6bwA
cjE6U6aq6XtFCKrKbdgtmxYhhSTmWdFbQrC2vOolz59CtOqYYGQdLgSJbUJ+
yofz0XOfVPaACCHIUTUsPMqaj6SfvC9/EyUEy/8TNGvqyUPG78r39McKwbpV
0WHdgDxUMNG/YJwkBMMsLvqJJOYiZ1vRGL08IWjMGGEFFj6gYcs3nUfaheAt
olbKreVM5O2SHmLcifODQOfY7EzEc++j+u0uIfgWmEYpX8lEGq8r0/L7hODf
tWbx7sYM9Kpv4Pb+MSGYY+8cU5yWji6ZiglILgpB07cOXyddU1GDYZIxtyAF
WqZ7H6/6mogc67AqdgoFauUzH5MxTkTcR/fLMolQYK1y+14b2ht0UbyMvryL
Ak2miVPnWN+g/r/tScNSFEham7XeMkxAM0/Y/hSqUCDLvGhg8444FLGlb5Sj
SoGmqbQt0/cv0WH36IosNQp8yOC0VazzEt35T+L5Gw0KFPQPez/9MBYxyR5T
Dj5NgWynjseoMr9AAvUu/peNKXC7tbzjgnwUKj36aeKSKQWWd0ixcXRGossf
6Ibn/6PA25Hl8i/uRaLE6OA9pywp0NDWw35vRwSStUlpOHCdAjNaFO8Y8D5H
6tudJEZPCixs7Cx5fPgJGnTf6bfuRYEfbnPf/S8pFD0Yv4ot3aFAHYpdTQdX
KPr6fabkjx8F1skHESumg9H51wTLzkAKVOtYsLSofoxs1NQzM2MokKzdNzRa
/ADd+Pt96VcsPr7YWzsBxQfINeXKKVI8BQ7OzHCavruPAoiPRvze4Hrwlebs
OBOAXvc27ryYQYGEo4RhaOuLUp+ZOwZnUaBneF1nvKcPyj7xt7QimwInRsmx
fE/vorJ3vCYyuRSo0fet1rLWG3XdNY3YKKFA7dCvCY42t9HAgd+Dhz5TYMHo
/e/lMR6INuJzwL6CAq/+872wv9kdLZxNbGqrpsD36VaRTqfcEK8IxpLZSIHs
XwedG1udkdAPL6NfzRS4PyX+vq23ExJ7SHjL9wP3V/jjcSDtiA78kQd+7RTI
0z0k++m1PdL97H73Qj+un5DSjaRpG3ThJktD0AAF/qQ4nThOvYrMJGMFK4Yo
8POJwjqhr9bI/klZoTSNAjctwpKlH1iim5p6TBbjFHhgulTleow58pr/ZRg5
iT+/R1NOuvQ/FHiZYWZ9mgLDqkw1Dx4yQU95otQP/aNA42pH7G7QJRRVu+ep
3RwFirEP0rWnL6IUOR3ZtiUKzN6nS+xtP4+yBnu82FYpUKGUg/CU8RzKi3Ks
P06nwOKWHzu0L+iiKvozm4wtCux7nuDWa3Qa1eeI5/czCMPpzsnfBWwnUYtN
/jYvszD8uPkz82UfRJ2CJw20WYVhru8tH94kDdTf1Pnal10Y3s808+FoVkOj
/nZTeURhKIg1eIupHUF/lFbVxjmFIYPV0zfnfiuhufHQkJ08whDrSjmkVKKI
Vl/t7DbkFYZpSy5KLBJyaNvgo1QQSRg+00tIIXNLITZmcLucXxgWKKlN2vhL
IO6S1tpZQWEoNPV6cnKCgvidbPikhYUhiKtye9HPh3aKL14x3ykM74BvJcm5
7Eiy43FOxC5hmJ3zKNDHfb1KLlhos15cGJqSUvbLGU5UHTr+Tm99tzC8GSyO
1eYwVf0P7ejetg==
         "]]},
       Annotation[#, "Charting`Private`Tag$876325#1"]& ]}, {}, {}},
    AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
    Axes->{True, True},
    AxesLabel->{None, None},
    AxesOrigin->{0, 0},
    DisplayFunction->Identity,
    Frame->{{False, False}, {False, False}},
    FrameLabel->{{None, None}, {None, None}},
    FrameTicks->{{Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}},
    GridLines->{None, None},
    GridLinesStyle->Directive[
      GrayLevel[0.5, 0.4]],
    ImagePadding->All,
    Method->{
     "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
      AbsolutePointSize[6], "ScalingFunctions" -> None, 
      "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& ), "CopiedValueFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& )}},
    PlotRange->
     NCache[{{0, 2 Pi}, {-0.6807879918364875, 0.3427129903206282}}, {{
       0, 6.283185307179586}, {-0.6807879918364875, 0.3427129903206282}}],
    PlotRangeClipping->True,
    PlotRangePadding->{{
       Scaled[0.02], 
       Scaled[0.02]}, {
       Scaled[0.05], 
       Scaled[0.05]}},
    Ticks->{Automatic, Automatic}], ",", 
   GraphicsBox[{{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJw9l3c81W/4/48R5xCdg8MhVKSh0iCVdd+prEqlaNlJEhmfSGYks4RQlD0r
e4+4EYXscY49zsGRPRvm9/19PH6/r388no/3+7rv63pd432uPaY2OvdZcTic
DRsO97//2RBXM+Okpwru//1J/apQFH9rBP4/r2YfNesJtfk/HhOw+FsY6vF/
/Eq4jiez7MX/ccoB/Vhcpi9IOGD8aJp/q6Ly1Pyx63GB4IBdTu1DFQHUq+Zd
nRz6Bkyme6THLu9Bq9LbpV2rQwCbGkM21fQQ8j1BtU8uewvcZq3vv5eTRQIK
cWXN+eHgqffUGK/3GRR/1pL9b8Y7ELrDb7z9CkQymnKX96RGAqnoVxNJ586j
squb4VpxH4B8vjfVPloDadyqG/wvMhrMCVoTjxtcRl1GofujQ2PBxWQr14dn
r6EcWZnTjNo48Pmln7MU5w309NTb/qfV8eAAy3v52i09xA4MpOLLEoHxTVUY
7G+AGlSr608WJwHv5k4V3dfGKERt/+OG/GRw2sv08DkbUySuvVC4lJEK2hYf
ZBzxMEej1/Tu+n1OA5TWRbQrwwJ90S3bEk39BN7di4g1WLZEpw1eaqjFfQG7
Isd2JDfbIB0r4d73oVlgpWHEjXbLAW0/DVCEejY4pO9XltLpiOZL7Egjtdmg
hsT8hzNwQoXl1AKH6hzA9cQ1T+2tCzpXG7cZU5YH1Kl2i61anmifesfVX2fy
AenCP/PPCl6Iq35bomxxPhC+ocFafOwFam20VK/LLwCLz7l77x99iQw65YIX
MopARS7p1zdHfwT1HtAVDxcDFzXBvt/vAtDe7kg5n8/F4C1KLddAgWiyb7Nb
JLUE5LQYRyjvDkJPGXV7zseVgQcvhfCB4SEoZFE/LyIUgQzF/TLtx96h1+Qo
pUqeSnBzikDfYHmP2nNe7jqpUwmsWHDvEqjv0d1J/THxvkowUp2163JAFHp8
l9tucboK2Is6yw2JxKAwZQu/SGINcHGx0L08l4B6eq4/4tWtASoH5e5v/y8R
iTsC7ReRNWDo2IFxu7+JKDVTUMBaohaYdedzSeOTUemu2lgo9x3MijQ+k1BM
RUMsEkXMm3VAce3VzqqZL2hvDE+U/sc6cC3k8XvJwHRkofDPtW24DpzOWIwk
SWegRftW1bKH9UAv5Hijn1UmYh91bw5yaQBvUnem0zyzkfT3vrGTsY3g0nbN
7jClPKRYq8Wl+l8rEP7RwCt9pgQRO1NaTr1tBVm30Dv6uxI0SmcNP5LXCgiv
Ei5K/ylBQbiyXcJLreCCQN2te0WliK546OS8fRsQgnuKjFW+Ir9cbuMY+3bA
2Sp+8dN9hKixjQWrdp3g2M9uMQrhG7J11jbJs+0G1ruG9j3TaUDWxfeDWCWH
gUOt/YrdaifWL9d2s8gNg/92smQsi3ah+R6l3K3zw0BB8aDiE9CFnqzxU9fN
h4F7wewZC+8u5KxSLfbn8zAofBJ1k8JLRT414hlTJ0aAFs/fRuXdNBTXSvvZ
qUoH5Knlc7VKPcjU5JtBx3U6MJTSFRA27EGSi5lzbWZ0UH35eJW+Rw9K4ffh
b/Ghg4sv7QVzq3pQhp7c3foGOjhO5fb0v9CLSvuDJ8t1GEAtfZ9xqGYf6pzQ
JKSajoIrxPmgy3AAzUUeuPDCfhQ48Z6xaNUfQFwXOT2NvEbBVcaXWwrPBhDM
qPkrlDAKbs71XajOGUBf7MCE38goGE443Bi5ZxB5rcr9sDQeA3yXcbx71wfR
0e27X8oYjoMmr6hCrcvDSKt8s5Lr8Tiwv+VRn3d7GJk9HlgfdxsHuNuCpdzm
wyiqNepJbPQ4+MsikerrPow4wslmxIFxILq+r/9L5jDqF+M+t3iXCTabbvv7
co+ggKO/cYW3J8C5BVpLQuEIcn6qeWjccgL80zXOdqkaQQ/RR11B1wlAOST7
XrFxBGlcOffZMWYC8JIfjtqMjCBOm+DrpxgTQGnP+sUULjryzpROKbb6BUTp
FwSlb9HRk99urRNuv8DUIdM1ExM6uqfStkp58wt8lTWPc7ekI9Xmp1ee5fwC
zJ9R8fdd6Yhltubvmd+/QHT4/c7xGDpyP2J0scxjEhzeEn0ZOUhHjx1yHaaC
J4GTUNc6zzgdGZRvi9uZMAkin26duD9DR0qX05ddaiaBxW5bs7I1Olqz+het
RJgCxEsj2WxCDDSZf6nOSmQK3DyrOBcgxkC967GLHw9NActZY85lSQYqfaWm
vnF5ClyY9rz16BgDOaW/nS8PnQLV6f7vptQZyGJ5XGQ2cQps031d5n6ZgW4q
KVwQL5gCl8KvTyzqMJB843CkO20K1Di8/+dmwEDLU0fOAdFp8Mtm6o6oHQON
ynpa2xyZBtM8RRIFDgzU4dL5LlZlGgxnOcjLODNQLrfL9JbxNLjb7/2r0ouB
7A7VhVcmTwPDILV95qGY/f3R53GF0+C625D9bDgD6cXhrJ7/mAbn0iTf3Ixk
IAXyGVX4axqELQjM/IxloPQrukd2r06D/SpfuWgJDCQeYEfBcc+AhDds2dXJ
DMS69WkWHZ4B+JbOS4pfGOjJme89scozgJrcLVGXwUDj/9FrPLRnQGCJrfLJ
bAa6lbmZZWQ0A8DLxFdeuQzUMCHyAdjOANdndbjsfAZSkjzls8tzBrw57BdQ
XshAmQbX7bZCZgAhoF/kSzEDhbYHqqO8GWCUXy2w7ysDsfOknYitmQGkxGGT
vHIGclSvEfPomgFne7647kYMNOE5jDcanwE6kUU3rSsZ6M7X9SWVPzOAzffX
WGQVAzX+pgyJ42eBTtszkZRqBlI5frJhkzILHO8uTYR8Y6DsR9cKBg/OAsEx
2QsGNQwkkWIdV6EwC3KOW0niaxkobNg/MObiLJA6utspFGOOnSmO7vqzwGbn
/tPrGDvpVpsYWs+CD5++3LrwHaufN4OXVNxnQZvU4RZrjPUbVk+Jv5kFoe+d
QxwwbmYXktyMnQWKv76G62MMgSzvYPYseL5l1S6Bce6zK//Kq2aB1ugtlTrs
/L35j0aj22dByEhC6yWMI2Z9W9wYs8Atvup5NuYv/mBSqcHyLEi+9FxtCYvH
+V5lsvK2OeCcuilMxng6uj9YTHAOdMvkLJGx+A27/7ps7JsD9bImjSuYXq18
5AcDp+aA2e/D0fmYnqqXj+uUa8yBQeRrpFPBQPm+l5Wjb8+B/LUzPC1YPvZV
PzzgZjkHdl01iZMuY6D36y/5DVzmgGeBLf/9EgbiOpWwqfRqDgi9PmbiXsRA
rnYVv0Sj50BLYYW3YwEDzX7p7VzPmAMqKUHPLudh9b2b/8vXljnwZdP6+7tM
BlK7ezTi4/AcoI8QhbjTGag4/KKn68IcUEoUlTb4xEDa26oPxRDmwYkAtVfR
iQzkK3F1Dy/fPMhJSKUGxjFQJRgUdBeZB8of//TcimagE86rLIaH5kG+T/S2
NxEM9Oid358m2XnwqIB56C/WL0n5gjPKSvOA4yaBE75hIMG5E91il+eB4pBt
3yNfBlq99yhz4PE8eEUk8Mg5MpAxacnpq+M88BCqj2Rg/Vpb4Xzug/s8YFYb
nnS0ZqBg4YDuW0Hz4CXJOP68GQPtb0lj7czA7ONMw32uMtAr12NNuQXzQEU9
Y0j7IgMtHix+F1I+D7pYdl1dvcBAFd4/Dl9pmgdauI28VQUGuqEwrtcwPQ++
qkUREv53/jCt93xangf2/FrbvLH5tCt8Zcp3fR5YU9q0Nf93fs2xe17YvgCQ
vNCoGxcDeSZLfqk8tADGcentb+boKJNoslVouQAeL7DYGOTTkUDFRH24/QIY
jNrMe5pOR88e2YY9cV4Aw/1iPxyS6Oj8d4+DJwIWQMvJgudiYXTU5xJzPePT
AigL+Wstak9HeGZfauLEAug6uOsh/QA2f8Pu2XvNLwCF1x9GVHbRUcfZKSWT
vwsgBMf/y4NMRzEfV9vE8Ysgo31oIIuFjk7qCK9H7l8EUu9OrdV1jyDTcr2r
weaLQOLTv/wZzxF0opDCJ2yzCJSpX9ySHEYQa1ZvR/zTReBut+240sMRlBBn
eDPPbxF8ZFzJ5706gugvzA2pnxfB9D+pOlVR7DwtRyvRuUWgv4eyYwP7vpnS
wv3Sni6B8VCH5GXXIXSi9abW8edLQKtNe4jVYAix1gtvL/VbArw6hYs9SkMo
oTT6zc/IJTCybJ3NWBtE9OjkiNmyJRC7L3De6ukgMjUrSJLDLYOWPpeSt2YD
6N5CZ2Wl3zKQjLHYpSDch4ydg28WBy8DhagX8q4rvUif7dJs1vtl0GmQ8iii
rRfpkb/tjE1bBrVRhqo3/HuRxpkcR7e6ZZD1JvXUxu8eJOP5+rACfgWo9Ggc
Vf3Zjf4S1d7n+a4A74Oi7VG3qejVsSKbFJ/f4LnIibu54m2ItbLolLfzX5C2
TX+1xaoYkaQvNg+9XwXrw4e0jwUWgnHOWFJd3Crw7PvyfeexIlA6tngjOw3j
WppZZGcRMIuP6ntevAp+RCm/uiBaAoopkxN7ulfBXN99vGJKGTDk9GczE1oD
R6UFcvKvVYL00e+nJ8LXQJmGmVuZaC3wrBZxbY1eA7eqewW2PagFenGPUXHy
GijIvezKmlMLcPqCav4FawDHc0hO8Px3cL3j/nXprjXw1zmw8I35D7Baxf7Y
SmAdMKkp5vBjPdCIVU2cD10HN3ZzjNz73gR861oHvD6sg06p/7ovcjeD7wuG
FHLSOnhNdez/dqUZXDjv8vpMwTpg6RhUNKE1A9VfeU+9aOuguObotPJIC1CS
lbokILoBwp7ffvXfXCtw0c/zSdm7AQ6tTY6p8LaB0pdnq04f2QCJDtfGGg+3
gTPdBvIGKhug+MXt4KMP24C8W8TuFOMN0Jc918Yx1AaO/eBcOZW8Af5NRhYX
V7QDm/nwow0ZG0BYwcrDv6cdZArvtdQv3AAezz8IU5bbwRErOPT8xwY447Mn
zfpAB5AmPauv/7UBjDaf2y6+7gB7705G35XZBHPk7qZ9lzrBPW+n7hn5TfBn
6+aNfNNOkJDBwf8cbIKXvLmOa06dYA9O0i/p6iY49sDsXUFSJxBPvms/Y78J
HIUC0tr/dgKhuUY1j6JN0LE7/6Tb2y7wPgSZ/EKb4Ofx/Rwnk7sARS7X9Xrd
JpBQMApOLsT42bvc/T2b4F+/FzG7uwuIsJmKt65uggEly/f3hakgKvnGGQW2
LWAoYcQ9eIAKdmqo30ji3gLSK1m1PKepQPT14QAn0S0wtz/kXqguFYgL/VnZ
rbKFxfUo69BrKogp+UUMVNsCs5W1325HUcEu/f5DK9pbwPkAG005lQp2x1eZ
1BttAfXeN6b7K6lA4tCrJjvPLRA/6cDrOksFCU3uE33+W2BhsZKs948KJG3t
2NRCt0Blp/rmKBsN7C3QOyOSuAW8tSkLuyk0sA/sSa6u2QISbu6/eZRpIGWE
v/Jw0xZ4NBP7IfAC9tyboy+iawuwpc0dLb9MA/vrp4iPxrF4Ig6VKhrQQNqj
wUNds1sg07ffxc+MBg7wtqlhyQD07gs1fo9oIOzg3QoCOw4mdtr7JTrRQKdz
508VPA4qr7zuKnenAYHGyz3/bcfB3BnipUBvGrgh9mM8jYiDQgEPqNwBNBD+
GC4PCOBgmkaOEXhDA1RUwsIvjIPfI6s6pMJoQJAku0NDDAftxXfvQu9pICJP
SjpXCgcjE60WCPHY++yxp5gHcfDB4vyNsiTsfT3KBVEZHFTwmtfYlUYDeqkh
OtdO4GCtjkv5yS808O4vl7GPPA6mjuslbmbQQLemt3WZAg52OGxOPcumAcqH
Ded5FRw8IJ7on5pLA7emHf2kzuHgE77DNr75NPBeeT78jjoOJvyQDxYuxOyD
Hia+uYiDp3cLjRgUYfbD9OyaKzjIG0HTu1GM2R/Xr/h3HQfNTc8w1zGO9Or6
KXMLB/cl5fvdLKGB3g7tnnv6OPjQ2+mEGcYiUnXj741xUP/2o34JjO84nl1u
MsPBUz+kPd5i9lE/SlnYHuJgwRC3QCF23/xw9rM5KxwMIGudjsf8UV9NWeiz
xUHrJxk2GgU0EMMf/bDuCQ62q7gVfsmjgeXDb0fynXCwpS6Q8DOHBrTU/G/H
u+Kglp+8eUoWDcQbebS9fo6D+L3mDcqYPn+dHDSdvXHwmWeQXMhnGtAOfVRl
7oeDzm2ZiXGpNJD0xeTM9Vc4OPXhmpAVpv9azc0cEIyDa/Wv3izF0YDO4OWD
h8NwMDyLhedUNFZPf87FU97j4G3+gWD5SBrYJCoIb/uIg2EnaaKLWL51pY+F
LMTioARVNu9hMA2wGoh6NqTiYP7BAa7XPpiejnz/Cr/g4OSGZcdpTxr2ncDb
JWbh4JWIgfw4FxrQr142di3CQe3Li9QkGxrI7ZvstijDQf+h45JnH9IAYWX4
qi7CQRehkNhIUxoo3N8EZX7goO0/5gWPGzSw/ey3YpGfODjK0+OxHeuHe3dK
jnG24OC/1/0EPaxfiK+Tdw9TsXgEueXET9LAo0U3XOg4DjLSiARvEg1Ucz9x
cp/EQZUbLcxlTqwepCznLWdx0AKXbCuxSQXfb+oNn/uNg0Zsk0V5v7B+L5ep
XNnGAq1e+FRkl1OBI3XvaTqBBTJr+wLxuVTQNCeS3czDAjVKOJ6LpVCBswRn
XCqZBZJ7vaQ8sHnS6TvkcVuKBU5M1XpI3qYC6fiuPxcOssD9WrlR5hepwLP0
p82JIyzQiV7u8kCZCmRmioy4T7LA7mvckwm7qcBPJxiUn2eBb/zucrGNdgEl
sbNbu81YIMtRYra7YRfQ0O+0N7Vggd5d4/F7L3aBGx8ejCdasUC+l4dj3E91
gUfCQU37HFig3MnrN67s6AJR5L4Ph31YYC/MkrUo6wR/eRxOn05jgcum17ei
t3UC9sv49GfpLDBFKjvh/UwHIL76sKssmwU2aC6JanV1gP1c1RwqJSzQbWOX
aXtiB9Dj2NF1rgGzF+z65abcAfI30+yuTLPAoq+PuezutYNKZaWx4HkWyO9k
Z7+o1g4aXVtutS+zQNZZ3eMy0u2AsboCdDdYYLDouy9dc22A748q711eVmjL
Iug+9bQN2M71fzY/xgrTisalrFxawaFh4qjbE1aY7EszM0psBNrM3dWOTqxQ
/UpzI9OwEdjOHouzcWWFRs08jzlEGkH++lV9E29W+G9G0x7/5idQEg7pPB/G
Cs3J8ZdknjSASzqkWu58VrjzvK6786E6YFVDSolcYoWbtl9+c+l8AxlpfBYF
9mzweIenZqpqDvDuQmfZnrLBTOqv4RX3bKDPar3zmgsbLDGMuqvjkAW49H80
T79gg1E5j+QCjNOB+Q7Xk3vD2SDfXRqTcDAF7HYcY31bxAaDV/GSG1xvQdi5
4mjbdTZYlWFp4PssDrkOGnQc8mGHjhrmHk9iq5Cd0umLMWHbINrXcv27LBXh
Hod5vHm/DSp+uYfkzagoKHY+7/nHbVDMhLUqNIyKPrN+Er2XtA0yHlduca9Q
Eb1OeHZ//jYoGJxV+zefhm7orgXndGyDIbJv282P9SD5xxXUGj4OaH7QN/kf
Xz+qiRXhLhTkgBxV5EUJ2I902hxBqggHHDzZ9O+QdT96LHc0zV+CA16YG2nq
+96PUldjn2mf4IAiRWsf05wGEMXXU7T7GgdMg+u6aV2DaDX2vOlUMAekyZk/
ti4ZRip/lEOPh3NA4wAL/rstw8hT+1T100gOGF328fbRsWGE3zgowZ7AAc/t
vvAvgDiCyHd20EXzOKDic5dhvPkIOibQa6LdxQGrxh8sLHLT0X0/G5McYWw1
spV5La3KQGlDFiF/xDihhy3fnlc3sP1S3rRKWYITSrmoXGk2ZyD7sRt7GqQ5
ocNo78+lAAZ6rqowQlfkhPZf3gRXtDFQ1Dq7iYAhJwz+zyZC5vYoarWNMn6a
wAnV6y8p3dMeQ9LBxLUnKZxwoVx0T9TdMfQiyyfc/jMn5JCptMyxGEPys/YN
1rmc8K2EEnruNYaiH12UM6vmhE83I91j88eQpfk65zU6J9RKaar4zjeOal7a
J2iPc8JRz9L0P2LjSCx5QunSJCeMWS234Dw4jloZnXbqi5yQ7eqn2WqVcXTK
JKNPmRUPz+WklaZajCN2fcMsaQk8nGtYcyMXjiNDl07NA/vwcI/wHNdBNI6K
orRGpaTxsO+ZpoBI3Tiy7JEX3nMCD4lqntfCesZRm96OF0KqePgsYNqZuTqO
pB1fipHV8NBdkzQ7yspE3uFrRXxaeJgnEXYFcTHR6U7mNI8OHlb5C33kEmGi
mGuVeuym2P2F/bH/5Jnor638Aos5dl5qaV2/MhPpBKcHbj3Ew4Nu8mYfzzPR
tpb3lat2eHhAI6k0+xoTWV2yk17wxENcCNPysCUTfX/ErJl9iYepysrpCzZM
tCvQwGjaHw9XLqz/DnVgovZ6zTBmCB563fv87OFzJjqjLrE5EI+Hwg8mV83C
mEguscHvRzIezrRJXF18z0QyW/b8OZ/wUEfzV5V+NBNJFH074J2Dhz+jUml1
yUwkym+VZ12Ah2UWoZYNn5hI0EZA5WYJdr72wp+0DCbi3o/9/q/CQweF3/T1
fCbieMEzyF+LxZNeddy+mIlwQwUWG3V4WH3okkFNGRMtR3C4t7bioR11/uhW
FRPNLmbiSzvx8GmnWQOzhokmtG++TezGw5y4iwcyfjDRAEdamuMwHlprHqVg
Ewx1m16VNR7FwymOr7EHW7B4K/6Wa07goZFlX5dJGxM1isRryE5j+UngKHTs
wPRy1OwQncdD3a4NRYsuJqpsXzDgWMbD/Qfr7snTmKhUJmpi7g8etjTmHRzp
ZqL8ANX/etbw0HFLN8Cil4myxic3qrew50KBXo19TPRZ9a1fOhsBxgQ47uAZ
YKKkGEX+CE4C7J73kpEexPK9yoj24CZAFh7TAckhJnqv9+rAwx0EKN66XWQV
49BcuTwdfgJcNfhGTx9mole8A8pKQgTYS+SRVxphIh/Ll3VSOwkwx7WYLw1j
z+9Hru/YRYApdn1P5zF2laAO/JUgQNHNdjMynYkc3d0t6PsIUMD2XI8gxra9
+5Z+ShMglx6laxl731K+xa1AhgDPeHjfysLYLPQpPvYEAWqssZipYWw4u+ut
nzwBPhTW+VuC+XNLq07MXoEAB8zdBLkwvp5im3ZXhQCn0zuQPOa/Nquw7AVV
zH52agVg8WoYVpXLqBFg/72kUilMD9XShxoULQJcc5IgTmB6KQnydbBoEyDf
j/a/Ppie8valBlPXCDCP+PApaw8THWs2nejUxeJRGg66g+VDyidvI9WAAG0M
K8PjsPztpt/1CzEhwBdRXq+CsPyKqLDzu9wnwATZbunbWP55f984oG1NgF0K
GzbeDUyE19nIPWVHgMnjBrvHsPphzUxW3uNAgLdvPgmUqGWiP/d/6yy7EqBr
xHyqLGKixaqYgYHnBChkIn6f8ysTTYupW/zwxvxrDVkuwup3uOudW9QrAlyS
YiRl5TBR/YUzafAjAV5aaBVYjWcipx/uMh5xBDjHyO3Qw/rngGZNfnkSATrk
fvgvAOsvn0tXqhQyCPBa/MQLxyAmOn/9fq9cBQGK4R3O3HPC+qHzs/F/1QTY
OvXrqpM9EyXqzY/nfCfAE+X4TDMrzP87LksyLQSYr4Tsvxlh9Wscsv3gMAHa
Rls2RWHzwpZODX0wSoA+B5iFj7B5sstMVDhlggAlZStMBbB54/4gVUpygQC3
a02ozezD8vG4XEWMlQuysvcrIg4mKnL5ZUvaywUT4z7RlyvGkfmGzJ8rB7ig
tEABTatgHAl6PHELOswFuTQetDz7Mo6eeG35c5/kgqMf/SXvvRtHJ/zJidvU
uGAE57X1+9bjKDPibNe/B1zwXM9ZH1VsfhtSfO+etuKCfAlfz8pxjCOeqMYR
R1suGNRF+YhbHUNW0Tdnl5y44LGW0riDI9j3I8mac9afCxZrZ0tWZoyh5JzI
M/QvXPBF+J91YYh9H9buSndncUGPeCWp0hNj6NgF8Z3NeVzwu1GfnaLUGPra
nbBeUsYFmZVk5ghhDHVsfUEhP7ngzl/fb1a1jSIW7YoLZ6e4oNJ0cwvj7igy
mKRfi5Pmhg9/+D1+c4eBJOSSVSNkuGHx43WTIE0GmnB7IPvqBDd83/fOzP40
Az0hTQs4KXDD9L2lkpNkBgo8vUK7osUNFTbLrD620FGJD6fh1kNuaE3VXzik
QkdkycOWhp+w9/USQ7bwIyg4NFDvUwY3VHmkdeb00jDiYp1SXc7hhltXh3r1
BoYRbuTTzoBSbkhVshtUyxlG0zH7mwp+ckMlQ1rhn5vDqEZE8hjPLDcs/3OX
9u7sEHrCJ/K77MR2aJ8oEed7rh+1s+C9hL9uhzZPTfe59HWhAvx/Ka31PHDu
w5sfCjOVaNzVc6aiiQemdhp/GvarRJSFN3IZbTxQgfDlb/jeSuTak17t38MD
3XNk/t0rrkDnP48Pqv7CWGeuIOB9GWq/eEeoAM8LNViZX3lMCtB80Fm/SHVe
aPjp0w0dxQR0mEy0NK3lhWjtmm54UzrgzypLT67nhVkrzUIKeplgVePB3EQT
L+TyVtEqHMgCda4VT2yovFC74Kzj2yc5wGzUysOVyQsZR5NSvznlg+i8+rAI
wg54ZvIBb5xpKeDVeYEaLu+AgSFZN9S7q8HKlAwrr84OePcPnd/36DfQ/7L3
/DW9HVCGK3OXqc838LnkeAPNcAc8dVYHzzxRA9T2DHeO2uyAs8KrNYK+tcB9
XmlyM2QHFInGvzDfWQfmgn4LyFJ3wHf1RYxvs41gel68+kjvDviqRz58v2wT
mNRRtzkwuAPuzhLq2/W0CYwKvq8XG98B/yZWCDRuNIHemDPu+N87YEuVGn4b
RwswvVmxNMtFhIzIwD21z1uBiZ+O1vR2IuyIN850i20FxiXjcb94iVAtnjV6
vbwVGOzcoT3KR4Tt+ucLL622gltDRmk9IkS48FUt9KhtG7j8gMWgRpoIfWyu
6LdcbQeX3ofnVR0mQiUemNPxqB1o1R/kQjJEmDM1shbl0w40DukUlZwgwo9X
B9PsytqB6lwCKUuBCDPOim2ckegAp5+e/x6pRYSnXvwW5xntAKfSukXfXSJC
Ni9vTem1DiDfY/VfmDYRfnNsfs7J1wlkFcN3v9EhwjxmjzUL6AQyW2PO3neI
MMjwW/1seCeQ9PU9amNJhMKRGm9eY/uVPVXO3dSKCH1L3n04o9UFKqXojbqP
ibCZ3XwqWb8LGNYoWSrZE6E/VaO74jm2f7EsJhGcidCPTe/dbG0X+HU1dmnd
hQjj9RWtIqld4HTcJdV5NyK0IAW8Y2N2AapK6hDVkwj7lxO9yZxUwOeqL5Lk
T4SNt971Hz5PBSY/CQ/fBWLxp/zxbtGhgmyRoqKA10Rof+Nb0TETKtAuIena
hWDx/fU2lnGjgoDf39+oRBIh6wujmwLY/tpz4cng8Q9EOK27ZyKtggoOhO85
IhVNhM72VbfWGrB9WNa1gTse0+PH9aAeOhWw2h3f1pNGhH/0rt9+toMGrlUO
Xm/8TISxcHJpVZgG4na8SkDpRMheEzB2Zi8NqGSOg5RsIvxn8/rt3CkaeL3x
NigyF/M/nSxnfpYG+i+dHXiVT4T6nUtLkVo04Dz1wfm/YiLcLkuNPa1PA/UK
mvXmpUS4tAyfJJph+3rAb6E7X4mwkNp8+qcVDRQeuFYAK4kQt33d4porDSze
55jnqSNC6Yvcpqfe08DZgjwVlgYsX+Lns/1iaCCY3eT18k8i5Oj4+SkiiQZk
Esuk+1qIkDru4TORRQNuixbPmtuIMKZlwPBwAQ00nhWsq+rA/M+daTxcSgOW
w7b302hEWGF8/869bzRQclQ8/0MPdn/Ruf8if9AA3uMn65s+IjTbWSXw6icN
3Gx2uuY1QIQlt2aVFFtoIEVsX5zDEBEWN5lNprTTwIpVx6zFCBHyPlATbemi
gfNfnyvrM4hQ5b1ye143Dbzllnl1ZYwI7193267bRwP0O329qkwi/L1V+TN3
gAa2K99MNf5FhIqFz6r/DNHAyV0d/7lPYf6p5mnMjNCAIcsV+HGGCANZ3U3e
MmjAj9GwvXSOCL8opBInR2kgt1ath7aAnc9ScmdhDMtHanXyyhIRlkutnv80
TgMcASr2/L+x+AUed5GYNHDUqlTl+F+s3is6+Q5hfEtbnvvKKhHOOHxmX8He
9zqWS7NaJ8JtHfvzbDBO55NJCtgkQr3gs/s/YOdTlz/ZpuFIMPdBmuVT7H4c
TUr5OysJLk7X+bJj/h0siSeMspPgP/Of/qqY/zofxKgsnCQ4kPLW9TgWn6tb
ZMIuAgmqDFRatfVjehqRbZS5SXDCOuOBRC8NtJwNUbzLQ4JhVY9d9tJo4K8k
D/7ZDhIcJXcV0jpoQILDvzOCRIIxFx8dVGqlgUsT7PH5/CTInefIuNZIAw4N
z63bySQYS/29QqmjgbogJw5eERIU1Vl5VFaB1ZftcvshURK8O8G/ElFCAzuv
28ZqipPgecjcszefBh4LWZx+KUGCgar0kutpNPD+H4M9cS8JKtTaX/0dTwPV
fUZtlftI8Juc/4jWBxogx920XJcmwXEqmw7Xa6xfvDrkRY6QYLpUiKnDSxqw
MLvCdvooCYpw/G4JdaOBrwfUP/4nS4JHXzoeGsPqf5zrm0XISRK0umR/VQTr
jx0zKiezTpHgHt1tkPUuDZjmyDdPKpJgXbTn8qQ6Dbx6mxuFVyHB1LU/K1PK
NFDgIPNgHyRBh/27FONksfpV2IczPY/pc/8uSVKcBjKrySd6LpHgKZO2qPA5
KuhOCtn8rU2COe6SpnIMKmDx5fkpcI0Ei4t4PH2pVHDj4jazq7okeIGrS//i
VypY61gO/2FAgn1HZjove1OBxljHv4LHJOgee/GiKQcVuB0KcsmyJcFLfQd4
n//uArl2Ghtp9iT4UW3mvOF4F9i5+RX30ZEE9Z7Qv+hj83JWMIXTy50EQ0xd
f97H5mmYhhNZO4gEK1uvXWue7AT1QSfeqQeT4INwchcHrRNsdk5TzoaSYJz2
rCvrt07wwMREVC4C0+PW1Mq5qE6g4KwlKRJDgnl5rB3LFzrB8BfR4+MZJIhb
UF46E9oByIvU3KEsEuShHu4pc+4AWqdD5HpySDD+fiYP+70OkFez7XRjAQnm
mwkoT5zoAD4Dsyq55ST4anp7/bbWdnB4R9Ul9yYSLP1Un/hjqw042d+3EJwh
wXsvzJV4TrUCv4uzDxTmSHBOczRRWbAVvNv79IHhAgm6Xtl8YT7ZAgqpfuYp
KyQ4M3fo36BZC1hSSDeT3yRBmSYt9kc3msFjtiVjvR188BzLXpNX4o3gXtjz
2xHH+WDuoLkmj9kPoF34QUvQkQ+yPu94boRKgMWS6RbZiQ/ek/vhtEEpAS+O
SeeTnfkg8bX6kpR9MSj+XCxGdueDwrYFUqKSRUAitmue34cPfjW8kFzrkQ9+
++54R4rgg9vVvpjGa2WCmNsvGNsL+WDv23SKeGYwmF23dGX7jdmvdWt0+Oai
eJpn/uP/+KHGp+1l9/Ja0HjtsaEgB36o8qS/mUBoRYfzhwhZT/mhOG3Q6fLp
VlQYrGw058IP3VmGz0ZEtKJGzX8EW29+qH+u5dsenTb0t8zWyC6cH47zmepV
17Qj7TgjridF/NC40N5GI7wLhQXxngwr4YdiHh73Nau6UI9ruVF+GT98h3/G
enCmC5nd3lmwjPhhY6Fditd5KnrGRzNyqOOHlzfMR8IWqCjJW7vAsYcfxvk9
axs+341WLZSNn63xw+4N/o8BTb1Iq0waOmzwQ7ODUUm8i73oAw9lt90WP/x4
Zlf0I8E+pJS7OGTBJgB7fsZezDLsQ+5rqUa3uAUgOiN/Qm+mD7EF8Rmd2ikA
LzXc45pgHUDb85gGKwoCcDY6WcqYZwjpb+tSXlASgNpeTzK9Dw+hjJvVYjMq
AnDxY5ee88UhpL3+cWBUVQBaRfgf6PYbQsEXdAw6tQSghYHY9BfcMBLo/qqf
d0cAPq689fbjk2EkthFy185ZAAoLD9WFaY4gDfrlkneuApA+v5+LdHMEPflB
EKpwF4ApNPaTj81GUEOIZzvXCwE4PPNhuMF9BDnus9NMChSAG6fE27lyR1Dr
launaB8FYOXcsRvBfHS0Jrc9fCNGAOY3Hu/aIU5H+0TqFiXjBWD49qWfjw/S
kesoyLRNFoASfoP1DYCODj47KsWVKQA/twmN/bWkI68EXn5lJABt//hkiZXR
UYZvg+29KgG4qUloDqyho24rn2b/bwIwpnVcvbuJjmRObfpTfwjAA86h/NuH
6ajv5wzOtlUAnu+Z8B1gZSCOnE+GEe0CkHhga02Jm4FORNz/+rVTACbZBRo7
8zOQn/HAU0KPAJTst2kMkmQguZXG2YQRAbg+H7UhqspARr1+l+oZAvBbN3eV
JbYfBWAbwtyYAOwYPbkn4ioDDft/va80KQBFXqOwl4YM9Fo8vb9rEXu++47F
dycGKmazUFhfFoCy/cUyOA8GGp2QfC/xRwB+zzTWFXrJQAp5H67brAlAy4hz
6QPBDDSuHtiAZyfDHl3aQn0KA5leNN4txUGGAl96XKhfGGhA+6TjWTwZrpon
eFVlM1Cn7tAe5+1kmB7rs3ailIGu3s5/GsFLhoV6axYlFQz0U9+/KZdIhvcY
E5ri3xio+p7ssykBMuR+Gm384icDqTzAt3AKkWG5mkPYixYGKrEc2LtXmAxZ
7JinTDsYKNvOt1VfjAzDFp42lfYyULRH3/5mKTL0IVy70MNkIMqLbLfJ/WRY
XZ2wrWmSgcJ8XnZwSJNhMIGiEjuD6fX6qAeQIcPoK7skhhaxfISwd909RoZ9
d+kdl1YYyDOsR9rpBBlaz0pzRf5hIKcPL6jZ8mQY0nRPoGmNgZZjbh1uOo3p
gTu1kL/BQDYJR7x+KZAh00jXzGWLgSaTWbu3KZNhnS7nf7tZRpH5J9oRCUCG
5ySnJJJYRxE9Pf2FylkyRD+Wn3GwY/tvtmfPnXNk6PrL2l192yjqztM7+vQC
GQ58XZF/wDGKbhQdevlWnQy7Jp9/vs85ilpLcX1ZmmTo5Xh7WBU/ii5WdB1r
vEiGa/Hiw5sYf6/67DNxmQwl22tyIgmjSLXWo5/9KhlS7cIN+LlGUUXdjRN7
dMiQo1Nm+hHGZxoP+infIMPDZG/zRIzzWzYHbuuRYdtaf1cJxsc6OmQdb5Gh
QnfXuSyM06lp/qF3yPDqPaPSFxjv73UbytQnw6iPFedOYZwwoHPypyEZOqS+
mqzH7hcf2R/INCZDl9TCEkWMo0bXh9nukSH9g9vXV5i/5Ik2+d33yfC+guPm
Vyy+4KmUV0oPyPBvnV9gIxb/9jkX+q2HZNj9e96mDNPHb/HqaYdHWDz5ysX+
mH5sv6WCQqzJMGn7NbtTbKPI/d8qI8OGjPWrUVoNpv/qesuZBjsy/CC7aHkc
N4occclvxv/D9Bemf3fbZKAFNucxVkcyjLPVb/i0zkATXHtDFJ0x/VJUWOP/
MpAZ77/xm65kCBKjblj9ZqAhUrPSE3cyjJm3CBBaZiAqxWki3YsM8S3fh3jm
GOia6GWVem/svkzqrtvTDNS4SyJszIcM44Vue3v/YqBv+xqBeCAZ3ik5uuXC
YCCtoCEtwddk+GNnzIXnwwzUurKoy/uGDFtT9+54OMBAgzXCVpuhmD8brA8m
qAz01/TB+8EoMrRjnaXF1jGQe4NzIvUjGaqqGzLqaxho24mgzOYYMtQWI811
VzIQH0tBTUUCVp/7oVZKMQMdjmVdiPlMhhPeuQtbWD/ncgiuR6ST4XFRn6Hr
CQx05vFBzjeZZFgmEXH7dTQDqSlfFfPIxfL1u1Qp5y0DmfRFaxqVkuHR3bSO
CncGihA6k7CrgQz/mSwms2LzR9TjUoZQI8ZgwqsYm0+J40bFO5rJ0KKZ96ju
Oax/C3yat9qw86nql8/JY/PgeufqUA9WfxzlQisiDIQLtrke94sMnz7Hq38b
oiOfP16G76fIUOJT3fA8jY54jCIeBs9g/Uq+8Y2jlY5EZMo9ni+Q4Ztx3bop
REcnG7nSjf+R4caLrsToGDqyxKew7yEIwi7vqPaim3T0k9Xv9TluQciiOOiV
cZmODm88FDTnEYTw0sRF/3N0NLdw5MAXkiCUvOJYsiRDR0/6Ci6eFBGEyyzB
+PfsdOSWWRuqeQiz3+nmAL6MoOAbY7vtLwvCm07vaWLMYbSg/eNT2BXseXbC
0z29w0hH89OJomuCsMfZ7TCxaRiRVazPr+sKQrxNlldW7jCK2r/ywMdQEMpK
KGeuug6jpFX2rCgbQWgrlMI1yTOMiuL2Kn8LEYQ0LtWe/vODaO/WeZeqt4JQ
yX6nl/GBQRRscL8EhQvC72F3K+u5B5GFSMrJr5GCcOvbtIJ2+wASCdsvUxAv
CLXsk0b3Gw4gV1/pXak5gtBZ9nqQmV0/Un18jCWwTRDG7vZ3gn69qFlRqVaH
JASfmnzlGvKnor+WH64bBwnBjR16tJtyTeiNcyKnR7AQjBRTefxpoRHtC/hS
Fh0qBGs1B8rFshqR7qcyyb4IIUjz5wq1k25Eucy+Jd1YISjXuqzSJ/ETWZmJ
hl3MFoL2pNuvDHfWo2GD6E75diFIYD+6X/1YLXpqneyn2ykEv30bFTVaqkG8
bhlKT6hC0PJ7sGBOYQ1S/lielNMrBM+6SHZtU6lBUb0DTw4zhOBSo0Poyyvf
0I2b4mSJZSFYfrBN4Yl7Faq7GqvLI0iBdkQ/pSN7viLLmrEKTgoF3tn2S+MB
KkM8pw8fYBWhQCk0tnbQsAzp7CpZ/S1GgQ89n0peji5FfTPtscNSFOjrVtmw
V6IEzQZwTOadpMDzoQf1mGqFKHjz0vXMUxSo+Nb5t8hsAZK1f/s17QwFNh+Z
aRoLL0BOt3e/iVamQJ7+4NDzzHzEekBBzvcCBXrtMD+qEpKHyLXW7nd0KbCU
7urRy5qDik/nM2/cpMCWwT8vx/Ky0Z301atXblPgtT/fW8LuZ6OYt76S5w0o
UNy9zvZ7eBY6YBpfd+Q+BQqIfv6YppCBlLY6+VgcKPDApaCj12vS0KD9Ttc1
RwqcUwzTi1dMQ57jJmMrTpg95wxNOz8VfW+aLZp0pcDwppl9TZ9S0JWPeINO
bwr8Jhs85JychEzPKKWmhlFgG/s1q8zhOPRwpmmlP4ICY8bV3t02iUO28Ubn
+SIpsBevVJr8IxZ5EF6MuEZT4MXJf56P+qPRx576nTopmB6vF8t/Ho9Cia/v
WvqmUWD0GGld92ok+nx2pvjrZ8w/h1izt/bvUcknot7+LCwf7FdunKiOQNRn
N4PXiygwre2Wq0bQWzRw5Nfg8VIKPLGkciS/IRSNjjgfefCVAn0G1YSYXKFo
STOmoa2SAqP8P1Z8iQhGRJEx9tR6ClS/+TynrygQCTU7Xu//SYHej3TVE6wD
kLgXPoHUTIGXv570Fj7kj45MHgKu7RQYojKvdrDcB2mV2j+71keBL1pVVmRV
vNC1x+x1PgMU6Bx/wXb8uCe6JREh+HWIAgclid/vH3mOHgSU5O0bpUDp5v7o
zKeu6LHKRVb9cUzP24K0mGhn5LjYfzVkggKtE3sewQ4n5H0HN7s2jdXjCTW5
DGsHFMgbqnR8jgIJnT55Xu3/odBqyUDzBQqs+WAqJa5uj+Kl1Q+0rVDgPRXd
Inv0GKUNdjty/KXAjYWvNS84rVB2qGWt4iqm73v1ynM2D1HF6mvTlE0K1OSL
uTwXbYZqM3fl9OGEoUiBkvYvM1PUaJqzRWQThlFbziJZF41Rp+A5bbVtwnB/
avTo6yR91NfQ+dGFUxjmnNp733rhFqK7m09lE4ThuMnoJ053XTR54u+ZcW5h
OFdTOWSupYMWxv39dvIKw7rHvL0HMrTR36idtKtEYdjpkU4/fkETbWlnSPnw
CUP2bELgX/3ziIMNPCkTEIZXgo7Pkn8BxFPUWj0vKAz/8+p4NiV0Ggk8MiXt
ExaGKjd46n3VjqOdu5aN7u4UhvLl12PXfu1DEh0vM4PFhGHh4Z63Dq9FkLSv
0EbtLmFoc6VANOo1Bzqu+Oni2h5hOEkn2H7nk6j4HwjXU+E=
         "]]},
       Annotation[#, "Charting`Private`Tag$876369#1"]& ]}, {}, {}},
    AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
    Axes->{True, True},
    AxesLabel->{None, None},
    AxesOrigin->{0, 0},
    DisplayFunction->Identity,
    Frame->{{False, False}, {False, False}},
    FrameLabel->{{None, None}, {None, None}},
    FrameTicks->{{Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}},
    GridLines->{None, None},
    GridLinesStyle->Directive[
      GrayLevel[0.5, 0.4]],
    ImagePadding->All,
    Method->{
     "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
      AbsolutePointSize[6], "ScalingFunctions" -> None, 
      "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& ), "CopiedValueFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& )}},
    PlotRange->
     NCache[{{0, 2 Pi}, {-0.7463487857271525, 0.3734843906304168}}, {{
       0, 6.283185307179586}, {-0.7463487857271525, 0.3734843906304168}}],
    PlotRangeClipping->True,
    PlotRangePadding->{{
       Scaled[0.02], 
       Scaled[0.02]}, {
       Scaled[0.05], 
       Scaled[0.05]}},
    Ticks->{Automatic, Automatic}], ",", 
   GraphicsBox[{{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJw9tnk81N/7/y87JWYwM5ZQslSUrVC4ThJZkhIqFFK2rIVCEqVdluy7SNJi
X0JHskcpEtnNmHoVIjvF9/m+3X6/z/zxnNv9ds7znOs8rsd1Pc9mB89j51hZ
WFhuEo///bNh3g/03dd1WP6/n8x/b/ZJxJyB/5+XC3Y59kZ7/h+PCTkvlkVf
+z++L9LM97Iq7P/4ibxNOsvLW4DVpV49rvuoU6s+pWSecQ+EIzT87r1n6HzT
v1GXE/0QTkyaXuC9+0dnefuG7UF1UdBh0b83zO2vzi2Vbp+cqhgodpw9p8/N
DkJ7M6o+lMQCz7uVLox5IXO/K/vii3gI9H3M93wXCXYaqh3enJsIXWkLT8R8
haHKbDXWKCMZdjemv3qrLAqHTjQPXkxMhXXrWtYrOUnClzPRcqnR6UCdvf8u
8N4WKFTdqUFvyICr4jE/zDNlwF89pt+/LhPirv2p/j4nDzr7FkM24Cy46im0
8O+OArCDrUxm1WOorOS03blLCVp161p2V2SD2PhEC0NTBaL05TxaS3LgYUln
zY1zaiBhOl028yIXDGZin/Fs1gTGUUvr28+ewm6/PJGiwn2Qb1G1Jp6bB+W1
7sutjjqgYXvzkH5GPkQ9EXz/Q0QXVu1+jn9LeQ4H3DK3oqcHoMHxSJRn4guQ
unTsZeShg3Dsgsi3hOhXwO1YGdaPD8EGDcBxBgUQub5+dU+8EUxVepNGGgqg
oL6Oy/C6CZTVdJf61hVCJ8edrfIGZpAIPNy1qAjm1qnnfzpxFILe7jvFi4tA
NuRau9qlY3CgIWM1raoY2JvE6m80HwdZg06z/zRLQObB5o/p6yyBt4XjsWpF
CQSaV4926FlBR5urQXNJKRw4tSu647+TYNulFjn9ohzwy8rvtQVnAFk6je5T
qIAjYU5fw9LsYGtPolr4swoYvuHvt8HVHn72rfaI5laCV4pC0hXSWfCnN2/W
y6gCsX2nafNZ5+HUuZWLERLVEPqW74hfrBPofFds7EmphhSX431FD5yB/Ve0
q3tiDTg+lbx776ErRP2xKY6LxhCUOn5+tNcDHggnadXy1YLe65ciXLxe8Lnw
puTuY7VwUHQp9cUzL7D+aTMm0VcLNmIPbjbOe4OH9XrvP+Nv4XT7Wlvf60tQ
vDBv7qRcB+nVh+nbfX1hMWZ0T79vHajc0MzQVvWD622VfxtW66DV4ufnxtf+
8Ejb+XaiQD28zvTdJ7gYAL295m4bLerBTyZkYvunQJDwA9OwxHpw0Tz0ds/L
IMh9SRFy39IA3x0ddXUvBcNryYZ0pNYIUjIpWudo12Fo3Zby71bNcMCre1ZP
4yZsTeNLsklphqaY0kONDTfBee9S0KfhZhgsYBf0sAiHPz4dulUuLRDGqLHZ
F3gL2BnBHyICW2Fpy6cywZE7sL2xb2x3ehuI3nzsM0R9CJ4OTc3P6G3wjLXF
Zbr5IZSsFuVLyreDnrq0qrxOJGir3/XmKWwHxmY2ibXZSDB7qvGvv+4D7DpG
bXnmFQ37Gox4dS92AD1e7+tQRiwIdD35qB7TAarRe3RkOeKAMcoaq1jcAXrG
h70DXeMggqVKUmSmAw7sTNpmoREPo/t27J7y+QSOQoVdvcMJcLtovV2az2dg
4zuoi8+lgO1bJ5lH0Z/BU6djVa8vBVQ63v28U/QZ5NvthRlHU6F/ItDX989n
OFcKiV66abBr2/hdE59OUHV1UnGhZUB3elvpsncXzIkLfojdnAX5L+UDp6II
vjGxpnsuC0JqbiBmYRcE5ox/CcnLgm19Wu8/TXdBkZCk9ivVxxBEeTH81PsL
MEkp5v2m2SD94MF6K+9uiJ8PWEvNeAJeAab2xV49sFfy0u+8onww8tB6yh3Z
A5JRm6zZuZ6DtMP237aveuBt4/WgCevn0G3IeZV7sgdM6M5/3bhegLbImwRb
t15gK3BzEz33EngrFDu4zn+D1/6GxmmpBZAztwFsrPthvzKiyOUVQ/B/yzcL
A/qhrO/apr39xXBi4EcbZ1I/8CwfK4zhL4ENDfWnCnv6gVygqSnqXwI+j4L8
OC0HIHPJAPsYlgJSHX9RYDYIeXYcirBSBu4V5yJYpYdBh+fw2MPo10R/PSq1
Tm0YTGZe2pzsfA1TvVpFa3rDcI8u0+ktVAWXVgS7/54fBu5D21ci4qsgQKdu
08KzYTAaSyu/lFkN4fUSL36pjMCs6l1e/O4NZHR8fd+lOwpPVBSLLDLeQtcP
Q55cBwYoyS/6coQ3wu9E+YNhPgyQXOnXuvm8EXiNua6fCWVAtvWMqWRnI6AX
9YvULAZc8J8rzJRqgnxv+HF7hAE61/umvlU3QeiyWpOr3RiczwyxXV5qhl0b
pG7uPM0EB+fBx+U338PdXfMsZSd/wI0TLk6nLnSA947m2Nqccdh/pWRzWUg3
LJ91ezngMQX8ZSllxU79YEeauVztNwX6V9h8uoL6oeFNwIHk4CnYp0Gp+BXV
D5Eid3tOREwB/Xux87+qfpD7+JS168UUkJdJ//L4B+D4XqZl6/gUsDlwSPCV
DsBLAfu1MtdpWEl/wSq3MAhCb360xPpMQ5OM/fpB3iG44ub16FLANFgcDO25
JTEEeo3XtqncnYbfnNNGxQeHoC8wzfxF3jQ4WZ+1Ox8zBNzf+3If/5gm/DLJ
SqYMg0ONpVnk+T/QJ9buKpY3DCplNLKI5x8YuRn6qaloGFhffevM9P8DQHax
da0ehqyM01bFt//AjpWrGVEfh2E07Pzp7md/gL9g/9L5OWI9I78L4r//wLFp
bCUBI+DwNfb2U/8Z4Es+5G/aOkLUn5WRcsgMkD/sL07/PAKsLSIbXt+egTeP
fAXHvo1A1uvUh+8TZ0C06DyHya8RGE3NiZusmgHBiDHhO+tHwcGxNFuNZRZ+
7JqwLzYchbPTXbW1t2dheUVo2Ll6FOwCIq0qImfht8nvG4V1o2DDZjL5KmEW
PG442Ew2j4Kl8Dux9KezcFySJUjvyygc0iz0u9o8C9PVMV72E6Ow8/oDhb3c
c7AvZ27NQowO23kN3ykLzAH7Fh+7ASk6yMWwn9pGmwOe73lnTsjSQSon8BZN
fg48Rcap25TpINjiMjpvMAd/n5wJU9Onw6KAfkLxrTlITbskPe5Oh9nEdbvy
H87BupqFyDwfOkxvqWnIip8D76A08VP+dPippjYTlTsHD6/89HoUQoeBE1tM
vZrm4NqPuQPu0XSoy1hlU+Sah+qQ6rt5xXR4s+110lb+eVis5IiOLadDVZGv
sjh1HlLauYwvVdGhpH789Hq5edDon6rkeUeH3B/fKv/TnwemBhcz5hMd7iuV
ez4Jn4e6SdFoqXE6tKp8ZfpGzIOBslWuzG86cO9esD0YNw+vW/+lbfpDhxua
6iaMnHlgU30mzVygQ6BuufyWhnkQd63i6WVlQIXe1/TptnmYUF2puMTBgHn9
BcrbrnlQa/olz8bNAB9jdQ47BhHvseB9ixsY4Hq8fCSVbQES+2TrdSkMeGr5
9YT7+gVwyZEgXaIxgHli4aOW4AJEcgmwJogywN5W/U3flgXgmWN9Vy7BgLQz
Vrvzty/Azerm+8VSDOi3938eoLIAFSkCk5lbGHDifHmSiC4xf+CBm6UsA+Kc
v5L+M1yAbXe3hYrLM6DLdeF2xdEFqJsSlPmyjQFmnur+VvYL4NAXwiGhyACD
y+XmsSHE+t4j43yqDLgR8LXV8fYCXBG+t9lZjQF1QQv71SIXoOzI+MeS3f/r
K+pKnekLoLyyIVpGgwFBYVa5WbkLcHirWJahJgNe3/SX8Hm1AMcsg7fa7WXA
nrvlG0h4AYzw2yuOWgy4dP9r2HDjAgTPxRge02ZAUcTC8qsPC9Cqln5TRYcB
U5FUn2vdhD53HypwAAMUY9T/Mx1cgNNeMtqtBLvGWtlJMBegqH132XVE6Bnv
/3VigtBvPi9ux35Cz8R405q5Bdj9/NpAE8FbU8ob7v9bgBrxljuWugxwSPuq
ZcOxCDWjz5O+EpyRsVC8g28RytnP8xsfYMBgFnXHitAiBMwe+F5AsHiOelar
+CJsi3i2mUePAadyrUSSti7CESr5jTnBCXn+kS4Ki/BCm6s6kuDu/HguTbVF
6BXeIVZLsODL8mBurUVoEFHtGyHYrODr3NcDi7BV7r+VOYIjihYu5BovQvjV
HVf/EtxWQmX4mS/Cuerx0/8b5y1Xt9a3XgT+EdXMYYIPVVp9Fj67CMZ64npv
CA6v8jccc12EH/FfDCIIrq+Jry3xWYQdjMfPzQhmrS1XvxGwCDw3h3w4CH7G
R37yL3QRmhp7LmQT5ztmfUHI/94ifAOfM6oELz9tDJ2KWQTz7+7GRYQ+j+el
pl1SFiG5PU95M8EmeoFn6NmLwE6VEgwh9J2N+tJu82IRts8zpz4Q+dBTvPvs
yJtFiDpr/ESLyNd4AIPW0rgIzU6ZIdZEfh8169zS/bgId4YrbFyJ/I+dnTm7
e5hYP7xD6sQ+Qo/Cw59f/liEfYnOG9T/55+1XCQ/TcTXN83CSfgrPMlmkxjr
EpgNa/Jd2sOAnT/K7j3iXYK97J7bBAl/du8mLfMJLgF1fahVFuFnuU8N3Sxb
l2Cf+pmZSCUGdEhIHQxQIObLXz0/tZMBly8EFM+oLYFriNQ8EPXQzLUrknlw
CWL4P/i9IurFy+LO6hnTJQg6xO3xUY4BtMf0C72WS3Dwq1zEkAwDXHQSDNuc
lsAAZVt2bCby5buOrejOEiiGnXefJeq7qM7ae0c08X77EH8cUf+nBMqGspOW
4OFL8SU5IQbk57tWx+cvAS93Z4E4P6HvaKfv1fYleE4uUDZiJ/TdtZMx/2UJ
muo5wz3WMSDl6u1jXoNLQPtVczBslQ4TNO1dZ38vAdeLOLHgRTpEHHny4xBp
GQr8Kksqf9FhTyqLVZ3IMkgsSi3Y/6DD4M9TDfu2LEOniXzHMoMOu8L5s3aq
LkNLcq7t6gAdOqovWwtZLMO21fVRhz/SgbTd+MNQwjLY6nxgJb+iA5MrndSc
sQwcvPWvBPLp8Hrsz/GCp8sQ9fPiYdZcOjhmJvWFVCzDjdsGF9+k0aGC9vPH
5p5l0JfFlOQHdDjNdYfNkboCs949Lr+c6aA61q9vIrkCmQoim8UdiX78Tumu
mtwKxPBFqKAzdCgK7hHgUF+BEPXfIQ4WdOCYl5N4YrkCsbHuiQL76fCc0ajx
I3YFTBc77EUodLheJxrUkboCnnn2umECdLDM8MAVOStw3Sh4cZiXDiw2FP07
pStQs+E25+W1UTDvPGe+/QsxP6YiVvjHKCy/Zfe4IPQXLmiLClPKie9nuu7j
qei/4HD879ULJqNwq7ljIDT5L1xbVQpp1RuFxunTNOHsv1C8zzVGTHsUDuoF
PtAs/QspT3g2xyuOgu5/xf6hX/9CpR+f0RzfKGipypgIif8DPtWgxJNtI6DU
xDWnnvMPjhqICGxFI0D93aZ/rXwVvtBrw9m5hiEhCtv/h1fhITVJJXppCGhq
RUHmzasg3M7xlGeM4CvxRXK9qxDGcZG1vmoIRNkcJDqWV8HUseO3i/MQSFAX
5qR01uB8fERp3JtBkIXNOXX1azDwJMrE6swAJDW9XsfmwoKCHS+/K7nyDbQ2
7V+TclyHTrXTck0qu2DHsADj6iVWFKrYc92yrg1ePCU7l/qwIbpCWnVdyFsI
GrTt3BHOjlhfJk01WDyDWQPt1HN32VGP3VUxcnkeXCgUd0qPYEce7qn1jbQ8
sL7Rt0yOZ0c08cebXQdyYe/2k1uWc9lRhLv40xS3HFjyPe7d3MKOQp4zjrUk
Z4AfnzH/OT4OFHR7RGx85wPw1tIwTnvEgabaJoOcFR9hFo9H1x4mcCDfgDVJ
m62xOCJ9qjgkhQPtYqowUkTj8DPWPPGz2Rwo9tKCEBdHAh5tFpmUK+FAo2vN
Z299SsbHLVYiCzs50Jj2YUp9dibe4/Gmu57Mibwf+Bft6X+Kl9P1HH5FciK5
pP6ozB3FWGdBO1o5lhPxLOQkuJ0txtdN1ev8EzmRvv5NzvXJxZj737Yt7Fmc
qJ83juM9TwkWPsU/Kl7Miah7shQtx0qwktA3e9MvnCg1a++dj4/K8LnbnvaF
IlzofpbdXtrXStzhlWTnn8WFJiwsN0YtYqxpsGV1IJMbOW5hKhxBDbjloOZT
lMKDFo0OCBw82IpzChM1R/N5kZCPlm3Cxg4sLK3gejpvPVo6mX/l3/ZO/Hkd
d6hI9QYUfJrj62n+r7iU++KTjhY+tJLHPeht3YcVhAVcHRo2IkrsLtyYPIR/
R8wLqXbzI3GTDRfZhEbx+JREneI3fqTwK33PTvFR/POYgaf8ID86ek/2qsXW
UcygJLRsYvIj3sk39olqo/hbmmYw9zw/+hnps2vu+Ch2sHozM8krgFY+V3L/
eTSKDzuts63fLoCW8oYa80h0bJIQW/xWQQA5x6ozTUTo2KhlGy/eKYCiFo6t
/JSi40M7jpVXqgggepLYF6oSHev+ziK92iuAvstu/aViSsca/nqNiUYCSPnk
7qazd+hY+tatXZ6uAujAydrBkyt07NOtFuxwQQCd3V26IMXKwLUyo20WHgJo
WDT8zQg3A5+u13LV8hFAVgPy8mYUBk5a9yebJ0AA3Qo/qTikxMDkIBvR7DsC
qFkwZlO1IwOzeitz9D4VQFIyk7HSTQx8tHbQvO2ZAEoPli/mbWfgDP77Wfi5
AJL/4y328zMD67xkwpMCATRV7JcVNcjAAb+SAy5WCKAq43d2oXMM/Occ5xRf
swAKuX/ObWTzGB491fdN97sAehFA+3T40hjeoG2Va/efAHrAl/0q5soY3i3Z
eTH4lwDauii6/CF4DN+mt254/VsAnXJry5S5M4Z3XXito7wogCQ/OP84mjqG
g64mZknykNAve3FG8rsx/OSMsKf2ehKKtGwaUWkewx/3R+2z5iMh3RBKxJu2
MbyF805XHImEzCYuNBR/GcPNEZc5N4qSkEqsuJ/Q9zEsnGHl+nc7CbFIeB/J
4WJindDOPaKKJLQnM50vdD0TOzseYdPYRUK81/e1H+Nn4mp5g5SLqiR0MPCj
XAeFiR0K93z4uY+E9r+a3qUpw8T3Y4qSuHVI6IRu/GCJPBOX+u50kkUk1B3/
SmerAhNz75VlcdAjIaW01PEBFSZ+WSes0mtCIvqT3VwVMHFPdtTqvClxvlxe
mS5dJl53i++90FESwmFjr4cOMvFxYw5HMwsSunaQZa7DmIlXOmdjm2xJyL98
LE/UiokPjXUulXqQ0I2XfVEH3Jj46o6IwFdeJNT0pVrJyZ2Ji7wP/XvqQ0JG
Df/sr3kysdhqNUuKHwnVPXB0eHCRiScpT7hCg0loaqVaUTKQiaVt7e4EhpBQ
zcSWImYQE594LLreN5SESvdvSsgKZuK3ux5udA4nIcObj70XrjPxo0OXhU0j
SCiJ16da6DYTt0SoxBtEkhCPRejL4DtMvNo1TtsfTehrFv6l/y4TO9nbi6vF
kZDiWZsXvg+YOCVXLE0xgRgndZoXRDBxx8QXSbkkEjpeGS469JCJ9wYYSYum
kdBcz/NZWjQTe9ay5whmkFDAsOmqdAwTZ3NiWb4sEpre9EBI+hET88Wobmd5
QkLCfqvm/2KZWLd3In8pl4Q+d09c7IljYj/Jp4ozeUT8Qn8e5MYz8XC+uDLz
BQm9F/VN25zIxMJ/uouGXpHQlrs699sINtKIUustJCENrSIHtyQmDg42Lvtc
TEIjzlmbVgkurufQaCslobPmMxU3kpn4B29tZUM5CX2wTFZaI1j8aMA+XElC
K9NZIe4pTHw0Xq2mooqEYju35nwkOHxgUqeohohv7HCSTCoTV0nn1eZjEtpx
ytbGk+Apl7O6OW9JqL0k+/tzgrcWbKpPe0dCrD23NAYJPjn/9WBCAwmd4b9q
zJbGxBFa0U1RTSTUKcq3SZzgulATw3stJBSj2ZknT/BCM+f7G+9J6Jva1cn/
sQL/W5PgdhIKd33//X/z7S0CP/h/JKG7XyOj2QmOS95t5v2J8PuW4l9DxH6t
I78/uXYS8Wdvn3pJ8JrcM3PHLyS05sua6k2wmofjF9uvhB6+J37LEexcImFl
1UtCb47ZMj4R501d7ukx6yOhkoTjVzwJ/oxiThkNkNBi6aN8FoK5bh3uPzBE
Qrc5b12+Sein1c51WnuE8GN/7sA/Qm8vwbqhPXQinrvhPS4E55wMslcaI6E2
54bzLUS+vqXvoW/7TkKfRkXCxQneyJxylP6PhMYCTdTOJjDxZZ9zzpQJEtrc
0q72nvDDbeNJp72/ScimJkyzl/BL/FZ/p9PTxPpvz8X2EH4q6759/skcCU2o
+BjmE/5reEU637pAQv32+1ODo5i463bSucklwo8a4XAgkoln9j533LNKQgZX
h50yCH+zCu12PMVCRrbNXj817jMxeaLmbDArGd29faGkjqgHpfSPDo2cZHRn
4R3t2S0m9mCbsbPkJ6NP7U/XmkKIeu4PtAsgkZFIc0TiT6L+7pey26UJkpEQ
6zfvVaI+850oZ5hUMlpVlepbukz4772mrZ8UGfk5Dt8ALyZezK6zSdpCRtbi
uRzjRD/gCja2ebOVjJhN8XV3iH4hq2RrzbmNjNZfGJ5OOM/EZx+FnIxTJiPJ
Gk8VThsm9nHnOflalYzCXOUytE4ycah+9InB3WTEsZRt6mjJxBmLj61k95JR
yQ1Dn1tmTDxo02RRpktGbdWPF8X1mHhit5nFNz0yqnfbHTmAmPjvxt7jq/pk
9J8q96WH2kQ/evvTXN+YjI4vaso37SH6jczGY93mZCQzZn4tnOinzqtxR5ct
yMiNt7jk01Ym9v8qeVTiBBk9jC08tHEz4c87ymbnbYj9rjfFW4swcefEcdN5
RzJCHAuyBjxMbFqWbETxI6NHi9nc55hj2HnGYU34Mhk1uZnYUUbHcJjS9hLh
ADJi1LHIVw2M4YpnFZuEg8moPwFRhrqI70n6lynBcDKqXrrSG1k3hudv8ceT
4siomy0lMT15DAs0dBuTEsio5mHSlztxY3gHaxoLKYmMJooaGpyixrBdkIKr
QBoZPRC+0r0aPoZbvQy1+Z+Q0VGFF4eSvMdw2skw+oYyMjrGk6N4T4/YP84o
YUMFGenqlm0r1hnDnZ2kwxtekxHtaojKR40xzG2aUbb+DRm9sNoSRlcYwz66
NXd4G8kosVq686rQGNbfsbCLu5uMTA+R2a6NMPDkX9cgtnkyOtelY8D0YeCT
J4y3LS2QkaxlI1udGwPXF+/onlwio5ADybIxxP0g0XV817d/ZHRp8AkWtmRg
vR53egGHIGqqOySbqsHAycVeRqeFBVHmZJZaNXE/OeTqR6tUE0RsB7jkAgLo
OPPr9RKPi4KIeaB/TsJmFDMblIYifAXRxoqE7o9HR7FCyRDPK39BpDHbquJv
MIrLIrXP/A4URIZveZLyVEZxm+ESj9cNQWQ1+PdDGvcoXqzyOuMdK4hed3XZ
9xSNYNOMM7yXygXRZxFnxgLLCF521ra7siKIKGZ2C/L8Q3jTvyhr7wAhpKFB
bf+xrRczDe61crMLo/sXqxoU8j/jd7JtIHFPGFl+c5MtudyCjSKGjCgPhBGt
M/mc7ZEW3DH3x2LjQ2F0xz7CXkq2BQ/Wi1xYjRZGWjw2qrOdzXjRwSlhMEkY
1ZXtefZVqRkrpLNOpz0TRq3KI54FvxtxHFUzS7JVGGW3v1ctCK3H4tdMXlDb
hNH483ehfTb1+DHzTAX/B2F0+5XeThn1elxQGv5h7ZMwShEcERcaf4ffm3ct
D/UKo+Tj9d1nrd5hlkhP84z/hNGBDf9qx1TrsCv3E/bNPBS0LS3YToujFr9n
vf3gwHoKOuLZarPtLcYK/1wo5/koyNj2w+NtlzH+Pa0on0+ioAfD9RKi39/g
S32lxrtFKWhtjJv7bXMNvvqyIdpwBwUdD/HqsUiuwpHHx6R8DlNQu6RZ1vfg
cjxt2pT36AgF3dJycfoN5fiYYZ5K+VEKMqvd+8BkXTkW1nHX+2tBQdZdTD+u
m2U4SW7OKfw0BUn1RcyERpbi7GX2V0meFPTtLL9CbnkxLs/Yqv0uioK4vuzR
2W5QgLeu6QW+jaGgGgPf9CBKAY60PVeJYynItSw9cU/bK+ws+mR3dSIFNTHC
ZQM3vMKij+R2lmZSUMuKbuWORy9w0K3tkrmFFCSyaVpUrvYZ/jFmZJtTTEHl
78DguPMzbKHnlvy4lILc22WfipCe4Z3rnlMzKimowOHfmwuOeXgoQFEg8S0F
6escWlAmP8W6Hkrr7n2ioIBZ50+b7+bgl21mcKeTgrSz/aca9+VgsR3eV299
oSC9Zy/CWCez8ez3wqWwXgpKfqTspHU8G+fYq/4JGiHWC/4XWCb3GHNb7KF7
TFPQs1vbDJJ+ZuAP+7QajpGoKLtItUC8JhHfrfBd0RekEveDBqFIg0RssPuV
8j5hKiqP7ztA6UzAtTu3pEmLUJGzqVjSiYl4XLyZ239WioqO3OPreLwzDidw
dcnHKVHR5cOh5Dp6NLa4wXfmrgoVZc34qghfi8bkdQaxwWpUtDWR76KXWDS+
t1y5zkmDihLy5qy1rKJw8ET6N3VERZ0mh6fkyx7is51u93uPUFHJL8Z3Odm7
WMosp679KBXN1Iq7Uq7ewf1tg4tvzanowZcJ65Uvt7Fl49Fzz6yoKHePdGjm
g1v4UKW6TuAZKjqRrB6nIHQT70xnmxL3pKKKydyNaUEh+KeYtizJm4rcNp3v
3cm4hnMT/Gw4LlKRi5ig8efPwVgq+r/mCT8qUtiyJitWFYQFb37MfBNMRZEt
r37tL7+MF12Tze0iqEg0x/uk6ilv/DDgMde1SCqa9nQOK3XzwrJ386tSo6no
n7/zhmBWT2yRVyXdF0dF/NzKc5lebrjoe9+MRToVNUb9MUzMd8SGC/TcS5lU
FPrXo/k931k8zDluHfOYihyN9gW/9LfH/LJ/33XkUlH4DpuS3HobfMFR/JFx
ARXFq/oY3JU+itkubT3kWkRFyRcLvoXWmuCkMIW/t0uoyKBc5+FWKwPckqXl
2FhBRSEBjxvcyjSx7IiN2v63VGTltzpkJb4DaqYcf5x5R0XbXoVutO/cA8dZ
3FOCG6joQPOguecIQJjkVfbqFioyriiaGH5+CIZtU7v2fKYidVL587fex8Hf
Pee2RRcVSXC/RkVWlrDx6gutS91UtFdgpqdQ6QRop9RkF34j9qd8Hzeptoau
/IYTHf2EXiVzW6U8beFCVfuG34NU5N1t1ustfAaSvg1cUqAT573yqjFEwh6U
f47JG48Rfj29Q1pryB6alyb6Xb5T0fOSZr/fiQ6wQFvVy/1FRQy+A+f8/52F
41YSwltmCX+ROY5mBp+HX+dlW9A8FQ0yn48XCDhBqN/Oq2cWCX+WG75TSHWC
wlidsZS/VLT66tLx7BxnOJSjn1i1SkWPJLe/r9nsAsMlpoe/sdCQuXj3bHCC
C2zsOl1K5aAhBRMvDzU/V8ihn3fZw0VDpXR5pv6AK2jNeGyy4KGh0EYJO21w
gy5W/08X19PQBdssAbkUN3AjX7sZzUdDxgKlGoKzbsC65ZZmIT8NkWpGZEgG
FyBJ+eHERxINbVAWVlGKvQDNZukWfBQa2u+vY8Iv7Q6u9WNvuGjEuL0vy4fT
7sCnoSDPKkpDiu9C/4bFuUNBvk/UXzEaGtCf+KHd6g7HJCuX5zfRUPr4uB/3
sjvMRq+dnZakodprpnZTMh4Qx6nf/mszDUn2WTqxmnqARsD9PUxpGnLXfGJp
7OMBfROf04dliHi9xNb1RHtAsL0IT58cDb1cK9tf+soDpL6c8fmyjYaEtwWv
/G72gLpDT/o+7iDG/9iKxA56gGP1uF6rIjEusDGicMoDOJVUX9bvoiGj4hnz
Y2sekPf4ChUr0xDHaXb3G7yeYEKtDalUpaE7eU9Hj5M9YfIu58/i3TSk2pPW
0EHxhMhVE/OX6jQk7jugyEHzBFWfmOqnmkS8k75714Q9oXusV+bxPhpaN+9H
ahPwhMsnpR6matPQ1cUPQ5e5PUG0/fxiPNDQcuDOXvI/D6hBL+yj99MQ+fxP
9dxJD7ArmWm9f4CGeHcPah0Y8ABW+b1qtw7SULiq6a5Z4nzZySGp1w1oKD+r
xaap0AMM+Js4gwxpyKFig0hrvAf8DOXz8jOmoZDC0RqeQA+4P2/e63WY2K+T
Kz/G2gN2uibpuh2hIX1Yr++r6QEdA8P5547SkFyv0nKVkAcIN7gHn7KgoWKW
Br+X79yhQqPk+3ErGhKwVFX2jXeHU8+XzY6cpKEHH8N/dDu7Q1rMLWk9WxpS
Pz/rmsnuDojrw32dMzQkxbJkJvPxAtADhOY17GlotSvAJSD+Asg7ZDYrnqMh
ZseT5A+bL0Drl+/K8k5Evnos537Q3cDdcGfyFhca+hJ1bN36bDcoUqpyp7rT
kK/mp8EUCTfQWusir/OloYtr9KP+qy4w6CMWtOJHQ3/MZ96aFrnAdab92Nxl
GtIOG36c7OgCje2T5T+DaMjtc9NgW70zHEnhtu26QUPSKWa1/L5O4KCplZv7
iIZMHXNpapmO4DLRPtcfR0PX4zIne9UdwSvzjB45kYbm7qYce9Z+Fq7xhI0E
pdKQHrusGm3BAVJ6W8SOPSHis1mu0T9gD91XrCL/ltNQ76/Tqi2PrWFA8b9B
5ddEfRpWtEc2nwLGSICiUzUNDQr3p4pMnoQZw7TWT7VEPg5PNPBonQAB0TH2
3BYa+qns+Dz9+3Eweu1z5Wgf4V9nVelencNw1IO9OXyAhvYBZbW53hhObImj
VA/R0Nrq0B2+o0bgdLeyWJZBQ5V9NtubrhrAjVMskyvjNPTamrbtpJIu3NsY
raX8m4aOvu7KNZ5GEF0nfe/8NA2xxJVy//IFyNxuIP9pjoaS3iRWaXfuhTfL
DxyerBJ6kPksN11RgYaXkoV9LCLomWxTPWe8ErQ5FK4JsImgX2fv2O2MU4Qu
ygFTfQ4R5FD2Xahlx3boa+1KCeQSQRNZXWw8U7IwGnz+VwGPCBo1v8oiwpCG
nyqLmsz1IkjhW6OZNM9mWEwS+2omIIKkLh+3TvpChTXTFzLhZBGUmiLdFmQu
CJxscKlKSAQ1XXmyYfeVjcBX3lE3RRFBQO2zuSHLBUJuDiRZERE06GlkGlrM
AmKSs2esxURQYYP+jgjDeZ0tnTdfRm4SQe9WYuUyV/7T2X6L+q9BUgRt+u5b
1ajaraO8L894ZbMIemQYnkhiTmv/P/3JXac=
         "]]},
       Annotation[#, "Charting`Private`Tag$876413#1"]& ]}, {}, {}},
    AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
    Axes->{True, True},
    AxesLabel->{None, None},
    AxesOrigin->{0, 0},
    DisplayFunction->Identity,
    Frame->{{False, False}, {False, False}},
    FrameLabel->{{None, None}, {None, None}},
    FrameTicks->{{Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
       Charting`ScaledFrameTicks[{Identity, Identity}]}},
    GridLines->{None, None},
    GridLinesStyle->Directive[
      GrayLevel[0.5, 0.4]],
    ImagePadding->All,
    Method->{
     "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
      AbsolutePointSize[6], "ScalingFunctions" -> None, 
      "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& ), "CopiedValueFunction" -> ({
          (Identity[#]& )[
           Part[#, 1]], 
          (Identity[#]& )[
           Part[#, 2]]}& )}},
    PlotRange->
     NCache[{{0, 2 Pi}, {-0.7068066270818705, 0.6334211773609446}}, {{
       0, 6.283185307179586}, {-0.7068066270818705, 0.6334211773609446}}],
    PlotRangeClipping->True,
    PlotRangePadding->{{
       Scaled[0.02], 
       Scaled[0.02]}, {
       Scaled[0.05], 
       Scaled[0.05]}},
    Ticks->{Automatic, Automatic}]}], "}"}]], "Output",
 CellChangeTimes->{{3.750625593998267*^9, 3.750625629545697*^9}, 
   3.750625690394271*^9, 3.7506266718438*^9, 3.750626954827797*^9, 
   3.7506273496749897`*^9, 3.7506274025513573`*^9, 3.750627504890024*^9, 
   3.750631341823943*^9, 3.750631504093338*^9, 3.7506316010803137`*^9, 
   3.7506852520825567`*^9, 3.750688656068368*^9, {3.750688726822229*^9, 
   3.750688740871169*^9}},
 CellLabel->"Out[87]=",ExpressionUUID->"74bfe69a-2993-4aa7-8d76-5daa20dba809"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 Cell[BoxData[
  FormBox["K", TraditionalForm]],ExpressionUUID->
  "d3c787e6-91e6-4a71-834b-8d49760493f4"],
 " matrix"
}], "Section",
 CellChangeTimes->{{3.725814385678084*^9, 3.7258145068856297`*^9}, {
  3.735058611652068*^9, 3.735058614481969*^9}, {3.73505865613306*^9, 
  3.735058656731237*^9}, {3.7396232985548763`*^9, 3.739623300142106*^9}, {
  3.741693899945702*^9, 
  3.741693902268856*^9}},ExpressionUUID->"9c8b182e-87c3-4ec9-8ce3-\
6172bbb3e292"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"K", "[", 
   RowBox[{"\[Theta]s_", ",", "\[Alpha]_", ",", "n_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"K", "[", 
    RowBox[{"\[Theta]s", ",", "\[Alpha]", ",", "n"}], "]"}], "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"NIntegrate", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"\[Psi]", "[", 
          RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
         RowBox[{"[", "i", "]"}], "]"}], 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"\[Psi]", "[", 
              RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
             RowBox[{"[", "j", "]"}], "]"}], ",", 
            RowBox[{"{", 
             RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "-", 
             SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
           RowBox[{
            RowBox[{"\[Psi]", "[", 
             RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
            RowBox[{"[", "j", "]"}], "]"}]}]}], ")"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Theta]", ",", "0", ",", 
         RowBox[{"2", "\[Pi]"}]}], "}"}], ",", 
       RowBox[{"AccuracyGoal", "\[Rule]", "10"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], 
    "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"K", "[", 
   RowBox[{"\[Theta]s", ",", "\[Alpha]", ",", "n"}], "]"}], "//", 
  "MatrixForm"}]}], "Input",
 CellChangeTimes->{{3.729015273124943*^9, 3.729015279775338*^9}, 
   3.7290153229023523`*^9, 3.729873251318688*^9, {3.729873425293119*^9, 
   3.729873470035809*^9}, {3.729873823245585*^9, 3.7298738569687157`*^9}, {
   3.7298739246431026`*^9, 3.729873951659837*^9}, {3.7298745255661287`*^9, 
   3.729874527372711*^9}, 3.729874622767412*^9, {3.7298756429757633`*^9, 
   3.729875656015667*^9}, {3.7298757028599377`*^9, 3.729875726503663*^9}, {
   3.729876101035668*^9, 3.729876139373316*^9}, {3.7304636445709333`*^9, 
   3.7304636465048723`*^9}, {3.7304637645426273`*^9, 
   3.7304637661030703`*^9}, {3.730464375055768*^9, 3.730464377343965*^9}, {
   3.730751289407498*^9, 3.7307513039142523`*^9}, {3.730751358202427*^9, 
   3.730751360392844*^9}, {3.730805050140264*^9, 3.730805085878983*^9}, {
   3.7309114818139668`*^9, 3.730911548056965*^9}, {3.731679918345953*^9, 
   3.7316799191841183`*^9}, {3.733064702144858*^9, 3.733064718044208*^9}, {
   3.7330648718772717`*^9, 3.733064909969494*^9}, {3.7330651362026043`*^9, 
   3.7330651427427263`*^9}, {3.733065719317165*^9, 3.733065719720889*^9}, {
   3.734627777334146*^9, 3.7346278097029257`*^9}, {3.735600097530471*^9, 
   3.7356001145078487`*^9}, {3.735600254068118*^9, 3.7356002900915327`*^9}, {
   3.735600476410337*^9, 3.735600595059184*^9}, {3.735600630150917*^9, 
   3.735600640215461*^9}, {3.735600778676319*^9, 3.7356007838339252`*^9}, {
   3.735600972552698*^9, 3.7356009784020643`*^9}, {3.7373946872590427`*^9, 
   3.73739469903185*^9}, {3.7396230385962563`*^9, 3.739623162891141*^9}, {
   3.739623508646983*^9, 3.739623513056185*^9}, {3.741657381572112*^9, 
   3.741657438474419*^9}, {3.741657481881212*^9, 3.741657491476029*^9}, {
   3.7416575342440977`*^9, 3.741657555517068*^9}, {3.741658480613*^9, 
   3.7416584813058577`*^9}, {3.741693910835019*^9, 3.741693913678966*^9}, {
   3.742049037856204*^9, 3.7420490616265917`*^9}, {3.742305832222436*^9, 
   3.742305911790238*^9}, {3.74231030508823*^9, 3.7423103083974257`*^9}, {
   3.7426395220734673`*^9, 3.7426395296170197`*^9}, {3.742689527983073*^9, 
   3.742689528028006*^9}, {3.742689588028388*^9, 3.7426895881957483`*^9}, {
   3.745687388610865*^9, 3.7456873904979353`*^9}, {3.7506301178697147`*^9, 
   3.7506301345963993`*^9}, {3.750630173455534*^9, 3.75063025869551*^9}, 
   3.750630447528161*^9, 3.7506305981671953`*^9, {3.7506309064349213`*^9, 
   3.750630926955102*^9}, {3.750631773523589*^9, 3.7506317740898314`*^9}, {
   3.750631893325206*^9, 3.750631896645576*^9}},
 CellLabel->"In[63]:=",ExpressionUUID->"2581c4b6-997e-4d60-9846-58c4feea8940"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1.8188779876130903`"}], "1.9349026422959088`", 
      "0.16462935152004432`"},
     {"1.9349026422959095`", 
      RowBox[{"-", "5.415791529368171`"}], "0.`"},
     {"0.16462935152004388`", "0.`", "0.`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.750630946778408*^9, 3.750631342067815*^9, 3.750631442975045*^9, 
   3.750631601998992*^9, 3.7506316741404753`*^9, 3.750631770908246*^9, {
   3.750631887633603*^9, 3.7506318973846483`*^9}, 3.75068283999972*^9, 
   3.750685286132913*^9},
 CellLabel->
  "Out[64]//MatrixForm=",ExpressionUUID->"e31e44fc-f50f-489e-bf6e-\
5c9310952d77"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"i", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"j", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"\[Psi]", "[", 
       RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
      RowBox[{"[", "i", "]"}], "]"}], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"\[Psi]", "[", 
           RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
          RowBox[{"[", "j", "]"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", 
          SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
        RowBox[{
         RowBox[{"\[Psi]", "[", 
          RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
         RowBox[{"[", "j", "]"}], "]"}]}]}], ")"}]}], "/.", 
    RowBox[{"\[Theta]", "\[Rule]", "\[Theta]\[Theta]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]\[Theta]", ",", "0.2", ",", "0.8"}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"NIntegrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"\[Psi]", "[", 
      RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
     RowBox[{"[", "i", "]"}], "]"}], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"\[Psi]", "[", 
          RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
         RowBox[{"[", "j", "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", 
         SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
       RowBox[{
        RowBox[{"\[Psi]", "[", 
         RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
        RowBox[{"[", "j", "]"}], "]"}]}]}], ")"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]", ",", "0.2", ",", "0.8"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.750683039979045*^9, 3.750683042608193*^9}, {
  3.750683080730702*^9, 3.750683132583828*^9}, {3.7506853018758287`*^9, 
  3.750685389268567*^9}},
 CellLabel->"In[80]:=",ExpressionUUID->"783c1bc2-4811-4229-b3f2-3a2b13ab8b38"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwV13k4VV8XB3ClKCJjJFOplEo0U/lSkkJRqCSJokmUEsrcgGSqlIQoPxRC
VIaK654riiRDKcK959x7Q5ShVPTu96/zfJ7zPOfsvfbea60928Vzx6GJIiIi
ChNERP7/rLGiHycm1mDk+KhHuOhl+KelOaal1mCe+oaiMrnLWDy8b3J2Rg2S
1XIPDs++jNg7LfbPCmpwP/2h+Hnjy9jdUz3aUlMDLwV/r5aAy2DCH0Dhdw1e
Sy5Lrhi9DFG25+toh1rcUVtY7TUSjvWGv7vDZr3G5mmT606KXkH98t0qvOdv
IO7ATnV+FwOZx2ZFkab1qOm1629wjsd8Ojf+ctVbMH/vGEyZdAOL5R6tCyho
wHNhudlVjVtIPr29a4/3OwRvto98apKEO/sUIk9uakS+i3zfVkEKxO4q2hXJ
v8d00cu6FfvTYD5vv41K03s45HhZm26+B2X9Opuy+00Y+i6r+d4nA5VBYWFB
7s0I7Wk95ZKciSZtn9YT81uw39ckPawxG2fsry5NHG7BP5uLE5/vyEF/de8z
5kkrevo6DH6/zMPPl15mnuEfoNnj7lxQno8jsvojc+w/4pJV7x9R9UIUZRi0
jiq2oaA20X6Z12NcrXArl2TaELJiPleGXwQj9c+Nxlmf4NjZy8y3eIKAiT59
+f6fIVc8clOq8Sm0jS11LLe0Y5J1joeaXQmyK47ekpDowPYzfuGuL0tRUp0r
MelTB1St4sVPbS5HuqxBnmbKF6i7aMxqL3yOcw4ihkW6nXhT7XOLknuJ454j
KWcSOnE6L/rhE9MK1GTE3jH/2Yk/rz081for4KGxXtXAuQsTgrMaX6RVopoS
X4LqLkheWsp9vIWFq1d9Wuy1u3HO5vny+A4Wbt/6OSEqphu6pWqv1/hUQZ5S
5TQOdGM4ujzCdbQKJurOs5Zac1H7MfWIsicbT0qLf2Q95SJI5a+dUxcbeWln
NBeVcLHOJSilt5uNjPCV23JLubDeLS7rz2Pjun1xdsFzLqyOi3vH89k49aPI
ubSKC3O/n3NK+tjQ1Smqf/2Wi+U1ykPto2xkJhY++MbnwnJRj3uYLIVbvo9c
VijzsHR1y0RdIwp3ss6Ktc7kYbnOF8c0UEj7gAd+s3ioFw39T96EwoM1Dd9f
qPNg2bhYfHgjhee/BoK3zOMhqWLt0IMtFLp9l6U6LeOBb+Cjy7WlsNiv+FO4
JQ+GhfPyQo9S0MsOCFq0jYdVzZpW749RWPlxk1b9dvL+i3/DHA8KRgYtR+V3
8hBWTIe+8KRgMzoymryHh8Jp9VGC0xR8/NYoP3bj4e3nEkYsiEKlX5ltezAP
mxxYG47GU7A0ZIwOhvKgG7meG3qNQstvmYU9YTwc7HgWevs6ha/n3MZ+XeIh
WTrmOSeBgkKg7H8KV3koPtakIp9EwS3U/adFIg9T9R6uOHePwrQr8rdLC3jI
9uD2XXpMIcHC6ILJYx4snD50BxZR0Jx25MSrIh7OhjXUnSkm87n6YkPLUx5M
TXMuOD+l4BRzpOf7cx6GmnFiYRmFgmsv1y2s5UHc2fFHOIuC/Z1jX2528yCS
2Pj6fAOFa81tGko8HrYc5XG2vaPQIL3VOYHmIeNXRYlGI4WtIQu6rgt4MIp+
e/HFewpwo7vjvvEwoezXo/4WCgv09jNXfvNwQ91YTrOdwp+qHX1BcjSiKDmj
pwIKa8Yql/yTp3GoqrzAXUjiuUr/RKAijRcefAWlrxQGsqb3n1emUb3yadap
Hgq8q28G/NRplE2V+qf6jUL9LrMhbx0asbcllxgNUrjbY/DHfQONj/VXqvL/
kvh6FBSIm9LQHKRTt45RWDGw4HDmJhpz8UefS1wzpNTMmNNodpP0kflH4cfY
cN6h7TS0WvNP753AgalMoYurI43dbVmijydzEB2/cKaoEw13H1phjRgHHxTS
3qbvpzHbuLCljNhjZuy6bhca4hKPj74U5yBhzokZB47QaLq84/WTqRwIVujU
OvnQyGCv0HaT4mDZ07Tg8bM08rN+qLQQnzeYuTrFj8YSC81SU2kOZI2m3G8/
T2PO4gujatM5MNzMBDiG0Sg5EepdIcPBhZp9K/5eoKEi0VgyT5aDeovmr0mX
aIxvmpwdQexqzd71KYJG58F3wVZyHETtSdd3iCXxVFeXouQ5aPk0UzAaR+P2
9qw4DQUONJ3iUhKvkXg0RpX7Ehe7BEt+TKDhMnH2JG1FDjqOOdG7kmk4CJNt
vGdwoB+okmj3gMYJ9e2ptcociMQN3fj4kMaNHb4fJGZy0HCvPn5fLo29U3Q7
txB71YRGHcqnkaBs6MAmzlfoCTrzhIZvVtH3ByocBGtT50ef0hi6NPtWF7G1
YapfQAmZn4KcpNIsDgb27/S+WE5jV5C/XSCx3sMytxssGjMvWJ8xUiX/f3HD
dSabxDOE5XuMuKHB0zmZopGXXWB2k9hrRMsh4xUNC8Ekkx5i46nj9gtraZjV
84/Lq3Ego/phZ+5rGlYBm1wMifNNoqyK62nsXO+ZdYE42NZtq0EDDe8jpb8y
ia3djTc/f0dj9HqgbC3xQNSQMdVEwz/CK0lCnYOK1Pr1W1poKJaqTVtIHFuY
ZVjXSuJRssHUjFjvg+OK5jYaEvsGpwYQi/Ss0t/zmcaaJ+MJCcQNYzK67e00
EjeHfM0jvivTo3PgC40PomsncYi9tChtupOGLb+H+4nYeFXq3CPdNAx/7YwY
IJbZ4je7j0vjjNnyIVENMp8TS2YNMzSqJHfP1CYODhFX9hOQ/fle7f0qYuvr
XQpjQhobzu212USsmVkmG9JDw2jju2s7iAdKbkhP7qPR6oPbTsQVbzwlI76R
+B9zO3iEOPbLlilSA2T9AsV+nCJ2/qE1Oe47me+szg3niJdOHp+gOEhj5ePH
9qHE40ofxm8N0ThbZbIknLhep/CP6ggNhfdbOVHEKeujft39SUPeO2Z2LPEJ
a7fhuaM0rhg1GMUTG7ka/8j6TSO3rU7tGrG0j0r/4r80KkbMX/z/fUf4UE/+
GA0Dy4maccR5SfWCFf9oPEwp3RBNHJiXRT8TYfC+bNn8SOJtlaHd6yYyyDy3
qPYCsXqT45cKUQaOFYd0A4m/Mas+m05mYOuTvsOH+MWozMdXYgyWX09f5UEc
Pa2n2XIKg34V7TYXYicNqrFhKgO24vd1u4h1l6W+tZVkYHK70mkrsWR0yIbR
aQzJC5YT1xILhK7FydIM5n3btliHmNpktmCDDAOvDr/PSsTpaQuSGFkGUYsD
5CYRB41JSF+RZ5ArNbX+G1lvxz19wUsVGbhRddIfiGfIFrr5KjOYOb1K/T/i
wePXP6qqMIj2su2L/P/+euVjWTmLQcLmTnKGOYgMXrtMQoPB2hvZnnrE7p/V
MvI0GdRsN7ecTmy6RkR55xwGb8XvP+wl+/1fP3ssaR6DsES77+nEny0yvYy1
GYQvlmo9R1yaGcHlLWDQsXST0U7i0/u31SxZzODVnd8xY+T8Cepbr79cxsCY
blttTEzplIofXMFg6PCwphRx+qU7/lNWMfgloRLxgZx3R6MDzjYGDC4kb24+
StyQ+3URFwxiB+/aBpB8kTu1LuWyCYPrP32XGBFHHnoku3gjg0eRgtAxkm9M
1U7/PG3GoGK1ffpZ4pKoMZaYFYnfH/rBHpK/EgRfVj7cxkDMcWi/NPFpU1bW
dmsGgdGLYiqVyPr9vRR9aydD+tIOBy3i9GMyDjoODMrNP6c0kfwYVP3jTf1e
Bkezn6kGEDtqNcN7H4M+kZ1yc4lnfEqcV+7MYLbNVVcPkl8jt2p9t3JnYKWY
d6OP5GfvhavCT55moHgwZPgeyfcTxx1tpH1I/OPmtK0kjm8MU3l4loFkpIwx
h9SHwnPvcnj+DHx6C6O7Sf34UXf83a4Qsl/NVy6dOo2DU6fuqRhFM1jWHGwk
PYV8b3Mtty2GwUlL6xnhpD7Fz/qeczaOgUrF4FER4kK2EQquM+gNux3RS+rd
9xltrnOTGMDgv+uFoiTflMnkSmQx+Jv65d7wOAWR2NU+mdkMtu9ja+8ijj3o
BNOHDGLeys99SuprvlTOu4A8Bj8ftF04Repx//7NIwNFDNr/c/nSMkrhxKRA
tFYySH/rWmY3RMHDqufdvU9kPMmVOjqkPwgtPN01t52sR1W5wWE+6Q+Vxgcy
OhikaPn13mcosLtkZLK6GMTV7V6nQlNQ8Vm1LYdPzotqbtBQFwVOSkht0SCD
4nMpB90/kf5oQKmKkuTD3+HHwk11FFbZpTWaSfHh91xp8643FCxKdbqrpfmY
X6s34fBr0o9cWD+hVpaPmynuFhdrKNQpu6JeiQ9fm+nLHlEU/E3yylq0+Djf
/HBmxXMK7+NNi/hr+Yi9KN7ZlENhxLjFIGI9H+Emk63TH1JQ7nd/qQM+Mv/W
zTjxgPRjFldqPTbw0SH2LFMki4JAtLFz0JyPkMqnb2VJvzd+er+UiD0fojID
49xbFHR2+bspe/GxtuGmnFkYBavJ03pLTvLhPmzT3h5Cwetx8sm93ny0Ti0Z
8Q6mUDy98nyKDx+uJ2IUkgIoGL8SvzY3gI9dMceGm85SsDO48VIvko95nxIU
fpL+N3hWvvKW+3xw7t+QG9xB+jeH+7z6DD5G4343etpQeHfrVr5tJh9sDUau
ZzuFOTNCzJ0f8JG7JVKm3ZJClcwO37P5fPAib4/lmVEQExtqzSjnI4Wb4y1h
SCFqYPXNic186I+f9UvRoJBIvVQsnyyAUii89Bg2hs7/Wm0oLkDpnn+aduQ+
Yb1C3+HZFAEOJ1Tb+pL7hti99JQiSQFEzlCmpR3kfhF8UTtXVoDPxdUmi1vZ
sFxrYZCiJgBTlhb5qpoNkYIWx+CVAmzNeKPXlcmG+52+e6aHBHBZ+udNhCsb
DXd013HdBIjV371D7QAba5I9m0IOCxBopy/5yImNKSnfJ708JkDPmbgZdXvY
eJA6fMjglADd9nejB7ez0Zs+pq0bJMC9hf8+CNeS8WRL5SglCtD6lxPmL8fG
+SdLir6+EeBg/yP+7NIqjH+X3Ze+RIjZdxOTz06ogoAtNNixVIgnql1c0XEW
Gm9WzpioLwRX3E8j+jcLmetPNhxYIURuXP+ZO4MsWEc0bJy9VojCkZ1Jd2kW
0jVjdO6aCzEjbcM9sxoWzLZP+5XsKsRSEcvwKTEs6M3hNVkdEuJTwMYFO6+w
oDJcVjDmJsSabSefJV1mof/2sWP7jgphxmrKmh/Mwi2mtkP1pBBHGsIorZMs
fA2IoJIChWjbqldXuoOFJusD6RbBQrxJTHXibWPhhZZB0J8QIdK5TLukBQvX
agRr9l4U4vbsoxm2G1lYp2ieoxIlxFvzjbZ1y1mYL9CIqL0qxNVrJ62ZpSzI
lP085B8jhORFW4OxRSzQzpkabfFCDEplN8ydy0LD8qC/EdeF6LR6H7hSk4VS
sV0fDRKEuHFRa6apKgv3P+o+Ed4UQnqPT6q1MgvROWLXEhOFWB+VJbtXgYX/
AcABobw=
       "]]},
     Annotation[#, "Charting`Private`Tag$714404#1"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0.2, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0.2, 0.8}, {0., 0.1728385825660629}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.750682774404161*^9, 3.750683043420875*^9, {3.750683083476033*^9, 
   3.750683128513287*^9}, 3.750683175434534*^9, 3.750685292319397*^9, {
   3.750685326955805*^9, 3.7506853903381433`*^9}},
 CellLabel->"Out[82]=",ExpressionUUID->"cfeefb63-6e9c-4653-877f-2c65a733bcba"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate","izero",
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"",2,83,36,
   26690378039402717963,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.7506853904393682`*^9},
 CellLabel->
  "During evaluation of \
In[80]:=",ExpressionUUID->"5e51c1d7-9687-4590-9f0f-b191bc932536"],

Cell[BoxData["0.`"], "Output",
 CellChangeTimes->{
  3.750682774404161*^9, 3.750683043420875*^9, {3.750683083476033*^9, 
   3.750683128513287*^9}, 3.750683175434534*^9, 3.750685292319397*^9, {
   3.750685326955805*^9, 3.750685390453143*^9}},
 CellLabel->"Out[83]=",ExpressionUUID->"80b43d68-897e-489f-b8ab-66ccf31907a6"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Psi]", "[", 
    RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
   RowBox[{"[", "i", "]"}], "]"}], 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"\[Psi]", "[", 
        RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
       RowBox[{"[", "j", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
     RowBox[{
      RowBox[{"\[Psi]", "[", 
       RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
      RowBox[{"[", "j", "]"}], "]"}]}]}], ")"}]}]], "Input",
 CellLabel->
  "In[121]:=",ExpressionUUID->"d8302f1d-7abf-4f8b-82be-4bc5dd2fefb2"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"j", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"i", "=", "3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"\[Psi]", "[", 
       RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
      RowBox[{"[", "i", "]"}], "]"}], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"\[Psi]", "[", 
           RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
          RowBox[{"[", "j", "]"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", 
          SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
        RowBox[{
         RowBox[{"\[Psi]", "[", 
          RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], "[", 
         RowBox[{"[", "j", "]"}], "]"}]}]}], ")"}]}], "/.", 
    RowBox[{"\[Theta]", "\[Rule]", "\[Theta]\[Theta]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"\[Theta]\[Theta]", ",", "0", ",", 
     RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.750631082830768*^9, 3.7506311144570513`*^9}, {
   3.7506312386703577`*^9, 3.750631263119049*^9}, {3.750682482775589*^9, 
   3.750682489589179*^9}, {3.750682869451111*^9, 3.750682889288278*^9}, 
   3.7506831650596037`*^9},
 CellLabel->
  "In[107]:=",ExpressionUUID->"7c351b89-51d3-4b3c-bfcd-0b9efe9cbfcd"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJw8l3c4lm/4/+09Hx6jJCklqg9Ci84L0SCzkJlRihIqIyOljIaszMrekT3j
ssmWlZnsB1H2iu/zPY7f7/v88xyv477u6zyv9/k+z/u+D5g/0L5FRUFBEUlD
QfG//9SYqXVM5tl5iv/3EyGVnxMMNoX/z5tZ/1n2BT34P57gvrNeEPT0//gN
fwNrZqnX/3GSqFE0RaYPXJcjKej7l56vOPVHQifmNUTZi/CKa3w/36/yoiox
6B0I95xkf/Ju5PymGIuYW1Ug7E9gPjOHSed9pHocEkuDgWffr0MtVovnuc/G
lLbmvQeviZdz+TPr52MVrGnWM8Jg9mX4QaGzFHDisvTVA8kR8OLM23xmXxoo
1dx5fyUmCgSXphIobjHCJf2G4YcRH6EmO31pvocNuk2DjnwMioZG2+l0kj0B
sk+eOD1WGwOdHiEOAw1EcDoVPOhUFQv/bE8b6j3kh/Pn1j1ZcBxcYm5svP9o
H9CAsUhsaTxYPvh2J3JKCBoVq77JFCWAk/GFXQulgxCocsS2MS8RtLiE+Q+1
iYCg+t+CpYxkyNocMexxE4dxLV1D37QUyBw5Or4SewLSr5fuCiSnQpP0UZS7
IwGnjV9eUolJB89MI63eTGnYuTkz1//hM6QrB0aM5shCraVG4IOIDLD0ZzC5
3n0atO/x94cHfQGXQwGdJc7ywHIacOjFLEC8n82ubQH8Kbbn/FWbBcflur07
NBSgoKwn/3FVNmhGOU8ZiF6ACGBkqEA5wLY2cyn4izK4VZ4zYMI5wLJ0u1/7
8kVQqo3Z+VSaCxcqn+nwFlyBwxc7NUln8uBAC2XOsp8aMH2jjT9ZlAfeaX4H
mBzUob3Z+mJDXj7MHl/guIS0wLhLOuBvRiFg66V2l47rgHStRs8dK4Lbv8Pe
yxP04NCPCGnvtCKIbTj7oMNEH2YGdn7sSS6GrMnHgUoChtBiLCl+S6QEhOQO
9/u/NYKsnxbuX+JLQGvUQPIFswk4jTUcuBBTCq5Y5PjBsJtgcGvrob/gV2Do
523Bd83g/NTxuh8fvgLXI96f+5TMgWY2yPp+RBkM2M8eesVkCYGLRrmhQRi6
wmnecy5YwVtipFwFawV4G73oec55F75nv9wvo10BB+uOHhuVswbDGaMJwYEK
+JfY4vMn8x7YGjLbL85VwvEfIhzLlPaQu7aqYyVZBefSDMNyi+xhPXhUdvBx
FVQkGh+4/8gBnjUXb9fuVMFqx+8U1q2HECJ/xzeCowaI10qli1WcoK9Px4bt
eg0IUEqq7xx2BkFHUPeKqIE2BfWyIhYXSM7k4b4vXAuxdHKPGqafQMn+2mgk
XQfZjFA0Me4BO6VZz/Od64C2N5Tu7uJTUNL/cEusrA74rRfukJ57QnPAQ3Fu
lXqw9WJwuBj9DH5SChdO6TVAywF25LvtBYc+sUYafWiAxDIntfCoF3Dn7IZb
x0gD+Kkz9rXIv4RFh3bF0rvfwElRz67itTfQjHu0+rs2QuNXxRxjdT+47Gmd
TV3RCNtcEqnS9K/AX0A3xJmmCfLumwwbVr8C3uvHDczfNgHBJHgqQ/kNiNUN
TMhEN0OJyKlfgw7v4IF5fUPaWDNMzVkE3TgYAHk7Oen7RVtgynFw+J9/AMif
emXPmN0C/db1PC42gaCZcvrfYFUrfP286WZ1JhjO1V5hUnzYDiuVniOqYmHA
0ZXUdiq4Heo/NPg8Tw6D8VGq98dz28HC6bXJlcPh4E9Rup9/qR1+BIYyvxON
gNFz4jJ/HDpgMKXj3SJEgW8O881PDt/hZL6sKW9JNBhXWomEBH0HC5mSMydF
Y0CqvXrGL+c7UOcAh/CjGBj87fr48eJ3MHKjprzGFgv/HZ17pebQCa/rn93Y
Uo2Dnujm/E37Lng2WLP2bTYB0jNFXf8EdsHfniRPIaVE8Cx7gSazu4DHbs/p
k1GJcHRArqnjbxfQN7tR+F1NAjeejJEU+25YO+Wl3FOSDAffvmXWs++B8VXL
nZ6CNLB7om6Wa/cDZsIWfCJDsuCKrVwKQwCZr3XJvPuWBQfNxRaMv/wAT2aF
GzU7WdBzmc6dYf4HFB/Od529mw3y/OXhxjZ90PRKw5FFKQeYio6309/uBxcB
zWteu7mQuMICRoaD8DQmqSwsrBA8SJsvs58MwqHWr3ls3wtBf2i6mS5yEAxb
RL6FsxQBS22NQfaPQeio+mL981kROIS4OdLpDoEXe5g37f1iQCfnMrI0h8Es
8JPo/gulcL/olj/VwRGQ91MTYZouJ893LSFK6RE46lMoLnAYw58+uZzdCyNQ
fllsc9MCw6Mtrp7t2yNwxFeqhesnhifnq/atpY3AlPMElrxL7uMawYxZqV+g
cu2zitv1Srh0nQlmlH6BcjBvEnKoBMbJlfbpa7/gNutPrS7/SnhD37I84fgL
blTUBljUV0LwFVe5kZJfcHAy++GaTBXEtPc2dSmOwvKrHO1o+mowN6s27tQZ
BaKSiEacUDUcXMxc6LAchVpDmgNGZ6ohicubq817FJ7NDzPO362GDF1pw2+N
o/BAQULWtL4abKf2/64fGAXTcxWiv4eqQcKZ+Wnd3CjMcgg+EFmuhryI0bhq
tjFwpxnN8t9fAyWDATNl2mNwrGqYl8+hBtzuu7l9tRiDsWNdxc0va0B+x4qt
9NEY8LqECe8nz42K/SBVFDoGtJmxNGHlNVBv/tslp38MvMtmlJJpasF36Qdz
9uwYcIhxXhMk1sKVFzUfv2yPQUrEzczDIrXQmhhV+VlwHEyHRNgmlGqha/oy
Y7L5OGwp30K/XWthIUJU2cthHEIrWidnfWqBSZX+menzcXBhu+n+NLgWUEbN
Om/cOAT/NvW7mVoLhibx0svZ40D00/JNz60FR/bndu2V45Ae9PKVd1ktpNvD
tO+vcWgVrSuYaquFOmHBQ7f+jsPp/UEONj9qYaRz21SBcgK8R05efzJSC1sv
BqL2cU5AvbCdGNd0LfDIlvRuCE0Ag5sKll6oBcmpcK4eiQlYj9/6NbRSC2rh
Tho5aAKqcw+as27XgtVl3df+mhOw1/IzRxVlHTzflK63vkm+//IU3qSrg4/p
XNQX7SaAgjlEqYC5DoqMFs8f9JyAfR9n7JbY6+A7a8eT3XcTsN8uUTCbqw5+
l38pGIiegM5YZpF5njpgsPNfLPwyAU9Ft4xT+etA+MD9EyF4Ag5uLib82lsH
8t9Vre3aJkDgv3297/fVgb6XWJLazwmY2H+js0awDh5KM46KLkxAhLvyk/v7
68B/Ymof7e4EDFleyfQnc2po3Y1fbJPwxzpe5wiZay4mvi8TnIRXCXc1Zcn3
/1z36og4MQnUSR/flAjUwWaqOavj+UnodpNczNlTB9yGCpe11SehNJZwX5Cv
Dv5jEXp5wmQS1MZOzm5x18GVsp0KJttJ+HkpzOgKZx1Y2g5tT7pPwvacXi4D
ax083f/1dPXbSfBtyx2TYaiDyPbIR9EfJ2HpE3G8naoO8p+5ZLlmTIKOs2DG
D7Le7VL6c3plkzD58csZtdVaoHtPtOQYmgSqYOsQh6laOKCyHD03NwlVnvMW
B3/Wgtza94GG7UkoFKjaONdTC/Y3AnSeCUzBFeGo9zFVtfCG6cE742NTkH3J
NWumsBaSSq82nZGbAmqzh6Xhn2thcB+z0qLhFOgx+U8JhNTCWivJo9VmCs69
PS3117sWCJ4NJWmuU/D2yr4KEZdauDT6UsoiagoC393h+WxYCxbBlraQPgWf
3MUEt1VrweOCUtre0ikQu7MwlXKuFnKTKIS7+qdgs9PXlMhP7he9n8ZZM1NQ
4pxu0UlXC9MM5RFvNqdgfdCjZ2apBgRtXDmV90xDok5St3BzDbz6b5Wi4MY0
XO+meRH7qAaeOF0Wn7SeBkFBxc/chjVwF3+4zuM2DV+upu18RTVwSUMpzfHT
NOT9uhGcylgD9A8CdE6NTYM+4WlYf2A1rBaMuVstTwNr1V7Wn/bVMLkjmxJG
S4IciXgGglY11PkPba8dIcEdJaZICtZqeJEpllR0jwSxMoNOc0+q4NGqe/u0
OwlG62PVtK5VgcX5jk2+dySQZgo5+/NYFSi2Omm4ZJOg8iZmeTRQCZTzNetn
Vskc22rkcKISPI6bqpY+nYHBz30juboYbB/nPJ4NmAFPN6tpsflyMC6jjdkb
NwMvRO+aTnqVg9zVz8uuNTPATa2QI5pWBlv3Nj7KMc6C80Y+Z/uvUnD+HPyn
LGgWJIcanrr8K4Q7y5N75uNnwY7FxJzHpRD05M4qC+bPwhX9w7dWFgtAtnkk
wqN3FoQ5bg8++ZUPy7PHlUBgDpZOmxVnpOaCvXjD+4rEOQh+OMJt0JoB47fG
PWMK5uCp2kmhoyGfQTeG4p5n/RzQ9cymftFNh7PEM4qINAd3/SsEZztSgGo3
dR4f+w2e5+PW9WLjIej764s49zdIDYVef6gWCDkuGhtllfOw78SmpZhSND6U
ZzP+8fs8uOhFLF2CWBw679PmPjYPZgQH776fcfiJRUWiPO0CCOpoJ7gfTcKK
VyW1yy4tAJFhwu7C8GfcKcSV/rVtAVgPOCe+sczHmxY2mUO2f+DWRDKzs2AV
vsm55PzV8Q+Mr7U+KAiswrXlT5SiPP7ApVRe37u01TiA/9UPff8/YKr+pu7I
fDU+0pZC1ZXxB/CLIK2+mlp87eykbuPcH3B9riz2+eU3nMlhtltg/Rc+rri4
xCW0Y/MyXc2A24vgb5dNd7foB5Yq4CPwP1gEEcdVw9O7PzDVl/7OWKdFkEho
1KBT6cNxMSZ6ub6LcDiqKDu8sw+Pet026UlbhLtzU+/s5/ux+RXHewILi+BR
uf64+sAQNu9975vitAQNKfeC00pGsFS73hVJzyWoNLT1TRscwVTf+FlKfJeg
+r3BWvDOCI4r+fiuKWIJfk07vRJS/IVHPyaGzpcuwd2/w76sDb+wuWV+gjTF
Mox+zWKMaR3FFn+7Kip8lyFHp+sHXdU4vvkkQK8oYBnuEyTi/gyNYyNqtfkv
4ctAUMoIyNsYx7rE6r3RKcvA5rI02CAxgS+dyXZ0b1gGGpLSdeLHCXzi2dtj
ZxlW4FbvDaaO+5NYjOlytSTHCmh7Dz4V9pnER4JpDI7yrUDFpvkdlZhJLJTo
6sMnugIujGWY+vsk5vp2d3T14goU8NEuXjw5hdc5VMJzfVbAnf6xZ8jcFF6O
oPwv/d0K6AnySv6gmsZ/hctq48JWYK5q+t4U3zSekZZeCkxegQsDlHO3L0zj
IX1hdbv6FajXO8rwM3waV8XsUB+nX4WjzrcYTE+RcPnRkshD7Kugo1364PgV
Ei7NeSwpwLsKqVefPKg0IuG8mjkT5iOrwBCYOnL6GQknT/cXk1RWgc8n28Ko
noTfSBQ+SPJehclpcJO4NIMbpXonH/uvQlPLit+U7gxmkFkzVg5dBaZ0qTaz
WzP4xZlTauOJqyCsNUsV/WwGuyoWigrXroL14c9mMoUzuOhCb/Tf5lWopiA6
7lbP4FWVNZ7KrlXyex9h4lX7DHZQPUV7c3wVPJc1N8enZ7D1tcJfH6nXQFrI
Wj+GZxan6Pbq32deAwnJiXcPDsziSf21NjmuNZifL3+3Iz6LzYxPlQ8Ir8GQ
p4isIZrFn0z1ZNLF1mA7wstb9MosHjRz+vxEag1edowv5unMYv3bhZH8imtQ
1ZPryXR7Fofe6eUkXV6D23tzrnfazuIu6zXfIq01EOHtyTBwmsWaD0456Zmt
gaGerFGCzyy+6Fyo895zDQrKpY4IJc7iF096Gy191+AS/vS8L30WV7mtKUgH
rIFce5m1Zs4sPv/slERn9Bo8ffPG3bF8Frt56SXHJa9B23iki0jNLC556STo
8GUNCOF7EgO/zWLZV4UsnHgNQi5TnYvtnMWP3vR6jdStwbm4Ac8LP2Zxjv/a
5pfWNbD8ypiRMDiL/wTwOjztWYMMGtusypFZfDz4FEl9eA16Zo+7R43PYuv3
ejcFJ8l6HMylk5om6xnm1Pv79xr8605XeDlL1jMiTL1sZQ36W5kOvZ+fxYc+
FNa++bcGZon0iWZ/Z7H5p145I9p1qO8VLJtZmsUxMWu54qzrcOsX522J1Vk8
HMcrvsW9DimXpiKl1mexQOKpuEaBdVhOmNH7uzGLDZL1+CMPrUPxh8+R97Zm
cXiqU8DdY+uQPUBlnrA9i3vSw+jPSK8DLSVDZvi/WcyVWejBILcOLS/3el7e
Ieuf1bvSq7QO0VG/+gvI7J+zdi9ZdR34qJWqRsjcnMc77qizDl/7752rJzNT
4SlDFcN16L3VrmZN5kvFet+JFuuQ1CS88428v3ep0+UJ63XwaDXSnyTHrykL
q8hzWIfLriW3MDk/qorCUy+erMOpjhLJ65uzOI2VkPTv+TowNX34z5N8Pm3D
e9xOr9fhL9v7A0bk82+m1D3/E7wOlK9v7CWR9YlfFfp798M6sKgaSx0h66d2
wdV0LGEdygva3HnI+i4HdrcYZazDteoonsKZWXzh+Ks0jfJ1kDsa+IFlbBbP
PRnn+1a3DiogVVM5PItDGs77KLatA3LdkynaP4snLJYsZEbWgQoHarK0k/XI
vvo9c3odEnzUtfwayf7ZTUaif9fBJSY9oYjsL+9Io317qTYgaktr8kTRLD4x
XfA6hGkD4t/eEHPMJusvw7nJyrUBsrqu9+zSZvGRjtoeikMbEK6xm+IRNYvb
BYWUnxzbAFLbC7rg4FnsfO9J7pL0BpSw3BE1fD2LG+j/C5hU3oDU5IpA+iez
2O66346p+gaU1765PGc3i/nix+716W7AasrNpJdWs/ju+fDLzVYboBOcp9BL
7kemx5TUOX4bwMUptbUtSvZ3laG9eNAGXJ5baXq7j+wfjoKfCZEbwMn6XLOb
cxanp1t/DUvfAJ9+jUi/jRmsNtr52L2FnN/dgjfvamfw8n8nxle7N8Dqdq7a
MHn+fHD31bYb3oCxcywnJ1Nn8G8++f8sFjZg7Wplg/DbGeyvkTR9iXMT5kqW
q/U1Z7DsRwq9Kv5NELVi13+CZvDwjEHtOeFN4DOUH9eTmMH/ebPHnTi5CT8/
NN+5xD6D2786G3JfJ3P1DDNTIwlziqm2/gzfhEM9AvdtTpLwJH00Z0PMJpwb
e3uoS4iESyYWr2WlbEKFmuQdejYStoyNHPAs2oRI3nuEtqlpXMQ3M33gxyZk
8YVG/CbPdxN6P2pL3i34Uq7mhf9O4c/jdaen32/Bc+utj7Gek/hZ1R639o9b
sBWk23foziTWjbHFRYlbwHJ7Yp+nxiSmMOJR8cvfAvOmnzEV+yaxTuctHbHu
LaApCKt0LJrAm5U0tve4t4F5+1WV8NQ4vhStGP8naBuG/KTl546NYZ+G9qHn
UdvA95Z/fJN1DNf9NeEjJmzDJq/S6d/zo1j5guvbM/nbcLlAx+9p1ihWJOU6
Pe/dhro93Er3pUax3EkRNW6Bf5ChGWJMKfULS9TTr5xK/AdS6xqpty4PY96F
ZpWnhTvQZamrtLvRjcMDsRkJ74DQJ+H9T4u7MZ90jptOww7MzYb843Uhs0tY
zpG+HWgyk5GJWO/Ce6jNBds3dyDtoJ3p+nonFuRdWxE6vwsfh7Q391B8x4fh
QGJVzS6kigxbXzBpxkm/uCqOtezCZijV8vupJnz4Bd1AaPcuOH9ilzSwb8JH
vs1y2Ezugj49Tb+CdyMOOWpYzkhDgTCR6Fqe34BDc0XEckQo0NENNc6xY7U4
sr6EkvouBdreWzIhwVqG/4xkuSzco0ChE6dOqWl+xRc3k/4O2FGgksnNI89D
SvHyseBfec4USOdI4ztloRKsHmRTeduXAr1yPHDLRLkQUxkLPGtMpkBs0Qz7
TOqzsb4jYaMgnQK1ipru2xDIxl/eMdjHf6FAYhpb/d8eZmGjquWbboUUaGPY
3/vQlUxccKQFnainQC9VNGkaj6Vim0V3iqBJCnQys0cwpesTrmJ+5OwxQ4HO
yll/n3/yAfOJWP+xnqdA3PJXHW4cjsR1erojSqsUiG9exCT1zXssVHaiYoWW
ErGGhIaI/eeHu3x+Pr0hQomuit88sW1gD2Kx3WvKRymRZbZU4oKDCzwraXog
dZwSFd6N52or8YQTvwtNmWUoEY2hD5vvR1/w1Q6AsguUKLn9nYXV32CQ26ew
K2RJiVylg6Zy98TCJaMuB/M7lCj3MO0kA20cXIuymoy/R4m6nEo6pOfiwIbf
v+XwY0pEyd7XxZmdAJHEgahj3pRo11W1eIo7BZKu2bLbvqJEn61+VmxOpEBO
MKXXF39KpBx48M1Ubip84xS1lgqjRMJvCde1VNJhnfXx6dMplEgvI/iqh3Im
0Fxl+OzymRJ9mld94fQvEzjeRO0vzaJE/hFDF7VzvsARpiq688WUKJPffEzf
Lwt06di7lRopUX+t6qZQbzbk7aTYa8xRouWO5RGeo3lQIS83EfCHEqnJnit6
fi8Pmt3a9L8vUyJ70twe28w8GNtcgev/KNGx2P+sGo7nA2FNkc2QjQox5qLj
X/YUgKBsz7MPBCrEs3pS9px2AYg9vrsyxEOFxHh+3tXyLQCFpXeDN/dToWeZ
uWeE/xSA3cJg2m0JKvR9Ws7bM7MQ3E7YCaZIU6HjlCJFvv2F4HufOoh0mgpt
NiykH6EtgujZoy73FKjQldXWbAndIkgTK5vLUKZCNq2sip/ciqDgrubNhctU
iDLOICA1tghapxwvOmhTIY2Vs1kxU0XQf5ipNFeXCgkXWJ4KYCiGiVsfT6wY
UKHDRp+mDokWw/ZYNdHZggrVMAc4KVsUA8NBPb9iKyqUhWMUuj2Kgdt8ZnvT
hgoFSFxkJEQUg/gIx7j7IyqETzxIj/xWDOpTQlWOzlSo69uUufNwMdjNS8Q8
cKNCG8WHWmQWiyFoBXnc8aRCnHa5sV00JZC3rWlk9oIKSY77f9DlKYEearOz
Br5UaCeH3aPscAmsM9nz6byhQjePDzIxy5aAHH9g14UQKvR7iahuplUCJkKx
OfLhVEhXuSbuvnEJeB7JDpD9QN6PV+in1Z0SiDtRaftfDBUiuVsMajmUQI1M
h5poAhV6kpjvdNy1BCblfokdSKFCSkF7graflwDDhb8Mez6Tz/vDnKPSrwTE
VCmnCFlUaF5174DbuxJQ0+asZc6jQt0JbsOSISVge+NAPE0RFVoRMmQaDSuB
gJuSz/6VUqE3e3i030WWQI6VgukqpkJ/Inw+n/tQAl22WvIL1VSobWSJi0Tm
1cdme6frydcX1l9GkJnP3X5jpIkcj96DWiOqBM6+eNbb10aFkEh8AHNECRi9
Dsz/3kmF7jhnSXe8LwGPoNjgpl4qFNZPuxkbWAIxEdn2NQNU6BAwzT99UwJV
MZUaZT+p0Ljw3gP3vUtgPLnjeMEYWf8N1ySbpyVA9+UX85cpKiQSphzi7lQC
ogV/ScmzVOhe9HXqRNsSuFJG2RCzQIUuzcqzj1uWwL0azqSIJSr0fOVkj5xB
Cfg3HXgRtEaFGJQpXLLVSyDru6T56y0q1P7j954LiiWwPKIl6E5Njbxyn600
kevJM222/ZieGgWdzPes5S2B0wv2/bbM1IhFizprgr4E3P4Fht7kokZs114V
xU0UwyeauEc3eKnRb/HWFaXOYqhgztHW3kuN4t8ciOauKAaaPd/ZLhykRj73
CrmuhhXDG1mCz5GT1KjlTeGGvlwxZMoL3xI6RY1iDT7u5B4shvYLUkr856jR
gOWGniYT2c862rtMStToeWw3v253EXx4EOQ0r02NznvXfuo2L4Jyx7jrU7rU
KKA2idL1QhGMuOecHDGgRodZT28/FCmCQ2++L3SYUyMVe+r5uxOFkJFCuJPv
QI1eHh2zSjIuhBfdWIHaiRrJ7qwNiJ0qBCOq+3u1XKmRAYx/32QvBCaj+tY5
L2q0Sni2xw8XwG12N5lD76mRTafCSDZPAQg5TlAFF1Ij77gasabkPFiLCxoc
KaVGK2arhdZOedDaBgUnKqjR0eemxJcqeeB+NPJuYwM1cnTNpZMdy4X+QfV2
yn5q5LlEsdbPkwshSkUf7bbJ+Yj9zSu/mw02drecyiloUBY/O4+1VDYofiRo
sdDSIM05vZbyzSz4u3qPJoWVBpUw5iQ+eJUF6mkHbH4K0iBRkYhMN7UvwMD5
+rS6Ag1aKW51oE1IB7dh405xbxq0N1ia12BPAixflP946xUNOqPwumQqLh7u
ZQtYRfvToKs37f8+FI8HwxcDm4QwGiRUvlsdKRcHZ8VuCG8m06AY170U8aYx
sPH4mn3DNxq0YE4nwX8jAhxZVdlvsdKi/Vx55TI7L2HeUazvEyct2qYWX/JW
eAG3RxjjfxBpEfUFy3KS93PQzW2QVROkRct5zUnnJDxA5oaKycn/aJGmI3dt
WLod+btHIYNSixZJp4nsf9Rrju3lTqt+CqFFh3aDWmiD3mAK25Cn78Jp0bmn
RxsrB95i/+g/uZ4faFH/qtOyifA7nEaVKmCRQEt+XxCrZesMwKMN/PNH8mgR
lTR3lI1SML52fSsgu5MWxSwJmbt1huFRb73auF5aRBNUQ51hHI7tinI3ggdo
kbJTjRbPZDh+s/eeueMYLTrNfm7F+28Erv01ePLcEi1ScH//9t1iFJa1Le+p
IdCh2oQVhcnIaFwTvYe5gIcOTRms5KdTx2DtDkdI3kOHTh/THpw9HYNtpf9L
8ROmQzxvSdEeMTE4eTPaRV2KDiUOsHzvsIjFsse2M0CWDr29cj9QLCgW1xjr
j0qcpUPUKVa8XjgWj1RwqHIp0iG2jCIOCd44zOfzTOCHFh2K848xsCiJw0lF
Q5rfrtOheUXeqviROCw9c+ZlyQ06tDqS0DxEG481ry7+/mBGh3yEdE7KXY3H
wx7qwv636JBAo6GO6YN4fC8rTffpXToUuvfMJ5eAeOzHZYHN7OlQv0+TZ2Bb
POZTxkvaj+nQGztONv/f8TjJca/oBRc6lPf8bPZTpgRc1dcZcPgZHZJbBCuk
kIA1mSVqeV/SIYVrh2LYDBPwsNybDQY/OnSUf/ps+8MEvBl9wXw2gA4dT712
8FhcAj6/Jh8k+Z4OiTJu8r4uTsDP1E9VOUXQob/qU1mf2xJwTaLEYtlHOiRb
EhsSMJ6AGf4dFaaJo0NaGW/tZdcTsOq1g9pXkujQ6Nbwv49Midg/XeB5QBod
qhD0W6vbm4i/U/Hk9GTSIQv1/w5liydiogH7qEAuHVL05tYyOpuI9bMZCBaF
dMjhAbNOw8VE/IGBSjG1lA5RnLpNt6KTiH+abtkvYDrUqD6pPmqSiIULl2Nl
aujQ4z2feQPuJOLbbPMdrg106Merwzo7dok49dYUZVUzHQpOpNo65pyI576O
SDJ00CGX3T0cAh6JWIK730y9m1zvQ1aRTc8T8UObzsCQPjqUsjSTcs47ERdU
NVf2D9EhmblErXu+iXiTv+6v0CgdSn3wPNvYLxHL2+MDVpN0iGbjwhw7mT0b
irQyZuiQ4yXXY14+ibhmf86zpXlyPj428YUvEjG9U3r2mSVyPrep/TI8E/GV
1oRfT9fokN2J+r23XMl6iXzirNuiQ2OdFTHDjxJxh1uYAgsFPdLetbXYcz8R
c3cF2GvT0KMhYl4+nyVZL/FXseEM9Mharnah70Yijnru1THMQo9uGvLbmagn
4uE+N0oRTnoUIKQflaRA1kvSUdKGSI/+UF4ezDmZiG/5PjDL5qdHUYpV8V6H
EnHKzzuBa/vo0b31sDv83GS9ZM0r5YXp0U+e18mPqBKxw8S1A41i9EiSlSHB
fSABF8ipa3H8R49eVux1PlaXgDeCLz7TPUmPfNvOL8V9ScCeimd/jZ6jR0/2
frtM8kjA1REnOY8iehSS7H6nxDIB0/09pvDgAj0Kzy/Ju345Ab+N3h+7pUaP
Li34hf1mS8CR2zRm3Cb06O0ZjtWC1+R+0dkJMDCnR1w8XY2/reLxgfS1ipjb
5Osdm93LivE45caM0PEH9OhVKPei40oczi9oG1H2pEdrT/9uWqjF4XXWbxxv
XtCjfQG8e2T2x2G5W1Xouy890is6kjH0JxZXceXHmATSI5NRH+HigFjcbhd5
0ymOHlGbvBvOr4rBYgEcW4+S6NFfK4FmqoAY7PXF+71DGj16V+p255pRDJad
d2i8n0OPnnPTRBLGo/FHG1Vpyyp6tJvl7Jl88BO2vr1NrzVKjyCFRsn3YiSu
eekQpz5Jj0ZGiio+VEfgfYnTcmoz9Gg/Rzpzv3wEbh/rsr+4SI/cfM697z0W
jk+ZZQzIUzEgVPv68t6V95jGyOSLmDADUi73ZSmRCMQmrl2XRQ8zoMp1xmHL
gABcGHllXESMAWnQ1aPiwXfYuk+W/4AUA9o3yMnncuUt7tBl9+JVZEDiqR+/
UO/44E9aFbo05gxInpimLRnsgs9cFN4ZimVAkMiNYr47gXR8o299IgOSY6RZ
v8voCid2HbiyUxlQsgGzpYeyBwgXVou+yGZAAue+tWTefg7MR8jf65UMyLvj
Y5YVqx8M0aWkOI4woDnFF47Nd4LhWd1xHfb9jCitLlqtkDoWvimfSUEfGFHm
cQ/Jaq8sSMyOODOazoRe+AYpWxZWgvWWodiPL0zoW9Xp4uj+SpBQFtzbmsuE
hL3D2lv/VcLXH3HbxaVMyGXIsoL2QhV07qbjwCYmRDtdaljfVgWU6uXKCrNM
aOjqP7l949VgPDOqFSPGjNjY6lg5t2tBWDpRMfQEM7qiIb2dIlAH0+5WJ99I
MaO4tUoOYfk6eMQ5x+18lhn1cXRI5bvXwevTK70aV5gR8/jlyav/6qDYm95k
9y4zMjvz+1bpSj0QDx6zNkllRgVdjXabdd8gIOi1bmoGMwpaNnxXM/4NmKhm
FZezmVGmvWuTFlUjUPxK3fuqhBk5VETw5cg1wtynIy35TcxI2fmkin1mI9Ts
OSjBOs+Mck5w9B542QTw6rmA/iIz8qLdQ5/7qQlKNn4xxK8yozL38O31wibI
+hHz6/QuM6LwPuwfTWqCD6GCwZYcLMgo0uCN/OVmeETYs1oqxYLuR3e9/73e
DAvPXEbpTrGgGXFTDSPWFrD++6NV6xwLelZFuWZ3oAVutoclTymxIJe7ihq+
l1pAzZ94g+s6C+p7z9Z2LagF6v89Uja5wYJOW71L+hHfAor3uyRTjVmQsFjD
+GZeC5xWC2aC2yzokLxsG09PCxxk4vxq7cSCJn6pzgdytcInF7uUfFcWlLIe
5j9/oBX4SG0hFJ4sKEasjGH0v1Zga/C/H+rLgu4yDdf4XWkFv1PzN369YUGs
zKqH5PRagSb5qsqxQBa0lfU0/JlFK2y+ZNlfFcGCpjf/bFU+aYXHKzbMrJ9Y
0Gc5NqbSl63wx7JpTS+OBWnWVFgrBbTCpNKr9t9pLCj21AvLnfhWMMslfT39
hQUd1UsX/C+jFQaFL6d65bKgU7Sx92fzW+E7JcNz/q8s6NWDvRzUda0wnxFh
olnBgiyjwxXOtbQCo4H4OZ8aFqR4d/PzbmcrHKIr4y1vYEFPHh9z0OhvBZSj
vrzczIKGB4Oyj460gpHJSLt4Bwsqs+Z483aiFZyYHDLMu1mQ1CzLwSczrRBU
QP0qoo8FPX3pGj4/3wqZ5u9vtw+xoJfCv9hmFluhke2IEv0oC6oYNki1WW2F
iZKi/ecnWdC6T8Jrh41WoLC6sv1ohgXpH97f/2+rFQS4Bn+kz7Og0IDeYpad
VjiN7+ePLrIgdea3mp92W+GazW4g/xoLGlssrMmmaAM73kBbzS0WpGMXIK9C
2Qavq4VVfXZZkMBp4pgBmZMe5B0pp2ZFITqE4Xny+qq9KjQr9KyIb/mIPiWZ
h+t7R8RZWNFRzQ7XQHK8zYd3y8w5WFG7n82jT9utwCO0FRHBzYoGPsnZHtls
BcnmN47tfKwo1uxJqPhaK6g5C+rQ72NFJ1219qYttcKdQ1n/nT/AiuYfz4nH
LbTCi3YFlscirCi9/ck6YbYVot06p9OPsiLOPadztsh6lojeqh09zooWuR6+
NiTr3dO1GssvxYpuD4cUSpLrsejp66Epy4oei4nauJPrxXp8j6HPWVbE+mhk
QL65FS68lOdeUWRFEg4vO/i/tsJNybY/4hdZkczTyfqTua3gOnSzxVyVFcUN
WZjg1FbIkXnh3a7Dil7qlBqJvm+F1l9EC3p9VuTx4BGiedUKpLfJcN6IFU0W
8HxV92gFocnG9XRLVrQOZ54fu9UKckFG3aN3WJH4uT865WR/652fz+a/z4qU
m1maSy+3gn8op43PY1YUdTWwdVW8FbZV9IfMvVmRTWG5xg1SC/AtkYojXrEi
KTUOd+neFjgZ7Rra7s+KPtRfLX5cQ+7vtU8a58NYkV3c+dwjH1vgR/JEJX8K
K1IBXctTF1tgWcfpk+ZnVqQuKZMnIdEC7BSMrj5ZrIhxfXjuFV8LqOgfk14p
YkUCko/STKeaIZ/hYVL7N1bEnxNZ7+LeDJNuz36Xt5DjVYzIals0A9/fd9IZ
HazoW+S/cx8uNYNb3+cqvz5WNEEaUPbibIYLaZPDiiRWVFXkcC/wQxM4Cq6I
SP5mJc/TjaAZ9yZICaK+v/8vK7oibLuRbdIEzK5CW1vrrCiyo/Sfn2ATfFc1
4M1nYEOeBmnB18Iagabijkk8CxsKzHxmRHzUCLLSTomBHGxI9XeAt4JmI0QK
hJy05WNDNbc/5fXTNYLZ7xaNI6JsyER3dYr+/jf446/gG3GRDS0NWhzg2NcA
wtSabT6qbIiU0lFk+acerjmZ8DhqsKH3aa+fC1XXQ5Gpa4K2Hjlemt9C8+16
eCqRX8FsxYbmPdK9l5PqgO276Ia7Nxt6mDuhY8ZUC0jlFLr3ig35zvwY4O2q
AYcSZR8DfzYk0vO95OTHGuiJMyeeCmVDB3qvm0Ydr4GPDz9I/kkkr7e9RyF2
sRqOETmszWvZEF3gN2NV40rg+lL6OfEbG7q7JVc2y1sJm5esFqZb2FCaXJpA
b0cFNLiVP3rQw4b4r71//06pAizH7z11m2JDFITd+j3y5fAx91tIKCM7Ytyb
7bJ4qxjYtL1w41V2dKxq3ea8WQaszJ6gYtNmR20dpzf7l9Jh8GX/BS1ddtSZ
fYKq3ysN0oolG3tN2FGBq3sJ04dkUDkw0jX+gB1dfBUr7ZcRCx5/5GZ2AtlR
hxLrOr/pS1jwX+U+2cOObKgaBZoFyO+bfwSrjvezo7RFr4i04iQ8o33xgegw
OzpivJtzSz8Fj/OEf9s3SY5f/DDCLCod938648Gwyo522yjPeuzNxuZ65Uvz
TBzo3XKYqetIIb5qRWlcI8aBGj6zabbsVuGDPj7/PbDmQEtC+y0Yd9uxQ4+0
h/k9DiR6mOXpnXMduEJktPm6LQdy+9w52uPUgU1q5KzlHDiQ9Fvna8N/OnAk
5WIC4xMOdP+1y4NnY98xwc1oT4IfB1q2vChl3tiFqewlaftSOFDruM3F3PAf
WKtiWKc5jQOdeKaz37z7B45hfxOHP3OgXV4OE1pCHz6fOQlJWRxoUTvttcjr
PvxkNurJwyIOFMN8dTHXrR8v3qL7w9rAgTzlzQ8oGg3iUYOBfsUpDuQxEf3B
jnoEs8jrJd8kcaB2N9v/4oVHsMz+zocesxyonkudtVxhBPuONbKULHCgDbte
07SnI/i/eyXnJdfJ90fZm1zdHMFu7hFx+xk5Ee/xyYvyM78wMUbPeluME8V8
KR/ryRrD5593yu45zonsGpJ3N5vG8B1LDerT/3Gil3krQSuTY/ir6MUPD09y
otggEcp7AuPYPFu2deYcJ7o7lKWi8XIcZ1YRpfrUONF4qc+Yq/oEvjTRuZFv
y4lSYpUKV0smsbu4v+sXO0709Ef6rF/rJM6xv/QvxYET+UUu3Jv9NYn37nyl
+ODIiSSH7RZ4GabwPE8S/XMPcj7Z17RYtadwyCVnoro/J2qSGGDbPziFv/lL
hV0M4ETK+a8ybeem8E7XHJ9CECcqIw4Z+2xPYSszMwHpUE6UTF1bTBCYxmef
XDm45xMnanwcJe+nP41H0gUkJzM4Ufhsk2xh3TQmLvbk/PzCiX6Yd0w+7pzG
V04HSvdlc6I9F+pUN35O49wa2tPN+ZzIvCgXrq5PY++h+fM5ZZyoT3ZY9+Bh
Ei49mFqRjjmRsXTKCUpJEv5z10IxsZITRb0wmIs9R8I3VnuVw2s50Rujz8ri
miR8jL1SzaOFE1Ucnh/RcyRhs+uurU5tnMjJIEjr4lMSDo2S0bTv4EQ/ozQE
1n1IePdImo5lNyf69qCS7WM4CUvbWnYb95L1W38zEhBDwnfyBPX0+sj6bsbe
VE4h4e8o2ODKECfSeMbWvFBAwvQ+VweVfnIiR+7b1hNlJCzXQm8i/4sTqfoy
7ouqIeHEG25mEhOcKNgjWlKlg4T7o2XHjk5xopXv0RrnekmYbfKP5UESJ2Kv
ZNxcGCRhZ4dbd3h+cyLF+TbnsCkS9lWdtzq7wInSmsXf/50j4bBDTlYmfzmR
SmXT7d2/JJz0b/f28yVO1H36YHPFKgkX9PjeTloh58+tmCizRcK1XzhvN65x
Ikzq7jTeJeEu38hb8xuciFFW/NwZ6hk8bnbwFmGbEy36nKiop5vBS2c/W8ru
cKJ/UfcvMzLNYCpuGUsDCgLSsOXD/1hmMOF3mYUHFQE9QVTMiewz+ECdikUc
DQG1B8fx7XLOYInoNvM6OgIa5uXpZuaewchZ33yGgYCU3+7/r4k4gzW0fpmx
MROQ1wdRXnneGWwiZm0mxUpAJ2MVn5rxzWBb6qWbuuwEJHv3iOFp/hnsPuh6
8wknAa1OnYmpJPObfJqbn7gIqHfIW/kfmT/4vzWtIhLQw4e1SrNkTrfiMZ3k
JaCdBr33r8lcgqJNmPYQkK+Lquwwef9GflGTEwIExEcwODRGjt+3mGWsLUhA
/W4kk0ieGTzddMbYUYiAZkxZpqjJ+a4nVBlFChOQlN1zvJ9rBtN7qBqVHyKg
scKiud8cM5hHr8tw9DABTbR/e2TLNoMPSxgb0h0loG6Ch2ks8wyWYZw0EBMn
IJZ81QJvhhl8YdTWQP04AQmM674QpJ3BFiGeN0IlCWgtQKD01j8SdrjPeKPk
JAGFXnqpSNwg4ecqQfrDMgRU2+3j4bRMwjHr8XqHzxLQY94TltdmSDir45je
FTkC+kn5+0zbOAnjtHxd2/MEJBHjKvBvmISHjeqvFygSkGBdmoPDdxL+LaN5
vf8CWY+kxv0lTSS8zdZ3bUeFgNa3CtXyyX7dWzmjo6JKQPpN2ul1+SSsL8Km
3aNDQIY8u+43Qsj9sBOqtXmdzD80hBRek7BT734tQX0yF8s5Dz0j95efpOZt
I7I/0o5Ry9iScOfva+qrluT68FevWyqR8Gjd0NU9VgQ0/vbbl7VTJLwYffvq
+bsE5FmjVYqOkTCHtrOa930C+s3NW77KRcLqBVFXeBzJ9X8opvG/8+LOkvku
0ZmALOLPUEuQ54mXhFge8QkBLQYZnxUlz5uitKJ9RA8Ckgu4LWCSPo2Fo7v/
cHkTUFfYisGkwzRe9WEP4wwloMO7p+diV6cwR22PKmc4ATVJVuezTU5hcapP
FJyRBPLz57WUQvcUvul2zJrjEwH1lV1yWc+Zwo12l+XZkwiItvfgSJHNFP50
w2uMpYCATtG35X77PomLQq+EsxSR68kh87uwfBJ3dnJeZSkhICtZFmH7tEnM
oB5TwFxOQPJrSUvKzyaxg2KZH1MduZ66IfFexyaxivjafww9ZD+JkDITXSbw
/La1G/UqAVGxZVy+/XcM39BXPbqxRkCDITcebveO4Zpc8Z75DQJyzf3mcq98
DEdYz/3X/4+AbNq0RgZejeELP+6PZdFyIanLB9/SCI/hqFy7KyZELpTR+ZF5
7PIovmTtyFcszYU25XPPCniN4NjeZ3m2D7lQYBD7+uT9PjxZK/HT/zEXeuim
8vGMWB8+lveT8YsTFzKV1TzqNvkDFwTImy64cqFm7yiRHpMfuPnyBqPdCy7k
fH90jlKjF6+X2pnav+dCTyPFrt040Y3VY0yZHhVyocPxL+Lm/rbjEH82mZBi
LiTs+SwjP7Md97mVmeaVcqGxZpHNUJt2bHljb/4y5kLuPzQuqze1YRdCr+nj
Bi6U83OM+3xBC054oZ7v2MeFSO6HkyfCv+HNO/I3Xba4UPS2oQljbCW+UiqG
Hv/jQonbfAo3jlfiKFY+IftdLtRLYznFUVKB5XIWf96h5kYy+ymbXFMx9thK
NtVn5kbOWWOejwW/Ymp/gumpvdzohV3MJXHpfMySO2W8cpYbXa6iO8yqHION
aLvl/8pxozM7jRxCNh9xhl7Vvt/nudFowtv8es1IrL79YWhckRtRqd4ROCAW
ggOUtY27rnCjXe85LZK4J+b+8dUo14AbuSmorPNc9Yd9/wIN7Z9wI32R+lzp
4c9wafRqcZgbNzomfuPSA9lMeFTPyFvuQc73xVqSzesv0Bj47DuTFzdiGYsq
ayVmg+Nh+8sJr7mRjcpdhytxudCuoXmq9wM3MhkqYtOwKYItaZb3/z5xo4zr
Ld8Gtovg8J6GxYOx3IjO+eStgjfF4DYOmXaJ3GixU+2IdnIJHHX5T4Qpkxsd
XtwkcTR8hedxbFzymBu90s3q/GhYARk+jXYWldzItSjGOTu0An7c8271q+ZG
sTL7oqfbK+DEqR2/nnputMR+5HalYiUMNP2msGsn50P62b9PoAroslNNQr9z
I9aUPWbqmlUgFXrr69cubjSwHTEb51UFvjeHnBj7uNGb0LSvRVNVIL3SPB/3
ixuFUr7Q/5hYDab9vmrfxrhRm/jGr6Xv1fAKX0hbmOBGmbdNzex2q2HE7+st
uRluxKdFmTh4vQbeCn4e7F7kRl5du4ffrtdAEfWds9vL3GjZ76GQ4P5aGJ8+
GC68xo0o94pZ/LhQC2dzo3QebHEj/8cZnY1va2Hy4utGBhoiUqnZWSgg1oG5
6k0hEToias7bF8ImUwdD6jKOCgxERK+1WhylUwdd138eeMJCREX/0fHd8K8D
zRt5TqFsRPTwtkbcq7Q6aDLya8nhIKLcbgvt5do6qLI46TLLTUQH6es2HDbr
4LwVQxs9LxFdg71Wflz1UGw9dOgQPxFdvHRqsl+8HrLsfdqN9hFRccKIjJRB
PXx8OnCkVYSIwlgTWxsy6oHPK8t95ggRMTKLmbZW1EOI98tOOjEiEjKxTCZ2
1sOrt/89hRNERLhGT3F/pR7oAmm6DSWIaPzD8IwHbQM8C+kTc5YiIpf3Uoo9
3A3gHOXVkyVLRG0rncQbkg2w/En/WMtpInrsaPja+3wDPIg7/px0loiQcJnp
vysNMJNI9YNWnohUxQfCy3Ub4HZq73FhIKIApdOoyawBRj9/9jqvQERmf2fN
D91rAOOsZ30GSkS01RLK2fS4Aa4Vir8MvkhEboUlt2l9GqC9hGLgy2Ui2mNc
CFH+DaBa3i3RrEpEmlcjGb3eN0BdZZr39FUiWr64M1EW1QCKtU8HaTSJKHy0
fEMttgHKG65JHdAmItmzc/ckkxrgTPNRX/lrRHQ+qNz6XloD5LXtDN3QJaLC
X34HdjMaQKKz86SjPhFZa2RXzWQ1wOeeFL8gAyKK+RD5XCK3AY70u//MNCKi
0qex4R15DRA3pC3TZEJE6/NM0JbfAIK/jryeuklE7IbrOWIFDRA5vj1CbUFE
fwcMpCfI14nTHbJCt4jo5dVJWgoyB8wmvZGzIqLvpjT33cj7syy4jurfJaI7
4fXZhtkN4LuoefqxDRHZXhwWjM5sAOpVEf/A+0QU8eAl5aX0BvDY2BzLeEDO
N841VS+5ATa328402hPR0/alO21xDeBIkfhu8iF5PZtvSNbHBvhL/WSCypGI
MjWPeW2HNcA9eo1z+52JqIb/gH9eYANMMx0KPPeEXA/hSY7+Vw1gybYxqedG
REPrx+7aejXAT85WuUceRKQW0LTq4NoAPXzO05+fE9HJQy3nvt9tAC2Bq+e/
vSCiD8N8Zf/dbIDm/cIhE95ENNnoQE9N9kP14WYQfE1EL4TEY8hfmXDF/+cV
nrdERIo+eExAhlzvlcXrbO+ISGMpVCvvaAMM1/Df2wkiItHH4zzvOBtg3dwq
fDiSiDwKrN6xD9aDR+OT+J4PZP+W3B6NaakHWin/zNZPRMSzEHSpobweCJT5
NeVxRHSic//+4eh6OBZN9fdTGhHRcN0duGhcDzl0PNuhn8n1ev9fcrpqPZyx
PUr/LpOIPqbdlv9yph5U5DX3Pc0h+21gKCeSux7MBj5eNi0hopKHWSoidXUQ
ynsmbn8jEc3YnS8X5K8DgadqGbzNZH81cUVyUdVB/KRpEXsrEZ1Vu7/n60wt
ZOV7t+52EFHSdjb98dJaaNLp2vzZR0TcD6XeZuvXAkXAA50YEtm/PsOxRj41
YM2QRHOAkQel9iVOjZVWQROV71slZh40dEpS7khUFRz7d5fnNisPys1OFoUn
VbDw97hoOicPij213LF4qgoeDeSryuzhQQEEt6vB2ZXgnlkbdFmcBxlaD2ZT
RVdAwLUJIYerPCjnDf0/hr6v8Fe9PjVEgwel/Yk/cSXoK2hfTpUq1OJBqiFF
3W9VvwLx/P0L29d5UP9zu8eHy0sh8siKlbcJD6JaYzL9FV8CCZs0XyIf8KD7
qkHq/g+LoDDmkHx1IA96Wx+vyHIyDw7tXnCtDOZB73osCaYruRBgfKsYv+dB
3FZFVxcKc+HOniSZrxE86IqzfFiLfC7sCTlyIj+WB3Uz6ygmq+WAm4/Y/uRs
HpThdpVjzikLpieuGCfm8qC9o7yR2eez4PoFm6j4fB7EdaHapoE2C05QfuaN
KeZBbqe+KA8e/QI/nxzniKjkQbwsw5/v2WaAoq0E5esOsn571YJ5+dMgs1kT
/Dp50JFLWTdys1Nhr7i9u083Od+ugTQ21VRYnsre8OrjQaNTyjmMnimQaHZy
0e0XD6LhFTU+tZQEDNdlx2z/8qBf1V56Hrvx0HpOrlabkxfN73e12LD5CK+K
Hm+pcPGi6Pr8Lwe6PsBFmS+S54i8iHIiwmpO/gNUnBD+dJCfF6WHMKpa80RB
7gEGp2UhXiTFnm3CMhAO4fRdoqESvMjayEy1OycErr9gNX0lxYskAmkvHZIJ
Ifv94nsPaV60qJ+SSygOhtebxZRWp8n73bHva6wNAo/f0f2nEC/SzZOhCVgM
AItOmzd9GrzohsWBI0sSr0FIM7GqRYsXsRqYvnpz5xUMNg+vV+rwoqb21EsQ
7we6dVq30vR40TWNgrCxA75wqfjUeVdTXsTsKctdq/QSaGXtHz0w40XmCtxx
9cEvoConLc3CghcJh8Xb1016wbnP+3jVrMj87bVDyvvncCKa+o/AA1700nxH
zJvHE2b2yh/mtOdFjyrPKUxteUByuKMR7UNelMLH7fao3R2EgkgNvx150dqz
rJdPX7kC18u22HIPXqRpI/ZMUcEJ2ikZf+R48qKDLBuZUXKO8MZDkS35OS9q
yBa8qyv/GOic85+88+b9nwbpPBrq/Y/juJDlRJbw/YwlDCFFzc1VLs3nbUu3
BaPsM2MJRYUkLVdSrj1C0iZ+ky2JrCnbp41rSQiVpUWWJEvSYq3f/P8+53XO
4/2EtBgHVr31UTzrd32PeyIFybeaRXPdA/CFU1miZ5IoaNcVLGJFHcZacQXV
N1IoKBbdkfai8iC2z6/W6EujQDR/fGff1gO49GPfjH0mBUc7PSTyX3jgv38O
5gXzKCjfzOjh7XHH70XGXS9m8f3DihbH6zhYSmvxSXseBcwIxm9PfSd8yEs5
dWcxBQtTpxbtJHdhoWDN7X6lFHj0GyywuFb4WsT6xZhyCs6Eiwlff2mGm26a
eDXcp+CB6mRu5zED7F5qgYarKDCNijc/dlkHzz7a9UyoloKcls6W0rsyWGuA
bWD6iIKwSqvJbV81Se0Xr1G3J3yPwcZPKeEMslfgcHpYPQXjuZKJMmZbSITq
aeGaJr5H8rufbf+YEcUNkfd7W/jeCrWJ0yGWpJSZcGiulQIVl2Ala9Z28p5z
o8vwBQWyUUa+kzm7yfHDOTH2XRTsOuQ0nCttS1aeLjQJfknBs4I3eVKYRbam
12aX9FKgmR3N0wnYQ7oK6p3a+yl4LgzC+4/sJYeqWyWm3vJ7Uj7ZcMrbnlzr
fRO8fpACBN4XmJqOhDE2rLNzmN/rJpmo0ClH0jg30e/7kYLejx9KCwqdiLv4
j6SYTxRc7prqlnVzJj/RL4u8zxQw6NlJs0IuJFFHZK5+goLdkfoeWekuRGvL
ysKhKQq2d0zfcF7rSmqtFDyEvlIw4LMY6HHblex1XC1P/0aBe1tkp4Aam3z2
0WqCHxTsCw+KPR3HJudC9E+7zVJQBW+VVcbYRDHKkBE2z9+3LOzRxBxScok5
nL5Igb28ZlNPLIdsz9l2tfoXBcrZUYXhzRzyvtx6d68AgtDyuiXOMi45/tRh
2ZwggpJK8YEKfS5Z2cWtoJYj+ItbpPvMjktyBn18DUURyItW2Q8d5hKTGX8V
e3EELJeGcNNwLukSPN5xdAUC86UGX604LjkoeyYyRRJB46d3vt3xXCJIjzYq
kULgYGlxJCmSS64xLky0ySCIMVgtdPU4lzBML/MmVyHwizjD2OjJJY22mfaS
CghGea5X0yy5xO/pcJ0oQiBkLD59U5NLJLes1xFURPAkZe3K10scUlwQlLyo
hACt097h28EhdqoP5n+oICiLpl2LzeSQbym/902rIoiskDx/8ACHpIlsa/2s
jqAn2N7CdD2H9E28yHy/BoHpVOyJxBw2CfOgifdpI+jPui53wJVN1LrdgrrX
IsjYpFg6J8EmXjXjFs16CHLuUmuFvVzJZJzIWNmfCGw6Lp3sHnEmSb927Sna
jOCftsTHyWecyaagizW3jBBcchfnHZd3Jiec1S7c2Irg/KohsQAjJyKo85dB
tCWC/JLJ+xU+DkS+/nCYiz2Cu9THyGWn7Mj9LeUf9zry7+Hi2exZFnG5M29r
44zgf5qU/2gIi2RcjNaw4PD9rdLf9AbZEh1PXqOeN4J7N8/Wf3u1i5j87pJd
dgzB5bFvURXzFuRtkFLoQgiC4Q9fjFeIW5CzIx7D308g8Hn95NNWRXPS0DpZ
ORaKYKVgd+mEpSmxSRfjdP3L74E4V69KNCZfpWwankchWGFolNQ0YERSIy5t
aIpB0HKn3Kly6xby2o8uWHee739sWwxT6k/iaWSSl5eKYOOqVItGpj7xnWj9
3p+GIKLk9BvtNetIIM/NQvYqAlve4BEbTR1yRjxiIPQGgrZMD60cfzqJqpXb
WJqJ4E7xlR76MVWSeCQ37CMPgUT3wkCEqSJJ72lSsstFsNy1rH/upDTJSnD1
i76FYFa5uuNFqyi5bTpxv+Y2AmysdsBF51fdg3xpB+27CGbOi61aFv6I+ZBz
M5tdguDt4/h9+nLTzEYZg5nkMn4PhpEzWSq/mS9POiYtVvL3q3qXpFgS+I3e
p7eMKv7/efIFZaIyeGjglN7+GgSUuYPJ8pJVePySRGh6HYL9BdW7hzQV8Mzf
Gc0dDxHIpDPynu1DeH5pA030CQKNvohN1v6KWLD00X6TegRXLKW/S1gqY3Gf
PfeO/IfgFFt37uGoCpZWHBbOa+L35MvYYeiiiqnnIXv6W/j+/60P1Nqthlef
E7sp8xzB1e2p2xe71LCW4bUv29r5fVzsDFOzVMd6Y+tw6AsEP9qziFqGOjbI
qE0o6UIgdu5RbOp7dWxiZ9M/8pK/F6igoSJJx+YiA7rKPQjmxSHBW5OOd1QF
nWT1IdDRDTpJ16Zjlr9wY9QbBNaBlcayCnTsRE9TqHmHYGrFhdrFaXXs/lLb
e3oAwVfvTwFl1ep4f9yDMq0hBJavIkxkgtWxP3OnIHsEgdna7JB5JXUc8rXf
NnkUwesWgxZ2uRo+neuf2TCGIMnVbOMGphr+10VgcmEcgezIvbzBYFUcvzLF
hDGFwHFtrhx772qc8lgj3mcaQXjc8+PjyiqYp2ul0/EdQYK9saVZkCK+9fZ1
iMgsgnOuvjW3hWi4OMWv3ngeQdSu9unAaArXzSd45v5CoPg5YVg5UA7XF6mW
9AnQoIaIpW/slcXPPEt+SwvRQPVsc84SyOC+5q70f0RpIPSdZV6sKoE/hPl8
LhanQW9qRjRjQQyP/TFrNLKCBvHSo77TI8vx7DWlV7bSNJh8J1XvLCqAf1sX
romSpUGltoZIa9I8U0QIB1fL0aC0bVJlxus7U7Ky/fEXBRrEBsiL0K98Ycod
9JTRotFgiRZhRg8YZSqpfnNzVaKBkl+bnhfvPZPeGVmUpEID81q7y+RpJ1M3
mlqqV6VB+53x1IDNj5gM4/ydC+r8+z8qOO+s+rf+H6NuZYc=
       "]]},
     Annotation[#, "Charting`Private`Tag$1036748#1"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->
   NCache[{{0, 2 Pi}, {-0.8748630327391735, 0.24222885215789425`}}, {{
     0, 6.283185307179586}, {-0.8748630327391735, 0.24222885215789425`}}],
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.750631106372221*^9, 3.750631133247858*^9}, {
   3.750631190808517*^9, 3.750631264857257*^9}, 3.750631343136437*^9, 
   3.75063170129851*^9, 3.750631767865642*^9, 3.7506824929800587`*^9, 
   3.750682891390669*^9, 3.750683167703678*^9},
 CellLabel->
  "Out[109]=",ExpressionUUID->"0ddf8f9f-57e5-471e-97f9-028805bf8b23"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{Full, Full},
WindowMargins->{{236, Automatic}, {Automatic, 50}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 896, 18, 387, "Chapter",ExpressionUUID->"ad850729-9755-4498-84de-cb71834c15f5"],
Cell[1479, 42, 214, 4, 46, "Input",ExpressionUUID->"653d9c4a-289d-4333-86f4-de334ba4329b"],
Cell[1696, 48, 212, 4, 46, "Input",ExpressionUUID->"a84267fe-8547-4235-8779-ea70f56be419"],
Cell[1911, 54, 213, 4, 46, "Input",ExpressionUUID->"26f78377-a4c3-42d7-8ce1-ac8760bcab6b"],
Cell[2127, 60, 209, 4, 46, "Input",ExpressionUUID->"8e7fb54c-ee22-4c34-9bbf-77f7e8b4ef02"],
Cell[2339, 66, 402, 10, 109, "Input",ExpressionUUID->"3b14f5c5-ea64-4841-a7ab-3bf3afa6d1d5"],
Cell[CellGroupData[{
Cell[2766, 80, 292, 4, 101, "Section",ExpressionUUID->"fdbf016a-a481-49e4-a4a1-7fef444f2077"],
Cell[CellGroupData[{
Cell[3083, 88, 356, 6, 67, "Subsubsection",ExpressionUUID->"d797f85f-1f44-4aa6-bb33-319d3a419081"],
Cell[3442, 96, 2132, 56, 78, "Input",ExpressionUUID->"6ad48be1-82f4-4c5e-848e-3955567c1543"],
Cell[5577, 154, 777, 17, 46, "Input",ExpressionUUID->"5c3f117f-edd1-4683-ace3-bed538a9f32b"],
Cell[6357, 173, 757, 16, 46, "Input",ExpressionUUID->"698cc67f-7248-4dd9-a3ac-3bd95fdd14dd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7151, 194, 173, 3, 81, "Subsection",ExpressionUUID->"1c31843d-c3a2-404d-b4d8-27c8c57f9dbe"],
Cell[7327, 199, 1582, 36, 109, "Input",ExpressionUUID->"9af92159-c78a-49e0-9fe6-46bac46cfd52"],
Cell[CellGroupData[{
Cell[8934, 239, 453, 11, 46, "Input",ExpressionUUID->"b2b8f0ba-19f7-451b-b670-3aed0e09fedd"],
Cell[9390, 252, 12656, 227, 346, "Output",ExpressionUUID->"17986e9e-d1ea-4bba-8002-8d5443cda6fe"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[22095, 485, 223, 4, 81, "Subsection",ExpressionUUID->"ea08ebac-e4fa-4fb3-a8a1-1284274fd018"],
Cell[22321, 491, 1870, 35, 46, "Input",ExpressionUUID->"06f9ecf8-085d-4573-ae5f-34752dd38caa"],
Cell[24194, 528, 841, 19, 46, "Input",ExpressionUUID->"46f7fef6-fb20-4fad-b080-d1185973ae5f"],
Cell[CellGroupData[{
Cell[25060, 551, 661, 12, 46, "Input",ExpressionUUID->"85373ec1-21a2-4953-9ad6-303bee858a0f"],
Cell[25724, 565, 50871, 885, 188, "Output",ExpressionUUID->"74bfe69a-2993-4aa7-8d76-5daa20dba809"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[76656, 1457, 471, 11, 102, "Section",ExpressionUUID->"9c8b182e-87c3-4ec9-8ce3-6172bbb3e292"],
Cell[CellGroupData[{
Cell[77152, 1472, 4233, 81, 178, "Input",ExpressionUUID->"2581c4b6-997e-4d60-9846-58c4feea8940"],
Cell[81388, 1555, 1139, 29, 116, "Output",ExpressionUUID->"e31e44fc-f50f-489e-bf6e-5c9310952d77"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82564, 1589, 2348, 65, 216, "Input",ExpressionUUID->"783c1bc2-4811-4229-b3f2-3a2b13ab8b38"],
Cell[84915, 1656, 7248, 137, 352, "Output",ExpressionUUID->"cfeefb63-6e9c-4653-877f-2c65a733bcba"],
Cell[92166, 1795, 500, 11, 59, "Message",ExpressionUUID->"5e51c1d7-9687-4590-9f0f-b191bc932536"],
Cell[92669, 1808, 320, 5, 52, "Output",ExpressionUUID->"80b43d68-897e-489f-b8ab-66ccf31907a6"]
}, Open  ]],
Cell[93004, 1816, 837, 25, 49, "Input",ExpressionUUID->"d8302f1d-7abf-4f8b-82be-4bc5dd2fefb2"],
Cell[CellGroupData[{
Cell[93866, 1845, 1518, 40, 146, "Input",ExpressionUUID->"7c351b89-51d3-4b3c-bfcd-0b9efe9cbfcd"],
Cell[95387, 1887, 26474, 452, 346, "Output",ExpressionUUID->"0ddf8f9f-57e5-471e-97f9-028805bf8b23"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)