(* 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[     94175,       1826]
NotebookOptionsPosition[     90650,       1762]
NotebookOutlinePosition[     91007,       1778]
CellTagsIndexPosition[     90964,       1775]
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}},
 CellLabel->"In[5]:=",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[8]:=",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[10]:=",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[11]:=",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[12]:=",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[13]:=",ExpressionUUID->"b2b8f0ba-19f7-451b-b670-3aed0e09fedd"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], 
       LineBox[{{1.2822827157509358`*^-7, 0.}, {0.0019271655319089223`, 0.}, {
        0.0038542028355462695`, 0.}, {0.007708277442820964, 0.}, {
        0.007996433256639162, 0.}}], LineBox[CompressedData["
1:eJwVmHc4le8fx23nGOUc6yiVtCSVysy4b76FKHtlJ1nZlR1pIEpGCZVRNsne
uSVKksjMyDySJI7zhGP9nt8/XK/rec5nvd+f+77O2WvvaXiFhYmJKQX/8///
dhOt142eRAK3zrfyxhJYw6DG3aas+EfgTmNCeZI1GTEkeSSDm+JAwTzL/f9Y
RVHEyT6frLoEsP2eH3ta1j4kcDq9rqP8CZjdd5v2cYckylBzZVt5/RQ4h+uc
ii2WRsfOyVzYm5MMji0I7zWtl0F1+ptPtNOfgUsPEiuqpBWRlnnrj2vJL4Cq
ueV0Ja8q6rWNP/QiPg24lX7ozc9TQyWnjilMtqSD+PTZrbZbZ5CffMKwX1MG
0NA6sKD8QBOxAesDGXWvQOiVGfEdUBe1qTd9kq3OBOmje8v/o+ujOI1DHm3l
WUBy8KVy1oQh2q27WLn0OgfkOmZy/uY1Q1MGppaR+blgru6xXqf2RVRgUrcl
mpMH9IhHIicyLJGC9T0tjfQCUPLx8h9NPTtk6CYymBT/BhSYdqpv+3UF8SgA
lKhZDCILPv1U8nJCCzXepPGWYvCc8Ipym80FVb7tq7jRVAJ8rk8fvmjshv5r
Sd9MrSsDrZh58fNWb3RQs1v/l2I5EKduWCU3+iCuT+yvTlWXg8b6z5b70DXU
2e6q2VpeAZZ7y3tHv95A1j0ysYuvq4DETWC3CQIRNHWaUJKqBmcGBv9zcQtC
+weSZcLzq0HJebc79mnBaHZoc2BHTg3YGHJOZhIJRX6TrXvPpNeBkzfvmPod
vY3iaFZlifEIPK/9TunniUQPBVOUG3kbgdhWU03NTCT6VnJvj6xhI5hfrmr3
a72PLGetqLuHGgGhXig+PCEaeVhye9Pm3gHC+MKDBbNH6LGKc2QyXzMw1Oh3
YfmSgL5/N7q6zaQZqHSlaXZffIx2+wLdO8nN4CmjxSZ+5jHKKRIScBdvAcHe
oYZD3Imodk9LGpT5ACxU/JeSXJLQKLN41U+zVpBqdMvyY9gLtD+VN8XqeSuQ
Dk5RCTmaipxPrwZ3jbWCV0f50vYNpSKaT6d6ncsnoGStvyS6Ox2xTYV0xAS1
gWCXIf22qAwk+WGIKpvWDiRrZsTD4jORUos2l/q1TrBjzFXRwLIA8fVkf5VP
6AR21wKGL7QWoKkJlidHyzpBWIn7x71yhSiGqW6PyFIn2Gt+11pP4DWaUDoi
u+DTBfYoBvc3DhahyFJuu1Sfb8BYqvfMmmYJ6ktrr2B494A8KDdRfqECFRRJ
BC3E9YDAsjLztLgKdOvtXThd0gPsnSa+0Hsr0OEh5c9diz2Aj8Xv8ZhtJQoW
ej2W690LzmQU59b6V6F9Dx9ym3n3AQPfT8teDTXIK1D3UpnXADhu1mX1LLkB
aXso5xJiB0AbWgtS/NWA9tlL/rV+MwAinbJFVBQR6jvHcZMwPwCcNG88vT6I
kIpIQ5L11e8gXlL0VldsI+KqPtrJ6TgILkm8/ek3/g5lYTzAynIYtAsVrznm
NiP36isxLPvGgAZTj5RzWxu+3wZizDJjwCkqq4aT9TNa+K5cunVmDIizOCWO
Kn1G19f4+9Ydx0DcO69f8kWfUaBq067l/DGgmqQUHfmkHYU37379++Q4+JDb
GvPEowOld/Z/7lGfAHI0V9/+/V2oZ+YcMcd+CsRtS9au1u1FUcf/MVVenAGP
P8zkafP9QN5HWp80Zs0B1lcJmYdNphDj8tWiEY8FcFk25+mTrl/I/q2pfqwj
DZxj83LeNTKPTlZSyCKeNDAku3LckDaPWN4Mdmf40QBxQ3MghPMvepluY1YW
SQPmlUyG6MRfNHHH0aYvnwYUHfY//X7vL7LX9nUT/UsDqw8kLKWEFpB9/5PI
XL8l4Fe5LtxbvYAuL/Y0NkbSAfjeoeH8exHZBcaaVcfSQU+7bPLs0iKyYj0/
/yaJDtQkc1kc1xeRqeD7nWm5dAB3cM9o8tKQlmKJ781WOpCOepZSeYyGjoU9
lDpNwECdzWC8jBcNrfBpJJVFYCB86MOoxR8aoiczHy94hAEVG7Hin3QaWhR/
2/LyKQYmz3tzuK3T0KyMzFJcDgaaVlP5LbmX0Ii5uK7XRww0WD4y65JYQk3p
m6xHOf+BLPaJDn77JfRAusozO/wfKLA1bXj2ZQm1neyfvhHzD7i9y0wZ715C
BNll67OJ/8Bw89yPXYNL6K6i/PmprH+gTczmWsD0EgpSr5IQb/kHSElfkuo2
l5CrcdX4C9ZlMJYuaGR/lI5yTfvN3bmXwbOZxnjiKTqaNl/+qsy/DHS/ij/M
VqCjS9byDUPiy4AjIC/i3X90ZO5YlSKivgxm2avdAi/SkaZ/ldGTW8uAmede
//xtOrob2N/mELkMUnXGB+ci6KgpeFlNJnYZpHH/oo0+oCPVMHnp7jQ8/5t4
haxEOpKLquIhoWWQ1G/MXJBHR/ufV7U82MCfuzeMpHyhI/vUfmUr9hXgM4HC
qF10lJ6+XHaEdwXUDXU+399HR6JZ8i/bRFfALOsW6dYPOuIvqgohKK8AxRr0
N3+ejlgaq+TvBq6AkHIDo1YeDOXzkrM3bq8AF8p/bS/5MGRo6SbgF70CCrW/
bHgIYOjVP7FFl+d4Pg9N9+GdGDpzNCpfr2EFaE5dOcUliaG5wCnKpw8rILpS
y9j7KIYet6pGqH9dAS+b7w98ksYQ9fLSZdmxFbAkXzqqK4+h8BSrXTtZVgEP
hVaX8x+Gjs1URj/mWgUfdxW1JWlgqE+WxODlXwXq7+/+DDyHoUNdLX1M+1dB
4FkRcbIehlo5j8dOn10F42yVwdEXMeRlcn/TVncVcObmjbBaYYjyatLtu+kq
+JyaseZmgyEX1aRz7U6r4I2rlSbpMoa4bjCzlt5fBd+9+LZnXMVQaZOl95H4
VSAWh7RS3TFkwVc5mpmyCloUEC3aE0MFBa71Twvw/C1FPgrXMGS82ixJKl8F
JQrn5ZavY2hdY09yVP0quHr7ydlsXwydn+i+cfPLKghXM97XG4Ah+vFjU/96
V4GbOu+YcRCGnt+MNPT6sQpW7Vebm4Mx9Ieicvzy31Xw89/7qeuhGEp0fPpi
eHkVjAqEbqu4hSHV8kVuUyYGWB04qEoNw1CMXvaMFokBMlo+PxS8iyG5F0xm
TSIMIPBF7YXgPQz9mLVoURJnABaXxDj2cHzeChWnKiQZ4P2zKnsqzsfDt788
dooBjJxyCBURGBroduHLVWKAs9SgoOuRGLq1tzlk7xkGuG1jW7bvPoYkPHf/
STnPAJh9bG4Tzp31/pYCJgwgNX7JyCAKQ/5c3Z8eWjNArYpS4VecxcyPKnA6
MkCaa1KBSjSuV1ZE9i0PBri/SdR+hrPX0rgAw5cBxnRkb83gTFFTvnMtFGe/
Ho0DDzDUGJO4OBfBAJ5xlYkGODsPL9g6xjIAPQc5u+NMktTpGE1igMXxmQp/
nKc500it6QyQN+Fvch3nWirNuDiXAew5YiIu4fzovUZSUjEejz1LA+LskJEy
dKuaAQL6tB/y4awYOr/bpZEBmBtETb/h9WyzVrc3aMXZpvXZPZwnTydmKXYy
gEUL90VJnKspszN7BxhANrQishHv9+E/FSmuMQYweZl5QAtn+544T9pPBtAa
a5V8h89LrpRaOviXAR4+Fk84gjN3rOK/pmUGyOxuuhiBz3vM/aFiwRYDbAR8
DujB9SjXGQ9O4FwDw6c958k423DeZ3UQXgP7HXKaLuP6nqIOa5zfswY2fmEU
X1x/wnvpKJlDa8Blwqgk4A7u55ABPnb5NfDB1iPRGPdLhJWU8R/VNSDJ7DF4
BPeT1elbT3s11sANV1Y9egiG2P8d2p1tugZWlWRqzHE/DnYHXYqxWQO94S+y
6YEYKir5munruAaSAxXKw3A/m7v7HdH0XQOsmK2YJ+73wqkPCjNP8Pw/Br8l
4vsS1rQjuPPFGqDdrjerw/fJNN0DVWetgQyGyWgnvm9MVkIa9yvWQBxD41eH
E4aMuq8YSfauAZ6uUzFCtrjfSmoSySNrQHS3S0cTvs/rMbyDjKk1wPhnsWJr
gaFs7Qq7z/Q1IFhnsOlmgiHGOzYPN4F1sD31VWSXNoa+ppmXGIuuAyaySTub
FoYybxbSlfevAwcnW6rEWQzpKhoH8cqsgw7uNzH/AQxlFL+KLDJaB22N+Xeb
TmJIK0391UL8OuDh++pJFcLn19o5cvvZOvjS86enm4yhD4s2FMHMdWBa8Uyr
eBuGzp4JeqhYsQ6yXz0yPcWBIfVfZX63+9cBRAf+XqLTkfKpA+cFRDdAqptK
iB1+fgdZlYVn798Az5Iarx5vpaPae2rvFI5ugEuUZtu5JjpSHLCWs1bdANsy
2K2lq/H74GaiWLbdBjiKxTeMZdCR9EdOTD5rA+j18XBNedOR58KT422vN4D8
dEd7wlU6KhLZ72pViccLiM87cYWOjrrB0VsfN4DTwe4aeXM6kiQFfPr0awME
Jlv6u6vg94vl7AvLY5sgaFvcZBgbHQn/bdcIrdoEm5T3DKd7SygpDl36hTaB
Z3SdLfHmEqLIlAYbtW4CyVnh9eTrOAc8LT30fRNUd9YuRF5eQjtY7Xd3MjYB
1bvaKRUuod3Cy5iY6hbIxn70FC7T0EGwN6upeQvobiHXiIs0lD3O3yj1ZQv8
s/qmIqqHP7/LMZTYuwXaKbKHM87Q0KFPv/muTm+BRX9Wt6DjNPT4sGUDkY0J
7m2SfeHDRkOJZQckSw8wQe6jUmPXChZRysdaZlYXJnjfxWSYML+AFsaKA/66
MUHreN+pyIkFpMnIXhzyYoItkyfk1vsWEF0qYbzcnwkyDHxH6tEC0o2/+s4x
kgleGXfUz4ldQCzWomFtOUxwq1+6+LD0ArpKu8kUP80Ez55CdJv2eaS8S21L
zIEZnlHmscy79BtpWfX42DszQ8G3tLA49d/I+JnT9Cs3Zlhsttxlv+83uioS
8+XgDWbIf7885N3ULEoRHHomFY4/z1TYds1xFq3w3lBQyGWGP90+vPl45Rcq
38z11ptjhkEFasJnTX6iI2N8Uzevs8CWJv35u0yT6HUu2bnChxW+KTemL8/2
oOAf1t1Hwtngy+KZDEv4EtE1VV5ciWKDCy3PjDrZ05FbiahTWgwbdPCTC9h8
+gxZ3h1ikJ+yweea/s7SFfHotORFcUYOG/xQuJXXtFcfrd4w9m79xAZldgV8
UuHKAL68Otuv8LJDPe0Gft7KGuCtrKCT+pgdtktvHZ2W+gKYPB6HPkpih8Sb
Cll3Mr6AmLSFslvP2eEe312jHYIdIJ8lT/RyJjs0iRcxdNvoABOtIvOHytkh
NT0366JDJzA2WYst6WaHDwp46gIOfQNyHg19zWQOqJNC1pGO7gXNaTu4K4U4
4BGmssNdqBcYdvmCnB0cMBNWHID0XuAhczz3vjgH3K2VfjXCqg/kMNICdE9y
wMZQY3qZZD+gRISJDhhwwGWoO6vzdgAw0s7Y/47lgO+uePyoqxkCqssq8See
cEAJ0U2m0N9DIExXvskvmQO+OcMENkWHAWHjsDjbSw7YfLxg8njoMBC02D4h
WsYBdypYvzODI0BaYPCSbi8H3KEfkRVY8QNcu9od9/g7zoa+V0cmfoDKpvZ3
gyMcsKiwpHt5+yhQ8UZ7naY54N/71EBDl1Gg3ZE5HrrMAS3uqlF1eMbAlUjP
SyUinPCc6HjY4/IxkDvqHLe8ixPasMRcX2odA3Ny9u9UxDlhs3yOtNTIGPCh
Gu9tk+SE/oIfoiTZx8Et9dPjE0qcsCS8KdDMaBykrLNdErDhhBuXHHLMZ8ZB
p1eKnd9LTljOGSYvT58AkrF8a9ezOWFwoGroY5ZJcOdN+BOffE6Y/Un63he+
SSA379PmXoo/d3C/9U1qEry4qiPj0ITHl+8+fvDyJHB1XOc0mOCEcVkif6Za
J0HzPZ+XutOcUNMrLQ/rmQS7smaUz89ywq1nB3nHxyZB52SPtyYNj7dbNEpz
ZRLIX3o9pMJCgKR9ZYUGB6YAm5XNG0lxAjwT2tl4wX8K2AT1nJM4SIA2TTJ/
mW9PgaoU7akDkgSYqNX241H0FHD9Liey9yQBdtgxOhRSp0CX6fY7wuoEOJD3
5vnbd1NA0vfeLkENAky+GxSs2jYF7j5ZqyJrE3B/Rr9N+jYFFHp+zvEaEuC9
ID357xNTINWg0ZTNngCN4p6ORjJTwYqX3CKzIwFqr+ikdhOowDC2MHrLhQAj
5E9dXd5OBexfkxoZ3gRYqj9LbNtFBW7nvSUXwwhwvPCRTYo8FXy4+rN5/h7+
vgJbXKEKFeyJtradu0+AnvZPo5/+RwXfPp17/DOOAB/1L9Qs61KBoqb45kgG
AYZVdUiNXaYCmVdtkR+z8PiPr35Pc6aCY1s+/CV5BCi2Z+mFojsViFe9l7hb
QsD3Pc5s7QYViPK7lblXEGB1YpTavkAqEPIUUDWrwePFnDh6IIQKuA/h9+87
AmTEO2KV96iA4w7vD/4WAvQ9bt2qcZ8KmEYrnDdaCVB8b/f94gdUQE/kCOns
JMCkVxIfSQlUME8rItT2EODHyJXTxEQqmNE1S3g1QIAS/hIPB5OoYIQjN9d3
DK+vLLOBK5UKBuz1T9lNESB5yOWpSzreb8PK23MzBFj/9Zbay5dU0L4jQ+vU
HAF6fxOrLc/E5+V7rlt0gQATLpRupmdTQeO3RWsOOgH6D3TwOOVSQe2xlJm/
ywRIVFKb4MingvIo9Wvf1whw0Yv95q0CKngzPbvRtEWAbwu1hroLqSBfPSGy
kJUIr1m4rDEVUUFmqhJ/IicRbu5dHuJ8QwWpjMkXodxEqPerNGgG5yTTBxIu
24lQvC/m+8tiKogvlSkz5CfCm3pFdLkSKniwbURFWZgIf25eac/EOdz1XuuB
nURo3qVkN4dz2IejRtv3EGHL8vQb3lIqCBbvG1kRJ8Lgn3p1BJx9Q0KcJw4S
4UGL2Nsj+PtegweXPksSYeCFJfYHOLvKfb1ZcYwIC1mwswI4O8T7EdJOEmHE
hvp/vng9NvN7EiLliFCbtX6zBK/XXLt1l89pvF4Rr4CPeH9G2V65lqpE6H47
vbzyNRXosoicOqtOhFkNqCQUn4eWzbu3xzSIkD3C20ccn5d6rYsWRZsIv3+O
WknNowJlIXI3sy4RPiohaWA5VCDnU2v924AIvR6qWB/A9ZDusJ/pMSFCgQlf
eBzX60B42UaONRGKfc670ZlGBWITlpFxl4hQeJ/rF6cXVLBDlY0/6AoRvvVu
WR9IoYJt/4wldN2J0DPl31+tJ1RAMNwolfcmwsH5nXWa8VTAUpSlsvcGER7a
F+Qm9ogKlq/8M6QHE6Gpdka8bSQV0N6ljozcwvWsAHwf7lLB3C5N5493iZBl
8sM9jjAqGOt9ejPlARGu3jB3Ewyggk9nFXPhcyJMuVlRf9WJCvw/hhwLTSfC
y70F3xXsqUDiXHP520x8/odPu3y3wvU8r/fu9Gt8niEXmJ8YUMEZoyuDMg34
PP+c41lWwPehJ9/uWhMRWtvvK7M7SQWvTBemSz4QYdnv94KpR/D6LYKWjn0l
wsccs/KP8fOg0S6O5/AY3v/whb2eW1PAa6Iv3mmKCAWPb+/zW54CexxERbJn
cD0XA9KN/06BEKecA/sWiTDo0eHcqB9TQNnjreouFi5oQPea/1WLn39Bv7xI
+7lgmap4VNjVKeC4cWxZT4ILliYnrXbZTQGh0Os3Y6S4oED47OVlkylw/fbW
fW5ZLjidanj1HZgCJ+8LvmLX4IKZc6nixqQpUJSo1rvqxAWZ0g8t3CmcBDaU
CEsFNy7IbZtqJJU2CXhT2sd9vbigVFh6f37cJHB7YTa/5M8F4z2NEiR9J4Fk
pjvn/H0ueDf08gF5lUmQVZKsOFGA5z+3OvL0/QRwXbOUHHjDBT9bhc29LZ8A
0md37+wo44Jbbq+qG7MmQP3Ay/WaOi6YzVY5YBUxAbq3ClDcZy6o9mrNdVFr
AjDrNpxV+80FE0YVWFJbxoH17IRBuiQ3jN7OZpiWNQbEZbLUE49xQ4e3wqnh
8WNg5qbTqQcnuaEoa861C6Fj4DppTsD/NDfMFe3dE2g+BqIVsH49bW7IFZpx
1pA4BmrCOW22XLihpQzH4cNtP4DgPilXmzz8fdlEN8A9DGLjo03zXnNDIb1Y
0eDvQ4CL5bc6vYQbCp4Uy7ieMwSYxvN2RtVywzeW+y+mqA+BudRDXyo+c8Mi
x9XFxhuDoHnHPmneeW5YN4y5oK8D4Dp5x7+6kzwwg73eqca4F3xjJtwWqeeB
OQpBTX9cP4P518k2+o088OfB8KHPc22AaHFEKaKZB9LYxkz0PdoALNWl09t5
IAiUPavn9QkU2T9x7BzhgYzIdjZFr48g+r24TsQWD8ySSOa1tnsPztxTEcDU
eaHnXI4wx7ZKUEG4lt35iRdGhXRfMs7IQ1KCfK72LdvgYKtxcOapHvQ35p/A
qb7t0ETvn2a9yiTaFxFx3NOVD/4RujHqsu038umTCbF344OVb5Sifkv+Ro0H
JtpNPPggk10OxULzN7JpVnZV9uGDcXdO+myE/kYpzLRMYiAfHL3rv9Ni4Tci
B1vtyLzPBy15+duNvswhFu8T7N9z+eDzzZKkBL95ZND4w6g9nw9Oxg+/DYmb
R+nbH7xEhXyw4UtbglHBPFItmgbZxXxwmRhOr/kxjwJ/Pwu8Vs0Hfwkrs/qd
+YtoVzgWeFv5oG3oufNf2RfQhMXQoPpPPmg8HBvL0byAeFTMcux+4fE96ryD
OheQ7J7uayG/+eDvSqH08eEFFDnZxlP7lw/mh1wpD8MW0HG3WtUTK3zQinLG
e/vBRRR8M/nlHiIJf19YzezuIhJMN3NdlyTBjOT0jguyNKR6u1tux1ESFG+9
wdUPaMjZQY9V4TgJygM+DX1tGqqX0Hx+7RQJfsibSSbb0pB9iVzHrBIJni+3
/NgbQUNFTYInv58nwbL8/iN6fTSkRe1erfAgwXHfVrUhxyV080hM0BsvEhSe
5uts9lhCpd5aG7k+JGhb71iY5ruEdm7WMz33JUGLy3yqh8OX0LxQNuftEBKk
ZVUYRGcuocda/oK6MSS46pk6FvFjCX2KOflUM5YE7V84+d2lLqHNnjmKWjwJ
Rj37KeYzt4ScLl0SlUkkQebXbOlijCV0OlB7345UErS5srL5SYCOxgpET0y/
JsHkf0Ns5zXpSJDWVzr6hgSJltcjzS/QkbZCnMz3EhKUfTixbGJER2XN7Art
FSS4rK2bsM+WjsJH5lVL35Jg6ZOjsURfOpLa/u58yBcS/KiucXsrnY4umQR1
+H0lwdnaINEd2XSU+ExW37uLBBMq+x4dKKCjrUP5Rg69JOhK5WLaXkFH32CC
hfYICbpTd20cwb9v+/tccRb6g+vFy7CN+U1HkTrzTqf/4nqFqClvX6Cjp/v9
nGwWSfBaXYxOCP79vbIv0jEbI8Hm/zg5D23Q0dLpQge5TRIsv9FPesWLIRYB
WQcLJjIMVVnNekLCEPnP28shLGTYHu7q5yuIIem0r/YfOMhw4KXxd5ZdGPJg
XbIz3U6GT+ryOzWOYOjmcJBdIIkMid70ewXHMPSggs0ulZ8MN7cu3988gaEC
JyHbaWEyHJrmuuepgKGZz4rWvmJkGM7m9uLmGQytZDZZpYiTYa/qfpKJJoY4
Q3SsGvaTocDM2dkd2hg6KG1tyXGYDIsCtAa99TB0+fGti4knyFDN+KLghYsY
8nEnXqw9RYYZW7XCSZYYuq0Rb/5DlgxFewVedFljKH3lldnB02R4uVtXjc8e
Q8VdUmbaymR4/2Cds6ADhlB+hamHKhn+6JLfT3TE0A+rjyaV6mRYbR/uUu+C
oT+y+iaDZ8hweCR7Pegqhta3fTfe1CDDvPLrcpLuGNr5btZIQ4cM/YUaSgy9
MCSZcs3I9QIZUg9obbR7Y+j0tXXDGD0ynPB6zip3DUPmB7YZ9hmR4fGqFpOe
Gxhy3kw0YJiQoW6iYyHRD0N+/XsMdpuT4Q1dl+6j/hgKL87VV7cgwzAzj04Y
gKHE+yf0Ha3IMPfq0Tz1QAxl2dfqRdmQ4Zr3jcungjBUrvSfXpEdGWY5S7Hz
B2OoWaBd95s9Ge7T03s8jnP3H2Pdfw5keI6PJpR2E0MTH0Yu7HAiQ3E3kKAT
giFamuMFVRcyROP+23/izBzw97z9VTJsS/oX5x2KIT5D//Ph7mQopHNo72+c
xY4wn8/3JMO/R2PfG93C0HG2KJ0ObzI0H2oPLfj/778jZB3aNTI8pvTIioaz
buUzbSFfMqS4pdlJhOH9L9lvCfqTIW+yUoYOznekJcsFA/F833eq2OGc6r7o
LBhMhtverSo64lydX71LMATX71baC2ucv/0M/SZwC++//6e7Bs5/9mtGCNwm
Qye2hmIxnAn225QF7uL9vdnlNofnF0/rXeAPJ8MPw+TsPJyVh59n8Ufifhmu
cTDH2UzEwYI/igwxGVC0gvfnY3pkO/8D3L+Jo5FROD9IoL0nx+D1ffzJ4MU5
u7PGnxxLhlXx5Vx38Hm94w07So4nw92VpQ2/8PkOa2tNkB6T4bVvNrvUcf4X
sf0pKZEM3c3PHn2A68HX0qdDSiJDFz8W7COu3xGWVCZSChlGcKsHYbi+dsFS
rnypeL8rZ38cwPUPrFnazZdOhqmtb1olcH88+Vfbvf0lGUoMTD/Yhfunzeuc
yvZsMrQsKW4Yvo6hqdd8tG25ZMgcefl0Du6/rdn+7G35eD0J+wov+2BI5ooj
37Yi3P86JTGlnvi8L96Z5Kkkw6d6j3bl4f6vTtRO4qkmw1pyhxCbM+6fbtIF
nloy5OB3UjPA94Wgm17J3YDv5+zlc034fvmov73P9YEMd3Z0bG5Z4PMKvavK
1YrvZzqy+mOOoZx6nSViGxnyx9Mlv5ji85EbtCJ2kOEMIV/BwRBDGkeWjxP6
8H36U7+1oIX379wwxTlAhoMj0y0+GhgKyrqXzDlIht9OV6j+/A/f5z0CrJw/
8P3i9nDNVcXQDoETfezTZPjo/PXnD05haH7dNZj1HxnSXbcdnNqBoYvmOodX
l8mQsdNgLksY34+yI33zq2TYSVfysRTAULLr3PHBDTLc8ejUWAZ+Pp4ZcJ8s
ZueHWcXHVV5s0dGzMi9tG0F+yOYQO2EyTkcc2w2WjYT5YUWw64DZCB15uZ7I
PCfCD9WHLu/R/U5HGntpGzK7+GFugP8H4S46osVcK+E+wA+fO2TfEWukIy1X
X0qNDD+Undp9lO85fp+0mLYUyfFD+v28CsundLR7r7xPpgI/1H2eL/k8Hj+/
+5c/P1Lmh984xFnWI+ko9WxAqOMZfpjTDcMUbtARXSx4mt+YHy6RFIte69BR
Rn9Yucc1fmjVeLK9a2EJTbdIj8bc4IeawW8p52eXkFT5KPGNHz9cTxyYqZ1c
QpWxKrZ/g/hhpWGqjm/fEmo/t0r0ussPw77oZ16pW0IrdV623k/44dO/Hayn
7iwh3XRbrutV/PDLCxrbH64lxHBWsQtY44eRvc8eC24uIu06SXhjA//87//W
6rBF9IyXIua9xQ/lBDsTLP8sIuVS2qgzqwCkbtCiw4cXUchajq05twCs/SZJ
e1C7iFhjyLbyOwXgO5t+R8cbi4in7Kc1dloA8gmAWy+nF9CujThL70ABSCPN
+kxe+4umNaPbCGyCkGm0+IgS2yx6f7Ad7I4WhAtKiopN8VPIlZDNtpcoBEmO
hjZKPqOoKn2/yvs4IagoG1DgsNyL9m+dCXqXIAQ/lYU/e9LXi2Ktr9SgJ0Lw
Wynp4lBFL3LekS1bnywEMzW3hWRd70U7Hh86VpEhBA+qlpyUW+pBwRGSe3JK
hOCXoallU3o3UveQZo7uEoI/zb0lNzi+oQ4l5RZDkjDkFWnPcx1vR1HVN9Y0
+IWhu0wrl1hhO9KUfXNCSVAYqnkeiWD1a0eNx8RT94kIw3Y1AYoJbzsq20vw
o4sJw/qbYoH7lT+jJM4eiURpYYgeOhXefvkJXe6++uC7njCs/PQwNCf+AxLT
z2r6YiAMTybUrate/oCG23+svDMShkNmVdc4ZD4g0w8GV/LNhGGsm/+tA30t
SKtGXjXIVhi26pY6Zom2oGNprAuinsLwylp0cGXpe7Ti+szILkYYOji5XHdi
fYceBb7iDI0VhpLByRwJfY3oYFRB3Yt4YficVytxJa8RmeTV7RtKFIYDlkwD
3gaNqPTn0JJJmjA8nJrGoqyBkJuD6GOdYmF47E9tRtvLejRm/aJH7pswJLcl
X+gvq0J+7lmRJj14frMnk53OVWjbzdfK1/vwep08IuHuKqTy/G1myaAwdDx5
Xvp2VCVKGRy5LjUpDOU/3dGMda5Axma7BcXpwtAiU8di/HQZatVPM+EVokAZ
/FaX8y5Crs3UBk4KBX586LI4yFSEeBWkJFh2UGBeYY+6XtxrZLinhvFvFwXu
2pm0lV1eiIb+fEsbO0CBz58VsgmwF6D5KI7ZMlkKdE2sq5xrz0Gxm+eNiuQp
sG51yzLWJQed8kmoz1WkQJ0SujIbZw7yvyj26IUKBVpYenZZnclGLBKnZSLO
4vEvelK+tWYiwRb3EAsT/PMGX9W6VzJQtUL5T2MzCvzdrcsW8jIDWRQy9PUu
UqBDX7Xn6vkMlJoQse+MNQX+uilobp+ZjiTsM1qPXqHAar4HfIyaVKS81UNm
vkGBP2Kcqspnk9EPn53Ba74UaNnTr77fPBmFTV+iYv4UePHjQST7MQl9+DJf
NRtMgQ8W+ApV858ivecE6567FKgvTVD6decJsldUzsl5TIG+vJsTOffikMuf
L9hwIgUyvz/UbbQVi7wybM+Qkylwbt8vCWJQLAol3hkPfoH317XHWGwxBj3/
/mmnYTYFnrcw3hCjRKNXDy1dI3IpcPPnNWmSVBTKV/tTXZ9PgUH0+nqBM/dR
TR6f6aE3FLgy4bLf9FYE6gswi12vosCFeUOW5sN30cjRXz9O1FKgeCHrN32b
O2hqPPCoUz0FKu6khC4m3kZL51Lbuhop8IzfeImHQBji20Fly/lEgfX9Z26U
3wpGwh2+RsOfKVBAH7pvPxKEdt8mvCR1UOAeIrdYwkgAOjp7BAR/o0CrojxZ
QUs/pF3rE2AwRIGhB+npnNU+yMCDrTV8BI//8cLfglfeyFw8Uah+lAIb3S1X
A556IaeomrKDUxTotmvp3qHj7shDVYfFapoCD6aejkSaV5EvbVg/boYCR5z/
IxR6uaC7Fkzza3MUGOxRGn1n6QqK3havfOIvBRpW7xb0vOCA4pv2RTsuUuB0
sb1EcZU9ypDUlOjCKJDO7chwCrRFuT8GfDlWcL/NxMfOF1uh4njXFiUGBXad
zLUTJ1igBsZD++xNCrxX/GI0gWKCWor2lAwxicCE2ey8b12GqN2+ZIuPVQQe
ebU82Vigj3qE/tPVYBeB68eaVBrOX0BDbT3PgzhFoHDN637LO+fQRIjj72Ki
CLxy42et5dpZNHtyRXGaWwR6nlzfhsrV0eL0/cid20Rgk/ebNV5bgFZSdvbr
84nAHqenTi79imhL9/WBcLII7FA+guK/ySAOVnC9TkAEVgZ+dHigL414qzqb
FoREIJO1Q2jJSwkkcNWedFAEzy+i7k/P24t27qHbWu4UgQ8a/0l9cRdG4t33
imJ3iUAvOTWRUwe5kWSE8EbLHhHYGO3x13b7csMJpTydtb0icH7Ii9AfSG74
H+fZ3vE=
        "]]},
      Annotation[#, 
       "Charting`Private`Tag$2176#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}, {-1.3467948060844004`, 0.67798504215147}}, {{
     0, 6.283185307179586}, {-1.3467948060844004`, 0.67798504215147}}],
  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, 3.750690749580523*^9, 
   3.750693909482127*^9},
 CellLabel->"Out[13]=",ExpressionUUID->"7bd03a54-6a9d-44eb-b570-3bc3fa9da345"]
}, 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[14]:=",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[15]:=",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[16]:=",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$6632#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$62625#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$135792#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}, 3.7506939742150917`*^9},
 CellLabel->"Out[16]=",ExpressionUUID->"92c8e5ea-b88a-40ec-a259-0579ee48f768"]
}, 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[BoxData[
 RowBox[{
  RowBox[{"integrand", "=", 
   RowBox[{"TensorProduct", "[", 
    RowBox[{
     RowBox[{"\[Psi]", "[", 
      RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], ",", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"\[Psi]", "[", 
         RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"1", "-", 
         SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
       RowBox[{"\[Psi]", "[", 
        RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}]}]}]}], 
    "]"}]}], ";"}]], "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}, 3.750690773035734*^9, {
   3.750692468554084*^9, 3.750692468646369*^9}, {3.750693013309111*^9, 
   3.750693013625947*^9}, {3.7506934767102003`*^9, 3.750693480731285*^9}, {
   3.750693813632524*^9, 3.750693823125329*^9}, {3.750694217437934*^9, 
   3.7506942255930157`*^9}},
 CellLabel->"In[23]:=",ExpressionUUID->"2581c4b6-997e-4d60-9846-58c4feea8940"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NIntegrate", "[", 
   RowBox[{"integrand", ",", 
    RowBox[{"{", 
     RowBox[{"\[Theta]", ",", "0", ",", 
      RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}], "//", "MatrixForm"}]], "Input",\

 CellLabel->"In[24]:=",ExpressionUUID->"44f971d5-4008-4e4c-a563-e0deae41cfc9"],

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,24,5,
   26690440541181691210,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.750694243925428*^9},
 CellLabel->
  "During evaluation of \
In[24]:=",ExpressionUUID->"825b10e8-cd8a-48fb-88e7-21321fda864c"],

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,24,6,
   26690440541181691210,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.75069424406444*^9},
 CellLabel->
  "During evaluation of \
In[24]:=",ExpressionUUID->"f3a54d85-731a-442b-9515-db1aad0485c1"],

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,24,7,
   26690440541181691210,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.7506942442355957`*^9},
 CellLabel->
  "During evaluation of \
In[24]:=",ExpressionUUID->"a2cf12b1-34e7-43c2-8e51-5fc93a2b016b"],

Cell[BoxData[
 TemplateBox[{
  "General","stop",
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"",2,24,8,26690440541181691210,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.750694244367633*^9},
 CellLabel->
  "During evaluation of \
In[24]:=",ExpressionUUID->"6cba59aa-d1ba-4a46-b73f-f3798d4627be"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1.8188779876130903`"}], "1.9349026422959088`", 
      "0.16462935152004432`", "0.3162145013786962`", "0.`"},
     {"1.9349026422959095`", 
      RowBox[{"-", "5.415791529368171`"}], "0.`", "0.89176108502928`", 
      "0.`"},
     {"0.16462935152004388`", "0.`", "0.`", 
      RowBox[{"-", "0.48786230182394485`"}], "0.`"},
     {"0.`", "0.`", "0.`", 
      RowBox[{"-", "3.355084167982151`"}], "0.`"},
     {"0.`", "0.`", "0.`", "1.9772849730968007`", 
      RowBox[{"-", "2.383753868702591`"}]}
    },
    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.750694245048985*^9},
 CellLabel->
  "Out[24]//MatrixForm=",ExpressionUUID->"c79e1352-a8de-4c5d-9d8a-\
98e75faa73be"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"integrandexpand", "=", 
   RowBox[{
    RowBox[{"TensorProduct", "[", 
     RowBox[{
      RowBox[{"\[Psi]", "[", 
       RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], ",", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"\[Psi]", "[", 
          RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}], "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", 
          SuperscriptBox["\[Alpha]", "2"]}], ")"}], 
        RowBox[{"\[Psi]", "[", 
         RowBox[{"\[Theta]", ",", "\[Theta]s", ",", "n"}], "]"}]}]}]}], "]"}],
     "//", "Expand"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.750693723922892*^9, 3.750693726839962*^9}, {
  3.750693841896378*^9, 3.7506938747778473`*^9}, {3.7506941944177*^9, 
  3.750694196416432*^9}, {3.75069424838072*^9, 3.750694249922039*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"76ff4425-c8b8-4171-99ab-18f9c6c87e86"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"NIntegrate", "[", 
   RowBox[{"integrandexpand", ",", 
    RowBox[{"{", 
     RowBox[{"\[Theta]", ",", "0", ",", 
      RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}], "//", "MatrixForm"}]], "Input",\

 CellChangeTimes->{{3.750693732636116*^9, 3.7506937577856407`*^9}, {
  3.7506938531951637`*^9, 3.750693853325164*^9}, {3.7506941990489407`*^9, 
  3.750694200080118*^9}, {3.750694252357121*^9, 3.750694253033738*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"3c67200d-8dfc-4a31-9b99-d00019929dad"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"-", "1.8188779876130907`"}], "1.9349026422959092`", 
      "0.16462935152004426`", "0.3162145013786962`", "0.4324669682753961`"},
     {"1.9349026422959092`", 
      RowBox[{"-", "5.415791529368175`"}], "0.4659745057268326`", 
      "0.8917610850292804`", "1.2115726682821257`"},
     {"0.16462935152004377`", "0.4659745057268332`", 
      RowBox[{"-", "0.8265765410839354`"}], 
      RowBox[{"-", "0.4878623018239439`"}], "0.3368557539613917`"},
     {"0.3162145013786959`", "0.8917610850292801`", 
      RowBox[{"-", "0.4878623018239449`"}], 
      RowBox[{"-", "3.3550841679821515`"}], "1.9772849730968018`"},
     {"0.43246696827539655`", "1.2115726682821264`", "0.3368557539613917`", 
      "1.9772849730968014`", 
      RowBox[{"-", "2.383753868702591`"}]}
    },
    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.750693751842778*^9, 3.750693760474391*^9}, 
   3.750693856356254*^9, 3.750693978833709*^9, 3.750694204779399*^9, 
   3.750694255993764*^9},
 CellLabel->
  "Out[26]//MatrixForm=",ExpressionUUID->"35ffca7c-859e-4e79-9e9a-\
98892c98d50d"]
}, 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, 420, 9, 109, "Input",ExpressionUUID->"3b14f5c5-ea64-4841-a7ab-3bf3afa6d1d5"],
Cell[CellGroupData[{
Cell[2784, 79, 292, 4, 101, "Section",ExpressionUUID->"fdbf016a-a481-49e4-a4a1-7fef444f2077"],
Cell[CellGroupData[{
Cell[3101, 87, 356, 6, 67, "Subsubsection",ExpressionUUID->"d797f85f-1f44-4aa6-bb33-319d3a419081"],
Cell[3460, 95, 2131, 56, 78, "Input",ExpressionUUID->"6ad48be1-82f4-4c5e-848e-3955567c1543"],
Cell[5594, 153, 777, 17, 46, "Input",ExpressionUUID->"5c3f117f-edd1-4683-ace3-bed538a9f32b"],
Cell[6374, 172, 757, 16, 46, "Input",ExpressionUUID->"698cc67f-7248-4dd9-a3ac-3bd95fdd14dd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7168, 193, 173, 3, 81, "Subsection",ExpressionUUID->"1c31843d-c3a2-404d-b4d8-27c8c57f9dbe"],
Cell[7344, 198, 1582, 36, 109, "Input",ExpressionUUID->"9af92159-c78a-49e0-9fe6-46bac46cfd52"],
Cell[CellGroupData[{
Cell[8951, 238, 453, 11, 46, "Input",ExpressionUUID->"b2b8f0ba-19f7-451b-b670-3aed0e09fedd"],
Cell[9407, 251, 15879, 279, 346, "Output",ExpressionUUID->"7bd03a54-6a9d-44eb-b570-3bc3fa9da345"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[25335, 536, 223, 4, 81, "Subsection",ExpressionUUID->"ea08ebac-e4fa-4fb3-a8a1-1284274fd018"],
Cell[25561, 542, 1870, 35, 46, "Input",ExpressionUUID->"06f9ecf8-085d-4573-ae5f-34752dd38caa"],
Cell[27434, 579, 841, 19, 46, "Input",ExpressionUUID->"46f7fef6-fb20-4fad-b080-d1185973ae5f"],
Cell[CellGroupData[{
Cell[28300, 602, 661, 12, 46, "Input",ExpressionUUID->"85373ec1-21a2-4953-9ad6-303bee858a0f"],
Cell[28964, 616, 50892, 885, 188, "Output",ExpressionUUID->"92c8e5ea-b88a-40ec-a259-0579ee48f768"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[79917, 1508, 471, 11, 102, "Section",ExpressionUUID->"9c8b182e-87c3-4ec9-8ce3-6172bbb3e292"],
Cell[80391, 1521, 3530, 59, 49, "Input",ExpressionUUID->"2581c4b6-997e-4d60-9846-58c4feea8940"],
Cell[CellGroupData[{
Cell[83946, 1584, 303, 8, 46, "Input",ExpressionUUID->"44f971d5-4008-4e4c-a563-e0deae41cfc9"],
Cell[84252, 1594, 497, 11, 59, "Message",ExpressionUUID->"825b10e8-cd8a-48fb-88e7-21321fda864c"],
Cell[84752, 1607, 496, 11, 59, "Message",ExpressionUUID->"f3a54d85-731a-442b-9515-db1aad0485c1"],
Cell[85251, 1620, 499, 11, 59, "Message",ExpressionUUID->"a2cf12b1-34e7-43c2-8e51-5fc93a2b016b"],
Cell[85753, 1633, 448, 10, 34, "Message",ExpressionUUID->"6cba59aa-d1ba-4a46-b73f-f3798d4627be"],
Cell[86204, 1645, 1213, 31, 168, "Output",ExpressionUUID->"c79e1352-a8de-4c5d-9d8a-98e75faa73be"]
}, Open  ]],
Cell[87432, 1679, 1030, 25, 49, "Input",ExpressionUUID->"76ff4425-c8b8-4171-99ab-18f9c6c87e86"],
Cell[CellGroupData[{
Cell[88487, 1708, 525, 11, 46, "Input",ExpressionUUID->"3c67200d-8dfc-4a31-9b99-d00019929dad"],
Cell[89015, 1721, 1595, 36, 168, "Output",ExpressionUUID->"35ffca7c-859e-4e79-9e9a-98892c98d50d"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)