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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    129801,       2504]
NotebookOptionsPosition[    125735,       2420]
NotebookOutlinePosition[    126797,       2449]
CellTagsIndexPosition[    126717,       2444]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 TogglerBox[1, {1->
  OverlayBox[{
   TagBox[GridBox[{
      {
       GraphicsBox[RasterBox[CompressedData["
1:eJxTTMoPSmNiYGAo5gASQYnljkVFiZXBAkBOaF5xZnpeaopnXklqemqRRRIz
UFAKiEGSLED8//9/gvjRo0fsN27cELp69arY+fPnpZExPn137tzhBqrRBGJ/
ID4CxMfRMTZ9L168YL5w4YIaUD4AiC8C8Qcg/gvE/9AxFreC7PQD4gtA/BGq
7j8ujKz3+fPnLFC9b6H24dSHTf/FixfVofYSpRdZPzSsAqBuJkovsn5oOF8k
5F9s+h88eMBBjt0w/devXxcGso+S4m9k/ZcuXRKHpgeS3A7TD0wrkpToPw9J
ywOmn1L3Uxp+lMYfpemH0vRLjfwDjQey8y9a+fGOWDPwlF+gsPhEKDywlZ8g
d0D9EniexPITGd++fZv3PBnlN6764zxa3QHCAIDESDs=
         "], {{0, 26}, {16, 
         0}}, {0, 255},
         ColorFunction->RGBColor],
        BaseStyle->"ImageGraphics",
        ImageSize->Magnification[1],
        ImageSizeRaw->{16, 26},
        PlotRange->{{0, 16}, {0, 26}}], 
       GraphicsBox[RasterBox[CompressedData["
1:eJzt2V1vVFUUBuBGTfTSSygthYYKNHyVinrrpbcYfwBEJN5gUkyMfw9KKaVw
KCkQKCW1RNIQKg1gEE3QtdN10pMRZCaBzNQ+kzxt5nzstfeZfbPOu/vkD8e/
e6+vr+/sR/Hn+ImfvpyYOPHz1x/Hl2/OnP3+9JlT33515sdTp09NfHHy/TjY
H6592Nf3QfyvqmpHBwbCYNr5GoMNAw2d1AEAADa3/obtLba1aJ7rb4xR9x1D
YVfYHYbDnjASPgn7wv4wGg6Eg+FQOByOhLEwnj4Nx8Jn6fM21dcfyzHq8cay
xuGseTDnMJpz2pdzHMk5D+caduWa6v6p+czafU71s6p1+/cGAADYqt52/9va
++6tXt331j3v0TTeoj5e966v6lv3Vv/uWfWr9LJm9tzMpN+UW3eaVV9s00yY
CzfC7bAQFhvuhjt57maYD9er9c+VcDlcCtMd1AQAADafqXQhTabz6VyL+vhk
Xl/uLX1D6UFmw9VwrVrvL0o/cqta7ztKD3Iv/BLuh1/DSngYHoXfwlp4Ep6F
38Pz8Ed4Ef5Mf71Gff5F3vM8x3iWY65ljUdZcyXncD/ndC/neDvnfCPXcC3X
NJtrnG48r3afU/2s6uc11QO/OwAAwFbzLvrfksWVjK30kSV3KxncUlgOD6qN
vnc1PK42+t7iaaq/r+U1q9VG3/ogx1rKse9krfmsrV+l15R9WDLmkjWXzLl8
yvuVsmdLJl3eu5R9XN7BNLPrhTxX3seUvT3TQc2Xbfr7LWq3JgAAwH95m33K
u9Dt5wMAAMDmp2dlK+vGnu/2+yQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAYPN42aZu
1AQAAAAAAICtrNsZIXRTN/b8xTbNhLlwI9wOC2Gx4W64k+duhvlwvVr/XAmX
w6Uw3UFNAAAAAAAA6EVT6UKaTOfTuRb18cm8vtxbcrOSwc2Gq9V6FlcytlvV
eu5WMrilsBwehJXwMKyGx2EtPElPU/19La9ZzXtWcozlHHMxa9zKmnM5h9mc
03Rjfe2uq15bvb6pHvid2NzKPiwZc8maS+ZcPiWDLnu2ZNIlmy77uGTVzex6
Ic+VbLvs7ZkOau7owEAYTDtfY7BhoKGTOgAAAAAAANDr+hu2t9jWonmuvzFG
nbsNhd1hOOwJI2Fv2B9Gw4FwMBwKR8JYOJrGW9THx/LaQ3nvgRxrf449krWG
s/ZQtZH3NdfY7rrqtdW6/fvw/9LMnpuZ9Jty646y6n8AhCP7tg==
         "], {{0, 26}, {
         2000, 0}}, {0, 255},
         ColorFunction->RGBColor],
        BaseStyle->"ImageGraphics",
        ImageSize->Magnification[1],
        ImageSizeRaw->{2000, 26},
        PlotRange->{{0, 2000}, {0, 26}}], 
       GraphicsBox[RasterBox[CompressedData["
1:eJy9lktLQkEUxyVdCBIELaxctYq+QduWbcKMPoCSRRsDC6JPp+Lj4vuq+EBc
iBtFfIYrpU7nPzQiN++9vmrgp/dx5j8z55w5c099z56HPYvF8mLnH4/37TIY
9L7fHvDNXeDl6THgv78KvPof/cELn5UfnjBHjI3J5/OudVBV9bhUKjmbzaad
iCxmcB9lAxKMmzmv1+v7JvpfG/DJjJkic8NrOut0OlYdfdoCjPXBqMx1q9Vy
7FhfgvX0MUa73bYZ6edyuTnrvPsZQ4Wv9PTRL5vNroTOGIiJu9FoOLT6UjuT
yVAqlaJkMin+0+m0YPEZbHTGQDwKyKtl+ugHjdFoRIVCgeLxOCUSCQGua7Ua
9Xo9YQNbozXI/aHVZ//RcDik2WxG5XKZwuGwANdo0JfrNIhDvFqtHi7zD+YW
i8VoMBgIvWKxKEDrdrtiHbAxiAF8pHAf5zJ9+Bca0WiUxuMxydbv98U64CvY
wNZgTyioI9r8QR/EEvrQUhRlrh8Kheb6sDHTZ1xG84ePMP/pdPpr/tvoL/of
Mcb9Jv7X+mcxP2X+QCMSiYh5VyoVMQbXMbP81I2vzM/JZCLyX+a+BPmPvbFN
fmr36pr713B/rVp7THy/tD6Y1cjF5zrahvVtB5jW5y21B/kVzpc1kecjfK53
Pu7kfNfO+x+/T/70++ob7NUEtQ==
         "], {{0, 26}, {24, 0}}, {0, 255},
         ColorFunction->RGBColor],
        BaseStyle->"ImageGraphics",
        ImageSize->Magnification[1],
        ImageSizeRaw->{24, 26},
        PlotRange->{{0, 24}, {0, 26}}]}
     },
     AutoDelete->False,
     GridBoxItemSize->{"Columns" -> {Automatic, 
         Scaled[0.6], Automatic}},
     GridBoxSpacings->{"Columns" -> {{0}}, "Rows" -> {{0}}}],
    "Grid"], 
   PaneBox[
    StyleBox[
     RowBox[{
     "Background", " ", "images", " ", "and", " ", "text", " ", "colors", " ",
       "have", " ", "been", " ", "adjusted", " ", "for", " ", "editing", " ", 
      
      RowBox[{"purposes", "."}]}],
     LineBreakWithin->False,
     FontFamily->"Helvetica",
     FontSize->12,
     FontColor->GrayLevel[0.4]],
    Alignment->Center,
    ImageSize->Scaled[0.5],
    ScrollPosition->{0., 0.}]},
   Alignment->{Center, Center}], 2->""}, "1"]], "SlideShowNavigationBar", \
"FirstSlide",
 CellMargins->{{0, 0}, {0, 0}},
 CellFrameLabelMargins->2,
 CellSize->{Inherited, 30},
 TextAlignment->Center,
 CellTags->
  "SlideShowHeader",ExpressionUUID->"c443756b-bc88-4e3b-a9dd-ebc17ad287f3"],

Cell[BoxData[
 RowBox[{"\[AliasDelimiter]", 
  RowBox[{"\[AliasDelimiter]", "\[AliasDelimiter]"}]}]], "Input",
 CellChangeTimes->{{3.7645623933871098`*^9, 
  3.7645623946405478`*^9}},ExpressionUUID->"825030e4-751a-4878-b036-\
d7b8185725f5"],

Cell[CellGroupData[{

Cell["\<\
MS-C1650 - Numeerinen analyysi
Lecture 2: A Gentle Introduction to Interpolation\
\>", "Title",
 CellChangeTimes->{
  3.559948400406288*^9, {3.764562310742819*^9, 
   3.7645623509502993`*^9}},ExpressionUUID->"db09f6dc-cd51-4f60-b989-\
d0bd041c5d56"],

Cell["Harri Hakula, 2019", "Section",
 CellChangeTimes->{
  3.764562359512897*^9},ExpressionUUID->"12ad6237-6c23-4157-9523-\
c690d4971f36"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["", "SlideShowNavigationBar",
 CellTags->
  "SlideShowHeader",ExpressionUUID->"65af3d05-da9c-4201-bc7e-fc9881fad03b"],

Cell[CellGroupData[{

Cell["Interpolation", "Section",
 CellChangeTimes->{{3.76456247696148*^9, 
  3.764562482407902*^9}},ExpressionUUID->"9f6092a6-ba30-46a9-ae7a-\
28eaf8d86d3b"],

Cell["\<\
Let us assume that we have a set of data points and we want to fit a \
polynomial which is exact at the data points. Consider a sine function.\
\>", "Text",
 CellChangeTimes->{{3.764562500761406*^9, 
  3.764562573960884*^9}},ExpressionUUID->"a6fb594d-3a5b-45a9-9ab7-\
4c59081a4f4d"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Subdivide", "[", 
    RowBox[{"0", ",", 
     RowBox[{"2", "Pi"}], ",", "10"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"Sin", "[", "x", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Sin", "[", "x", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", 
       RowBox[{"2", "Pi"}]}], "}"}]}], "]"}], ",", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
     RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}]}], "]"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.7645625770288982`*^9, 3.7645626693124657`*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"e2944af6-5975-43fe-82c5-734786dd8e16"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwlmXk4lN/7x61ZEjPPjJAU7SsppYXuUwlFSokWIWkhSURIKpWQZK9kS4X4
RHZJTpSikKiEopAsYx7rmJksv/O9fn+d63Vd82zv9/u+73Ou0bA/u/e4mIiI
SJaoiMj/VnEsW9ux9upmXun/rwt7SjfNibSFu/mNw7PJKnyu5dAUcRbORecZ
fCHrH/YpfkHEZbjytMwnmazNhtfLn0TcAbWrSdX7yPrVNmJxfEQirNs+w2J/
eCLsdVZpvheRBWux8JX37SwIH7LOjYnAsF7xqYGaH4ZNFTtlt7rXwY2SV9tl
t9eBq4/Z0VzX71B44Y1OSVMjnCk6Hio2/xcsypjzyZjTCl+6d8ik2ndCZ1fT
wQrpDgjW4okUHOwGDal00Z3MLji3vDL69RMOCNzfr83gd0Pn8c4rSQUcOOdZ
NH+teA9YJok4X3nPgedXZqgVzeiBjYobtqIeDtTb6Sx6Oq8HxKaecvGKftja
M3OLrmkPRNTfMsK5/WB/MpktmtADOd67Ba/KuGDW1Oihq9cLwmOnM3+6DEDG
zO7au/Z9YMcc9irxHICSqcNzXzj1QUWpz7YHfgMQ4RBZ0+DWB2Eqwd8PhA5A
kYnnplH/Plj8KU3sy7MBWGcWvKg7uQ8sNnZZfuAMQPaxhbdCfvdBJuPoVIHT
IKjEWJwpP8gBdml3VbTbIKSUiVVNP8oB79OuUed9BsHkR7Oh2SkOGLy7vHR1
8CC4N4gnl3pyoOViwr5nTwfhTo7bE9NIDkj/bUl91D0IHw3k7WZ/5ID9K8s9
YSeGIGA/5vmt7YfVBcqUytkhcOPND3bX6wexrOaGhxeGwHveqmrbbf2QnGRj
lRs4BJoP7SzmmvdD+7UTNt/Sh0A3OXzzRmei005P59n0ELxe71vzM4lwY3Rg
2oVhYJze2F8rzoXVdVY7ta8Mw5HOxpx8GS6IVanIFQcOQ6BDISdSgQvJxfF3
Pt4fhp7D4QXrVLnQHv8khvtyGA7k3TqvuJoL9g75j3VERuDXG9dV62y4cGzw
y+vXgSPg6O/YJsjigp1PmFVR2Ai05sw+5pPHBWtxU27WvRHw1K3ZMVzEBUvF
N6qJaSMQGle0qZb4aLwh2/NS5Qi0Gc2PXNLABc2rt1dslB4FmYnPGi4jXFgm
u+ONNmMU1GKdLc7zubA4UuLQUuVReFu0tMltnAvqTy7eVF4yCrdB6esBCRpY
VY7tPKNR2Bzme7GSRQOfYXgv9+YorMuIy1yrTcPIfVGtjDuj0PYxN+y1Dg2D
815VJN8dBX9RYf629TT06ugMh6eOQlHd/XergYafB+aZub4fBdFg/e5cUxrK
kybFV0rxgMtRPr/lBA2lS4tjFyjwoMrjy4I9jjS8zPHQnq3Eg/6ygulWzjTk
veXYTF/MA1py0GuXGw2p3c0vegx5cG9+TsmPSzSErCo8mxLAg8NaZ02uR9Lw
YXVjl0coDwZ3JC+XiaFBeu3Yke0xPDDmiNwIuEfD9Q26pp1PyP1ZHUlO8TRc
3Fq4ZF4FD0avKJz5nkJDkUFj4mA1D/a2npim8ZQGnuHYzLIvPJghfYRtn0GD
m4mupF0nD5pviIbUZNHgZFH4O158DN7tjGjZUURDmmXjgTPTx2BlzMbplsU0
dB0Y+6THGoNSkVMZh0poOHpEt7Rl3hjUR5zctOs1DQdOFMaqbB0DB92RkMb3
NBh5Fe6LvjIGOnPFo658Je/r0/jBIXAMlp/aGMNqJPr5jm3RCRuDwkg1scTv
NGy+qruqIXEMnvu4Kia20LAuuFCOicfg9VJndYPfNJwPabz2690YPLmqciK2
nYac0DFhVu0YrIrOpXo6aFgZqdtj1joGnFLXxLNdNCyIK6wImRiDvHWVfhp9
NNgnNOpZS/LhWtPN4Y0cGpKSxnKXz+DDzIHuatN+GmY/0U3+MJsPOVtD/jtI
03Ao1UoldgEf1J3NcywGaLj39EKY4wo+bOx8pWo8SPKVWegnrceHptspysrD
NOx53jjauI0PA5p5GSOEQ3PGnFNN+DA8ceth1QgNsoW6hw0P86FWUJdnzaPB
+IVVveIxPkyf1fp51hgNAS8v7PjjRN5n9+iez4TFXhfqXvfhg/nX6POLBTSk
z6BSJvz5ECLxe88XwnsPO7Mv3OLDz7RGHS8hDcK0d/4DkXwIrzRRZ/2j4RFP
fdAxjg9aXCbrCWFTg4u2HY/5ELBEdIbmOMl/+Nca62d8yJcdnpFJ2GBlcPru
Uj5MpJkviJyggePTqVz1jg9mGn4b+ISjKjff3PqJD+MvVuy3mKRBb+b90ZeN
fLCK23YhlfCfY8PH1v7iwzrPjMRBwqHZu+ozu/kw5b27dvUU8XMqFS0Z5MMz
GzXx04RbTcWyHgr40GAruTmWcECstZqqmACe6Av9XhPW7C64FSUrgDG9wYqf
hL+tZQpnsASwbXE7NUTY79rpUzdVBZC+8/3xCcKLP1d8E1kggBHv+NIpwnVz
1Lf7rBBAtJyDmoCwl7NP7rCOADJaZ/v3EFYv/qJxRl8A43uqOXWEK6W0wrq2
C4Bn7H0ki7Dr/qBJWzMBZK1a+eU6YeVHHc5NlgJo4w3vMSf8ekC/Za8tYfkf
DWzCjpvv7ag+KYBHrjzrT+T7qZChwu2uAnhxz5JzhXBxk+ki7CWAzYGU/1LC
9otTo9ZfFcCPt1vnVBF9ZT1ExXOCBFB+QKzMlnBO+eFzyyMEkFDv4tRP/DnE
KGh7HCuAO9vTZrkRFrNhmM15JID5Oz5/7if+ZmQ4ldzNEICJQD7MjrCF4O0y
Zp4AdvEC938g+Rg3nHs/uEQArxwc5y0nbNre4HGpRgDnX1ENn/kkD1qanbyv
5Ptv6BUqEY67FLjXtVUAyDf0kQXJY7+yvtYxWgCXDwdG544Sf3endBszhbDy
ocsfySHib7yIVbmKEDLG7JkSpD5aew9VbJonhKM5F40EpJ60AhSSNdcIgdYx
bqgk9VZX4nWYvV8ItSaqbuPdxB/ZhqrbR4SgLk2P5fwl/hxYuV7qhBC09K8H
25F6dh3+zRZ6CkHc4lxjEql35jKT2rZ7QtgdZala8JP0J6lEZmWSEPrOuGmh
H0TvP0MWz9OEkBDpe7G8mQaHh7EtV4qEkBNeX5JF+k+Rcm+3xnchHD6qRa/4
TMNtnv4K2V9CYFgeqvH6RPz5En526K8QFsdEdJXU0DA9bAOvfEwI3KDWjys/
0GAjFSTuoPQPOOd0XXLKaVjz54eh6dx/4JlOab8m/VD6zapgncX/oDYhcvW7
UuKn33eGpO4/aBnxHC4m/VSSt3hOiuU/kBos9THKoeG/znfru6P/QZHwuGJ/
Ig1Xy2f51sX/g89+mSKXSX+3THLBRU/+QY3shjUyD2gQsZ5pGJT/D8pPXHKd
RubDvobj+5Z9/QcfB+SOPL5F+kGZhIszexySeOWhLh6k3yRufTQQMQ7vAuU2
bdlGw83Kup/+D8ahSNVjYRSZb+8GbZQVH5Pfb1+l3LqJhu0GF29vyB+HuVc7
ZlmR+bi1J/eCf+M4/LyenftzAan/NQtN2bMnICNoTHMzmber3kuN6j6ZgKsx
HTkWJVw4OxCt9eHZBKDnN3qCC7mQqbLAybpgAt5aXd9dlMOFlc6o7cr7CZCR
nZ/FTyPzneldVdUzAQW3t3cqRHNhweHe+MOak+CccdTx/mkuKNHVhpcLJ6Fc
st13PsWFe+H4aA+ehNCrxko2clxQ1snx3Vc5CS9Tl9qHTSPsfTdncdMknO6o
Uvgh7IdZ4vZz6oSTsFJLgq/Q0Q9zlMZG1TdPwWj78Nm+7H5YBBpPyt9OQYrX
GPPSjn6IfV8sKu4ogg6zfy10t+PAwK/n3rSzCLJ/HpemSfZ/RsKUwRZXEXQi
PfZghzkHRlZE/s7zEkEL1I7uWr+NA2YRp8tOBIogReekZ2ELOSB2ZPbVD6ki
qLg+oH28uw9OD10SiegSQb6NIiJ1ZH+qp7ZlSt1BFNXdDDxYc6gXjK2/uNmf
EkXf2950aZj3gsWDk12PnEXRIcEN0bNGvXBaJbRmkYcoavsX2zu0phdiFVse
rAgQRZ+Kjn9/KNcL/Bke69eniSLfUsOohJIeyJtMO7ebI4puRhgGLlbqgeW/
GJ2XzouhQ7rDnWqFf+FZGnUq300cyfGkPrm2dYJv65GG5QES6Ofe40uvufyG
c3rrTRKiJJFvnVHUGr8WECYa2PeFTUPNXaYPKzW/QJ1rrN2FZCmkJue9xO3L
e9hgNG/y50Np9Pmb3IlJkedQtX1DGoqTQRvqbOZkhWTiJ9n3N7RnyKLPLHW9
JxcqseL8FU42T6cji8U5Cwz/NOB6UWl/lRI5JPEh8alOcwvOl3ZPqauagaK8
V+W/ifmNVygynOwr5NGKsEcHrF92YjqUx17zTQGtus3sWFT1F3MG5pSvbFZA
holviqsa/uLevUZnl7QqoLf0vmcnWv/izpn3qtS6CI/f+R46/Bc3J2zwk+Yp
oI0RppN31bqxvVXpMFeWgdSX2W+95dqNd50UPfJ2GQNti3do7GD24Pk3b2qd
dWIg6SpL/4Ytvdjtm46fvTMDrSsQj325sxe/Xthevd+FgZa0njgXv68X27zV
c9JzY6DwllZLs+O9OFZ06LGMDwOd7vtc5hDYiylf61mPgxjIyyzM0L6mF4ud
05ZsSmMgx6PjdeF7+7D569Z91ekMtDJI8c+zQ304SSEkGf/HQL4xel5v7Pvw
5swuSHnOQD+yJdktbn3Yp++Bj3sRAwX9ZXFeRvThoePTBmZUMlDdp+ILWz/3
4fZDLc1b/zKQ9vYtsy13cLCcvlWqXQ8DJbxq+aRnzsFr5za4+/UxEFdeIUb1
IAcHdnyQK6YZaJalm3nZKQ7Wci7erM1nIKfCaUbPbnKw76X7yXNlmCjdLtJE
8i0HKyZZOY0vY6K90iEaqev68Wb/hnWzVjJRycl6+RD9fnzKYbf4ei0mMlwS
XH3KoB+XLDGKc1/DRHq9Sdkye/uxffa62t5NTPSmTru/27kfZ5Yrrm4yZSIN
9eSbqg/7sfGfBkG+CxN1zjfNNhDn4kvLQy9muTKRUUnU8EtpLs45ZzyR5sZE
l8LVO5fLc7HqZIlInCcTvW8r+DCqzMXcmSlS/n5MNBmmE6StycVRxl6KZqFM
VL9S+j8jKy7+lTFbu+sZE9UkKnD+POJixaFvOW1ZTLT6uWGgaxoX71wfrtOU
zUSp579mDP/HxblvJddX5zNR5p/o6x35XBzwk7s55xUTLWeY2Di+J0dZhTJT
vxomEuRWVcj1cLGX2/FTM/uZaF+yvW3GIhoHmnBPbqSZaEvVqmkRy2h8d8GF
kzaDTHSSI/b6rCaNC74FnkgZZaLgPqNTjHU0Ht74n8O6SSbqcB917jKgsYv4
sJ2lAoVCHE2tO+1ofOnHRTsfJoUWjc3M2+RA45B8CbsEFoUM/x6VDz5J44yT
M227lChE8bLfTXehcffHDUc81SnUb3JoxSMfGh+LunIwRptCDDmQ1IigsdsZ
mYPFayjUc1dRpTuKxv6GEQda11Jo/8RFjZS7NE7iP7JatJFCkbnf5KTiadxq
/X5/wVYKzQ9xXG+dSuP+tXv2NxtQKNF/b0vzUxqPyzdZTBpSaNPFj857/qOx
alnvPkMTCp39lnBAPZvGBxbK7/22j0Ka7gVoXzGNT03GmAv3Uyg4hG8QVELj
C41zzeccoNBSyXOrC0ppHBOkveeENYWWPw96N1JO44Z+CzOeA4UcZULKJT7S
2Kzgwc6ZnhTaqDp862wTud+w/ZSiF4W2yj3Yw2+m8bVVy/IUfSikN9qi7v2D
xkXpRWqKfhRyM01vt2+j8bzErwOsAAotfOL++1snjfV+xD1hBRK9xOv4c7to
bKXicIgVTCG7rePKR/8SvSOH3lChFHp5/4bbxx4a824q3GXGUIiFWf6IS2NG
xTcT5j1y/ab7jTtoGi8XSxBhxhI/ah6vMRmgsZ3vCidGAoUeq+yeXDtEY58X
w3MYSRTyMHR3njtM42hecYNCMmGnJb9ERmj8wXWHvkIKhRIcTjRmjNK48xlj
SD6NQhZx1w958Gg81duYIp9O9A21aV83RmOd4ycY8pkUmm0eMPWQT/RJXlkx
4zl5vod3lImAxo5tI94zcihkJHRb2U844eC1DrkCCp181eSk/I/oE7PznlwR
hYa3jMsnE25oYO6SK6YQO/9OwfxxGnMZTaJyJRRyKXA/Gk9Y2iypYHophdZ2
BikoTBA9b508Pf01hTa35L/2Iqxfqak+vZxC4ro/zzcTPiDJ+yL7lkKp4X+X
60yS/G19FST7jkIT10r/3CAccvn6ZtlKCkU82ZVM9ts4tcRkWOYDhf4L8T/K
nKJxuYBKk6mm0N5a0/kmhH+sa7aWqaXQnsn4v76Eee4PmTJ1FNqQ5ZiZQpiZ
feqddD2FzriRaUR4BVfrovQXCil9N97WTthw+ZiW9DcKrTm7izlC2O5UaafU
dwqFVWT+JucPfPHJjftSzRSqPemUR84fOKbd1EzqB8nD78tB5PyBn89li0u1
kvosGbDrI/zBuqVw2i/yPP/XG78S/nM/2XlaO8lD8cjMQsIijY4a0zopVMWJ
Gg0jPIut/U2yi1wfn/7NnrCOOT9YsptCt823FK8gvDsUg2Qvuf7WiSQu0cPp
Y8CIBIfoU7s0KI3wdWmzpxJcCi2OjTt/iHDCdkUbiQEKMV3q7SUIv/D/QUkM
Uehc5a99KUR/7riTrziPQqFifzfXE/8OHjBZKhij0Gl7kXgbwm9zl3/jCijU
WFGyq4P4f9+Jo9U8QSGR+pjUr0IaS7yr/vFpiuRnzeiebYRdNJ4FVYiykIaM
166nJE8G3890PJdkIZ/3V5ANyV/WGrOwFCkWslhoqfuM5HPWHU39OBkWsmvn
Xh0h+R3YTsfcnMFC01PVWM4k3w9yXXfaKBLOebVsapDG0xTMx/YpsdBBX+5b
BmFXJ+3HO1RYaOynV94sUk+GGkMTOmos9OxiYJQSqb+hUPfs6QtZyLPA1zaf
1OeR3n02ootZ6D3XMiOom8aV23Wmjy1hoYTiGDNLUs/x4yMO7SvI9cKzcd9J
/Rs7eSq/0GGhrBvyNldJf8itsKzIXMdC/Z23zyi30niOhq7b4/UsFOo18i2V
9JPhxrGPd/RYSH+ouSir6X/6e18+YcBChT0r3xxroLHMw4MrrA1ZqPujMOPV
ZxqfH9/QZG7MQldEUwYV6mi8I1e4Wt+UhUxq10bHV9N4RN23i2XBQrNseMZH
Kmhs62sdKWPJQnkeDgt83pB8NeqhKSui5x41uztlNE4Mnbjfe5iFdu02Fia/
ovHOcb9dZcdYaMmf7+9c8mn8sPFqnos7C8XMuqN9/xGNuypWtYV6sJAQFRRK
PyT1kNcmk3WBhWaAwe6zCWQehenb0hdZ6AYnIF/jPo2rdwhkXK+z0AabBnfV
UBrzX7ranosm9zfXdDzmSePN6XODw++ykFbNoWMG7qTf3qvNy75P/K1rvjjb
lcbyHitkh+JZyCvo1PQcRxrP1+zOc0thoUt9pqVu1qQ/JdnKni9kIba4qLEJ
onFUqPzaqBcs1GxWfzdBj8ZNvq9s816yUE3jccWe9TR2OKiaP4JZaCpzWsIx
bRp7U422HpUsFLBaXZzWoPHj62b5nk0s1MWxZduJ0lh4St/O+x8L1dcER+S+
IPP95TLkMcFCb7LTryIyzx/MUFY/N8VCzmckOyqec7FezlDbKXE2avD7KlKc
ysV+/1JtD0xnI84lzUqjaC4WD6VsdVXZ6FmWSYDQhYvlcv8eGd3IRmX/7V9Y
MZuLrSW/6g/qsdG26jnJ9kpc/MyqXK1/MxvVUPxCPpOLzcbjfnZuZaPYiFV5
8lJcHLZ975EvO9koNcGD3T7Qj9nfS6xzD7FRd+XCtMVv+rHaRPjhcz6E/8aq
rLcn+6P2XS/u+rJRe22ulsPhfnz+vYxSqR8bFQzUf7hp0Y8/hF+tl73GRvIN
rHv5hv3Yc9G5HY9vsdHPa0XvMpb147rde3Qb49joqoS5vsQgB/sny7P0MRs1
FynuP3+Bg7uMbn2QllBE1rb9az6d7sP2JnbqC6cpois7fLSVjvXhn2ZrPbdI
K6Jxwe9QK7Jf/bK/TcNHThF1aYm5lxn34fJja7z72IrojxJauHFhH46/3LK4
dqEiMh9xsSn70YstCpffiDRSRA/rDo+tN+zFbxZVw5xbiqjx2FsNe/Ee7CSd
IqEhMxN5bpfdWWPbhQuTFui/CZ+JpF6tdCrR6MC1m/Qq9jKVUIBbccWbu22Y
7/Rgn12oEtIqtCh+H/kdV+5J3D9jpjKaaP3+Xj20Dttv0EtNjVJGcZ7WaWpx
Zdixv2b0R4wyKpG001SzLMOuD20NqPvKqLrn42klRhm+LHPtt2+8MmJN22Q/
fu01jmuqUt2booyCz79dpNhYir95W4WNFyqjhu6ZntHuxXhnsZu3eYsyqndO
mPErJRubu0hUBvxURhu/N7a0qWTjA/NiZpa0KaP51766N4c8xyeDX+Qu6lRG
x/ypQ1f/ZOLrh0S4/zjK6L2P5kCGbzouFd62T5lURh5PoxdFPn2IKzLnZreI
qCCe01lT9p0kXG2fPcUQV0EtYXnV6j0JuOXDl7iLUiros9H9+FG/WMyPVW3c
w1BBv7clFm9LD8dTZs8WBlAqyGHRzk2M3lA8TRzOv2SroNXD9tF3Lgdj9ml7
5iIVFVT39e/JC+lXsercEdvDqipoAdU2r+KyD57XcCMzTE0FvUjxsI9MP4eX
3VSaqJirgtbHnblclH4Ua296avJPQwUt1L79v/+3Sv8Ph2LZdw==
        "]]},
      Annotation[#, "Charting`Private`Tag$12553#1"]& ]}, {}}, {{}, {
     {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
      1.6], GeometricTransformationBox[InsetBox[
        FormBox[
         StyleBox[
          GraphicsBox[
           {EdgeForm[None], DiskBox[{0, 0}]}],
          StripOnInput->False,
          GraphicsBoxOptions->{DefaultBaseStyle->Directive[
            PointSize[0.012833333333333334`], 
            AbsoluteThickness[1.6], 
            RGBColor[1, 0, 0]]}],
         TraditionalForm], {0., 0.}, Automatic, Offset[6.272]], {{{0., 
       0.}}, {{0.6283185307179586, 0.5877852522924731}}, {{
       1.2566370614359172`, 0.9510565162951535}}, {{1.8849555921538759`, 
       0.9510565162951535}}, {{2.5132741228718345`, 0.5877852522924731}}, {{
       3.141592653589793, 0.}}, {{
       3.7699111843077517`, -0.5877852522924731}}, {{
       4.39822971502571, -0.9510565162951535}}, {{
       5.026548245743669, -0.9510565162951535}}, {{
       5.654866776461628, -0.5877852522924731}}, {{6.283185307179586, 
       0.}}}]}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}}, {{}, {}}}},
  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,
  ImageSize->Large,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "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.9999998592131705, 0.9999998782112116}}, {{
     0, 6.283185307179586}, {-0.9999998592131705, 0.9999998782112116}}],
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.764562649920072*^9, 3.764562670420122*^9}},
 CellLabel->"Out[22]=",ExpressionUUID->"7ad71d8b-f25f-4735-89df-6af91bbc5bd5"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"poly", "=", 
   RowBox[{"InterpolatingPolynomial", "[", 
    RowBox[{
     RowBox[{"N", "@", 
      RowBox[{"Transpose", "[", 
       RowBox[{"{", 
        RowBox[{"x", ",", "y"}], "}"}], "]"}]}], ",", "t"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.764562676462694*^9, 3.764562718475502*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"c750fc2a-6f2f-496a-b487-242d7e074886"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Sin", "[", "x", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", 
       RowBox[{"2", "Pi"}]}], "}"}]}], "]"}], ",", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
     RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}]}], "]"}], ",", 
   RowBox[{"Plot", "[", 
    RowBox[{"poly", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0", ",", 
       RowBox[{"2", "Pi"}]}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7645627332954073`*^9, 3.764562753830508*^9}},
 CellLabel->"In[27]:=",ExpressionUUID->"a2082dcf-a84f-4748-8969-810cc7ad8e0b"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwlmXk4lN/7x61ZEjPPjJAU7SsppYXuUwlFSokWIWkhSURIKpWQZK9kS4X4
RHZJTpSikKiEopAsYx7rmJksv/O9fn+d63Vd82zv9/u+73Ou0bA/u/e4mIiI
SJaoiMj/VnEsW9ux9upmXun/rwt7SjfNibSFu/mNw7PJKnyu5dAUcRbORecZ
fCHrH/YpfkHEZbjytMwnmazNhtfLn0TcAbWrSdX7yPrVNmJxfEQirNs+w2J/
eCLsdVZpvheRBWux8JX37SwIH7LOjYnAsF7xqYGaH4ZNFTtlt7rXwY2SV9tl
t9eBq4/Z0VzX71B44Y1OSVMjnCk6Hio2/xcsypjzyZjTCl+6d8ik2ndCZ1fT
wQrpDgjW4okUHOwGDal00Z3MLji3vDL69RMOCNzfr83gd0Pn8c4rSQUcOOdZ
NH+teA9YJok4X3nPgedXZqgVzeiBjYobtqIeDtTb6Sx6Oq8HxKaecvGKftja
M3OLrmkPRNTfMsK5/WB/MpktmtADOd67Ba/KuGDW1Oihq9cLwmOnM3+6DEDG
zO7au/Z9YMcc9irxHICSqcNzXzj1QUWpz7YHfgMQ4RBZ0+DWB2Eqwd8PhA5A
kYnnplH/Plj8KU3sy7MBWGcWvKg7uQ8sNnZZfuAMQPaxhbdCfvdBJuPoVIHT
IKjEWJwpP8gBdml3VbTbIKSUiVVNP8oB79OuUed9BsHkR7Oh2SkOGLy7vHR1
8CC4N4gnl3pyoOViwr5nTwfhTo7bE9NIDkj/bUl91D0IHw3k7WZ/5ID9K8s9
YSeGIGA/5vmt7YfVBcqUytkhcOPND3bX6wexrOaGhxeGwHveqmrbbf2QnGRj
lRs4BJoP7SzmmvdD+7UTNt/Sh0A3OXzzRmei005P59n0ELxe71vzM4lwY3Rg
2oVhYJze2F8rzoXVdVY7ta8Mw5HOxpx8GS6IVanIFQcOQ6BDISdSgQvJxfF3
Pt4fhp7D4QXrVLnQHv8khvtyGA7k3TqvuJoL9g75j3VERuDXG9dV62y4cGzw
y+vXgSPg6O/YJsjigp1PmFVR2Ai05sw+5pPHBWtxU27WvRHw1K3ZMVzEBUvF
N6qJaSMQGle0qZb4aLwh2/NS5Qi0Gc2PXNLABc2rt1dslB4FmYnPGi4jXFgm
u+ONNmMU1GKdLc7zubA4UuLQUuVReFu0tMltnAvqTy7eVF4yCrdB6esBCRpY
VY7tPKNR2Bzme7GSRQOfYXgv9+YorMuIy1yrTcPIfVGtjDuj0PYxN+y1Dg2D
815VJN8dBX9RYf629TT06ugMh6eOQlHd/XergYafB+aZub4fBdFg/e5cUxrK
kybFV0rxgMtRPr/lBA2lS4tjFyjwoMrjy4I9jjS8zPHQnq3Eg/6ygulWzjTk
veXYTF/MA1py0GuXGw2p3c0vegx5cG9+TsmPSzSErCo8mxLAg8NaZ02uR9Lw
YXVjl0coDwZ3JC+XiaFBeu3Yke0xPDDmiNwIuEfD9Q26pp1PyP1ZHUlO8TRc
3Fq4ZF4FD0avKJz5nkJDkUFj4mA1D/a2npim8ZQGnuHYzLIvPJghfYRtn0GD
m4mupF0nD5pviIbUZNHgZFH4O158DN7tjGjZUURDmmXjgTPTx2BlzMbplsU0
dB0Y+6THGoNSkVMZh0poOHpEt7Rl3hjUR5zctOs1DQdOFMaqbB0DB92RkMb3
NBh5Fe6LvjIGOnPFo658Je/r0/jBIXAMlp/aGMNqJPr5jm3RCRuDwkg1scTv
NGy+qruqIXEMnvu4Kia20LAuuFCOicfg9VJndYPfNJwPabz2690YPLmqciK2
nYac0DFhVu0YrIrOpXo6aFgZqdtj1joGnFLXxLNdNCyIK6wImRiDvHWVfhp9
NNgnNOpZS/LhWtPN4Y0cGpKSxnKXz+DDzIHuatN+GmY/0U3+MJsPOVtD/jtI
03Ao1UoldgEf1J3NcywGaLj39EKY4wo+bOx8pWo8SPKVWegnrceHptspysrD
NOx53jjauI0PA5p5GSOEQ3PGnFNN+DA8ceth1QgNsoW6hw0P86FWUJdnzaPB
+IVVveIxPkyf1fp51hgNAS8v7PjjRN5n9+iez4TFXhfqXvfhg/nX6POLBTSk
z6BSJvz5ECLxe88XwnsPO7Mv3OLDz7RGHS8hDcK0d/4DkXwIrzRRZ/2j4RFP
fdAxjg9aXCbrCWFTg4u2HY/5ELBEdIbmOMl/+Nca62d8yJcdnpFJ2GBlcPru
Uj5MpJkviJyggePTqVz1jg9mGn4b+ISjKjff3PqJD+MvVuy3mKRBb+b90ZeN
fLCK23YhlfCfY8PH1v7iwzrPjMRBwqHZu+ozu/kw5b27dvUU8XMqFS0Z5MMz
GzXx04RbTcWyHgr40GAruTmWcECstZqqmACe6Av9XhPW7C64FSUrgDG9wYqf
hL+tZQpnsASwbXE7NUTY79rpUzdVBZC+8/3xCcKLP1d8E1kggBHv+NIpwnVz
1Lf7rBBAtJyDmoCwl7NP7rCOADJaZ/v3EFYv/qJxRl8A43uqOXWEK6W0wrq2
C4Bn7H0ki7Dr/qBJWzMBZK1a+eU6YeVHHc5NlgJo4w3vMSf8ekC/Za8tYfkf
DWzCjpvv7ag+KYBHrjzrT+T7qZChwu2uAnhxz5JzhXBxk+ki7CWAzYGU/1LC
9otTo9ZfFcCPt1vnVBF9ZT1ExXOCBFB+QKzMlnBO+eFzyyMEkFDv4tRP/DnE
KGh7HCuAO9vTZrkRFrNhmM15JID5Oz5/7if+ZmQ4ldzNEICJQD7MjrCF4O0y
Zp4AdvEC938g+Rg3nHs/uEQArxwc5y0nbNre4HGpRgDnX1ENn/kkD1qanbyv
5Ptv6BUqEY67FLjXtVUAyDf0kQXJY7+yvtYxWgCXDwdG544Sf3endBszhbDy
ocsfySHib7yIVbmKEDLG7JkSpD5aew9VbJonhKM5F40EpJ60AhSSNdcIgdYx
bqgk9VZX4nWYvV8ItSaqbuPdxB/ZhqrbR4SgLk2P5fwl/hxYuV7qhBC09K8H
25F6dh3+zRZ6CkHc4lxjEql35jKT2rZ7QtgdZala8JP0J6lEZmWSEPrOuGmh
H0TvP0MWz9OEkBDpe7G8mQaHh7EtV4qEkBNeX5JF+k+Rcm+3xnchHD6qRa/4
TMNtnv4K2V9CYFgeqvH6RPz5En526K8QFsdEdJXU0DA9bAOvfEwI3KDWjys/
0GAjFSTuoPQPOOd0XXLKaVjz54eh6dx/4JlOab8m/VD6zapgncX/oDYhcvW7
UuKn33eGpO4/aBnxHC4m/VSSt3hOiuU/kBos9THKoeG/znfru6P/QZHwuGJ/
Ig1Xy2f51sX/g89+mSKXSX+3THLBRU/+QY3shjUyD2gQsZ5pGJT/D8pPXHKd
RubDvobj+5Z9/QcfB+SOPL5F+kGZhIszexySeOWhLh6k3yRufTQQMQ7vAuU2
bdlGw83Kup/+D8ahSNVjYRSZb+8GbZQVH5Pfb1+l3LqJhu0GF29vyB+HuVc7
ZlmR+bi1J/eCf+M4/LyenftzAan/NQtN2bMnICNoTHMzmber3kuN6j6ZgKsx
HTkWJVw4OxCt9eHZBKDnN3qCC7mQqbLAybpgAt5aXd9dlMOFlc6o7cr7CZCR
nZ/FTyPzneldVdUzAQW3t3cqRHNhweHe+MOak+CccdTx/mkuKNHVhpcLJ6Fc
st13PsWFe+H4aA+ehNCrxko2clxQ1snx3Vc5CS9Tl9qHTSPsfTdncdMknO6o
Uvgh7IdZ4vZz6oSTsFJLgq/Q0Q9zlMZG1TdPwWj78Nm+7H5YBBpPyt9OQYrX
GPPSjn6IfV8sKu4ogg6zfy10t+PAwK/n3rSzCLJ/HpemSfZ/RsKUwRZXEXQi
PfZghzkHRlZE/s7zEkEL1I7uWr+NA2YRp8tOBIogReekZ2ELOSB2ZPbVD6ki
qLg+oH28uw9OD10SiegSQb6NIiJ1ZH+qp7ZlSt1BFNXdDDxYc6gXjK2/uNmf
EkXf2950aZj3gsWDk12PnEXRIcEN0bNGvXBaJbRmkYcoavsX2zu0phdiFVse
rAgQRZ+Kjn9/KNcL/Bke69eniSLfUsOohJIeyJtMO7ebI4puRhgGLlbqgeW/
GJ2XzouhQ7rDnWqFf+FZGnUq300cyfGkPrm2dYJv65GG5QES6Ofe40uvufyG
c3rrTRKiJJFvnVHUGr8WECYa2PeFTUPNXaYPKzW/QJ1rrN2FZCmkJue9xO3L
e9hgNG/y50Np9Pmb3IlJkedQtX1DGoqTQRvqbOZkhWTiJ9n3N7RnyKLPLHW9
JxcqseL8FU42T6cji8U5Cwz/NOB6UWl/lRI5JPEh8alOcwvOl3ZPqauagaK8
V+W/ifmNVygynOwr5NGKsEcHrF92YjqUx17zTQGtus3sWFT1F3MG5pSvbFZA
holviqsa/uLevUZnl7QqoLf0vmcnWv/izpn3qtS6CI/f+R46/Bc3J2zwk+Yp
oI0RppN31bqxvVXpMFeWgdSX2W+95dqNd50UPfJ2GQNti3do7GD24Pk3b2qd
dWIg6SpL/4Ytvdjtm46fvTMDrSsQj325sxe/Xthevd+FgZa0njgXv68X27zV
c9JzY6DwllZLs+O9OFZ06LGMDwOd7vtc5hDYiylf61mPgxjIyyzM0L6mF4ud
05ZsSmMgx6PjdeF7+7D569Z91ekMtDJI8c+zQ304SSEkGf/HQL4xel5v7Pvw
5swuSHnOQD+yJdktbn3Yp++Bj3sRAwX9ZXFeRvThoePTBmZUMlDdp+ILWz/3
4fZDLc1b/zKQ9vYtsy13cLCcvlWqXQ8DJbxq+aRnzsFr5za4+/UxEFdeIUb1
IAcHdnyQK6YZaJalm3nZKQ7Wci7erM1nIKfCaUbPbnKw76X7yXNlmCjdLtJE
8i0HKyZZOY0vY6K90iEaqev68Wb/hnWzVjJRycl6+RD9fnzKYbf4ei0mMlwS
XH3KoB+XLDGKc1/DRHq9Sdkye/uxffa62t5NTPSmTru/27kfZ5Yrrm4yZSIN
9eSbqg/7sfGfBkG+CxN1zjfNNhDn4kvLQy9muTKRUUnU8EtpLs45ZzyR5sZE
l8LVO5fLc7HqZIlInCcTvW8r+DCqzMXcmSlS/n5MNBmmE6StycVRxl6KZqFM
VL9S+j8jKy7+lTFbu+sZE9UkKnD+POJixaFvOW1ZTLT6uWGgaxoX71wfrtOU
zUSp579mDP/HxblvJddX5zNR5p/o6x35XBzwk7s55xUTLWeY2Di+J0dZhTJT
vxomEuRWVcj1cLGX2/FTM/uZaF+yvW3GIhoHmnBPbqSZaEvVqmkRy2h8d8GF
kzaDTHSSI/b6rCaNC74FnkgZZaLgPqNTjHU0Ht74n8O6SSbqcB917jKgsYv4
sJ2lAoVCHE2tO+1ofOnHRTsfJoUWjc3M2+RA45B8CbsEFoUM/x6VDz5J44yT
M227lChE8bLfTXehcffHDUc81SnUb3JoxSMfGh+LunIwRptCDDmQ1IigsdsZ
mYPFayjUc1dRpTuKxv6GEQda11Jo/8RFjZS7NE7iP7JatJFCkbnf5KTiadxq
/X5/wVYKzQ9xXG+dSuP+tXv2NxtQKNF/b0vzUxqPyzdZTBpSaNPFj857/qOx
alnvPkMTCp39lnBAPZvGBxbK7/22j0Ka7gVoXzGNT03GmAv3Uyg4hG8QVELj
C41zzeccoNBSyXOrC0ppHBOkveeENYWWPw96N1JO44Z+CzOeA4UcZULKJT7S
2Kzgwc6ZnhTaqDp862wTud+w/ZSiF4W2yj3Yw2+m8bVVy/IUfSikN9qi7v2D
xkXpRWqKfhRyM01vt2+j8bzErwOsAAotfOL++1snjfV+xD1hBRK9xOv4c7to
bKXicIgVTCG7rePKR/8SvSOH3lChFHp5/4bbxx4a824q3GXGUIiFWf6IS2NG
xTcT5j1y/ab7jTtoGi8XSxBhxhI/ah6vMRmgsZ3vCidGAoUeq+yeXDtEY58X
w3MYSRTyMHR3njtM42hecYNCMmGnJb9ERmj8wXWHvkIKhRIcTjRmjNK48xlj
SD6NQhZx1w958Gg81duYIp9O9A21aV83RmOd4ycY8pkUmm0eMPWQT/RJXlkx
4zl5vod3lImAxo5tI94zcihkJHRb2U844eC1DrkCCp181eSk/I/oE7PznlwR
hYa3jMsnE25oYO6SK6YQO/9OwfxxGnMZTaJyJRRyKXA/Gk9Y2iypYHophdZ2
BikoTBA9b508Pf01hTa35L/2Iqxfqak+vZxC4ro/zzcTPiDJ+yL7lkKp4X+X
60yS/G19FST7jkIT10r/3CAccvn6ZtlKCkU82ZVM9ts4tcRkWOYDhf4L8T/K
nKJxuYBKk6mm0N5a0/kmhH+sa7aWqaXQnsn4v76Eee4PmTJ1FNqQ5ZiZQpiZ
feqddD2FzriRaUR4BVfrovQXCil9N97WTthw+ZiW9DcKrTm7izlC2O5UaafU
dwqFVWT+JucPfPHJjftSzRSqPemUR84fOKbd1EzqB8nD78tB5PyBn89li0u1
kvosGbDrI/zBuqVw2i/yPP/XG78S/nM/2XlaO8lD8cjMQsIijY4a0zopVMWJ
Gg0jPIut/U2yi1wfn/7NnrCOOT9YsptCt823FK8gvDsUg2Qvuf7WiSQu0cPp
Y8CIBIfoU7s0KI3wdWmzpxJcCi2OjTt/iHDCdkUbiQEKMV3q7SUIv/D/QUkM
Uehc5a99KUR/7riTrziPQqFifzfXE/8OHjBZKhij0Gl7kXgbwm9zl3/jCijU
WFGyq4P4f9+Jo9U8QSGR+pjUr0IaS7yr/vFpiuRnzeiebYRdNJ4FVYiykIaM
166nJE8G3890PJdkIZ/3V5ANyV/WGrOwFCkWslhoqfuM5HPWHU39OBkWsmvn
Xh0h+R3YTsfcnMFC01PVWM4k3w9yXXfaKBLOebVsapDG0xTMx/YpsdBBX+5b
BmFXJ+3HO1RYaOynV94sUk+GGkMTOmos9OxiYJQSqb+hUPfs6QtZyLPA1zaf
1OeR3n02ootZ6D3XMiOom8aV23Wmjy1hoYTiGDNLUs/x4yMO7SvI9cKzcd9J
/Rs7eSq/0GGhrBvyNldJf8itsKzIXMdC/Z23zyi30niOhq7b4/UsFOo18i2V
9JPhxrGPd/RYSH+ouSir6X/6e18+YcBChT0r3xxroLHMw4MrrA1ZqPujMOPV
ZxqfH9/QZG7MQldEUwYV6mi8I1e4Wt+UhUxq10bHV9N4RN23i2XBQrNseMZH
Kmhs62sdKWPJQnkeDgt83pB8NeqhKSui5x41uztlNE4Mnbjfe5iFdu02Fia/
ovHOcb9dZcdYaMmf7+9c8mn8sPFqnos7C8XMuqN9/xGNuypWtYV6sJAQFRRK
PyT1kNcmk3WBhWaAwe6zCWQehenb0hdZ6AYnIF/jPo2rdwhkXK+z0AabBnfV
UBrzX7ranosm9zfXdDzmSePN6XODw++ykFbNoWMG7qTf3qvNy75P/K1rvjjb
lcbyHitkh+JZyCvo1PQcRxrP1+zOc0thoUt9pqVu1qQ/JdnKni9kIba4qLEJ
onFUqPzaqBcs1GxWfzdBj8ZNvq9s816yUE3jccWe9TR2OKiaP4JZaCpzWsIx
bRp7U422HpUsFLBaXZzWoPHj62b5nk0s1MWxZduJ0lh4St/O+x8L1dcER+S+
IPP95TLkMcFCb7LTryIyzx/MUFY/N8VCzmckOyqec7FezlDbKXE2avD7KlKc
ysV+/1JtD0xnI84lzUqjaC4WD6VsdVXZ6FmWSYDQhYvlcv8eGd3IRmX/7V9Y
MZuLrSW/6g/qsdG26jnJ9kpc/MyqXK1/MxvVUPxCPpOLzcbjfnZuZaPYiFV5
8lJcHLZ975EvO9koNcGD3T7Qj9nfS6xzD7FRd+XCtMVv+rHaRPjhcz6E/8aq
rLcn+6P2XS/u+rJRe22ulsPhfnz+vYxSqR8bFQzUf7hp0Y8/hF+tl73GRvIN
rHv5hv3Yc9G5HY9vsdHPa0XvMpb147rde3Qb49joqoS5vsQgB/sny7P0MRs1
FynuP3+Bg7uMbn2QllBE1rb9az6d7sP2JnbqC6cpois7fLSVjvXhn2ZrPbdI
K6Jxwe9QK7Jf/bK/TcNHThF1aYm5lxn34fJja7z72IrojxJauHFhH46/3LK4
dqEiMh9xsSn70YstCpffiDRSRA/rDo+tN+zFbxZVw5xbiqjx2FsNe/Ee7CSd
IqEhMxN5bpfdWWPbhQuTFui/CZ+JpF6tdCrR6MC1m/Qq9jKVUIBbccWbu22Y
7/Rgn12oEtIqtCh+H/kdV+5J3D9jpjKaaP3+Xj20Dttv0EtNjVJGcZ7WaWpx
Zdixv2b0R4wyKpG001SzLMOuD20NqPvKqLrn42klRhm+LHPtt2+8MmJN22Q/
fu01jmuqUt2booyCz79dpNhYir95W4WNFyqjhu6ZntHuxXhnsZu3eYsyqndO
mPErJRubu0hUBvxURhu/N7a0qWTjA/NiZpa0KaP51766N4c8xyeDX+Qu6lRG
x/ypQ1f/ZOLrh0S4/zjK6L2P5kCGbzouFd62T5lURh5PoxdFPn2IKzLnZreI
qCCe01lT9p0kXG2fPcUQV0EtYXnV6j0JuOXDl7iLUiros9H9+FG/WMyPVW3c
w1BBv7clFm9LD8dTZs8WBlAqyGHRzk2M3lA8TRzOv2SroNXD9tF3Lgdj9ml7
5iIVFVT39e/JC+lXsercEdvDqipoAdU2r+KyD57XcCMzTE0FvUjxsI9MP4eX
3VSaqJirgtbHnblclH4Ua296avJPQwUt1L79v/+3Sv8Ph2LZdw==
        "]]},
      Annotation[#, "Charting`Private`Tag$14806#1"]& ]}, {}}, {{}, {
     {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
      1.6], GeometricTransformationBox[InsetBox[
        FormBox[
         StyleBox[
          GraphicsBox[
           {EdgeForm[None], DiskBox[{0, 0}]}],
          StripOnInput->False,
          GraphicsBoxOptions->{DefaultBaseStyle->Directive[
            PointSize[0.012833333333333334`], 
            AbsoluteThickness[1.6], 
            RGBColor[1, 0, 0]]}],
         TraditionalForm], {0., 0.}, Automatic, Offset[6.272]], {{{0., 
       0.}}, {{0.6283185307179586, 0.5877852522924731}}, {{
       1.2566370614359172`, 0.9510565162951535}}, {{1.8849555921538759`, 
       0.9510565162951535}}, {{2.5132741228718345`, 0.5877852522924731}}, {{
       3.141592653589793, 0.}}, {{
       3.7699111843077517`, -0.5877852522924731}}, {{
       4.39822971502571, -0.9510565162951535}}, {{
       5.026548245743669, -0.9510565162951535}}, {{
       5.654866776461628, -0.5877852522924731}}, {{6.283185307179586, 
       0.}}}]}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}}, {{}, {}}}, {{{}, {}, 
     TagBox[
      {GrayLevel[0], AbsoluteThickness[1.6], Opacity[1.], 
       LineBox[CompressedData["
1:eJwV2Xc8ld8fAPBrlJFy73PNlGgnJS0NfE4yEqkQDSstSVYhI5WWVJJSJHu2
bInkRClaX5JkREkyrvvYd1i/8/vL6/26PM+5n3XOeVF39rA8IspgMLJFGIz/
/xTD0l/+rLugf+BRXW/56gv6i7rLNqvecQSxGKtm15uOIMzROtwY6QGPrTpX
aN30gL9yLvznkecgofSijOLNc9BkfKkiLfIWOP7+cU144xbUO0YuiYtMgCSz
5yy7Gwlg6abcFB2ZDS8s7BauupYNtwft8u9FYtDYMpLGD8CwuXK7tMGpGggo
1PpuaFgDngEWB/M9f4C8h2cZNDbAyRdHwkUX/AKPZxX75Dmt8K3LVCrDuQPk
LGXT30n+gTCtUcbzfV1QNCNaypbVCV7Lq6Jep3HgemV5SzG/CzqOdJxPfM6B
tBSl8q1i3WCTyHA7/54D0e4+Oe9mdsMm+Y0GqJsDIpH3Uovmd4Po1CMu1uyD
SwObbxuZd0Pk1+smOL8PxObG8WfFd0Oe/07Bq3IuKB5ZZ2ao2wPCQyeyfrr3
g3rF0r4k515wYg2dKfXth/I7ZkveuPZCZVnA1tjgfkiNmot/evdChHLYj73h
/bDo/djUZEgvLPkvU/Tbs37YL1wdMZDcC9abOm0+cPohwPWIbtTvXshiHpx6
7joAKtTzjOp9HJAr66qO8h6A+X6tvdRBDvif8Lx7OmAAVkZqW+5x4YDhu3PL
VocNgPjSPcGVvhxoDoy3evZoALy8mhdY3eGA5L/mjJSuAbBtylgw/yMHnF/Z
7Io4OghPi03ML67rg9XPlShlj0EQD3F476fbB6LZTXVJfoMQb71q4NDWPkhO
dLDNDx0E+9ylpgt290H7xaMO3x8PwmhDd52eWx84b/d1m0MPgk7zwcxficQN
UaGZfkOw5vvAsRoxLqyusd2ufX4IzM0ilYqkuCBarSxTEjoEPslRelGyXEgu
ibv1MWYIJBVVWnVUuNAel3aP+3II8ibFzBRWc8H5cGHqWsYwzF3vfGK9AxcO
DXx7/Tp0GLYO2k3xs7ngFBBh+yJiGCRdy0/5F3DBTsycmx09DH5ixmaDL7hg
I/9GJSFzGDYdM5j3meRx28Zc37NVw6C+wEd3SR0XVl64qblJcgRqZR3t3Ya5
oCFt+kabOQL8gPYQbz4XltwR379MaQTO+T4Veo5zQS0t8KrS0hFgs9V+2ojT
wK4+3j5qMgIJsw0t37Fp4DONo/OvjsCcz7nL1mjTMBwjovXk1gi8OTHKKFtL
w8D8V5XJ90fg99VwlS0baOhZu3bodsYI5NuA4iqg4efe+Rae70dgbnAxM9ec
horESbEVEqPAkZtlpX+UhrJlJQ8Wyo5CmN7mMfPjNLzM89GeozgKIqr5jVZu
NBS85TjMWDIKL7Uvrd/uTUNGV1Nxt/Eo3O2cd6DxLA03VhV5pF8ZhUVP1ued
v0PDh9UNnT7h5PmBq29Mu0eD5DqevdG9UZhvEtgWEk3DpY065h1po+BRFkUf
jaMh0KBo6fzKUajVfJ32LZ2GF4YNCQOfRgEMuzfOfUTDqDFPofzbKPj9qDVw
eEKDt5nONKeOUbhRW5pTnU2Dq3XR7zgxHvjOdeAavaAh06Zh78kZPNiQ1KRs
WUJD517ef7psHvxkjRfalNJw0F6nrHk+D4biDfRNX9Ow92jRA2UDHqhdrL9Q
954GkzNFVlHneRClWOgdVE/WG9Dw4XAoD+zXZ3nLNpD4BfG2rI3gwZYez5YH
P2jQv6Czqi6BB08ON/c+aKZhfViRDAvzoHre4X74TcPpGw0Xf73jwcmy0rX3
2mnIC+cJs7/wwC5a9VfHHxpW3NHptmjlwZ8Vr11PdNKw8GFR5Y0JHry9qmUy
t5cG5/gGXbtpfLCd8n+5nkNDYiIvf/lMPhjMv3B3Wx8Nc9J0kj/M4UN8WPmp
PTQN+zNslR8s5EPrsbP+u/ppiH7kF3Fckw+b1Mbatw6Q+soqCpbU5cOdLXt+
sYdo2JXTMNKwlQ8nKipO9ROH5/HcMsz4oGkh51o5TIN0kc4B4wN8MM+MCd47
SsO2Ytuv8ofIepgnHyrwaLjy0s/0rysf3jw/NOczsejrIp1LAXw41znTYIGA
hsczqfSJED54HombV0NsecBNzu86H27NGRE/LaRBmPkupP8OH0a+dw7MGqMh
ZVRt4PhDPhT8Me5MJDY3DHT8k8qHFnnO72XjpP5v13+2e8aHq2+zOx4RG64I
e7yzjA/eoeJT4RM0cAI6lKrf8SH2l7XCMPHdKv2rBv/xQWezwdpdkzToKsSM
vGzgg9WlvL0pxH8PDR1a94sPFprhF/uIw3N3fM3qIu9b8rpg5RTJ51QGWjrA
B2b/Os5R4lZz0ewkAR+cQ9qW3SO+8sBuroqoAAqdsk+WEq/sen79rrQA7G8/
KGok/r6OJZzJFsDChGgJLnHwxRMuV1UEsDMmxV5AvKS28jtjoQAS2C+Kx4hr
VNWMAjTJ54vrVUaIz7gF5A+tFYCa+/jFv8RqJd/UT+oJwE5m9eAn4ioJrYhO
IwG8Sw86+oTYc8+1SUcLAWQc+/vrPLFSyh+3RhsBVE/4HdxB/Lpfr9nSUQCB
Dw26mMTH9aNNPx0TwMWxPT4fyfenbgwWGXkKYKltpeRZ4pJG88X4DHm/Z3ry
ImLnJRl3N1wQQPR11tZKEl9pHxGxvGsCWH5KvecAcV7FAa/lkQKoPz10v5vk
Zz/zeVvqAwEInDPM3YlFHZgWqink+911lOwh+X3yxLX0/hMB/PyDPtoRWwve
arAKBMDL8Y56R+pj3HheTFipAFa0qR5bQmzeXudz9rMA5EMeqH/hk3rQWtkx
Wi8AbPlBWo744dlQS89WAVjorhvfReqxT0lP6xBN1lsWzc8eIfndmd61jSWE
MSUNB9FBkt84hm2FshCSnxXcYpD+aO3ZX7l5vhBuXj39cYT0k9YV2eSVa4Tw
qfO1cyXpt5rSMwfk9ghBUdZKWdBF8iNdV33TXgiSA7pp2f9Ifvau2CBxVAhV
C3h69qSfPYd+ywl9hfC8ZN/9ONLvLA2zL23RQhji1wzn/STzSSKBVZUohHfn
RGX1Wki8/w5a52QKofj8ioO4iYbDSQ+az78QQjT3b9ZTMn9eKPV0qf8QQuzX
S8JltTTcHNXTlP4lhMdrlNt8/iP5+XbbY/CfEP66v50s/kzDjIiNoxU8IVzS
oQY0PtDgIHFN7LDiGPzL/PEsu4KGNX9bjM3njUHT0sN+r8g8lHyzKmztkjEo
2B937k0ZyWfwD+Y0nTEY2ipp8YLM02mjS1TTbcZAb9vzfsM8Gp52vNvQFTUG
tcf2velJoOFCxeygmrgxeLaqpDCIzHebRHf8Im0MfrcdpKfH0sCwUzC+VjgG
1vbcTeJkf7CqO2KlUT8G7+yRYfJ1Mg/Kxd3d5MbBemh2jpsPmTcJBin9keMw
uLWvFrbScLWq5mdI7Di8dbZ4G0n2t3cDDkryqePAWCmCWzbTYGQYeHNj4TjU
KifW7iH7o0F3vl9IwzgkXtIOa1lI+n/NInO5ORNgzgwx1CP77ar3EiM6aRMw
ULllq3UpFzz6o7Q+PJuA7gzHC2FFXMhSXuhq93wCWEZ13KI8LqxwQ23n30+A
/Htje34m2d9Z/tXV3RNwWbEkWzaKCwsP9MQdWDkJ57IGT8acIOc8+pPxuaJJ
UJOO5yyguBB9Gx/sxpMgG41uOshwQWltXpBV1SR4yFi1Rkwn9r+ft6RxEiS1
1kW3CPtgtpizao1wEtxTT0Yz//SBqiJvRE1/ClL7VuhzcvtgMainVbydgqrr
4auDTfvgwfsSEbHjDLRwSVHFaScO9P/K8afdGOhIyMZtq8j5z0SYPtDsyUDV
Bf1jHbs5MKx553fBGQbSGBMb2biVAxaRJ8qPhjJQ0E5jx8hFHBC1n3PhQwYD
VaY0pE129cKJwbOMyE4GOuiy9WctOZ/qzt0ypXZYBBlPLTT+b38PbLP75u3s
IoIc+Oy6Bbt7wDr2WGeKmwiqORfU52nSAyeUwz8v9hFBRqJF9cNreuCBfHOs
5hUR9F5auzxFpgf4M302bMgUQXP8mL6Jpd1QMJnptZMjgs72yFgtU+yG5b+Y
HWdPi6KRJBvneUX/4Fkm5VLoLYaGnjzx8mrrgKBW+7rlV8RRjsKBqhD33+Cl
u8Es/u40lJ9gc1A7uBmECYbOvRHTkSc7dLxy5Teo8Xzg5JcsgfKomunu397D
RpP5kz+TJJE3b5XpGCMHqo02ZqKHUsi6wG75oxtZOC03ZmP7E2lU8EVtepJf
FZZfoOnq8GgG+lh44fzWv3X4q4hkiHKpDNoc/TNkdVMzLpQ8lV5TPRMl59xx
eHPvN9aUZ7o6V85CL//U03YvOzAdPiq35rss6ppQO7K0+h/m9KtWrGiSRXNG
bsDHun+4x9LEY2mrLIpj1Kx1af2HOxSiq+d2yiKv3ANHIob+4ab4jcGSo7Jo
45R9VMzcLuxsWzbElWai6a1XODc8u/COYyL2bzWYqLrbMOkvqxsvuHpVy8OV
icSesN3rt/Rg7+9rg53dmGj+ktlXXm3vwa8XtX/a485Edv99tE+w6sEOb3Vd
db2ZKM0bb911pAc/EBlMlQpgoj2LF2QfDe3BVJDd7NRrTCT3K3714c89WNRL
e1pjJhPxDApz71j24t2vW60+PWaiW2L5b7P39+JE2RvJ+CkTXb+QZ1Pp3Iv1
szohPYeJtvRL9Ld49+KA3tiAUy+YSFL6YcWryF48eGR6/8wqJkpXZGwxrO3F
7fubmwz+MVFLWmG5rSkHy+jZZjh1M9GTzR7e+rs5eN28ulPBvUwUxc5fN3cf
B4f++SBTQjPRWNaM3goXDtZyK9HX5jPRtA0tdVlXOTjobEzyPCkW0jq9vHT6
Ww6WT7R1HddgofFsiZ2Z6/uwfkjd+tkrWOhqRrrOTb0+7HJ4p9gGLRbS8H3U
f9ywD5cuNXl4ag0LMQdWfJK27MPOueu/9GxmocKTwqFutz6cVSG/utGchVab
eR2ek9SHt/2tExS6sxA9a2CVkRgXn10eHpjtyUKP4z1OlkpycZ7XtolMb/L7
v5bs1ZzFxSqTpYyHviwUWiZcParExVyFdImQYBby7w5p1F7JxXe3nZG3CGeh
3peVH0xsufjXkznanc9YqF3+Y9LfFC6WH/ye15bNQqo6OjqemVy8fcPttY25
LMSrH7YaesrF+W+nbfhUyEIh3xnz/hRy8ZWfXP28VyzUMVPyp8t7cpWVLTcP
/sxCBZfXyMh0c/EZ7yMuCn0s1HkDVj9eTONQM+6xTTQLretVq43QoPH9hX7H
HAZYiNIVv+i+ksbPv4ceTR9hIRPjoyqy62k8tOnp4fWTLJRsoc/4a0hjd7Eh
JxtZCvUeaMpud6Lx2ZZApwAWhdS1c6Y2HqbxjUJxp3g2hWqLJCxDj9H4yTEF
x05FClWW+wil3Gnc9XGjva8ahS4Hn0JJATQ+dPf8vnvaFNIMiW1SjaSx90mp
fSVrKOQCzzr/3qVxiHHk3tZ1FLr0wepfyn0aJ/JTbBdvotDdWP034nE0brV7
v+e5AYXk7uvU7Mugcd+6XXuaDCkkLht24scjGo/ParSeNKZQhsupkR1PaaxS
3mNlbEaha6p6X+fm0njvolmW360oBAXxD3eV0Nhl8t5u4R4KLVseGnOllMZ+
DfN2q+6lkGRTzrn8Mhrfu6a966gdhRoCjaYNVtC4rs/aYvQwhZR4BR0iH2ls
8Tx2u4IvhY5Mivm5NZLnDTlPyZ8hfy/wWDfSROOLqzQK5AMo9Gv5+SmfFhq/
ePxirnwwhSKyZQsd22g8P6G+n32FQkdzemPrOmis2/IwjR1KIbvpbs/mdNLY
VvnwfnYYhWRXVr+3/0fifWfwDRVOoVWSLfJV3TQevSp7n3WPQjJ96rJ6XBoz
K7+bsaIptKBng6MxTePlovEM1gMK4eSOQpN+GjsFaboy4yk0Q3bcZ/UgjQOK
h1SZiRQqSeprVRmicdRoSZ1sMoUkss7tmCD+4GmqJ5tOoZ3N6ZsyR2jc8Yw5
OCuTQhc+1r7wGqXxVE9D+qzHFHpbV7V5DY/Ga48cZc7KopBYkd3OeD6JT/KK
ypk5FOqxTGgzEdD4eNuw/8w8Ui+hglPdxPH7Lv6ReU4h59U2mXJjJD73tkfL
vKDQ6L5Z2+KJ6+pYO2RKKHTmomnvvHEac5mNIjKlJP+X/92OIZa0SHw+o4xC
/ms+bJ4xQeJ5/diJGa8pRB+s7zpNrFe1Um1GBYVM27nR34n3Thv9Jv2WQsJL
k2arJkn9Gby6Jv2OQpQyVySE+Ma5S/rSVRTa6vaohJy3cUap2ZDUBwr9267q
O3OKxhUCKlPqE4Wmnddda0Lcsr7JTuoLhXbUjwyfIR49lcSSqqHQYrb+i2Ri
Vq7LO8mvFEpkzAh+S6zJ1QqU/Eahq4esjVuJjZfztCS/U8iKIcfqJ3ZyKeuQ
+EGhdQ+MWsn9AwemXY6RaKKQwlDnM3L/wPfazS0kWigU+pF3jtw/cM48OTGJ
Vgq1tPtb/yP+YNdcNP0XhabzPJfXEv+NSXab3k6h52/bxfOJGQ3H1ad3UMhA
ouHXDeLZctrfp3VSKO2MLXYgXrubHzati/RLrUfiUuKd4Rim9VBIq3rFpR4S
D9ePV4bFORQKkE9yTSW+JGnxSJxLod+uzVY2xPFG8g7i/RQ6FdQPDOLikBZK
fJBC5WMjK5NJ/LnjrkFioxQKFF2j8IXkb99es2UCHoX2jt313k/8Nn/5d66A
QrcKctXbSP5jXDlaTRMUqne6e6ZWSGPxd59a/pui0KRrrBoQu6s/u1YpwkaP
D7WppJF6Mvxx8k/ONDbaH1Egs5/UX/Yai4h0CTZKvhcl+ojU5+xbK/UeSrFR
vZaVeT+p334j+t7VmWxkN6XYcozUd2y+53YHeTZa6TR7UDhA4+myu3lWimxk
5bLohgyxp6t2qqkyG9WmpAYqkH4yVh+cWDuXjR6d4R5gk/4bDD+VO2MRGx3s
iVuaS/rTvsfKQWQJGzVOLPe53EXjKqO1M3hL2Wjs/mZlS9LPcePDh9s12ah5
SPL4N9L/21x9lYrXslGwRqfWWTIf8ittKrPWsxHlN4HkWmmsqq7jnbqBjRxW
Kj9NIfNkqIH38ZYuG51ZUXH7SeP/4+9/7qghG/ntnkp1rKOxVNI+TTtjsp6p
3RHFtTQ+Pb6xcfc2Noo5lPRjRg2NTfOFq/XM2Wi7Uf2FmE80HlYL6mRbs9HJ
7F0G+ypp7Bhkd0fKho3Cgzer+74h9dWgi6Zs2ej1lQin6+U0TgifiOk5wEbO
by1FE17RePt48I7yQ2w063tEm2shjZMaLhS4n2Ij/+Eir6gUGndWrmoL92Gj
VMemCfEk0g8FbVLZfmxkWB0WeyKe7EcReo50IBsV6ofOUo2h8SdTgZTnJfL8
LYu6FMNpzH/p6egVRd6fl6ri5Etj/cfzwm7fZyOJFjPVLafIvI3+UpAbQz43
9QBlTxrP8tGUHoxjo+l7y35lHafxgpVdBd7pbBSrySr2sCPzKdFR+nQRG5Wc
tYvahmh8N3zWurvFbGTfsqk/VpfGjUGvHAteknj2iPt3bqDx4X0qhcOYjd6v
uKV0UJvG/lSDo08VG2ktOPmMo07j1EsWhb6NbHR7g2CRgwiNhS56Tv5jbFSh
feZ2XjHZ319qIJ8JNrr8NOsWkP08dqaSmtcUWW+YAaMyh4t18wbbXMTk0FI5
jm5xBhcHj2U47p0hh4wv3d5qEsXFYuGUo46KHEIDt64I3blYJv+f/cgmOVSp
sc/k3RwutptWrzegK4f+7UMNhxS5+Jltxdw+feIia4aQxcUW4w9/dhjIoZUX
9y1gSnBxhJGl/bftcigvIPTrn/4+LPej1C5/vxzawpeT0njTh+dO3D7gFUDW
I6v/dbMzOR+17yi+HySH3IsOMlwO9OHT76UUy4LJ80qS8q9b9+EPty98lb4o
h2bkH84tNu7Dvou9TFOvyyFKh2GQo9GHa3bu0ml4KId4MfXnpAY4OCR5FlsP
y6F0XXmtQD8O7jS5/kFSXB69zGvd0HiiFzubOaktmi6PhC4OFuqHevFPi3W+
WyTl0cJXiW8OkvPqtz1t6gEy8shopkbDx229uOLQGv9eOXnEWk53Gi7qxXHn
mpd8WSSPDsYurf/U0oOti5ZfvmMijxofTis2Mu7BbxZ/AtXr8iisSFjmKdaN
XSXTxdWlFJBqtVtPk2MnLkpcqPfmtgIKiEjfUqb+B3/ZrFtpyVJEbS86ve/d
b8N811grp3BFJO36Rd3qzg9ctSthz0wFJZQWJ/C0v1mDnTfqZmTcVUJi9qfe
BT0ox8f7Po+03FNCwybBqsHW5dgzydGQilFCgSnHTANmleNzUhd/B8UpoeI9
fPdDIa/xw8ZqFct0JXL/2Vw0WleGv/vbRowXKaEPa0xzHnmW4O0l3v67m5WQ
95b98bdSc/Fud/GqKz+V0P5QmbpYxVy8d/49hdI2JeTYcnFZSlgOPhZWnL+4
QwmZH5nZaPI7C1/az+COcZTQgRmnDf6eeYzLhDed0yfJepPVZjakJ+HKrHm5
zQxlVBEbrHX2RiL+5Jw7xRRTRic2Turrd8Tj5g/fHgZKKCML2cXR2QEPMP+B
SsMupjLi/exeZJ5xG09ZPFt0hVJG9t1evrUd4Xi6GJx+KaeMar00+eKBYVju
hDNrsbIyWrA4/I5TxgWsMm/Y8YCKMrqcEdOnEBiA59ddzoqYq4zotqehbhle
WOOq4kTlPGVk2Dzt8rWMg1h78yOzMXVldLW36v//3yr7H2U5yCU=
        "]]},
      Annotation[#, "Charting`Private`Tag$14890#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,
  ImageSize->Large,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "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.9999998592131705, 0.9999998782112116}}, {{
     0, 6.283185307179586}, {-0.9999998592131705, 0.9999998782112116}}],
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.764562755582024*^9},
 CellLabel->"Out[27]=",ExpressionUUID->"a185c92e-4c1f-4a66-95ba-d51bd36ac4c9"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"err", "=", 
  RowBox[{"Expand", "@", 
   RowBox[{"N", "[", 
    RowBox[{"poly", "-", 
     RowBox[{"Normal", "@", 
      RowBox[{"Series", "[", 
       RowBox[{
        RowBox[{"Sin", "[", "t", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"t", ",", "Pi", ",", "9"}], "}"}]}], "]"}]}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.764562856391041*^9, 3.7645628580115423`*^9}, {
  3.764562894495171*^9, 3.7645629138619127`*^9}},
 CellLabel->"In[36]:=",ExpressionUUID->"a3adec70-cba8-4948-99ec-b09fadfceb22"],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", "0.006925270707504705`"}], "+", 
  RowBox[{"0.023266839083691426`", " ", "t"}], "-", 
  RowBox[{"0.03432493464974673`", " ", 
   SuperscriptBox["t", "2"]}], "+", 
  RowBox[{"0.02903417761216412`", " ", 
   SuperscriptBox["t", "3"]}], "-", 
  RowBox[{"0.015466335279239435`", " ", 
   SuperscriptBox["t", "4"]}], "+", 
  RowBox[{"0.005371368307109145`", " ", 
   SuperscriptBox["t", "5"]}], "-", 
  RowBox[{"0.0012155682188676733`", " ", 
   SuperscriptBox["t", "6"]}], "+", 
  RowBox[{"0.00017291247487945136`", " ", 
   SuperscriptBox["t", "7"]}], "-", 
  RowBox[{"0.000014041898606899012`", " ", 
   SuperscriptBox["t", "8"]}], "+", 
  RowBox[{"4.966305719678991`*^-7", " ", 
   SuperscriptBox["t", "9"]}], "-", 
  RowBox[{"2.318610723662395`*^-18", " ", 
   SuperscriptBox["t", "10"]}]}]], "Output",
 CellChangeTimes->{
  3.764562859159561*^9, {3.764562896801812*^9, 3.764562914821203*^9}},
 CellLabel->"Out[36]=",ExpressionUUID->"90025396-220b-4f1d-95ec-175998ea7f79"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{"err", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0", ",", 
       RowBox[{"2", "Pi"}]}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], ",", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Thread", "[", 
      RowBox[{"List", "[", 
       RowBox[{"x", ",", "0"}], "]"}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
     RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}]}], "]"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.764562772859906*^9, 3.7645628806667233`*^9}, {
  3.764562927741186*^9, 3.764562979635928*^9}},
 CellLabel->"In[38]:=",ExpressionUUID->"2f94ef69-4336-4f58-9a12-5021d868d4fd"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwl23c8ld8fAPCLslKRrJStJCEhSudIGkaiUsoISZIyKjISKqSShKTse+17
s2mdRyIzKTLSwBchyVZUfvec31+93q/PPOd5LtcfyTq6HzjJSaPRnnLRaPhf
Loq/6T+tYJDJ1HfOtQ2gFIfQNql7x2HcMZffdaoB1Fy+mlNntDt01ShlqNAC
qP6VLr9Ko6/Am0qR08/T/amPu69VMqLvwOyxzTq0UT/qw/HodYnRybBe8b+E
2khf6oCbxMf46McwaPN6iY0L3pSADqTi9uTDMMvVPuod3tTYE0+hnup8WB1q
OX62wJsqfdFWcrGyAPK5vROMc/Kmdlan/Et6VgRdN48PRr+5SNm2akaNM8tg
2+2zpxeyL1B3J2yK4qIpyHvQZ6g0wou6LZKgV7G0Avod1fmITnlR7wuuS2sd
qIADeV5NtF1elPWwTb9UVwXsKrnlvYnmRZ2zXuI5MfIS6i2EvNXz9aRitruE
PxCsghxb57RKLnpQXznkyr4dqYWid4c+SKicpRSSlibYPKqF7zYG212adKNc
tv4OeNddC1+4GjK3PHOjJryaDZ6droPSP9ecbDNxoxb1BTZF+tdD/q3zH695
nqGUX3f1ayU3woBLTfMqTaepbdXG/Abnm+Ev/6jzi+ucKcHWjLdb7jXDfLHz
1y8lOVN9vZyxG4uaoeibwoRNF5ypSNozaYnJZrjhSNTOCzLOVO+2DVpjXu/g
9pldbeoBJ6nwwiX2SV7voXdeRqmGvhPVltxYMufZCpdISsGmWQcql6XkP3a3
FVbrpvgZf3Cggl5c0x8oaIXc97k3cxQ5UOu79BrejbfC80eEX692d6ACRJnd
WZ4foFr6Bz2uYXtK/vbtJUc822AQn9TkjdbjlIefmUORRwdc63xZuFPBhjI+
p5fFG9UB1U91HjIdtabkHZV/2j7ugGE6mpxryq2pNiPuy7yjHbDO4LB/q6k1
tV0Cxdue6YRKDnDxbb9jFH/5xmYe54/Q0O3ZjccDVhRjWgDaWH+C67ys9DT/
WVKBQ3PXC/w+QU+t7tHD7y0pq8+DjdwJn6CbnKxYTYYlJVBddayg4xM0R10T
t80tKa+YAG/uw5/hVychqezMQ5T+5hFmvvkX6PSv3LrgxEHqbPnJSE75bnjo
25rbPIIWlKqRhQyHZje0cens4x82p8Y69QoXDLuhoEa8xdEqc+rCvHDbH+du
+DiV3+2DrznlByrXzOZ0w/P2pkW+zfup0Cop5neNHjhLSxF4Qt9H7bXkh8M7
e+ATlQuFTId9FN/AdPPgoR44MWKjsV56H3WL581Uv3cPbA2qlix4aErdM/bX
637aA0fNmmy/J5hQKc3tDa0GvXAy95OcYYER5ejwyrblYC+U0k16s83biJKf
YP1859QLWZ2mlnnbjKgM4VDht6G98H3muqypmr0U87CmdV19L/SJWnR0ycAe
6umnqOEXB/6DxcYjFY+27KZaB434Mh37oPfBPXnflXdSPx8o7brq1Qd54j87
BNJ2UvwmPMHHQ/qgIL/DhpNtBpQ+s+qXWFofDHbe+FzlqgGV6wkHw3v64F5w
eH9M9w4qZE6zxtW+H67j2pZmkqNPJeYKc+3x6IeV6ODk52v6VLnNBJAP6ode
E7Rnz47rUz/Q49Ku5H4Ykngp9YioPmV1VTnD9Gs/XBHUsmXxMkipCchcV7Ub
gGLg4Cm373qU8Yt/FfznBqByyJrRsGw9yunc5z8DlwdgEvWcoeyiRyU0J1xI
ThyALvXLhjm+baO4Y0WcBD8PQKUXSwc+D2+lPq1ZsnPC+htUqaKLBPPrUhFq
M7TSo4PwGK/1c6d0LcrPx2jDgOsgFJTUG7vopkWdph5ZigYMQhujhozF2lrU
3v07c7yTBuHEpou5j+o1KR73qINb/huES1WLP97+vZm6xlLOKHcbgi8GDMyz
zmpQF2YuNw9eHoLaljdaVYEGdQK8mxO/MwSXbnvyXmG5BmXQ5LPft2AIih2Z
z7cu3ERxjFb90p0ZgpHF3y9++qtOBW48bvLsyjBUFjGe+2CuRp27WHjxe9Qw
NBWvuc+1oErZvlicIpk2DDmNCqSmmaqU3r68Kf+qYXggbXHEmmWq1Lzb70Q9
vu+wMVOu80i7CnUp797Yi+jv0LWq3UjyljLluaE2toIxArsP6oqzLqyl+k72
BaWUjkC60cLELd211OEUmltQzQh0ipe3+vNPkdoqomugPzQCadp+P8JuKlKc
C9mjlMoPOMXf8VM9V4GKfn9zD1X0A4rP9ce85pKnFi3N0kiu+gFj1s4nyL6T
o7z3VK258uEH5M2ATXNJctSx538mwewPGM/ScpPXk6PkMs6moK2jkPbcr1It
QJYq9N3/+8XLUeixdPOSfStkKIXiM32J70dhswXt36omaSpuNOzt5f9GoQvP
Cn+fS9KU34kKxvbFP2FKmnXos7dSlMG+TQde7P0J1eWzv2wNX0O1yAjnPn/7
E+YHBfzzlZak5k6cYX0+NwYr4s4zX38UpeyFJi899x6DY/I33nqkiFLVyG/n
w8Ax2FyHtuQ4i1JREhEdVpFjUF9ff7x6SoRa9zaLs5U5BmlF/QpXxESoQ1sH
DtePjEGPYaHMsYvC1NNvZ2Wzp8ZgyrZjzjRDYUo6dvp72B92vgKj6tYKYWr4
56LgXQLj0MO79Z9owQoqmCGfW7FhHFZU/oAtk0IUS9BhodR1HI71Hd1GuyNI
rUSDdbFe47Dbbn7gmJMg5XvGI+aC3zjMV419KLtVkDJ8fWW9RgS73q3H/8ub
5VSXf9JBZjY7/6XW7p6IZRTvt67M9MFxqG+xOcZHWYByfHHYPMp5AtLuCB9L
nOShNErFV0i4T7B/X19cmpnPQ3E+/tiS6sN2Qn7Y2XM8VFqK3ZGicHY+32/x
9SPcVO9VZ7u2nAmo/3BfrOLYYsrR2Ntt9U92PFvrabz4Ikpjp44qY3YCVuwz
HL3TzUVx6s393EibhBUdQ9/WZnNRaRsDz0OhSaj/Isl7ox4X1SsY6ue4eRLS
buh1nzvNSTm2x4Zn+bA9qLpg+5VGaTQfMd4UxHaZIUdUPo3irJMQeBrOdoB+
8t4QGpX2NPFOw4NJGLS1vnvPOhrVm8iIG33G9gP79czRf8jRqYSuSZuCtC9G
nTTbP0jD1sf5BS/b9luahpb9QZyHdZV2C7Ft4n+9+OU8StvzIveIHNvmxtv6
18+jXuXqQj9DttevFHNfNIdOjLdWVISzzRFyfkX3LLL3izpSHsV2hIDTrvRZ
ZMNlOvo4nm3hxTmrT82iwyKvJJOz2M7p21YwNoP26hZ4X65l2+1IiMGSGaQa
fFtlK+80pI3s/x1tM4WU+Y1ebRJkm7lJ1WD9FFp3b9Gx9eJsf3n8wH9mEskw
/MPEldheOS0UFj2JhOtO987sYbvXXb+9eQL9EtwdXxTG9nTK1JqT42jqAYda
7h2cP5z2Snccjcu9qE67z/aRRpHPy8bRsKbm5N1Mtud05J2fjqHPVnJmHjVs
U9HXmaJjqDLlH9dGnhlI27tzjOP1D3RLvcw9I5Ttxga+rOtDqF6jfeBiJNvq
/znd2zuEeLVmbXfFsa1QJCkhMISu6W4x7WOw7VL8aCRmEPkblCnJVbPdB8eu
Mr8h10NlPYlcs+z4SAXn736Udbjd6uwStvOStzdU9aMBq9m3esJsq++TWnm3
HznYbkFdcmxrInnZDf3IyrksQcKA7U+j9QtOfWjPpbKDsUHYtp/jR3vRNb/2
eqdwtvVsM9HLXlQZMLtDM4rtijN/zWJ7EQjeot6SzPYFGeuc7b1IO6JMQIhi
W1CEQ/leD1J4VFZ96y/bdNWTRubdiLOibMs1v1/s+u1xqmKfUM7SFRl/Q9gW
PH9rTVMXOmDtttLnJtvi84Xl17tQ+ozM+OlHbDdnjgzPfESGGyNy9iO2Xb6r
LurpRKEJNmskOX+z37cxGeP37Uh1sPRmDD/bgjp5PPfaUZuW0NxSYbbHNPgO
HGpH695Vt9EU2A7fYxza0YZqedSiBnaxvTeycvu3D4j/IgdX4Q22Peq4rq5u
RYWV1p4botnWHzU27m9BxwRLv9IT2Kax7K6zWlBuruvz+7k4f5e7hUELMu1t
uXj5DdtR6aOS7u9R5P6Mwb1Cc+x4Ft3+UzMSUjZp+hrPdjenafJ8AxrgSRaq
TWE7qNXPMLcBPe2fOJSfheOiwd7WDcgpNaErqJxte/0y/xf1qFx8eFC2g+0U
V7eBsDpkx3ODy0lsnr2f04iyVg3K63utMxjLtj1Ha+dQJQquXBXQnIid5/cl
qxIdTjlHlTPYDmrYZOZSiWg2ortvlGCnPT4z+BIdbDl5UPkD27QdlTd+VKC5
l4vOua38w3aSoIsIQnuTDdLHotmu2LqmW7wMhdU2fw55yLa+CdMwoxS9HrcT
F6GzHcTD9NEsRbsM/W/rluD6jymGB0uQwVCRT0g7dpWibUIR0tusaLpy9V92
/tNf+fvykb9NUWiGAts0WhP/28fo6fUdL3U2Ym9/pLHoMdLtsNW2BdhqA7+8
mEj7cpxMhj126U9txxykXsMzvYWBvcZmbygduY/FqtUzsSkxnXfpiCWh4GpT
iq2YendNOtropv81qAbbtuLy01SkLORbVzeE7VKly5+EFKyHE61V/7GdsXPr
k3tI7Gfj7itl2DLRiu0WKP4u5TBEYdNuuQcuQ+KahQEHa7H1LZ7d3g/Efe8X
rusk+QYfHFzAKi5HqeY5bHM/2pYAICU2Oy0DFvB91P/SuA2SngwJ3tyNfW2u
eCESSNt82jBthr21vepEFJBJfelQdxx7pthAIRrIbbj1xjMYu/YFR3gsWAtl
GZVV2CuvhXQ9BBk9whUqb7BXOyo9fATWXuPuivuATVNws0oE6+q+C54ZIPM6
3zUkgZj11ohvEU2fRhvkkWxIAa1+rQ2AF1v32kXZVLCycV/neQFs0Tvp3qkg
9pz+1OeV2P8KAtakgbgiReVCReyLEpq26aBtUfKWb+uxuce25aQD0cPiu1ar
YrN0Vs+kg/u/+O1DtbG9c10i6CB++1jssT3YKWAngwESap5ycJ3GRsWx9Eww
1p3v+9MNOyn/elcm2DOXMd7lgT0lHCiUBaZU7vUUX8IO94/wywJm0WdeOodj
OwrI78gG9FwH3YO3sEd8Tnlkg/mqIwUwCrt2z5OkbJA1uzNVPB77Wn7Yr2zA
abs6uD4T+8Ce3pQcYOW94ndpLnafK702Bzy+w+uZ/hh71NL3Zw6wqZyyDyjD
zj9otzUXlK57o69agz3Avb46FwjseFW+qgH70PnY/lxw4tgTdZ632Le8Vy3O
A4K3GTLdbdj8xp478sCZicu06AFshepTj/NA5ZILlwKHse38WbV5QFzRdcx1
FLuZZ2VPHnh95HD3zhnstj6wnAlWe5paqc9h531crMgE5yMMmlf/I/3GR3WZ
QOaFasX0Yg7c74uyAxN4tyno9PJhs1yunGeCNz9X5Tctxe7z+32NCfzkeFIy
RbDTChwYTNC87a9YjAS2XrxFMROss5y8E7QGOzrCtZIJWsO+XjmqiD21XvIT
EyinfpjdtR77QGbxNyYIftrgrrERe7TZb4IJOlpefpPahO1679wfJlD9UXZ8
iRb26s9Ri1ngGjerfVYHWy6pdykLdEnT9/fpYYfUOomwgIZuQk2zPnb+fvHV
LBB+IAq+MMQ2WP9PlgW+nAkty96L7Wwqto4FtK4HqMWZYnMWOm1ggVtJXpkh
5thNh76pskBvmYu0+yEyb1PKJhbQeWd339oKW2dX1GYWuDN8aPleG7J/VLkm
CwxwmYRp2mM/FVmlxQJ6a3YsyDhhy7SXsON7bVq9HF2wBZ65susPPTw1kO6G
XR6gy+5v/3HuaL8HdutZefb8MxKRb9ZexNbuXa/MAt5WsjtcfLGV1M0UWSDk
fnFx9mVsw/y70iwQ2bZH6XsweR4vpsVZIEGk66FKKNm32F+IBTIOnVt+LgI7
dEqBjwUK73FcfRxJzvt1YoEJ0PuYmbFo7ObOgWkmqBNSctW4jx23l/M7+3ma
P/t8/iG2eczur0zQfcfMoiSZPM91Je+Z4NfSizo6WaRfu2gZEyzax5vnm4ed
1SqQxQSCtx5KP8vH7m7XiGe/L/yV3OAJtkLNSm8m0DCy9LvyAtstveUEE4Dw
wR8VL7GPOdeaM8Fh7uUfdtZj18cdVWICx11pe683YUct+r2CCc5d03r++j12
wNHWP3ngOqdNulEXmd9k2JQHiv9lee4fwQ4a2uKSByq26/VHjZF9rveY5oHG
gLdW76ewI0+2qOeB/+amoeVf8j7CuJlcsGLWYJn1Mk62945W+eYCKe224Ecr
sEuVP1vlAuWLp6c/i2IXehtvyQU7Ju98spfG5qxUnsgBHj8/5TirY6+QVLfP
AQGqHlJZmtifXkhtzQHhZ7mih3Swh5PdhXNA8vf1vm47sK1aQVU2aPrmvcfr
APYtpUqpbLChW7Dv8gXs1ZWd5ZnA7JtMpfclbNOqr7cygceoeop7ALaTwAX7
TFD8x9zG4Rp2nsQ+7kygJ3G31TAGu29u0DQDmB4Qql5SjD22j7eGDtyqhDIe
TGLTF9F8UwEza4VLiRcX2+LnLxrHgGsfqB1cPtjNWhte3wM2nGclLfyxU1Y9
3nEP8NvUNI1cxY7KVNOJBs7LA7QUYrHH+n3k2b/PvPs575VhyzADRsJAzM7y
RI8/2BW3FdedBWc8Tvog2iK2aaEb+U4Dg8QVFgKLiTeK/nAC4zNui7KWYgfZ
91w8CsxyZM98lSJxHbPzeohX6KaO2Q7silUlbz1RwBfblg2h2M0ztwzD0dSe
7YknI0i8rFPlBnIrWH0qORJbvWZeOAJZX+uaW3EfO6o698tNtFX5qNxcJraM
xKaTd9Dvi4c8a+uwx6rMJe4hj6+bt3E2kf5nFkLuocG9wov13mOnTE98v4c6
JN/HP/6Ibf4n+VkMKn+5vyLuO7b+kjaLOOS91GT5yaWL8f7uA0ceoFFv5c4k
Iez8wKHiB8i5my+9QwRbf62QUAI6XFSrbSqFPTb9+nUC0jq6226zGrb6aMK6
R2iKvoPJYYHdXeP6Igm5LZP12WZJ+p8YWpKM+nxoOy4exbbX9DqajFqMqdZB
B+LyR+PJqHBM789bL+woVbGbKchTT8ckKQZ7b2PU/VREOxdz5U48iZu/KE5F
kcljRUGPsF3kvd+lohzO7NUn6NgVnDf50lBvrcToumJs8UPaF9KQx9xFOYly
7KBvT26noX8b3h/mf45txXidkYYk70RQI6+wL233a09DhyznowpasLMK9Den
o97QI9Vp7aSf9SKjdORRXvT7XhfZJ8vMLh3dknRz9P4Pu1zSOywdVfd82rxt
ksw7zteajg4J67qozGKHi8z3s/sZxj5aM0/uO8xzNh39y9y3mMbFjft3fRen
I+1zqK1qBXa8s9BhOqpKXrWkVBQ7S1jBiY4OvPOGmauwZTY/9aSjc5pqWTfk
SPygx006+nPy5ie/tdhj0lH36Sji/jdBN2WSP62RTkeZc8m+ZhpknhvHE/Y8
lT9MqI1dPqLyij3P1qpXfSu2kmZXI3venWJROUDqdQTa6Ki7QtBE2AA7fIz6
wp4/4XZl0W4St5oZYM+XryuaNsK+dObxKB2JhwWv7rAg/eKy/9BRRvln8zpL
7KCiYU4G0hzWvf70KDbtBouXgSol457k2mKnSE8vZSDzfRM/Hjlg619FKxjo
S6CZXORJbKuX/GIM5Jafc/jKaez84c5VDDTXw33T4yw5j5CKFAPdED5BOXhi
D5pwyDKQ+C5q8sBF7Ob8o/IMlOEtqWToS+7jmJYie58sHxuty2R/p5i17H06
W6LWBpPzjnqvY++zRL1a7DrZX+Iz21/0bv3mvYFdy1vDttu5wY1zt7DNezXZ
nks2dPweRfZvWsvuB2a3R2+KJfsreSswULDZlkqfB8T6mnIMVMVQn3iRiM17
0EGagXj/rpdblEbuN+WvJAOZHJI/YJyB/cuVT5yBInNXh0TlkPP+vibMQO85
RQvbWNgdsR7LGEjk2PLe1UVkn1uN7Pu1KuBdcaKM3JfKA/b9P+LlNMh+hu1R
3zZPR1+Pz3v+pLDtn12foiO5sqlUrSqy39GMETpyXjb6zr+W5NN0++go++Q3
jspG0k8QdNGR+sqPDmYfyP5zKTV0dP5My92YTnK/aQvP6ai0svHlx8/Y3ebv
C+houycle2oAuyL4wwM6Cqott2AOY+t8577Nfh+lC4MnR4mVWVfoyLiJ3nNl
lrxP4vYn6ChSMUno9Tx21Oszh+joXcD9HQI0Hryv0U9DOrLaEJEaz4vNO3Jc
gY5Ohrs7FEhge6yd/ZSOsr663J1dQ+L9zXXpaETb8eV2OWx9rnWl6cir/5Bs
vTL2pY/Gt9NRkMHWnt5t2GMzIVrp6NWDzULr9bFpqRZS6Yh7XGWHuyHp/z2V
Ox3dTpZOnTcl+2jXf0hDCX8WOay0w84XMDuXhr4c/Bd1zJHkLwMH05Bs7mxF
ijNxcsmWNJR1dFhmozt2RRxYSEUlpW+7dwWR+rCm8FTU7JFg75OGbd+XHpqC
lKME5y9kkPzOcucUdPVxaKxXDnYz5bE7BWmPetWfLcROcWtclIISz5hoOlVi
d6+MnktCrs5/eCx6Sf1dkeeP0CIbu8fKcry4P5du3X1k599qpLQWW7/G/OJ9
VJZg3KeojC2j1i9zH7l2akvIahD73/SOQ+8OL78qZoAddI5rdSxKsqg4vMgR
O8XyoEU00t0j9+9zKjZt6stAONJMrw+vYZD8m/T14Uh1wUu4IBu7QnDoTBj7
fXyldK2AzH/zb/Q6WrLu5EHll6R+pvFnCPrMnZXl3f3/eRPIFwW/3nhwuTQf
jqe0emigALm2z7/kiGkPaBLAOzDQpXct8djCGACu2m8vl6gSawuuPAQOZnhk
WQPijvd2zkAxtOhvpi3x4jlvH1C3SzdL/xG2vg9vYxi4VBOoeiWFmOPug3Cg
ZFRV/IJO/DXQ+QYINd3/ciuT2N77XwQwPHjyoyYiTtyzNhJU2N8VWN9NLHdL
KRp49LZFn+ojdtp/IxpIO62WyBjErihNHIoGgacyFeXHiX++y7oH9M69AGs4
+XH+gqVsLCjzH/IQUsBOeUj9ug+c/6rO7lcicQEJi3ggeuXC5UgVYsuOrHhw
IWThxhIt4ltXDj8AGjdE0hfvxu7+ppGTAFhxOz78PoVNe2ayKRHYiYdZ67iR
/MR/gYlgaUJjj7cHdpDJdEMicEs8Mjp5icQXZ5xIAsr0szyjN0g99flGMmAU
PNDtzSUW+9ueAlznrZU7HmPb77aaSQHqu6Qkm4pI/dlDK1PB8460P0+ekfw3
KvtTQctCLnW3AbvC07giFcQbncsPe4stM/ygKxXY3VNPvdxC3HNtJhUMry0J
Od1F7qNAeUMa4DBDu3Z8J/ME6HfSwOv7QdpbfhInv8hMA7d6DNZtnCRmnaXS
gPjFGl6JebJPT+SPNKD2qLlhjH8Jnu9fbZgOpvqjnw8sw6ahkKPp4KmaJfPT
ChKfenc2Hex+1RlZuwpbX1YhNh3YDvdapCgTS974kg7kNBkGcaokP7bmZzoY
vHxq8y0NbJmFezQ6uCA0svLSVhJ/1SpDB1utWYvPAVK/Za86HdAYHjMnDIjr
9kA6uKkz3b7fGDvl7pA1HZiHlNXuMsO2H4s4TQeijb5Pth3A7q4o9aaDtOP/
EtYdww6SzrtDBy7ZFTfX2JF+gxce0oHqZEiAsCNx/OsMOngSymO3cJqcdyvX
czoIbK4zmz5L4r4/qunAcNUt+N2T9N9y/C0d8DuZqfdcJPvLWXXQQTNTULbd
l9QPtnfTQdzse6E3l8l5db4M0oHNjljOV8Gkvtp1jA6+tUr0sW6Q/TxF/tFB
ntSnVvptYtrmRQzg5ZJUnXCX+PEHPgbQKbQvjYol+2z/u4wB/s3LZYY+IP1c
s4UZoGpX//2ARJI/1ynGABF3MsO9Usn8wkhJBjDvPO3rwiD3d7RGigFE5FVc
7bJJ/NlVWQaIir55OJtJDEflGYCf87vBVAE5n9lrRQa47mmsBkvJfTOXrGMA
Wk+2ZMRT0k8Bse1vzsf7AZH8K+zvE2CGcpmSfkX63T7NtqdabbdrDdlX2nEt
A4wkrXtT0kDql9YoMMCpZWFPaM3EqvFyDNB7eYBh0krOe/yDNAPY/tgVHddB
4lcCVzNAhw0jsOcTOa/XPXEGONi46IxKD5kvtXIlAzRtczri00/qz9KWM4BR
7qudlUPk+Zhase+zapW8+tJREs+V5GIAGBGy2mqC7Bdo/IcOnv7u4U2fIflU
3xQdaJ3eMf1jjsw7OTxCB/kdKT06C+R8ntZ9dLBhL63pKpcA3u+/LV10kFF2
/GkTDzat8so7OngUJ3XPSZDEv1mw3y8x7sArj1diVxh9yqeD6Iufz8yJk/j+
93Q6CD/0yDBKljiJ4yb787Ji1cwzDWyZhx7mdPAz2LeXewvpv9tanw5cxzua
LLaR/MxCNTqwb76f+W0ntr1VlgAdmEaKHBW2JPPqSqh0UPP3wi67o9hBJwRy
04HB2dZN2bakX8Y79udZx/QeP3Qm9XbZp9OBPL/Qc1cfUk+7z58Oknw9skr8
sfXXuo2zf54MvY2hBRFXFbangWW1kWfjwolF49LSwNx1AenKB2S/yy83pYH3
HLwhEs9JP4FV2qlglPnAzryCnM/iK/vnI9+xDdvCqsh+cvKTKUC/0GxqqpHs
81aJlQJYjrHOzZ9J/pPtUing5is5k7AF0q9id1MiyHAvXoe4luLzaMKARFAp
uXvRNA92UL7x+kQwd/70C0dBbH2Z1KBHwEUhXw3IYldc8lV4CAyvb185bUDi
TCuDeGC/6e3Yhj2kn7Rp733g/9n+jaMJiaecC74PCrWuhTYfJHH63hdxQGag
/leuE3FU+IZY8Ge31WfHUJJvz9d3F5Twns9oriP72LhlXwcDAcE/0Bti1+Nv
rgHx8TuazHekfuvJsasgoDOv8kYniXM93RwCDHMGvhgMkfgqmv1l8N7kmFgJ
7zJsBZMl58BY5I7wB3uIZZp0zyIVEUFXx2riPMXP0Uj48bM8Rh12kOgV+3to
bu+pn4NviE1z/ruHagPQBfc2YlmXoRjk1Od2JeAb8SbpiTiUWFQXE8e3HPcT
c5pOQMsOXKXq92Hr2623TEXT31U5lx0gNrwalYo+Xf9oaHGYmHa7gf337pNN
9e12pP51ln4a2i3b3drnjh0k9XZtOgoc0xv+d5fkbzfsYX/fjRhUMYgj8fxX
ogxkqhDjfj2B1O9oM2GgVVbfp5akk/x7zGIGKqUeLIgXY1foTAZloJ+RMys3
txEHmXFmoZExqcqNH0m+v9zGLDR8YI+70hdiiStHslCfaHzdmgHS/2FjThb6
mKQbyDtD6l3ummSjDprjBq45Ele+6pWNPjhGdPz9S+I2E/HZyOXe2Qe/OQVx
v9s7+7KR4xE0OcqPXXEn4kIOcgg/YDwiIEjue8f9HGT/ZCBlaBnxNdcnOchW
crlZ3wri2WV/cpDV1+NZnauwacf8/HLRYcGpf22rST/D5PhcZLkj3LJVisyb
3VCaiw6k53O+lSP28PmZi/ad4rCtUib5CoXH8pBpfGzRSxXSr3zsfB4yrlvP
T6mSOOvWrTy0d8OBsicaxMz2Z3lot83A0jJNso/atfd5aNdtP6dibdKf48Zg
HjL4mSb0eCuJJ8StYKIdMltc8vRIfd+jtUykb9GAssH/zzejy0TbiybP0A1I
fU+SHRPp9YVVphqSfOXv7ky0VWS1RPJuEh+6HsREOj6Grx8YE4+lJjPRlqyO
1fdNSb8WGRYTaXe6nY8xI+dxHHrGRJr8HPV3zUncda6WiTZvi5W5c4C48sAH
JtJwW+9z6xCZpzfezUTqiS/e3DhMnN/6ncn+vt7vd+0Y6fft5AITqaj7vQu2
IfNthHhZaIPDMqUrdsTaf5ezkHJ0WmCAPbGFkhgLKb3S/uDrSBx5dw0LrZuq
3+DjRPo1b5FnobWKx0MuOJN9fksosZDC4ckOTxficW0VFpIPC1NzdyVOvK3G
Ql5tmoGObqTfqJAGC1Uo9jZaniP7Jn/YzELLL96RNPIg+XntmixkV6XnqudF
LCymxUJM4eFytQuk/kM0O/7H8T6PvPf/5xmx600KDQ+LXiL9zm7bxEIJHBN0
Pj+y7z5nVRYaMk+e/ONPfL9OmYV0UkwNxi6T+oNOa1ko7OfvqP+ukPgtLVkW
agOZX9uCiU31JFlIMfKQav1VMj/NZyULXfjMcfnFdVL/aECAhV6pPG7ID/v/
5yF8EQutCLBZRb9B/Mp6nokcGvhO379J+onZjzNR/qqysojbZF/rmAEmWjjt
xB14hzh39iMTmT0RsvS8S+pX32hiokReKt3pHpnfuOclE40ccZs4Ekvy2zWK
mChi5vUd8IA4ICKGiTp3Xfiy6SHJt5q+ykRKsbIbFROJRW54MdHrzQH1S1LJ
PrXr9jERp+emxZ1ZZF6a8q88ZFHx5WBjDsnnrO7OQynLb6VReaR/anBtHgKs
AZiRT+xuF5uH/L4/9DtfTvIPRinkobqtRnXOT4lL1nLnIfGIGbFjz8k8hc6B
XFSqZFGiX0Him+4wctHESe6xpbVkfuoG8Vy0o6QIcNST/N+3J3JQ1CKH21MN
xA+4G3OQavoz5a63JH9yS0AOcu32OJnVTvaRCX+XjXqPdX00+Ebyl5vuz0IC
249k2g8Rf4GiWUhLuuV84HeSv2/lp0wU/l+9wNOfpB+XjVMmUnN7Cjb9Ivlr
DpzKQFZm2kv2zxHLaytmoBD1wna3PyR/YqqHgdqmsj2yaEI4LkWzYqCAyw/S
pPmIH+hq01HGcRH37Uuwg550D6ejtzvubrNeSuLANikdyXHfaI0TIv7UwpGO
aiMvcS9bRfyU8TQViaQccf2jTExjDCQiENKivWoj8SN+5UTk4rSfS0eNOIzn
zCP0XGnPo/ObiZ9U9icgxwLtpuFtxAVS6D5iVYpodJoS7/b8HYU66Hf/zZgR
W1aJRCGOsKUNKy2IPXSPR6JDJoudzC2JVxaWRaD5lqnYGlviaxUqV9He/pbf
JeeIU64pH0WXN0T6P/YgriiT2IkKPff+zfL6//5H3+kAyX/PaY+8icsbBE+A
UdEMnpBA4vzDw74gZu8lEbNI4tBktUjQnbt60wCTWLZjcSoQmWgr/PqYeLtb
Ziow1rmr2VlAXKVunAaKqhbrNJaQ+7Y2u5cOQj+PgsIXJJ47uTkDqCx/aRr4
hsTPlTzOBZe8TrqI/sCuOL6vsQiEm4ye2vqTxPkf6BaD+wo+p+zGiafyMotB
aVu4c8Y0tn7wptASMLk1z0n7H6mPqzctA+e4Ju0PL1+B53HqcjwDJ2KCjsZt
wu5epRj+Enid5Tv6dDN2UGdw40sQsjva6osWtj3vUaFKkPIr/cjardgpkX6P
KsEXmxrLUgNs/a0Z5a+AleKyA20HSX15FN9rYFb60FjUm8Q/2ok2AJdJxwWR
S2Te4TjnBnBVXblYxI/kp9uUNYDynPI1IoFkv70PjzUCueQPY8KhxFp8WW/A
TNjy+0JxJD9xQa0ZCFa3mQjFk35CIzbNYANnEk0oAbs5RCWiGdgHqLgKJpH5
43n9zaDew2j78gzsiu7qlHcg6ejV/wRKsaN4Niu1gPI443iBctLP7v2xFtDS
IrRP4Cm2R8mH2y2A1yyldAkiHhKYbgFeBi9u8L8mXmFX2wp2b5hV420j+71L
utoGRv+4BnDNkPuqXnf4IzhqZbL+9yyZt+Rq3EdQVbShbfQ3OV/56faP4IHr
iNrHv6T+w9djXcCw4+x/+YuFcb5Vkssn8LDIw9hOBFtQaCzlC+BebjF7UAx7
bM/74S/Aw3UT3UgCu7vIRusr2C078VdzDbb9Bq/Gr2Ai8nzBEkXsWh+34G6w
19Vb/Ikm6UdtresBRdWHq1na2OZLRyZ7gJTsFi+6DnbKEKd0L5hsn224o0f2
kXL27gVJu3yvOBuS/osUlf8DUzIBA8KHSP4qRkYfSG0PLj53Hls/JODvNzBQ
rf418iLJNzmrNwhUir/yPfYhfrMkYBCURm0//tOf1D8/9ncQNBr95vO4ht38
L5B3GPx65nHcM5bEb48YjACQIx1x9z62B3hzewRcjW8qLniAXVES3DkCll1U
4Z9IJPbwOP8DyKsOFntlkPNx5heMArOU4/wXyrCzRoWZYyAmcplWzBNsmeM1
nWOgM+DF8eJn2PnH57nHgdNRyZIpClt873fHceC7ov34xVrseAUx2QlAv2ZW
4t1J5j0XLJ8Ew15/v8Z1kXhU1I9JoOaQx1/2Gbu8Pld+Cjzdzm8/20P2f/3s
7hR4O1PNf2mYPA9d74vTYM5lu73vPLbLal7bWWD8TFn/4l9yH7SO+FnwcKm4
jOcCuY8U4w+zQK9w4qsL10ry/KTNf4HA+czjVkuwZWxK9/0GXJErjm+RxNZJ
2+o6Dw51L4DNa7BrdwsWzAO6xg8pNWnsvf4Zv+bB7vaaL2vlsZvTVkf8AeEy
l+1EN2CHW5wu/QsEir7ZTm/FVnBUk6RBm8Ufto/rYXc4rN5Fg8wjlWt+AOyK
dj93GjT78+hznwG2Uk1UNQ1G7Tpg22qM/Txnhw8HXNnx3KboGLZ9Xc04J1zz
9661px92uaW3KTfc27vvyf0AbL2bP29ywws1fGIoEDulVLqRG9bfDX7Pf5U4
Jc+MB3qv9TSi38SmRT224YXN+823tD/CdvoilcQP5zUFYv8mkX1Pzg7xw7Wr
aifkU7GLWZnaS2BAH2R5MMj5reLfL4HrfdUU+VnYn378WrkUhqQtE95OkfPt
QtXLITOs3uPES+yJ+RWCgrDDLbTpxitsr8Dq7YJQdcu/G2012GM3cxMEYVfD
D5pHM/avouRjQlBzunE0rQf7y23zqRVwYM/Net5FImz7NefFi0JHE3sZRW7s
mb5N7aLws5mW9w5ebE51W1Ex2Gr5VdZPAJsemxcnBitPbPb9vhJbOfBYsjhM
vNK1rkkRe/WY9ftVUPxq/uXhddhK/7hWS8KY0Ost3MrYX2zcnCVhxG21K1AV
+z3X9b+S8NLDq2352tj13jxaa+Chsg3X7+0h8+UH66Vh81Na12Mj7B1ppYIy
0AR9UG80wV4h5KInAw2qr3xaZI599IJkrAxUb2nZ7G2FbajUbSILBX7691qd
xlZXeP1WDr5a2wilbmJ7nb6rsxYaR341Fr2NvUzqufta2Dw9YbnsDvbFZycz
18IvVRJu/6Kxc1i64uvgL8dT8V8SsIer2WeEKsmc40k52DEvJwfXw0Ju0T9x
eaT/rnAlZah7bj3PHRb2Mfc8F2W4e7v5miuF2NWyWcPK0KEr0ej4U+wm21vT
G2CcmG6adD32Ujk+JVW4+oopU6wR2ylcxlMVpg8cL1/ehI1Wf3yqCvNLQpsW
3mGPFOmYq8GGg61zXzux1yqc2KEOaVHuB1OGsA/5Sl/bBF15MxbJ8onin58l
v8w0YQNn+O2dS7ArOX+HaEKVv6dFnZdiIx6vMk34c3yjUq4Q9mGu97Ja8EJX
iYnWKuzFWU5/tOBlVnW00QbsBaWCui3wS1ampNtG7B2ChYt1oH76DXqkGja4
v81AB3LG7ytt2YydWj7yXAeGBn/otNXD7hXtL9eFUYf6Zbz2Yee6yTVug+Nm
Ndkx+7GV+dSE9OABo2yNMgvsJz0fDutBEXDW8I8ltuqPmf/0YMK66VOhdtj1
9pmLAaTPLXqc4I6d9/rbZn1YlqKw/dVdtoPUlvF1GECFBUP/l/dw/IzLcc6d
MMr25BMqFsebX/Fv3AldVmVoPX+A7VbveHUnXBWzTrUkFfuCtZKmIQwIU5bO
LMD1A56RGbvgYL+xLaMIO+Kmf9suaGl45mF6CbaObg33bqjKkSeW8gSb0RLg
sht+9dso+OAlNqveS30PNDinznHzHfZM/6qWvZDVaA5vtGBznwd8RlByg+fl
sA/YwbwN0AhOfSv4fbUTex2/HMsIMhw2TwT04P3qGj7eMYa8ltr/nRvHDogS
cTGFTdv0qg8IibF9Te+MjjmMKL84v1uYbdomQ+Fj5nCP1uNN20Sw2zf/9jeH
FapySfIS2IUHgirMYZEsr8+UDK63sA81tYDxPK1Kceps63+frvM4AE+0nLnV
uR/n33AW+HEIypgzKt9Y4PiiIR0xS/ip8cuvlwdx/aao/B2W8PBri5M5R3D+
gPjV+5Zw75MtwP84jifWvdx9GKomc42tdsdxO/2lxUfgL9eHB+0jsWVjdN4c
g3f80nmuRGHzb9rIbQ3XRuQ+S4zG9VH3bkFraJn9TL4rDsdtjkYWWcPCb12T
lsnYQ0P5KTbQzWl1jEk+2xXu6kuj7GC3bWKr9nu2U7pkRJ0coM9ZRrhlK9vd
U+tHwhzgsstMvQttbMu0buXKc4DbH72gF3zE+e58M1MOMOHj5wsq/+F44+3f
NxzhoSNSInJTuH5zeEvFCVhrnmy5VFScna8nUe/kDF2r+hGPONtBjxmRkc5w
qY6KEucqtmnMsvvlzvCA9JO5mTXYiaOXlp6CXT/eJ3cr4vxHFiFPT8HRCO7h
Ii22ZRbyeeROQ5Hqs4HHLHHcWT13vRss1yn+dugInnd9Y5WNGzyWN2e+/yjb
3WXpqlFuMOlemLyhLY5/efJ71g0qOabWbjyJ5634eqHhLNRbaF3BcZFt+z27
ToS4Q0ddvczMGDxv/1PWA094+seb6U9xeF6zmdNzT+iRetxwxQNcv3LB8qsn
vMJ3tScgke3mBvGfil7wUWed5IEMtgWHVnCWesE23yNRf8qw7enzPeeh8VMv
X4suXN8t2nnGG1qcW1Qb+hlbY7YmzhtaycWJPv/KdkXwwr+X3vBUxJOitX14
v+3r90j4wGvHaKPzIzhe42TS4APR3G3HjH9438Dtwdt9YTVLuqCLJsHutzFc
4awvbHQsWBDkYtv+3hPRRF/YVd/6yJ8HWzZ66K8v/JUg2W4uyHZ3V4jlKz+4
YMZUDF2B6w1cAqb9IDcXvPBsJbZx9fQ6f7jyjKPQWgm29T+U99z2h5LSU8et
Jdmu+KNr+dIfyrVcZ0WtYTtIkHfNtD9UDhP7Wy3Ntgx6oLs+AG7alm0yL4vz
yf/Hhf8DB1ir6g==
        "]]},
      Annotation[#, "Charting`Private`Tag$25185#1"]& ]}, {}}, {{}, {
     {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
      1.6], GeometricTransformationBox[InsetBox[
        FormBox[
         StyleBox[
          GraphicsBox[
           {EdgeForm[None], DiskBox[{0, 0}]}],
          StripOnInput->False,
          GraphicsBoxOptions->{DefaultBaseStyle->Directive[
            PointSize[0.012833333333333334`], 
            AbsoluteThickness[1.6], 
            RGBColor[1, 0, 0]]}],
         TraditionalForm], {0., 0.}, Automatic, Offset[6.272]], {{{0., 
       0.}}, {{0.6283185307179586, 0.}}, {{1.2566370614359172`, 0.}}, {{
       1.8849555921538759`, 0.}}, {{2.5132741228718345`, 0.}}, {{
       3.141592653589793, 0.}}, {{3.7699111843077517`, 0.}}, {{
       4.39822971502571, 0.}}, {{5.026548245743669, 0.}}, {{5.654866776461628,
        0.}}, {{6.283185307179586, 0.}}}]}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}}, {{}, {}}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{Automatic, 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,
  ImageSize->Large,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.76456287388422*^9, 3.764562881547727*^9}, 
   3.764562932530139*^9, 3.764562981529163*^9},
 CellLabel->"Out[38]=",ExpressionUUID->"66c96ffe-6c76-4900-8672-080a253bf775"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["", "SlideShowNavigationBar",
 CellTags->
  "SlideShowHeader",ExpressionUUID->"ced35eef-05a9-4c5b-beeb-93fbec39c1d7"],

Cell[CellGroupData[{

Cell["Runge Phenomenon", "Section",
 CellChangeTimes->{{3.7645630223709583`*^9, 3.7645630237789793`*^9}, {
  3.7645637902066517`*^9, 
  3.7645637930981417`*^9}},ExpressionUUID->"0827e0f6-3063-41ef-99a7-\
bf91ae511d9d"],

Cell["\<\
The example above is misleading. Let us consider another example. (Notice \
that it is misleading in its own way!) The idea is to examine what happens \
when we increase the number of equidistant points.\
\>", "Text",
 CellChangeTimes->{{3.764563046213896*^9, 3.7645631206906977`*^9}, {
  3.764563481872431*^9, 
  3.764563523093379*^9}},ExpressionUUID->"bf486cab-c44e-4920-9d6c-\
04f7f907337f"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"np", "=", "5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Subdivide", "[", 
    RowBox[{
     RowBox[{"-", "1"}], ",", "1", ",", "np"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"25", 
       RowBox[{"x", "^", "2"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"poly", "=", 
   RowBox[{"InterpolatingPolynomial", "[", 
    RowBox[{
     RowBox[{"N", "@", 
      RowBox[{"Transpose", "[", 
       RowBox[{"{", 
        RowBox[{"x", ",", "y"}], "}"}], "]"}]}], ",", "t"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"25", 
         RowBox[{"x", "^", "2"}]}]}], ")"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
     RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}]}], "]"}], ",", 
   RowBox[{"Plot", "[", 
    RowBox[{"poly", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.764563125983479*^9, 3.764563207264735*^9}, {
  3.764563260945092*^9, 3.764563290960785*^9}, {3.764563398071158*^9, 
  3.764563475925356*^9}},ExpressionUUID->"34c62287-5f88-4042-a6ae-\
06271f35e9c8"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJw1mHk01F/YwJEiIalIJVHIviRZ4qEQWULRokSh7JIl8bOvhYRkXxq0yRJh
hGsnS4yxDGMbW9b5hpIt3nnPed+/7vmcc8/dnvuc+3ku/31nY2sGOjq6O/R0
dP/bSgLX3M4Ohli2MySUi3DgvH4qfH0DQxDZMChZhoPCYhmhlT8YqjlMlhGs
woHsGf370/MYEjcQnOT5joPzB4IH2/sxdPSzRcgpCg7UJpZakwowtE+vkw9/
IBuMwzo+nb2HIcbDurzrbtng2Rns+LCaiq7ekYsQ08iBnp59c7WWi2g31nBL
TOYd8NY4Hzp/fAF1jIo+sFf+AF/tG3bvH55D+Ig9t6l+n8AiUubXQOAsct/S
W6k0zIe1WO9wBDPoh/iIDvFTIeDYR3gSKdNI7sYRT+m4L5B6RcaVGDGFjuK5
RukbSyDQwfKpiM4korzv0c0UL4ORaPPUikuTiEByc6KeLwOlIrMaXdVJNBIm
H6txqQyW/5gwO52dRE2aOhvMt8vA0u9K4hfeSaTC/OboSlgZqMWeLVdenkCD
ZlWC16bKYLts95p+ygQSNKi3f5pbDt67Pnq6Lo6j4Q8CdIIaFeCUZL7HYmYc
RYoIW7cZVsB9qYOv9SfGUbK1UMKzuxVwxczny5mBcXSH9a73ikcFHCk2WBxp
GEd0vK5Shh8roOT+yn3d1HF0+NwOVHJ+A2rtBYPTuuNIqcOsL3HuG1j6d57u
e09B6F2Cv2F5FZgOZKgO4yiIp0Sn52RzFVyRdbk5mU5Bu95q1m72VoHcJEfk
cjwF/cXHRRFWqoBZ23iFPYCCCh+K7PyTqob8/b21l29R0PM37lP3P1bDZvqg
OX4vBd2cSp3HvUPwumoqMdVuDFUr6TdJT9SA+dNk3/PWY6j2XJpG/0INnDl7
1ar73hiqejOfHrJaAxXvy6WYTcaQvE0h08beWhiJe9H0BMZQpJto1EWZWhC2
k/l95dAY0myWfXsroBZcJTbSl71Hkd1R+h6ycB1UjJ5hzRMYQWKJS7KyMfXg
c4H+c/KREVR/zTBEJ6UeVJIG9CPYR9CqF77UNrceaoyfv7TZGEYCDL0DDZX1
0Nw4z8nfPYymB3x5fs7WQ8+nfJ6EgGEUkjlgOaHZAFRPOWH/sSGUWT2tu8XY
CKc41C5eTycjJn7HxqGsJoiWpjyziyMjo3Tl7NTCJlgzDCz2Dyejcx7eM7ao
CdpfNZz+7EpGeikBQ+LDTeB2UIdpz2UyOnjx1yFTnmZo4jJuL8MG0evR7UOX
4pvBjtfK9Jj6IGLkSIv+E98CPSq7X8rIDyIH+48C/jktoGqe23xZbBANa/Lx
HC1tgYMZMwpuhweRAX7psU9/C1SedDzWMTuA+Msjy7eOfgf20x5jvrEDiD68
iWkz+zt8EQ23G58goTvTxMdOja3AR7DAiw2QkDyM/yzsb4VIT0Vm9x8k5Gqb
bMs41wrWDXM5THgSSn4tTh5hb4Mj5voU8WgSik9wj3hwqw3+e8V501OBhL7k
eS5WLrWB9lqqJmtUPwrp+nlcWboDRhuLTsrL9yE37EAV82wn+Fxb11AQ70N7
SPGFGdudcISiZqskQGM5lQOCLF1guNVZpMreh1Z8b88s8HVBjSz14uXpXsQj
Qrm4otsFmRmi1jcTehFLRmXARE4X3PfCfXi22oMKaNFMtiTA9p75Dp+FHiQ8
mvu3yJEAyfGyy77jPUhJrxR6vAjQXVCrFPSjB1X8sDqlFEuAS9OjbS9ye5BF
LbtzZD0BTl87Tk01pY1n9N/E9TPdMCX+WhaVEtHLgvkrjRvdIOqbQ535RETh
pAnjZ3uJ4NT59SNnFhGJRdWzqh4hwtrjPgGbF0S0xuJEXDpHBJYy7oNsFkTE
mPXrNqsrEaTVk1du7yWi72eDJo8uEsHbJKNk1awbcfg5nGTFeqAmt8DlpFE3
Cui3//x1pwcY15D4Fa1u1NuUMOXM0QvRyWPZadLdiO7M5UQu2V7IHOV/rbG7
GzEKkjokPXqhyTbbPTafgCz/OVck0vfBAd8P8pL0BPR15sxig2A/hJaeaeH4
24XoZXP3Wp7vhw3qu1srC10o9Bz7Yzadfhi/l+tTTupCFH0FkzcO/fBFHVd3
sbALXXA8Mu5S0g/Gu9MMbtzrQsNT359oXSaBMA9/B513J1rkadzu8RoAcXlM
PGelHUkPRW2cPDwENpJ+fvjhdqRlZde2V2IIMoT2d3c0t6PBnAHBvZpDwMkt
5bma0o7ajM7idN2H4O9fp9rLGu0oRYwwPdY3BHV4qslcfBv6uNvwoEn6MNxU
ofpKybei6xzRnCmqoxCksdiFf9qM+GaK7uXep0Cd9sJVuQfNiHKR31H5MQXo
9Od/5Os3o/WqWF+yHwV8TGbbcQLN6G2IXvyFNAp42Ey1RLU3oTKTptYAEgUe
hQ/XWJ5sQt7HuG0OGI2DXntH4d6WBuTtdtVqS2cCDl/Pj7nFXYeEDFVD/7sx
BV9WG/e3nK5AW2SG/C7LWYj++LZOdR2PUMHLKp0ns2Bv7uf+tQOP6Exf+1YF
z8LpJgXyWw88Kuj9bzXu3SwkvP6U+19LOSptz2/8vTALz+ReqZx1KEOaj8Ml
NDzm4KLrHfv04hLEgGuyw/zn4YSQ4onD4SVoeI6TQ+TVPGwOHCa8uFOCcN4P
rK9lzUOJeqe81+4S9Fr1b83T2nkQOnCR4fqNYlTuWvhyi24BWAqFk/ZuFqHN
gl6FZJ8FyCFPRc0G5yMJNRP6ogeLMOjDlZF2Nh8tf/fwtXdZBPYTlwuNxj8j
JmHdRvb/FsHD4n03Hj6jJPU/YewJi6A1Zcf9fOMToijnEeaaFmGGimWKunxA
vlp97/pOUeF4zMkvIyc+oNTlQPdNSSoYyhjVx3a8R2Xk4hBGJSqUP/kytSH6
HvFMXzCtu0qFiDU30bapXOSzN110yosKogwbxfZm2eiP6/fCyBYqmONEG/lY
stE0zoXQS6BCrIZZH7Ech37vOLLuDFJhI7RyTZkLh7wkzAd/LVChbZ+fCish
CxFu/VLQ3Y8B3ecCgxq/LKR3iPo8mhsDOYOxe26SWehLkt3RQj4MUmPUA4de
ZKJffN0cUVIYkESZo29FpSOmA5Xm7/UxUHO0LaydS0NLUowp/dcxeF/Q2i2i
nYZ4y3/bTJhh4CUXxb3BkIreHng0/dIWg6OqnFnJT5NQqGKa2vkgDAL9n9Qz
9CWi2eQbRZoRGMzX9UzZnU1EczwbilIvMai8/EZUmZqAWqdezcQmY3DX6HgJ
+UE8ujURoZ1YgMHCyh9F94Y4pO35reN+CQY+CZ2IXTAO2e6f/8SMp62fHNim
Pv0KXR0v+DVVi4G47x0jstYrFLxmdYS3CYNvJ+X73d7FIJMjd45JtWJQ0tbd
eJUQjdKe8qpNEDDo4QglBMZGIb7fC6IhvRismCgOf70WifgrWYPoBzCQHcv4
fawvAuma2h2JGMXASPAa3dU34WgupFMofhwDF7s9rIE3w5AUXb+G1xQGBb8d
Ts0MBiPR2bC7hDkMfiielDqWGoRu7yqQ0VjEgOpLVDK4G4hOuXq0xGAYsDWE
agWcCECtJy4dKFvC4CanW/Yj5Ic8b7/YKl3B4E3q4Lcf6v8htdV+75g/GPQJ
qRPlGp4hWe3GgEt/MThU9G4uWespEuZlZOhaw+CaMjsD/Xd3dPvN7aVzGxjE
NrrxPNR9gl6uZ2u6b2LQdZUs3dHhgt6b9/2J2sLgRElQjE2mA8LfXV/3+4dB
Soshks15iHSSpXV1tzE4MsxL3f5wHy2rlU1gNBab6FVMuGuG9rlQqh13MBBi
dXkj+dAI+UYeH2igsX/PmZ5j9SqoPBCT+kNj8f3DdrslTkPEykbtNo37lgIz
ONa0gVt/PZpC419XhZedXU1hpLP6TRqNveO0Wu3lLCDET3FElsbbvl4pGVLW
8HPHxTqLNr+/XZ4DUdQO/t4NlZmlrXeX6agKk5AztDYUX2ahcf3JiLA0DlfI
P6GXz0jbX3CCgMv5aDeweFvhQKKdhxZr5U3CPk8I1wt8EbKOAVOgibp9hBeY
aaocYqWdX8tfqshuJh8okE094LCKwXPHcM6MYF946mH9Muc3BrlLMf8dD/AH
dMstuXwZg4cPRW1KtwNBnf/sK3sqBtoCnreDVIKBKWY9Zd8CBmeG6w0MfUJg
T5C0a8gsBnPGd8/PrYfB3H5C655JDFrZPoqVKUSAc1mHGzsFg08tq3zBns+h
S8Fc8dcwBg6qMcwn/kTCqbUBh4v9GOivD23NnY2GIrf3+4qJGEiWiCyVub4E
3ZWO09udGKzFUxbOrMXABQkumdMtGOBZvHvr3OPgzkWjvi+lGGRKHCYWCMTD
BSFMLLYIg1DDgq7Uzni4bHfQQysPg+sJE20eIgmwt+t4kmQWLT78enWiw4kg
3+ukxEDL336NacT9PAlGHkroRgdiUP3Qv4rxfDKE8rYQf3lj8OJzSflITAp0
qb0UVHKi3Q9F3oJYjXTYrDfYnDDEgPVOWZ7fUjoMstxdFtGh5Zuv0UeH9Azo
5pqT1lbHwGK8zviUTyY8aeUOYJeh7e/MX/1/vVkw+Gw4+js7Bt+L72kURWQD
44WGF+fqqXChUbxeeCwbhLXxCXx4KhT2rauny+fASumM9nQ+FRLW4+DFRA4s
xomQWJKpYAUtStYq7+BraYwK3okKO63SMjxLH+AyPff+5QNUODdOfyLgZgF0
2LDdELy4CJIVN5akPxSA6tkb9bLnFkEo9nPD2HoBaD8ivhU9swj7svLOZeoX
Qoenyugg2yJsGd4c2qEWgpwMi64maQGGCwtEkMwX+Jbu3B9guwAZLuYNKmUl
0Mzxi/+13zyc+vVtXakOD7qfFQJNI2ZBBzfQYXCzDtglAttjTk/Br5B8m0bH
Ovi2NSLddnAKEh8F7SgH1YHJx9hS1l1TMCMhISOSXwcC5PTkHMokhJf7x+9i
rIe3H6ZIMhmT0NIhbFZeUA+VCVeF2o5NwuU19xl+pkYwSWPrzzg2AZr6nLv+
fG0G/HbYgLUiBaIt9y6tf28GvzS/5C0xCpDc6Ua3h5shKPPfeuoJCtilUyuY
97RAs589bmMXrT/W6nrctAWy8YH/RDrHgPQqaPzSagsMMHvwmz0cA7u+37Wx
8q0wZaG4yk/zl5f3+v2ly9ohwobVR/vGEIQuXIDp9nZ4MsfG9wSGwNfr7b+U
8XbIPx09VSw8BE5xjs+Y2Dsg+tNGkecaGQyad7mNWHVAV8QcB1MKGdilpB9G
cv4Az1W/VKXJQYjeDtefceyEGUeFmeTQAWg24GUz5CSAdJSHGRdzPxySpHC5
HicAfYj2StpqH1iw5ZyMFyLAP/GV4nNTfbDRJi5HUiJA0sO3z17X9YG4zgUz
iwcE4F73f3/Itw9iLpl9cCkhgCTPMmn3Wi/cVkjSfGXaDZxn7+w3/dsDi/yH
fbuTicC7IDjacpIICt5Zj2uziXBKIhg7f5AIgT0S1oX5RKBGHrEs3k0ErjAt
veg6IlgN9Ct0znWD6oLn0StzRLBZprgmfe2Gl6WDpbWKPXD3Fteco343yOhm
YoWkHmidVDrXHEIAdzcxy5dcfRDKVrDVz9UFXCxYpSl/HwyFydw/vq8LyjO+
HDkh3gcww3uTRKs/tloVu/LUaZx+w8OsqxMC+bWh1aEPIjjjyE5unRD9w4p3
d30fOKiw7Pla8wNyz6QPPHPuh/KEIwU8th3QP8RhbNNCgq3w0chHv1vhygtr
kls3CfjkRvjFJlqhUrHCPGiIBNLdYk1bhFbITHhgn/mLRItz03Z7QSs8MiwN
HjwyAAzt7IF8Dq2wUW9WZmA7ALLzigl+P7/D8U+5x8+zDIJv5qxZ40wLWHhd
+MmkRwbOb+a6VpzNIGrV0ks1JUPY3ti8q7uaYcXgekOvJRlw/kPV2r+bIPS0
QxbOkwx7GGIrHvQ1QV5nihngyPDU+KGjckoTrAtudnpskGHHf2bJSrgJYrvx
5VPvh8BScyxERacR6sXOP6/fPQLzET9LBQrrwVzqOm5p/whwnwjdHMyuhzXZ
x5V8R0eg7lPnrdSkehBX+kT1lhyBndenQCmoHuK1+a7J3RgB6w6tIoppPdhY
Mx3PeT8Cm6w8T/t36oA5o/9zqO4oaFpplAjQ8lxob+VMtvYYmC8ra3ierIXl
Zkm7t/pjMF6ltX6fuxaqQ7PmM4zHoPjuwVAT9lowZQyjJt8Zg/hyAbarWzUQ
smP0O8ZlDNicV3njSDUw+efnjk/iGJxS/dd8K7YG3o4f5DaZGYO8RPXnpSw1
wFtpr7k7nAK3/Sd3Bo5WwyGejQ2nKAqQX9T8ecpSDawe4YWkWAoI2U7mn96o
gi2pnGN5tLohzbFW4e1AFQziRpaMiykQ/0PUQC6pChJeGKVnjlAgWF7p+4ej
VcBudv6v8rlxiB3KOBkmVAl0GwwfXMfHQXUocin9XgVYTAys+P0ch/Nh9wU2
jSugpr1QNWphHBSdS9gttCrAP8O8593qOLzoZt+8JFEBOxoVO0MsE6BLEdV/
somH7ZePTS+fnYA8oqambRIetgTHdh0PnoDPajERfuRy+GtUZd4oOAnNYVcM
ej3KgCx375hg8BSMtjOXc3WWgMsNjgMCQ9Nge3VQyyHlC/DET95k5J0BrHBv
cRepEPLNMt+tPJ0FM96AK/M2+XBDxU+7pWoOCoFoO1v5CYrMHQWUtuaB5ckV
dzPXD7DGcmkX3mQR9j+6qXrr0jswJijL8yZSQe/woeu6njnAoryvtjKNCueb
l+eqnHOgNntQ9w6OCnET2U5yj3JA+qmXZQrtncrz65KRuJUD7HxlkTyNVNC2
H7mjr5wD3x1kJw4v096xvKF1TbocUGMWiWHXw4Be/YfFnqhsWv29djTfCIMg
s4+RSyHZUEBuztG/gYGG0di1Md9sOFFg8y3yPgZE568fW1yyYcske5rFC4Pe
lN/yDdezoQzHp8KUi8Hou8IyK95scGLDmnM/0TwV/8c7nCsbhDyrjbVontAD
U18L9mdD/JW7j0IqaZ6bmxzAzJANrkvJcbtoHmLI/V5hYhoHomZ2J3AkDAS7
fewExnBAaVD8cHEEg1eSgztWAzgwTCRVB9A8yKxaTe93Gw6YGd7r8NM8+M7j
/m21Rhwge8+eGppnodoDrK+qceDRq3XPguZtbNSKx5NlOPi//1r4///a/wFW
E3H2
        "]]},
      Annotation[#, "Charting`Private`Tag$31785#1"]& ]}, {}}, {{}, {
     {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
      1.6], GeometricTransformationBox[InsetBox[
        FormBox[
         StyleBox[
          GraphicsBox[
           {EdgeForm[None], DiskBox[{0, 0}]}],
          StripOnInput->False,
          GraphicsBoxOptions->{DefaultBaseStyle->Directive[
            PointSize[0.012833333333333334`], 
            AbsoluteThickness[1.6], 
            RGBColor[1, 0, 0]]}],
         TraditionalForm], {0., 0.}, Automatic, Offset[6.272]], {{{-1., 
       0.038461538461538464`}}, {{-0.6, 0.1}}, {{-0.2, 0.5}}, {{0.2, 0.5}}, {{
       0.6, 0.1}}, {{1., 0.038461538461538464`}}}]}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}}, {{}, {}}}, {{{}, {}, 
     TagBox[
      {GrayLevel[0], AbsoluteThickness[1.6], Opacity[1.], 
       LineBox[CompressedData["
1:eJw12nk0VP//OHBFiYSUtEiSFKkoFNFTkrRYkkqUrSxRyFuIsiXZs4dkD2Nu
lqiQPO1bSBKVdaYwZsaglKXS937O+f3+mTmPM3ee93lf9/l6vl7nnrvV2tnI
ZikXF5ck+fG/7z2wjvnv3ySKtb06pFqcBc7z24LnFybxfV3c86PHsqCoRFHm
x89JXOFLP/vpcybs26lnPcqaxMPuBoppPJlwYHXgl7beSdzaOLzhiWU6aH6d
bk0qnEQXlcanBleTwehBO3W/xSSKDxUKFY+FQ+yRvhO7TSexqjNc+KB3GHT/
ZozJnJtEaUq481vhUDjvwiO98dQkdhWNNa3RegBmFw+lcB0gz7825JrIG3+4
uosS3i44iVHhYiHChtfB413gDbsqDkpu1DEu8buN4Sz33pkyDu5UiYw5I3AH
M3mvHQko4eDBl2KJlrE+2AZ6ok/yODi3qNjHbvFHqSLRN10xHJwSUCDUCoOw
IypH4LAdB3NmBa8L8kbhNyLR/a0VB02WHnRdeBSF882hwyaXOGidY3/p7Y5o
lF7iXPrfGQ76qfIvl9COQa+bBy5RDnGwooy7Jd45DmWMmqlrhTmYOrU65nhM
Iqo7VYhm8nOw7bmSZ+3qJDwTSvjuXcbBT/ua+uSikvBObdTZkwsTKL5d8XBj
SDJ27bv42/fbBIZHXrFQcUnBu2uYJ1llE8hvsz5VdVk6dnevZNZYTWCUwv5m
PvUs9L/6oop6aQJd1jSn8dpn4Z4Z89j4CxN4p6J6R3xsFgavea5+TW8CU5o8
pCzHs1DDyOShsOoEvpzbXvr6YTbmvMtWshCewMCdfc/Lmp7iWQt9vhMrJ3Dp
Oc2/reynyDU5O7Bv+QR2ntxt6rI6B00FTz1Y/oeN1g8HpG5dzEEhvenPz8bY
WGap2Xn9Ww56tmr4/qliY8CQjf/R0Vzcbjp2drScjQr0bRvTefKwazxqZ2cp
Gw8PbJGIlcpDeb5v77Py2fjl6At1jct5SDseKn0qgY13/Hyyc9vzMLJ3/7xS
NBu1ZxLs3jDyUN1uoF0inI1xx/GwIzcFE4IUPL77szHw4s8/Dw9Q8GRDT0vS
DTaOVg75cT+m4KyxX+o9ezZ+xwxi6DkFs7/J/nfjChs11zJjzrVS8B/3XfEj
F9moVRtvJTRHwVKtbc4MbTYaUpc0Vevn4+Zq57UHxMn8/vlE6Y7nY5uIspua
GBtt7dXe+C/ko5fNwgcNETb67PpaY7aSij3892O1V5DXG7JD7McuKkadSxYx
mmGh/Uk8XWpPxcN5Fq7nJllYJGsbN+9OxYkF6S4TJguNPG8cogVS8WR6YbTF
MAtFbhQ7JKZRkZtVL3yjjYUzciUZ1E4qPtcIdXFpYiExc/ggXz8VLaMMOv+r
ZaFro0sJ/xgVK5W/PPQqY2Fye92Ghd9UdAhOm7xbwsJ+WrsNYxmBG/quGvgX
kPErDUy8hQh095sUDM5mYXaB/PeHWwmU/lDqFJbGQp9fcgur5Qj8sN2rIzKZ
hTxXxq6p7iNQ4S1PZHwUCzkH/I/cPkLg4ObWicQwFipIVvLF6BIY7vJQLyWI
hcK7XRYMDQgcF924KusOCyW/porWmhGYaD90PceDhc2XFNPdrAjUeZ3dRnFl
4ejj62s6bAmcWeUg/+wGCwdrdNI6HAnMtNwbXmTPQgG1Z/+5uxBoWDLDKrnC
QvGKwOoGNwL/Las49cqchVePUF5UehJYYOJLrbhIjs8GToDZHQIvUbVXVhmz
MKFns0mmL4H8i3yONQYspJ9buBodQGC54bvW+pMsNEX1zzvvE2ifFSfXfIyF
vwLjv9k/IHDdr4uhbzXJ8bAmXp4JIbBedwuz4xALBc0lfb+GEuj6+NuJLhUW
ljY9uCERTqAkh0L5qMjCIFWX4qURBL7TdOb7LE/+/sjcK4q0T6zStf4dLAwU
fjfbRFp+dL55SIqFOuLKFymkvxzEnV83s9Dk7u92JdLBYYHBo+tZKKcZeN+B
jK8yeIIxvoaFKq93Fp4MI3BEQUh3QpAcn4eqXh/J/OLudedO8bHQNmdRdGUw
gVo9SbwzPCxsVH5VzCKvb2qnhd3sPyYKVybd8bhHYJq3dNPCAhO5lnY/LfAj
UL9jXGbxJxODN4Zdir5L4B/JwqAl00xMt+57t8WLwPz/3EZ52Exc5/9pj6k7
gSaNqjorxpgobxeRp+lK4AvH+mVCA0yMqpbIE7lG4EFJJ8WcXiZ2mdKCwq8Q
WNEtdlmji4krulc+zL5MYLWGY6ljExMrz3ziLzIk8/2+dnhpLRNFZP78zDpB
YENO1crkSiba3g6bV9EisFVI5EpzMRNPRzsePETWo17960gLgomDeQNMQpbA
Tk+bil85ZP4jf7OqJAn8SCtbvT2Fieebl6+vW0XgcIlFtX8QEwUzmwPPf6Wi
tT0fe70/Ge9NJLOjl4oj4iViRd5MNL50THH8LRWZ95c7D7owcaN8tR5PCRV/
Xni2+ZAZEyt62Mp9d6joKXDhRNc5cvyOhK+75kTFhWquW9cMmaiUY1qYbEFF
LjnjtkfHmOjndnhmXJOK/H8WvGb2MnGPWvx++JePEUXZOWFyTLTndjs8PZGP
wjb6XVLbmRhzofKjbn8+inZkyJ7ZyER9rIkZf5WPkum6n57xMLH+b2D2rEM+
KmsnKNt9HsfA2/c5ntUULJsDq8UP42ipnfg7lUJBtWfj4fEd4+jJ7aPkFENB
zXUa3+rqxjE94ZcQ3xWyXzK+xkgWjKORvO1LFS4KWoYrTn0OGMey13+GsuTy
0Ks0Jqj97jgG7S5yMxPOw/j+H+I1nuT/leYDWn/mYqv8yxN5TuMonv+n/h3m
4v52tSx303F0t30y7m2Yi8sFtc+v2TeOXevoEwFXc5B4eL7yNI2Bmev/NoaZ
ZGPjq1dGmv0MtH0vfiFMORtpQ+vH9/cysMNkdK2uSDauU+gT3dTOwPPGce3z
rVno22nhPF7GQHFpOYmjqllotPqaVFAUA00oeW5xfJk4F+MdjMDAuLUjJR9m
UnFMvIapo8bAlHez36IfpWJPznK9DiUGSl7eoFqtmoolFdGrB2QZOFdkzOd8
5wk60XOTF9Yw8BJvuffW34/xq2I3ocwYQwHnthG1wUTs6JR/T40aQ9MXiZHa
h2MwS3BwQyJtFNM1NtH3sR0xxgx4VftHUfPMmYh2IXsMyEuf+dwzigqN6zW+
7b2K1lpX3m1qG8XhZW5Hpd+aoJQ7IzDt5Sj6PfFXbGpQg8yB75M5YaOY6Dw/
mlR5E9IJ3sYXyqP4pX6aWyAiBFJOKrp+CBlBn+1pBiuepELAdStP2RPfcMlb
SwfZgSLw5s73cJ2g41hHePFa2muw8nsn3ZNHw6VTja36DnUQ/2YkMcVhGOtG
BRWuzDZDxdBOAUJqEPWSbL2PC72DbcKaWsapfSjTJvBi3KQLnssFO9C/fsKC
V+1fjhl+hKGGYkkVlR6UfLTwhbb2E4zIx+/Dlx9w7nmBms/dL7Dah6KyZ8l7
7L8t0bxqqB/kVSbln/5ow4J7NlyxvYNgu8fXt3ygDbPlH+4OGRqENBmhrvam
NhwUGT7kNTYIImJ7PX49bsOy1Wk/Ts8OwuysU81x7TZcJ+kbX7puCGrLOeeY
cW9x52O73lDjITDR4PjsVWlFfUHV08fahuCe9kRnuWcTCmpcul13bhhqddkG
Slea0PSEUzz3xWHg0mN1FOg1YYXh5X+al4bhzrnxtiypJjyVklZeYj0M7rYj
zRFtjdjxwLvY03kY7IMHqq0kG/GA3d5/vsHDcLqtvYivuR43DR2o8ikbhtDO
tr2Bz+vxYcAJMbvXw9Dc/bZgMaUeRX35c05XDcOxgRbix816PKveaCVUPwyH
OQ15A5vqMVfG7/B/74ZBQbgqvdi5Dn9ban8vGRkGUeOCqItitXjXdzF2rTAN
nsuEVRzjrsVlPP/C3EVoYDhv901xsgZT/OncPWtpEJq69SB/Uw1GvDG6GbGB
BouMuKEK9xrc+WVeaVCKBqN3vfeI91Sjusix78uVaeBkEJ34RhZxd+PB+ZsX
aJA3c5IdvwRRbZMh172LNKAn8Wg6fa7C6/WXPaPNaHD+m8eYREgVFnOyynMs
aAC3zVX8GG8wamLa7bkdDYSzdnVr5Vai1Qnxm/oeZL6/GoSapStw44rDxfFx
NIjMz6w9PF+Ol+X2B95OoIGjue+tF+3leP6oko5pIg2kGw/2ZbqXo9xsYLxI
Cg0S4qk5d5vLcPrHC2PzLBp4KUVr7L/+CiOOhn23KibzYdyYomi+wrYHvX5r
S2iwP+VklqToK2xq3T9bV0oDNjcPn2DVS5Q455m8oYwG5h/cuxlCL3F8ergw
rYoGWq6XHFNLSrE8Xqlq3VsaSMioSogGl+KRn/Qnj9po8Puz6PuwS6V4dJmj
hWgHDUqPvFO5vawUz1Lig3jf00BmtdZS4wslaGRlbtbYQwP+oh1JfL+LsUwq
+v0zGg2e9o1EjAcWoFFLRuHWnzT4cmdd2pP9BfhJjGOx6xcNBCWOF52hP8P+
RIrBvlkauFvmdZXDMzxdHcU+ME8DnREHsdAFKhZY7Vsn9ZcGDM5kupwLBflK
jv2+yUMH8SjJ54MSFBStbw4xWEYHQ8UzdTHkPnmPQYix3HI6lP33fGSB7NsU
Dn9qLy8dQubc5N6O5OBfluwf4ZV0kFu6UOJolo3EgyhcFKaDeZZcwxb+bFxt
Gf+3cDUdYrTNej6UZeHAHvVb5iJ0WAiqnDu0Lgu/XxgQLllDh7crfTUE3mcg
6kX+VFtHB65nhfrVvhkYth//fCStpD9s4bYnA19k3FF3EqNDStSRgP6wdOQO
2myTsJ4On+RWRF6MSMXrRo8vNGykg+aNa0U1zCcYfLRzy8lNdMgrbO2S1X2C
sv0ZW9pI31aKEFtYmoLeqgUNTeJ02HhYJCPZMwldf8XkPZGgQ4Dff3VLexLR
2LRTTmALHVi13SMO+xNxg8jwVw/SlccfyR3iJKD5bfEFHUk6XD4jXtp3JQ5F
x0pK67bSgf3jp+qt+ljk+R2jslGKDncS3qHg9lhcH2LN70Q6pS/g7ZHRaDQd
fpi7ahsd5H0unenTicZah3aXi6RfS6r0uuVGYdV0RVgm6dK3XQ0G7yNxu5Pj
S3lpOnQLB70PiInAZi+Rxhukf5xTHXhxNhxvPlfcQZDeN5w2s6knBB805QxI
bqfDme1nuQweBaOCyMi+C6RdHJYLBJg8QCPr659DSRfOXN/G+BKIZV/lNzJJ
d6hK7t2Ucg9VPxmVisrQgePzQU3/cgDmmjxOO0x6VX2Qjr+EPwZuXey7StpE
xC3bHn1RwOu/a8GkH6V8ed1x5C6C0sxxCukemSMflOq9sJe47dFEem1xLjNZ
xxPPv/07Syd99pDg0iUtt/D3Eu+W36RjGtw22J36D2v12KzVO+jQadCn0N7u
gg4qehbbSUuU3ouyTb+Obh1J21RIP242xH1P7XBzScdBbdLrBzZzFinWaCMw
nmpAetfXj6oJl82QL2L0nAlpGQGXR3vszmC4a8Nlc9J+3Tu7N9Vp4K77gS+t
SMsLDTgs2y0NY2LbLK1J90wHpAnP6YLfigxTS9JTBju+O7ueh/9w8akZae9Y
nVZHJUvgUdDQMia96HP7cdpeG9j76fLek/+L70Bc/yDnAGXV1jc0SHOfH9Lg
lXEGibv6f3aTrpMMefBE2BUydkn3byIdmCDlciDSDY7eGhPgJa0jUGnyfqUH
1ESkRE6S48MbcO6IY8htyNiva/WRdPMsR3YZ7x2QKpoMLCMdeiNYJC3QB471
Js4nks6Zjror7u8Hrbwn0Z20nZ2c7cvFAGh8MawhR1pXysP0nkYg+OV8/PWP
vP87B+r0De/ch60Pprk+kGYaXT7AnH8AnJnaFTdJt67K3/XqYAiYnQpfrk6a
2vxrS6BHKKhJZ5/jIX39cNQKiZ/hsGmUPRZO1qPefP8f5v5IuMpR3qVHek+p
7PQr14eQ8+ov8pOei6Oxd85FwY8aGLxL1ns5v/fH2luxoOq/Y/IoOV/Sd4t+
KJSKgxZhXtsxcn4FGRZ2pryLg4Dj9188IG2c8PWtu2wCzN0oH6wi5+PU1tO1
cgOJQFvy7BIvOV97tUdRLDQJrrB3Po4m53eVnd8bngPJ8C1EsGQ96bBnpWWD
UY9BLTfUefNmsj5UNxfGaKeC1WyyxBKyXwhcekX4TqdCWH3JLleyv/zwOZN/
PTUNjl1SXDO8gQ6W9FqjbXfSoahi9OZzsh/t2Tmr9/djBsQyVFwVyf7VUmKh
XRySDT9Pt6Qnkv1RvUG+bsdwNnzmSNj2CtGhqGf+SKrKUyh7ZrFiDemE+VgI
+/oUvC+bf/NbRYer0Kxmo5ELGeFP/qnw0+Ffq4LihmkK6HItnRLgpoNr/9/C
h8fzIWkx+fGGpXQYnWjdszw1Hywf2+zYtoQObcI28j9OUkHikfew/D8aJJ9P
2tGeTYDW9PGhjb9poExfIuFvUgh+SXU6Vt9psKfiwrQCpRBEdRtrDKbJ9Snm
Wf3wfCH8Cq19cGiKBiszCOV0vSJI1In8yc+hwR9Dk/5/nCIIaxwyDxmnwUBR
oSwqPgfLoudPeIZpkOZiXq/xqhT2a6vF/WqlQcolvoL4oVJwufaBEdlCg0Td
0kcTvC9AvKYwSrqZBlGS/I5PTF5Ag6VV3/EGGvh1vlj9d/4FeD5fY+OENLBS
WGVRqf4K1thORB4i1+ttU6/n1WrLQWjblsBP5Pp/p/1o4MxEOShvzfb+Re4P
evLfChZsqADTvzpmq+PJ/Y1N37atNytA60VDt0Y0Daa+LOjxSr6GEgbnzc0Q
GlQ1qGZ+uFsJ+g61H4/dpoHp41cnr6si2F5Wp2aeJ9dvT/i43Qbh8+WQJ+bG
5Pp7vsliKAohssEpWMyIBjXCvW5GDIS1VTon/fRosPaot0HKqWpYPb8vcrc2
DbzLkm8HLFTDqqftSbMKNDiR9bld36QWNjh2/dNbQeZzv8C24UYt9L7ZEfpy
GTk+9vf+HbpXC11H+NTFucn1e/duRdmCWnArmFk/9HcYgsv84rh56qC8QKFF
6we5X2zfYVZWWAcbOVWxAQPDcHzuFmMrbwMI+Ed7MwrJ/aOeCPfPF03AGpjV
u3V6GCKt+KbnW5qg0XgWbU4Mw6dbXEOLA03AXyczfFZnGBxSORUrljeD8JFX
u2U1yeMnW13FzzdD7ztZ4bL95PHR9+hHfzWDTfTdYLuN5PE9MzUxKq2gtvuo
nDS5f35o0eun8KoNXIteV3tPDEKT/uZVhiLvYe7lmen3nf0wsVXUpyv5A+wb
OW9w6tEXuOW2y+rhuh7guXpB2kfpE/T2CxvZNn8CuvXhx543PoLlbfUx3tN9
4OP1bffXki6o23UgtG7ZIAiZ2T/yNnwHMnyVjGzdYYgwMgh4GNUCmysdjy0L
psGax670Q7/rgGthKcWVTvYJnxz9txqVMHvmjXnD9m9wPm/vSrHpYuhTsti0
PXAE8g8lV+qopYHtvLLQ/ZAReLX9gIvv01SYeiPAPRI5AvobrGnCq1OB53gF
MztpBKTW374F4ymw56JohXThCPj0uZR/fJIMfnffmkh/GYFbqT1ZqusTQKbx
QMI2hVFwCLAZ/JUXCi4XhFdL9Y/C8H3LYu1/dmidtzd5O430VL3yimRHPDev
v012dBQ070muK1d0RrXHEcoKU6OQd6HVkr3rP1w2xG96mGcMqsPEtK6CFybb
LssylR+DKZlSr1OrArH+1oJS7J0xqLrQFKuVE4UvGzdUJfiPAc/pn0+LBaKR
IqZ6PDloDAIfhK4vdY3GyHKPixlRY7BRbqjzpkYMXvw7c7cwewy0kxbsmtpj
kRM42fi2bQxK4/rroj8n4Ia4byY8mxnwqbpVzvhxCh68t2WdmxQDtEO3J/xa
SMHzrqYfvu5ggAAl7KOA6ROMMezUq9vHALa12qbYdakosKpSK0CXAVfluZsX
gtNw8X6sPJcbAwpym7QcZTNQ/FbHuPNtBpwUT7WsMc7AQ1f5cod8GCC4aDfX
Tu4TPbV8t2IIA55/lirj+5iB04sO63zSGNBGeWKz7E4mfvPQ4vrTygC/W5Gv
+19n4VK7u28cO8njzbZd9aNnoeT5Mq++jwwIfXOKv3dFNl5W2vOzYpgB102i
nFKMs7FnagPz9i8GzOW2OPkzsrHl2vSH2a3j4JOZ9S1gSQ4WmKXn/vAcB6MT
Gh8//spFpd5tC8t9xkF8i9764NV5+Noo9/TGe+Pg5dO5cnpXHjafKJjWjBiH
fiU93VGLPPx6oPJQZMY4aMZ2TZvU5+H6tZ/e7Wwdh+bTS6uUfSmY+tBMSv3d
OGg/2j6lHk/B7SuH3Ay6x0EkZWJwNp+C+5aObnAfHIfRPwpKwR8pqDf140rd
93GIaRl4sGpHPga0Cc1d3sSE0fRxXv/X+bhCN/bkTUkm9LhtyHnWlo8P69Y9
CdzOBAX7Ru6wgXxMqdisRd3LhEuZyzUuLObjq7xd4XNHmSBYsNlKQZ2Kh6UL
B1eeYELb/lsvBk5SsT5tn+IWfSYsR6Xjahep2JWg2nPMhAnsmr1WXLeoyA48
Lhl7nQlKwiOSr/Oo+B/XW9ecm0ywb4wQGSil4oK3fkO5OxM0lfe/IaqpyPvf
OYdhPya4/Qv3Mu6lYiTnU+WP+2R885e75OlUXOtwSYg3jAnV3o1dRWwqSlld
Kd0dzwRLOQmo5iIwr390+ZFkJsjNPcw8yU/gXhOHi8ZpTBhcVD/ov4ZAdYOb
i14UJsxYqKh8lyawrnXGMLKACcMK1uryuwk8qeOZlVHChLIdwUWrlQm8oOGr
2/yGCe5ZGyUGjxLY/3NUqLaWCY3Rb9oyTxJoXaDf+7qJCUTmgofoGQLHbF+m
vmgj49HNp5UvEHh9i4Rt4XvyfJVO/DyXCfzee383pYcJofNP43ytCfSImpjJ
7GPC9yoPqzw7AgOWvLmXMMIEp7PhK5bcJHBFhfSpKCaZj3bofYVbBEa4houE
TpLjx50UKnCbQJFdM5/vzZDjF8Y1HneHwMSvZhl358n7pcLl2uJLoERKnb3H
IhnvdccvSgCBWca7FG5ys2BQIVlX5T6BO1fFzjqsYEE9X/BBpwcEFjQsVF1d
xQJ7sReP/vd8WsnHOshchAV/gpUPDIcSWK7SqmcixoIo1lrGhnACD08qihqJ
s0BQ7rz3Iun63KT+U1tZcD7iR01EBIEnLJdkH5NhgXjlYPj/nk93rL/mCLtY
MFfB9YxK2vh95z5VBRa4slV/qJL+HHJwYZ8yC34JWyu4kvHMtdJr5NVY8FxH
V9k4jMBvC7whMsCCUb6SriEyv2slzoaS2ixIVLL8vDaYQI5jr9jGEyx4uciz
8Qd5fW7SMLRGnwXrFvVsvO4RON+fk7PqLAvcZIbC8v0I9IkXdOI1YYF5y1nz
e3cJ5NF3V15ymQWmCYef8XgRGLp88M+CFQtak2VP7HEnUAiP1c/YsmDmSMJS
LlcC4zyehXEcyXyN/1Ju3yBwo4LoWYYLC5gLHUKP7QlMY9zZSL9F5lPmynP1
CoHSGd9ofV4skLp53KSbrIe9IqUu7wJZULav2zyfrJ/S1k0HW0JYwO6ReC50
ikC1e/f+1UaS8an4dJM2gdozRpEvE1lgu6O1QIasz1ai4lzRExYsXPC/up2s
XwMbqc35mSyobJf93ULWt2nPNPUJwYJvT0NX8pD1Pxx58b9HxSyI8z+1PYmP
QJvjNWrRL1mwPF7Q8sM/KrqURbUEVrPgurmepwqLikHJe8ccu1kQ8C13svgN
FQXOPiqw+cyCqXPF7Q7FVIxe+e+WxSCZr4/mqaZsKqbc6eA5y2ABV2HEFYdQ
Khab35BS+8MCnvWGb86eoaL1i6nN6UvY0OpphR81yfks4LZhOS8bNoY5Hxjf
S0X3ci/hrtVs0Pq0NbpOgIpqa4MXr+1gQ1tVxmaoyUeWw8qFd/Js0BQpXKrz
jOxXNZE/lfexQeZQrWN/Yj4uOsWzl2iw4fqFTZ0vnfKxtiXzS5IRG3SGv/Jt
Es3H/yRlev5dYENBTdWBR38pKO1BeW9zmQ0LK6gz6SMUDJIuala0Z0P/CVPJ
S6UUPOX75mXLXTaIHBi75HCagt1Kn2Lm89hgnSW2r9A2D++HmUVaFLBBe5lh
3ejxPFShD4Y0lLChe3/ZVO7OPEx8OOIXXcWGua+VWxsZuXiZ+d1JtpsNPLvu
jDhcycXRtFWnLv5jw8H5l2sktHLw0a8oHeSZAP4ld/zeb8pBXb21Wtv5JyCO
kxjK//Mp5i9sVJ1eOwEp3f63BHKeovO5nTtC5Cbgk/myZN2lT3GO/yh3+bkJ
UOqChGtZWSjCEYyrNZ2AgtYA2ZueWbjr/RfpNosJ8LMXo/LrZaH5I1edoWsT
YClT+vT7TCbWbcsMWeYzAcxW5YIoyMRI9SXCRjkT0BZkn9yD6Zi3pT3djDoB
C8e6Vb5FpGPN0iRFm6IJUGM5O/uZpeNMs8JZz4oJqFxHr2pipqHpOcuE1A4y
31uO1RnHU3G7E4qzfk3A6Z6qsZ09yXj4TNizmd8TkMedYztyJRlNlC4cXuTi
gMmKILuyqSQMXeBYrF7JgeYl4pDOm4RTQRJZB7ZwQKJH7Kaj3COsTLsrG3ic
A53lYVzRR2PR6P0hlc2JHLheuGKPesB95D+0sqbyCQc+TR+4t0QzEGuyv5y6
lMUBdbHAd+a/A1DB87bV4wIOGB6f9e938EPBLa/CNzRwYGfmziv3L7pjy/V9
X0W/c4BL61K8fsJm1FwhGyV4ehICOhL7lgQ9gLmbcxsLzkwCj8eIaZxYCBT2
NT3VuzAJay3P5b+ghoJEoe3rcOtJUFE9lkHpi4A/57JH+W9PAndVRp767hh4
lbVFgzdnEkZdaqttmIngtGqyKYc6Cb9aovQnm5JAxqPKSKd4EoIi/SRbc5Ih
7uRl+/uVk9C5KXOlu30KuE4nx3J/mITWERmrst9pIGfmIJH1aRKqgirKODfT
gVavStEanAT5Q4ctiNF0MEz8VOU/PgnWQhE6cl0ZsGJp3omtk5PwfefaHSzd
TEBHj+7qmUnI6nG2ManOBPePOhaWC5MQfTLjjOvBLPh/73vA/3/f4/8AawIB
7w==
        "]]},
      Annotation[#, "Charting`Private`Tag$31869#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,
  ImageSize->Large,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{-1, 1}, {0., 0.9999994076421085}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.7645631983389606`*^9, 3.7645632077386293`*^9}, {
  3.764563276309475*^9, 3.764563291465564*^9}, {3.764563398772921*^9, 
  3.764563455482463*^9}},
 CellLabel->"Out[82]=",ExpressionUUID->"e1d77504-d688-43ce-9669-7d652054b834"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["", "SlideShowNavigationBar",
 CellTags->
  "SlideShowHeader",ExpressionUUID->"1f165dea-40d2-49ed-a901-a80aa8d15f1f"],

Cell[CellGroupData[{

Cell["Runge Phenomenon (Solution)", "Section",
 CellChangeTimes->{{3.7645630223709583`*^9, 3.7645630237789793`*^9}, {
  3.7645635520509768`*^9, 3.764563557173545*^9}, {3.764563797140655*^9, 
  3.764563800492836*^9}},ExpressionUUID->"1f0952eb-56fc-424e-b45b-\
c771242e2c61"],

Cell["\<\
Now repeat the process above. By choosing the interpolation points \
appropriately, one can control the oscillations (and thus the error).\
\>", "Text",
 CellChangeTimes->{{3.764563719380584*^9, 
  3.764563781890662*^9}},ExpressionUUID->"09eecaae-e115-4c10-991c-\
404b57942b51"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"np", "=", "5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"xp", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Pi", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "k"}], "-", "1"}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"2", "np"}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"k", ",", "1", ",", "np"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Cos", "[", "xp", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"25", 
       RowBox[{"x", "^", "2"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"poly", "=", 
   RowBox[{"InterpolatingPolynomial", "[", 
    RowBox[{
     RowBox[{"N", "@", 
      RowBox[{"Transpose", "[", 
       RowBox[{"{", 
        RowBox[{"x", ",", "y"}], "}"}], "]"}]}], ",", "t"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"25", 
         RowBox[{"x", "^", "2"}]}]}], ")"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "]"}], ",", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"x", ",", "y"}], "}"}], "]"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", 
     RowBox[{"PlotMarkers", "\[Rule]", "Automatic"}]}], "]"}], ",", 
   RowBox[{"Plot", "[", 
    RowBox[{"poly", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Black"}]}], "]"}], ",", 
   RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.7645636499224*^9, 
  3.764563710884986*^9}},ExpressionUUID->"25942195-732a-4c7a-83bb-\
aafce7829194"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJw1mHk01F/YwJEiIalIJVHIviRZ4qEQWULRokSh7JIl8bOvhYRkXxq0yRJh
hGsnS4yxDGMbW9b5hpIt3nnPed+/7vmcc8/dnvuc+3ku/31nY2sGOjq6O/R0
dP/bSgLX3M4Ohli2MySUi3DgvH4qfH0DQxDZMChZhoPCYhmhlT8YqjlMlhGs
woHsGf370/MYEjcQnOT5joPzB4IH2/sxdPSzRcgpCg7UJpZakwowtE+vkw9/
IBuMwzo+nb2HIcbDurzrbtng2Rns+LCaiq7ekYsQ08iBnp59c7WWi2g31nBL
TOYd8NY4Hzp/fAF1jIo+sFf+AF/tG3bvH55D+Ig9t6l+n8AiUubXQOAsct/S
W6k0zIe1WO9wBDPoh/iIDvFTIeDYR3gSKdNI7sYRT+m4L5B6RcaVGDGFjuK5
RukbSyDQwfKpiM4korzv0c0UL4ORaPPUikuTiEByc6KeLwOlIrMaXdVJNBIm
H6txqQyW/5gwO52dRE2aOhvMt8vA0u9K4hfeSaTC/OboSlgZqMWeLVdenkCD
ZlWC16bKYLts95p+ygQSNKi3f5pbDt67Pnq6Lo6j4Q8CdIIaFeCUZL7HYmYc
RYoIW7cZVsB9qYOv9SfGUbK1UMKzuxVwxczny5mBcXSH9a73ikcFHCk2WBxp
GEd0vK5Shh8roOT+yn3d1HF0+NwOVHJ+A2rtBYPTuuNIqcOsL3HuG1j6d57u
e09B6F2Cv2F5FZgOZKgO4yiIp0Sn52RzFVyRdbk5mU5Bu95q1m72VoHcJEfk
cjwF/cXHRRFWqoBZ23iFPYCCCh+K7PyTqob8/b21l29R0PM37lP3P1bDZvqg
OX4vBd2cSp3HvUPwumoqMdVuDFUr6TdJT9SA+dNk3/PWY6j2XJpG/0INnDl7
1ar73hiqejOfHrJaAxXvy6WYTcaQvE0h08beWhiJe9H0BMZQpJto1EWZWhC2
k/l95dAY0myWfXsroBZcJTbSl71Hkd1R+h6ycB1UjJ5hzRMYQWKJS7KyMfXg
c4H+c/KREVR/zTBEJ6UeVJIG9CPYR9CqF77UNrceaoyfv7TZGEYCDL0DDZX1
0Nw4z8nfPYymB3x5fs7WQ8+nfJ6EgGEUkjlgOaHZAFRPOWH/sSGUWT2tu8XY
CKc41C5eTycjJn7HxqGsJoiWpjyziyMjo3Tl7NTCJlgzDCz2Dyejcx7eM7ao
CdpfNZz+7EpGeikBQ+LDTeB2UIdpz2UyOnjx1yFTnmZo4jJuL8MG0evR7UOX
4pvBjtfK9Jj6IGLkSIv+E98CPSq7X8rIDyIH+48C/jktoGqe23xZbBANa/Lx
HC1tgYMZMwpuhweRAX7psU9/C1SedDzWMTuA+Msjy7eOfgf20x5jvrEDiD68
iWkz+zt8EQ23G58goTvTxMdOja3AR7DAiw2QkDyM/yzsb4VIT0Vm9x8k5Gqb
bMs41wrWDXM5THgSSn4tTh5hb4Mj5voU8WgSik9wj3hwqw3+e8V501OBhL7k
eS5WLrWB9lqqJmtUPwrp+nlcWboDRhuLTsrL9yE37EAV82wn+Fxb11AQ70N7
SPGFGdudcISiZqskQGM5lQOCLF1guNVZpMreh1Z8b88s8HVBjSz14uXpXsQj
Qrm4otsFmRmi1jcTehFLRmXARE4X3PfCfXi22oMKaNFMtiTA9p75Dp+FHiQ8
mvu3yJEAyfGyy77jPUhJrxR6vAjQXVCrFPSjB1X8sDqlFEuAS9OjbS9ye5BF
LbtzZD0BTl87Tk01pY1n9N/E9TPdMCX+WhaVEtHLgvkrjRvdIOqbQ535RETh
pAnjZ3uJ4NT59SNnFhGJRdWzqh4hwtrjPgGbF0S0xuJEXDpHBJYy7oNsFkTE
mPXrNqsrEaTVk1du7yWi72eDJo8uEsHbJKNk1awbcfg5nGTFeqAmt8DlpFE3
Cui3//x1pwcY15D4Fa1u1NuUMOXM0QvRyWPZadLdiO7M5UQu2V7IHOV/rbG7
GzEKkjokPXqhyTbbPTafgCz/OVck0vfBAd8P8pL0BPR15sxig2A/hJaeaeH4
24XoZXP3Wp7vhw3qu1srC10o9Bz7Yzadfhi/l+tTTupCFH0FkzcO/fBFHVd3
sbALXXA8Mu5S0g/Gu9MMbtzrQsNT359oXSaBMA9/B513J1rkadzu8RoAcXlM
PGelHUkPRW2cPDwENpJ+fvjhdqRlZde2V2IIMoT2d3c0t6PBnAHBvZpDwMkt
5bma0o7ajM7idN2H4O9fp9rLGu0oRYwwPdY3BHV4qslcfBv6uNvwoEn6MNxU
ofpKybei6xzRnCmqoxCksdiFf9qM+GaK7uXep0Cd9sJVuQfNiHKR31H5MQXo
9Od/5Os3o/WqWF+yHwV8TGbbcQLN6G2IXvyFNAp42Ey1RLU3oTKTptYAEgUe
hQ/XWJ5sQt7HuG0OGI2DXntH4d6WBuTtdtVqS2cCDl/Pj7nFXYeEDFVD/7sx
BV9WG/e3nK5AW2SG/C7LWYj++LZOdR2PUMHLKp0ns2Bv7uf+tQOP6Exf+1YF
z8LpJgXyWw88Kuj9bzXu3SwkvP6U+19LOSptz2/8vTALz+ReqZx1KEOaj8Ml
NDzm4KLrHfv04hLEgGuyw/zn4YSQ4onD4SVoeI6TQ+TVPGwOHCa8uFOCcN4P
rK9lzUOJeqe81+4S9Fr1b83T2nkQOnCR4fqNYlTuWvhyi24BWAqFk/ZuFqHN
gl6FZJ8FyCFPRc0G5yMJNRP6ogeLMOjDlZF2Nh8tf/fwtXdZBPYTlwuNxj8j
JmHdRvb/FsHD4n03Hj6jJPU/YewJi6A1Zcf9fOMToijnEeaaFmGGimWKunxA
vlp97/pOUeF4zMkvIyc+oNTlQPdNSSoYyhjVx3a8R2Xk4hBGJSqUP/kytSH6
HvFMXzCtu0qFiDU30bapXOSzN110yosKogwbxfZm2eiP6/fCyBYqmONEG/lY
stE0zoXQS6BCrIZZH7Ech37vOLLuDFJhI7RyTZkLh7wkzAd/LVChbZ+fCish
CxFu/VLQ3Y8B3ecCgxq/LKR3iPo8mhsDOYOxe26SWehLkt3RQj4MUmPUA4de
ZKJffN0cUVIYkESZo29FpSOmA5Xm7/UxUHO0LaydS0NLUowp/dcxeF/Q2i2i
nYZ4y3/bTJhh4CUXxb3BkIreHng0/dIWg6OqnFnJT5NQqGKa2vkgDAL9n9Qz
9CWi2eQbRZoRGMzX9UzZnU1EczwbilIvMai8/EZUmZqAWqdezcQmY3DX6HgJ
+UE8ujURoZ1YgMHCyh9F94Y4pO35reN+CQY+CZ2IXTAO2e6f/8SMp62fHNim
Pv0KXR0v+DVVi4G47x0jstYrFLxmdYS3CYNvJ+X73d7FIJMjd45JtWJQ0tbd
eJUQjdKe8qpNEDDo4QglBMZGIb7fC6IhvRismCgOf70WifgrWYPoBzCQHcv4
fawvAuma2h2JGMXASPAa3dU34WgupFMofhwDF7s9rIE3w5AUXb+G1xQGBb8d
Ts0MBiPR2bC7hDkMfiielDqWGoRu7yqQ0VjEgOpLVDK4G4hOuXq0xGAYsDWE
agWcCECtJy4dKFvC4CanW/Yj5Ic8b7/YKl3B4E3q4Lcf6v8htdV+75g/GPQJ
qRPlGp4hWe3GgEt/MThU9G4uWespEuZlZOhaw+CaMjsD/Xd3dPvN7aVzGxjE
NrrxPNR9gl6uZ2u6b2LQdZUs3dHhgt6b9/2J2sLgRElQjE2mA8LfXV/3+4dB
Soshks15iHSSpXV1tzE4MsxL3f5wHy2rlU1gNBab6FVMuGuG9rlQqh13MBBi
dXkj+dAI+UYeH2igsX/PmZ5j9SqoPBCT+kNj8f3DdrslTkPEykbtNo37lgIz
ONa0gVt/PZpC419XhZedXU1hpLP6TRqNveO0Wu3lLCDET3FElsbbvl4pGVLW
8HPHxTqLNr+/XZ4DUdQO/t4NlZmlrXeX6agKk5AztDYUX2ahcf3JiLA0DlfI
P6GXz0jbX3CCgMv5aDeweFvhQKKdhxZr5U3CPk8I1wt8EbKOAVOgibp9hBeY
aaocYqWdX8tfqshuJh8okE094LCKwXPHcM6MYF946mH9Muc3BrlLMf8dD/AH
dMstuXwZg4cPRW1KtwNBnf/sK3sqBtoCnreDVIKBKWY9Zd8CBmeG6w0MfUJg
T5C0a8gsBnPGd8/PrYfB3H5C655JDFrZPoqVKUSAc1mHGzsFg08tq3zBns+h
S8Fc8dcwBg6qMcwn/kTCqbUBh4v9GOivD23NnY2GIrf3+4qJGEiWiCyVub4E
3ZWO09udGKzFUxbOrMXABQkumdMtGOBZvHvr3OPgzkWjvi+lGGRKHCYWCMTD
BSFMLLYIg1DDgq7Uzni4bHfQQysPg+sJE20eIgmwt+t4kmQWLT78enWiw4kg
3+ukxEDL336NacT9PAlGHkroRgdiUP3Qv4rxfDKE8rYQf3lj8OJzSflITAp0
qb0UVHKi3Q9F3oJYjXTYrDfYnDDEgPVOWZ7fUjoMstxdFtGh5Zuv0UeH9Azo
5pqT1lbHwGK8zviUTyY8aeUOYJeh7e/MX/1/vVkw+Gw4+js7Bt+L72kURWQD
44WGF+fqqXChUbxeeCwbhLXxCXx4KhT2rauny+fASumM9nQ+FRLW4+DFRA4s
xomQWJKpYAUtStYq7+BraYwK3okKO63SMjxLH+AyPff+5QNUODdOfyLgZgF0
2LDdELy4CJIVN5akPxSA6tkb9bLnFkEo9nPD2HoBaD8ivhU9swj7svLOZeoX
Qoenyugg2yJsGd4c2qEWgpwMi64maQGGCwtEkMwX+Jbu3B9guwAZLuYNKmUl
0Mzxi/+13zyc+vVtXakOD7qfFQJNI2ZBBzfQYXCzDtglAttjTk/Br5B8m0bH
Ovi2NSLddnAKEh8F7SgH1YHJx9hS1l1TMCMhISOSXwcC5PTkHMokhJf7x+9i
rIe3H6ZIMhmT0NIhbFZeUA+VCVeF2o5NwuU19xl+pkYwSWPrzzg2AZr6nLv+
fG0G/HbYgLUiBaIt9y6tf28GvzS/5C0xCpDc6Ua3h5shKPPfeuoJCtilUyuY
97RAs589bmMXrT/W6nrctAWy8YH/RDrHgPQqaPzSagsMMHvwmz0cA7u+37Wx
8q0wZaG4yk/zl5f3+v2ly9ohwobVR/vGEIQuXIDp9nZ4MsfG9wSGwNfr7b+U
8XbIPx09VSw8BE5xjs+Y2Dsg+tNGkecaGQyad7mNWHVAV8QcB1MKGdilpB9G
cv4Az1W/VKXJQYjeDtefceyEGUeFmeTQAWg24GUz5CSAdJSHGRdzPxySpHC5
HicAfYj2StpqH1iw5ZyMFyLAP/GV4nNTfbDRJi5HUiJA0sO3z17X9YG4zgUz
iwcE4F73f3/Itw9iLpl9cCkhgCTPMmn3Wi/cVkjSfGXaDZxn7+w3/dsDi/yH
fbuTicC7IDjacpIICt5Zj2uziXBKIhg7f5AIgT0S1oX5RKBGHrEs3k0ErjAt
veg6IlgN9Ct0znWD6oLn0StzRLBZprgmfe2Gl6WDpbWKPXD3Fteco343yOhm
YoWkHmidVDrXHEIAdzcxy5dcfRDKVrDVz9UFXCxYpSl/HwyFydw/vq8LyjO+
HDkh3gcww3uTRKs/tloVu/LUaZx+w8OsqxMC+bWh1aEPIjjjyE5unRD9w4p3
d30fOKiw7Pla8wNyz6QPPHPuh/KEIwU8th3QP8RhbNNCgq3w0chHv1vhygtr
kls3CfjkRvjFJlqhUrHCPGiIBNLdYk1bhFbITHhgn/mLRItz03Z7QSs8MiwN
HjwyAAzt7IF8Dq2wUW9WZmA7ALLzigl+P7/D8U+5x8+zDIJv5qxZ40wLWHhd
+MmkRwbOb+a6VpzNIGrV0ks1JUPY3ti8q7uaYcXgekOvJRlw/kPV2r+bIPS0
QxbOkwx7GGIrHvQ1QV5nihngyPDU+KGjckoTrAtudnpskGHHf2bJSrgJYrvx
5VPvh8BScyxERacR6sXOP6/fPQLzET9LBQrrwVzqOm5p/whwnwjdHMyuhzXZ
x5V8R0eg7lPnrdSkehBX+kT1lhyBndenQCmoHuK1+a7J3RgB6w6tIoppPdhY
Mx3PeT8Cm6w8T/t36oA5o/9zqO4oaFpplAjQ8lxob+VMtvYYmC8ra3ierIXl
Zkm7t/pjMF6ltX6fuxaqQ7PmM4zHoPjuwVAT9lowZQyjJt8Zg/hyAbarWzUQ
smP0O8ZlDNicV3njSDUw+efnjk/iGJxS/dd8K7YG3o4f5DaZGYO8RPXnpSw1
wFtpr7k7nAK3/Sd3Bo5WwyGejQ2nKAqQX9T8ecpSDawe4YWkWAoI2U7mn96o
gi2pnGN5tLohzbFW4e1AFQziRpaMiykQ/0PUQC6pChJeGKVnjlAgWF7p+4ej
VcBudv6v8rlxiB3KOBkmVAl0GwwfXMfHQXUocin9XgVYTAys+P0ch/Nh9wU2
jSugpr1QNWphHBSdS9gttCrAP8O8593qOLzoZt+8JFEBOxoVO0MsE6BLEdV/
somH7ZePTS+fnYA8oqambRIetgTHdh0PnoDPajERfuRy+GtUZd4oOAnNYVcM
ej3KgCx375hg8BSMtjOXc3WWgMsNjgMCQ9Nge3VQyyHlC/DET95k5J0BrHBv
cRepEPLNMt+tPJ0FM96AK/M2+XBDxU+7pWoOCoFoO1v5CYrMHQWUtuaB5ckV
dzPXD7DGcmkX3mQR9j+6qXrr0jswJijL8yZSQe/woeu6njnAoryvtjKNCueb
l+eqnHOgNntQ9w6OCnET2U5yj3JA+qmXZQrtncrz65KRuJUD7HxlkTyNVNC2
H7mjr5wD3x1kJw4v096xvKF1TbocUGMWiWHXw4Be/YfFnqhsWv29djTfCIMg
s4+RSyHZUEBuztG/gYGG0di1Md9sOFFg8y3yPgZE568fW1yyYcske5rFC4Pe
lN/yDdezoQzHp8KUi8Hou8IyK95scGLDmnM/0TwV/8c7nCsbhDyrjbVontAD
U18L9mdD/JW7j0IqaZ6bmxzAzJANrkvJcbtoHmLI/V5hYhoHomZ2J3AkDAS7
fewExnBAaVD8cHEEg1eSgztWAzgwTCRVB9A8yKxaTe93Gw6YGd7r8NM8+M7j
/m21Rhwge8+eGppnodoDrK+qceDRq3XPguZtbNSKx5NlOPi//1r4///a/wFW
E3H2
        "]]},
      Annotation[#, "Charting`Private`Tag$33257#1"]& ]}, {}}, {{}, {
     {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
      1.6], GeometricTransformationBox[InsetBox[
        FormBox[
         StyleBox[
          GraphicsBox[
           {EdgeForm[None], DiskBox[{0, 0}]}],
          StripOnInput->False,
          GraphicsBoxOptions->{DefaultBaseStyle->Directive[
            PointSize[0.012833333333333334`], 
            AbsoluteThickness[1.6], 
            RGBColor[1, 0, 0]]}],
         TraditionalForm], {0., 0.}, Automatic, Offset[6.272]], {{{
       0.9945218953682733, 0.03886990416599716}}, {{0.9510565162951535, 
       0.04235006897143931}}, {{0.8660254037844386, 0.05063291139240506}}, {{
       0.7431448254773942, 0.06753742310194143}}, {{0.5877852522924731, 
       0.10376363605465243`}}, {{0.4067366430758002, 0.1947090760181396}}, {{
       0.20791169081775934`, 0.480611696544933}}, {{0., 
       1.}}, {{-0.20791169081775934`, 
       0.480611696544933}}, {{-0.4067366430758002, 
       0.1947090760181396}}, {{-0.5877852522924731, 
       0.10376363605465243`}}, {{-0.7431448254773942, 
       0.06753742310194143}}, {{-0.8660254037844386, 
       0.05063291139240506}}, {{-0.9510565162951535, 
       0.04235006897143931}}, {{-0.9945218953682733, 
       0.03886990416599716}}}]}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}, {
      {RGBColor[1, 0, 0], PointSize[0.012833333333333334`], AbsoluteThickness[
       1.6]}, {}}}, {{}, {}}}, {{{}, {}, 
     TagBox[
      {GrayLevel[0], AbsoluteThickness[1.6], Opacity[1.], 
       LineBox[CompressedData["
1:eJwUW3k4FF4XVhFJkkrqh2RPkqUFpSNLpAghKZQkhJAk2UJSJEt22bcMQ4w1
3Oxj3/d9mWHGjKks2fvm+2ue97nnnvue95x755znmTll/lzvyW4GBgZTRgaG
/39KAhf53z8a4tIpH5w3CIfn64IB6xs0pPeIeXfU4UjIL5QWWVqhIcsLZ5FL
QTTIiGmZExdoKD5Kx3Mn8iu8mDZmmCTSkNyFyzyXxhMBF2+VMDRFQ9aed25P
MSbDpUN+w60DNCR9QEPQkiUVXFtCXRu6aYjTR+rBxGoqlL9L5PrZRkNBzV/M
VgfT4MpGuV5hLQ3tkk9zNXudAR44/K+cKhpqEpW4aSqXCVX2/cEZZTQ0gDOw
06dkgtLM7+aYPBp6sFTjYX/+G7z9+s8qPJuGtkdNXWSavkHN3QPMn9JpiPf9
tfGz+tmwh/O/dP8kGgqInFrM6csG1VYxFe84GtqYy7Pc0MLAO/+LU68j6fwl
zILNKjHQoKTq9SKUhjhu8l4VEsoBlk1dXrsgGvqZUTUd7JMDN4rMfli+pyGs
xYjP7GAOfHxud++hDw1Vdo7XWonmQuvpN3/vedD9b7PQTOxz4cBsQMQdV3q8
yHyVEZsL2gmRslovaEiKK+31V2IufDZK67puT0NWUgdannBjoZOz4LmSNQ09
5ZZZjlXBwqE2dEDBgq6/RFm5qxUW9N63YWTNaGjUky1d/j0Wwq+N3DhrTEP9
F3Z2uJKw0Ls5PydiQOfD9nBcoxALR4tX3/Hr0JArW7Do3mosGDowCp24SUP1
MZRQvyYsRIlz1hy+Tl9fupez1oqFwdmTDw9co6HOy+FTES1YOJF4dmfvFRoK
F/3kHlyPhfv3LsczXKKhCrufOYd/YCH+8A2FDWka2jlk89EgBwvjbYaDSxI0
lJWGTvjGYOFkgIULVZSG3EOP6Tf7YOGhstOROQEaUolmEjCix5O85VUwyUtD
f94/szbSxMJ08SedYW4aSjb0ZFsVw4KgY9xiz2F6vvPQojEjFizOfAtqY6ch
i4OUtbzRXEgnFIs37qOhxMAxHt6CXCAm1uF/MtKQ0d3H95r8csH6yCQjbmMR
HXL7be8ikAvZ7dSU3JVFJMlmnPWTmgOUgE2lzF+LSFv1dc1ocQ7Yb3N5xBIX
UVZAd+xhlRzIKxH678vUIvp7zajIYW8O/HKUKfs0uoi4X8o5COEx4ETUWvHu
XkQ3MzlN9qtj4FWHn93TqkV0S1Yqbf/FbAhacBlYLl1En24q6D8gfYMUZutr
PoWL6NwhdbmGuG/QClpHv2Ytoh+OGp8f72SBQP7Ryu6wReRStz+CoSoT2kMy
2K4+XUR2/vd98xXTYTYn2qXl0SI6zdDn6URLg3X8x0mjB4tor1uU2HhyGgjt
eo57obuIap+0+2+wpIGb46UH3y4voiqetze4p1Mg5NPphosXF1FJtUPtdmwK
pH/7T6pOahFZqPSf2W+QAp1T//ZMCC+iq0HktM7WZBDRw2OOcCwiP3RE7kZr
ElyxLz+awrqI/uTvDiz6nAS6H3O8zjEtooDDw5+T9JPAvSbkjuYGFTkvcL4X
bUiEsHGfyoFlKlLwtHpteC8RMjecRZ/QqKgyOILxJy0BumXubXrNUtHud8eC
X4okwLz2zScHJqionTModrDhK2zbKHbEDlHRi37Bp2k2X0Es9VRqUTsVXb6K
umIr4kERHT6g0kRFFmISWX028aA3wvSqs5aKrlqvlMXwxIPHYbLmQikVwU+1
B30f4yDs3CjOtZCKbGp/rxiqx0HWzXa+vVgqaijeP/eSOQ4qn/4MCM+iItHL
/bYXWmKh27fgD38qFe0I73BFh8XCfGLaA+xXKhpKu9yaaBIL2z8iGy5HU1Fx
02S7lkQscA4GSDWFUZFmoUF1yk4MiC27xRp+oiLvdz7ZCb0xcJXDjnH2PRX5
PxHdUMLGwB0JM3tHHypawAh8CgiMASsN3cEddyrCeqd/f2EbA54WKspBr6go
69DdMwy6MfDF+0LOcSe6PhGWl87Ix8C3eFGuTFu63reDQn4LxUBV6XHv80+p
yLBi3en2kRjo7d1Prn5ERa0i8klKzDHw1qKoCvOAimKujXiq7USD5LJpeMRd
KpouuHPk4Vo0jPiyWHnpURFX2chI4Eo0BBwuuGKtRUUSK2Y7tXR8IfX+oTsa
dP12WAz2rEfDtAwT8YoKFe3LLHum+i8aQmqw5SJX6fEe8PV9yxIDinpGnznk
qSg/kj0CR+dDntplsSFLRbOUbb0BwRiIdsTIzUpSkf7JaxIzF2JAbZfBgfbT
VDS/ba7dpRkDf0J2pkqE6PG3tGTGmsdAIn9WcfJJKpKy2CBd8oiBW/m6gYEn
qCiHWhiYHBsD67Bp9vIoFZ2WFW/pKI+BjI6082YcVPTgJ3d++RhdbzPtfTf2
0/cPDx4w3hMLDLS/YzJ76f7tbN+ln4kFY/ab7/duURDvKQtDPr9YYE5Yvv9r
lYKqlLeZFXCxgDubIDX8m4LYj+a8HyfGwkGt30O5cxR06p5tX/ydOKgYjcVG
TVOQC49xceqnOLCxVfV9O0ZBy73ZwQea46AuMErCoIeCvo5Htk/fiAfXZkWv
rSoK6n/0KeK1xVcQNp67QyyjoFipCCmv/K/QTQoR68RR0CSj8KXIna8gsW+2
KzWbguJ/Ka63JSXAlPpHoZuRFBQSIxZycSsRNOv7m2LsKIhZmuXwN59k+Kvv
neBrRUH4dJmtb8XJkDZ7+oXdYwpabV1Lub+QDP/2ePBcu0c/fzdjbblhCuCU
BZ/Pq1KQO+7X95bzqcD78/mRSzwUFPFV/fn9fenQynnBWeEYBRn9+e9gvGo6
uD3Z6FHkpKDIbYZiae906Gd9F67KQo/PMntBfS0dQgxiOfWWF5D/U0O1LwsZ
sGehjsOudQHdoai7Wg9nAenoiQOp7guIylzQPVOaA9FWE7YZrxZQALx7mb6S
A9d/pLV+c1pAZxpXMiRkciHl4bmgfKsFtFgX2OGenQsPMKr7q/QXEFPYxnBV
FBY6lJ7vG5JYQNzh+lSCVD54hp+3HhVdQBEumLbkW/kgQVzHTwgsoGV0r07Y
Kh8CAv0CiNwLiPd08lXvhHxQ7o9hXmZcQNWztjtGLN+h6Fkd08ExMhLct+uQ
Z8d3kOO3l84YIKOrC5qCPKTvUN57zESxm4z8Xt1Rf7S7AH4qPsM9ayQjbt8v
rSOyBdB8kPMx/jsZOS0MnHUMKwCtuh/BZjlk1BqhQhb5VgCdrk/KVzPISK5t
Du+GCqBvqvSQcDwZnXLgOLdOLoDJQrOfb/3p9iPiH+UUC8Hcah+F+y0ZCcQr
/2vQKQQCT+Gx/DdklBajdYT6uBDI7/Y+H3cgIy69hZDtgEJYuZvLe/k+GS2q
sW6fbC8EV7a7N7oN6Hg50Oj+eCFs/GR4aa1DRvxj/GLii4XAIK7fGqVGRpp5
O6XebDhg3dpwWz5HRtm6Ein7VHHwKT8tI1CcjH4dknw4qIMDjifa3QLCZLRZ
X5cgZ4KDo+3Jp3VPkBFlU1At1xkH0T43DeaPkJF4rmxEtRcOTlxa8fY6SEa3
AS4+/YgD/iSNwVxGOp/9cV8sE3CQov9nj9o/EtIuM3L7mYkD4X3x50bXSWin
pnZXTj4OxJ1o/qw0EvIXK4mTq8ZBrkhMQTKJhIIoV44M43EgNaI8LjdLQt8Q
bBzoxMEF1cgLT4dIyOywioTkGA5K1+DRTg8JVTw7efy/GRwo5JKCItpJ6PfS
nqK4eRxUPgovlWgiodsZnOb5VBwocSnO1taSUEHWvQ/3/uCgppl48H4VCXF4
VMRGreJAzSvk8p9SEvrDPsb2bAMHeFmFpx8KSehJEqNJ7zYONOdnwvixJLQc
xXawl6EI2uI/VZVkkdA/KfVcmz1FoKN7iaydSkJbHFb4aKYi6GaaOkr8Sj+v
xvavKXMRGJR/vOYRTUL6x98NVLIUwaD9ebsj4SR03WNzf/G+IjAWHI/GfKL7
F+S9qsFaBGMD7+uUA0go47wA70s6fhgk/WvIh67PyFFTRTp2w4X5t3mQUDT5
v0t59P0Ro0s81a4klHZVPKeZ7j+f0aAQ94KEkveXmITQz2+WKL6RZU9CUvtN
fjHsLYJZ/WOTcdYkpGDWzMbLWAQ77q4uny1ISOh5nhlhVxFwpw+x+ZqRkDPK
jTH+hwPZNoVUF2MS+q/h+nOPLRxor8TJ2xiQUHN5c7T2Og6seLc7THRIyPSJ
enf7Cg581EwtdW+S0Beyy8jObxzE26Et1ev0/EdhX03Q81EcwR8ud42EWPLb
7BxJOOiqfHta4goJDV+UCCiaxcFedlXDwzL0+lHZenx3mF5fF9Mpe8+S0DXS
eHNxLw4um+713RAlISvyz6SWdhw4YPF5U7wk5KNjwny8Bgcf+09f7+MmoTCV
JMM7P3CQtvNxFH+YhNr9ukKv43AwpK21L38fnd/fY9b303Gw7IJNTGWk1xMm
e9jzKw7YEw9ejPo3jzaaGWdMI3CgTOsy91yZR+P6ot5GfjjI+WxYcWtqHnF0
aNFumeGgoaRET2l0HjGFf2F2MMDB1AQ3SXZgHnGvi7+7exMHXFIjR/9rm0cK
yY3Vry7iwKvT7DmpdB75fKsTO8GKA71D1gL+IfMo6awctiezEGzlm0tfB84j
SU7LmfLYQnj/6MxtO/95hMlp+PjqUyFUfKe+0XefR81btakvnQpBWM+pT8Bq
HmUd8zdNUyiEtbA3AQjm0fWztICmugKY46kmX1eYR92SnOnquALoz9ir1X5+
HrnPjsl+Sy2AwvLQQ2On55HpWLEvj08B2E9nxm4cnkcxurdJFMUCmJHuzbkw
P4e8lRQl3Qu+Q3unRBcmZA5ZsX9Ieh6TDxXGTrKygXOILPveUTwgHzCzJRHl
7+bQ3o2Dtqsu9Pd5Tdm4yW0OFQS9VWzUzwcV/nszRIs51CN/zNPxUD6UPn+3
ckp+Dlmosykc5c+DzPXmu99k59BHnNVHEkMeRPpylEtJzqG1NB6RfdNYeBEV
5wmCc8jN826HbDoWJNF3ZpMDc0gobDHB9RwWUtnHj0dPEdE0n3tFuUEuhN0H
ZvlRIrqyy+ZS65Vc8MlKWh7qJ6LhG8YBGkK5YK78uOO/ViKqJrHf+LCcAwIu
836JxUTU8HJW1iw2BzhrbzgpfSeirbuWrzb9cmA3B8ZsCkNEGWpFoocccmD6
m52CYDIRvTlYkZSnngMpY39oGYFEtGYVpTy3hYEwcf0xdX8iCuUcfOqxgAGf
V0XN895EdKtw9XfMMAbMD7mmi7sQ0dGyO9K25RjQMx0Ka3EgokHlmh4hDAaU
MQrets+IKGHPwyLreAycUts2xj4kohtaLXVfvTHAEWaqcfs+nV9ESNbnFxhg
mEAXfhkQka6M1aXDTzEw6erDIXOTiDRXn7zsvY2BzvqZ7W41IlKXfZdyRA0D
PznVFl4oEVGwQKLLggIGknKYG4ouEJHyz63qR2IYCFm3KjSUouvDnOC5hx8D
3tebk/6KE9FjGtVWiRsDDuFngqOFiUhIofXr8UMYeDgZ9Eaen4g+xSydCGHF
gM7ZRavhE0RUcvIvMY0RA0putw3fHCWiuX+wbfQvG6Qa81V4OIiII36/Xu5G
Npw8wildyUpExhPW3Smr2cD+6AWfKRMRuc1etlBYyoad3N79//4REPu8DNH5
VzYsblxYT9wgoMF3HefuLGbDuHoUUWmFgAqUHh3qp2RD+5e1nikaAWno79Xf
WciGqql71T5kAuLrxJS10jFW8gdWkEBA5M0D29fo9glveOLrJgjIek9Zlyk1
G4LxHh+eDBNQJnnfrCAtGzyPTrjs7SOgN/hzQ6G/s8HeXMkis4OAYjTO3v62
nA2mecm6Gs0ExO80PGqzlg3aW7uBVEdAEzjioYGtbLh6w0LiIyIgLYURl8Vd
GJCMrD9+ppyAEtx4P31nxgDvjAhzK46AHnD9aednxwCbVMCybR4B3f/tlXfl
KAa23ElTB7Lp/kN7Ypl4MUBp0uzAphGQZMX0yhthDIxy5VTcTiSgvfdD+JIk
MdD6+ED2rxgCSonOf+Ykh4GKfPuo0C8E9EPlwNklZQzEa0o79XwgIL9bgWK7
jDEgtF+zudOPgI6FFZt/f4IBTIu5QLsXAV07Wqa44YiB8ltfuvEvCahQfNPC
MIhenwdyxRscCMhN5u2CdSwGmtvqfWqfEZAei/N7rm8Y+nv7V6bqEQGphm8q
XKbP2w8PcgT+eEBA2vyGjwsHMTDfITZTepeAvjIU1tSRMLCmYxxeqEVAWP6U
g4UHcoDrTuVSljw9vzklIhUPciDhcP/NjPMEJHvo9fMtxxwQ7l1MTT1HQBUL
L6/2vc+B8wb8+gnCBPQz+0jp7aIc0Lvriws/REAbdlLM5sdyYfhYPFsoGwG9
xbeHPpTOhUeDOItgZgKKdrYvXrmZCw73iEc+bM+iClyT+7RPLny+f+OlJ2kW
qRUH1/z4mwtcPOZtb2Zn0Wr8sT1PuLCQMOom/HpiFmlPvatNvoAFrElO34u+
WfRiLnuy5yUWWs0OXrSunkXF6fwmbNtYYLHoW70TM4tSV48/Fj2cDyFCi9q6
X2aRZGyN7F+RfDhG2Jup/XkWlYaKMLFdzgcRS7m7N97NIq+EbzJ7LPJBzSqu
5KrjLHJ4VcnCW5YPPraPXE/fmEXDfkdZUp9/h/Fg0/hylVk0L/qxK/bDd1D4
fv/nzauzyCVlvmMk9Tv8WTFgsZedRUfM/dcvDn2HR16a0QW8s+jXxOhTJfUC
UAqTLb38ZwY1l2ntY5MqhJ0SpjWtuBnEucrx/t05er8yvJtnImIG0WzJvLd1
i6B46x84hMwg7oAj5ZovisD+2sb7sHczCPdgGNdbQu9nmhePDTyfQfncjw15
VIqhYmzg0kPVGcTaOba74lEJvNmT/cqJOo00b1W/4y8rg0feHUL9WVPoR9Si
5MLTajAcSrw6ljqFeDZtKVyvq0FTxsFoNmEKSW/m9F78WA3nZzmC/nyZQrd0
lw7o5VQDi4beEvvbKfSw3ShxkVYN2IN91er3ppAK35lD/s41sJkwbFq2bwp1
LpZc93SshV9/s11/Mk4hzeWkjELvWpjVeRPW+G8SBUsew/R9roX23TwNfSuT
6OPe3Inx3FpIsTQ5+2dqEkX2X+5lINWChuTkpnj5JGp6zW3y4H4dRFQSouNt
JpHAjknI03P1YOoa63npySQy2pSn/rpSD2Kyty26zSYR+9h+EVPNeijPKj3H
YjCJGuzFy2cs6mE8PLDhBUyi4Z65iOjoehC1kV7WPDKJ8hXFzSM36uGXEHGI
wD6J/E8HFentbYCyiVjkvW8SrRX+w/w+1ABR2SxcdtsTyKW13HRArAEmvbEq
qt0TKKRkReyaQQM4nd1I+PNmAkkwKDcHZTUAY+pjGzHHCcQjcLjUsaABIrnb
LphaTqCEmw2a8hUN8GNPUmuTzgSyevDuo0MH3X5YbSNJeAINml4ley/T7f1D
DW53jCNu0sScnHwjiG1t8L+rG0ftXS5tt681QrmjBaW8bByxkka/aN5ohIkH
l3xF0saRkHGfMsWIbi87lr/tOo5whmUFMy50+wkxthyBcXRL8ODb2txGcL+y
KzeWexyNZsw4VxQ2gmLMkNYH9nEUfMeENbGsEX7qffxsuTGGbKvzy07VN0Jj
/QLnqe4xZOlxpPjUSCMECNQVHsSPIVb2S4q2k42g6RWvv1M5hron+H8kExqh
XU4rcuTbGJK5QTNopzVCLwZ7PPLtGGrWe6MjuwcPkSzvy/1ejaGPaqudw8x4
uPvE7P4LuzHUWrrRbc+Gh2E+jq869+j+Yt5ImRzFw1SoI/9+6THkgtu4fU8E
D4uvzot6T46iHPacjF1qeMjvZcPb948ibqHSrkANPDhKE6xMWkcRB8shib23
8LBMjshWKB1FVn2Ozq16eNg0+Xt25fMowo9Fp2WY4aGivKN95t0owv4IjSgz
x4Pnsazn3W9GUV0ZoafqCR52dxl9z3s6ilitJvWTn+GBRbX8vA2MolWxsDvt
LnhoTgrrM7owiojbDJFvXuMhaNvGRf3MKOo9G/+c1x0PB0v+KxU6Noou7kQd
V3qLBy5xj8sTiyNI1dekbl8gHgQ5lJT1E0aQ1Rst7eV4PARLTbnZhI+gn7UM
nIqJePr77VPoHTCCRpeTzT2T8dAaWieU6zSClgNlzSjpeHA+fIN5r/oI+hh4
dewtFg8NXHqtJbRh5Kb1nrhShQepS0uM7bPDSN98+8DWTzzE3v2iODs0jOZd
CqZ31eLheXQ/lqNuGHFLf3Bkb8QD94n7odbRw+gW/sQduQ482PBaGP53bRh9
bDtxTWESD72KTJ+lLw4j7QUFDedpPFw1zWhUPzOM2ORu5RbN4uFw4ryc89Fh
ZP+iKUiHRNeX3+6/NtIQ2j3OGLn8Bw8i19j1Z8aHUMWFas9XK3gIeZQXtN4z
hLy33euZ1vDwJOX3tnDVEHJvX3K6vo0HdiGXSc+wIZRt9dPfkqUJXFWPHY94
P4Rsaq5G6u1vgmmLUl2M+xAqrtHg0mZvgpL0jZoByyGEV3V8+fJIEzwU9UqX
ujyErijnHftyqgkKxANspmcGkZHuRKQfNMHJrodlZ4YGkfbN6d1hKk0Q9Eqe
5WX7ILIKu6lWrN4ET+rI6cxlg4i1t11FTacJuE21piSCB9Fs7/DMR/Mm8GMU
kXLxHUQP3iqiB0+b4Hf2jidyHUREx1MmN2yboPlvHo/e40Hkhq73Obs0gUco
p9EruUHkLzVlLhDUBORLCxk/zw6iOo/8nYHQJjAcr13ZJziIuDwDVXKimkDq
zMvw+AODKIO5tbQglR5v3UB79fQAitXQGVKtbAKNtXg1tk8D6CeD03PSahMU
Jbz8YuAzgJyLVPVqd5pAQE17JuHVAOLCyvvV7m2GzdB/3tLmA8jlUTVO7lgz
YM+YVxheGkD6+le19yg0w3/dCvuTJAaQ0IBTAY9KMwS4HjYmnRpANqdNTzy8
1QyP6uv+vmGjr8uSGi3NmuGImahs8lQ/Ern1Rdn7fTO8ZWLwIff3I6Mv0g6r
oc2wiBnslG3tR/lfMhVi45uhce2DfUNxP+JjfWptW9AMbmGUbwuB/Ygn3/+p
+mQzTNR/5794sR+xEG4MXVJrAfc766pyEv1IOdagvlO3BbinlKwVBOj7NYSU
okxbQGer4/tV9n5E4XeSxrq2wE+ZRWV1Yh/i+NZy6EBeCyQlij8xiuxDfJqW
Xn1CrXDlrNMH46A+FKD8NsjufCsMlpflPvDpQ9OnrUXlVFuBo19j9aF9HzKq
NZCxtGiFt2xP31tf70Mslhdzzma2gvnr1G9uq72I5Z5wTNbFNtjZu9DmTulF
12UfGgTeaIPYLzJ/PKd7kfITR8b0B23QnVet4Nvei/Q5mRU9fdtAhTjREpjR
iyw/SO6097SB0B2exXjDXjQ8ejuCz6sdCBIRMqi4Bz37Lw67wdAJ4p7pi/OY
HtT3qrr23MFOsO8oyuZM7kEVHPjDH3g7Yc2xX8AysAdlaQY0ZSl0AmvJscMH
HvagJT0PyUSXTpC6FrtkvK8H1e/F/HNb6YSXYdl5fjvdaMe9/OEOUxeUz5Q/
wy51owDBmpuZXF2g4j8ys2u8G0Vvi+1Yy3WBYQtPb1ZBN/q7ZGyJ8+iCNwaJ
uNX7dPtUppqnB7vhZ0aeA79uN+qOLKu/IdANjGtIQvN6N2oqqN6jd6EbgmMn
075KdSP2P7aqffe7IWniVIQqUzfCrq1pZmZ3Q4N12sswbBfKfxc54na7Bw55
frsouasL6TWKF+CLesG/WAzP8bcT5Qx0WI209MLGYua9JUonSt18sMA+3QvT
ZhnupYOd6B1fOf8gex8UXEutUc7vRLQPhY2/rftAj+mr9l2zTlTxvkX2vEg/
iB4/1cbwpgNtqemDbfEASFykSaQvtSKR4Jxn5edGwFLSy6tsrBVdhiX92tsj
kChysLutsRVVRwiYkZ+PAOexc69W41rRfOHhldD8Efj7175aXbUVNY7WjyrJ
jkJN2aIB+UsLol36cJOkOAZGioue5y42I5kvl2xTLSfAV5XaWebaiN5XGrOu
bkzBUX1syL1jNSjw2IOnN9OJUCASWK62pwY9UXhc+baYCDrrT2eladXISonz
ZGEjET4mnJJjbaxGfQtvbfaSibAz/2Wi3KUaOWC4CTaSc0D0eCPJ0/8TaWIx
7/YVz4H97dDoytMIhSb8ZPyM5iFrWZMSsQuhU03WCyMd8zAdw6hkP1SFKEJh
bCcn58Fw9tUc34cqNCD/xPI9AwngtelF7/lKdP3GjTujSiTgSD3Tq5xZgfze
lPK5IBIUrNYfxAuVoy6xs9lPCskQnJ1Sc3W9DL05odz3vJYMz0y9Xha1laG2
aN4j1j1kEGqQG0lxKUNmgsPcYktkiIzAZHjgS1El95ZRj/QCuJ0PVZS1LUH2
z5YPlWIWQNnpwbOEQhz6KH8HBYdQgE9Enu9oAA5RTIpviiZQYHPoaFfgAxy6
buuV9A1DAdy1jouvmXDIsNSf+VkDBUQOKe/Wv1uIqtzVI0c3KMCaLxqzb/M7
kuPM679nRoX0EcInkh8WjTQpidYeXYRhd67Er7JYpG8mWxt3chHY+dTzdadz
UUzdn3MmpxfB5WFWdxnkouUqc+7YK4twnWBz7OMGBqV0Pn29/mgR5hdpSeIO
31DOHsYJ+axF4AnhLxjn+4bsmLKMvuQvgo60bm1YWxYqDRo511+6CKUvCggb
4lloi88s5WDTInxYcxZvIWSgyovOBnfnF0F890bhs/tpaIjslicmQAPTVPH6
k6xp6JUvYxxNjAZhqvf7e0pTUZGdMzb+HA02/CvWLnOlIrcEO5ZvV2jQst9L
ka0rGX04wVQSYEADhtw87Z9eyWiH8RK7130anNeeNHOWTEYixQOs9x/RID7k
ms9oYBLabjlZWmtLg0FxluB7nxLQzbn8m4Y+NFCys86vJn9FUybWF1zf0yAr
r7n7tMZXFHSLrOcVRIPX5z8d29gdj9wfvFu5EkmDE1c5k2NdY1DZ5IaHcxYN
fLxf1O7uj0b1b39nVeXQYKGml2AjG43aG+Qyqfk0qFCPEr+8GIka5rL+rJfS
wESXBzfy+AvStdnMaG6kAWVpRf5lXTiKQns0T7fQwD2yA7ELh6Okl6d97drp
/Ed8Wq4RQ5HdlndURi8NJDwf6I5cD0WCR6Qfxw3Q4Af/xQHnzBBke3Mz/eUw
DXAt3fW3u4KRa0ie8/AEDXo5/Lt8wj4h3mDRJYtpGiwZyI8V3QlCvkV9bb2z
NJCZTFz+r/8Dusnc5GVBooGu8B2G21EBSMkRe9NvgQYONnvZfIzeo7VMVpMP
VBrkLdsKzg/7oattKldUftOgXZ7/3H/xvkj+Qx/D2h8aLHr2KGib+CDhZvW1
0GUaHKjzv/6W7y1aHYnjO7RKAyNO5zQr5IVkZo86vPxLg6j44R/t1zxQE/UK
qWqNBv0i13rO17mhtaVhf+o6DY58zyTHXndF7pQttV2bNLhzmX33rqaXaLDz
j+AGHYfVOx9/evMFas+4zTe4RYPO2yNSbW0OiC2AJBO3TQM+nG+IZZIt2soa
MlPZoUEcXgfJpD9FjabuqZ10zD3Gu7jzzRy5CzluKf+jwZmZPvlIk/vILjHW
OpaORdgcoiSf6qKiFD5yPx1794r1/leriGakLT3W6Fji4JgN01kh6PfhFfxH
x/2/fRI51jSgbUd2hEzHv26L/nnuZAhflEfTf9Dxm/Drzc/OPwQSb5CvIx3v
eL6OSzz3BHaFbTmz/d+/TY5tj7gNaHWcfhNE57fHcEKRWeQ56MrcifpNj6eW
/8P7rxxOUAmbzZfp2C9SwOFSsDOkJs4de0aP/zpbhVHX/ldAOHbK24uuD7OP
wbVnH15Dy7slptcbNMD/XTzNxOwO5f1UzF26vh/tAjgT/TzhHMO08390/TN+
h3jwvPWG5tXZx3X0fD19Km5ZvOMD1eNl1fVLNNAQeGXsq+gHMfcY5fjo+RYb
q9XWcX8H/+K1Z+7/ogFZz+QSef09eNkXkP0pNGg+kH2mRO4DBBUa6ziTaYDB
r570e/UR5jj4d67P08D2aggL30oQHJ3cVI+foYHW+ugWWTYYBu4U0fimaCCJ
O/27xOkzHEJdq/7jNFj7MkURWwuB+FY9ReYhGpSxvumreRkOB4n1Dk2tNEg6
e7QnT+ALaJ3rCpBsooG/Tl5nfMcXEHL8EuJcTwP9yJkWl9ORIPNaJy61kp6f
U7dqxMeiocKnnhBMv78DqkR07GMM8JMVqrcyaVD11LuS8VIsRGrmUm+m0iAw
F1c6HhIHg569DsEx9PqQ580LU02A1YDvxel+NGB7UJLj9TsBhiL5sSe96PfN
UzfbNiERaM7C9m9e0+DhdI2eoHsSuDRc/zhoR49P7K/Wdl8yKF7R5TpDf8+a
Cs1Uv39IA9nR1cJrfDS4Ui9RKzqZBlLmc7mTXDTI71+/lnAxHVzSm5PMD9Ig
cj0cAmfSgY+ds+kgAw0sAK/wRDETXhTcjjs0uQj/mqWkj//+BvFclw/iYxfB
aXQ777N6Nhjhb79PDl0EIrVZcm9CNhzzaX5jFrAIrRxPJJY0MWAkmaX24eUi
xBrGiLal5cD02RazKO1FuDC9i++tUR5kVkb+KN6kgmT53d9S3/Lg5Ogr5sTf
VBAJy62bXM+Dc4wvfjnMUWF/cs6FJK180J+TYOjupsKWjtHov8V8iJD0oQVm
UWEsP+80ki4A/r+7vFZvUyHRwbROsQQHbGrfRGw/UyD+wT5sxAQOUtTYXCt9
KBCtgYuiMhdBZ6TZ+uZLCoTwsz77alQEwWbAeu0BBbw7iw5trxeBJUfL9KIo
BR5JHTCruFICleGTUys/FkDw1491hZoy4GhUfRg/QAb3NhW/ZWoZiLbEzkU0
kqE/u4Ude7wcdt/V1fUsIcPHJyOCpxzLIcFL05cnigy/hje0mPl/QPHRdLtw
AzJU1cun9HhUwJalsOyXDhIYx5Vo2sojON0xj5Usm4cbqUNt2kY1wPblj/gb
OyKoaXHuWSlqBJIkB2tVwxQEP9r3e72pEerVvugXF03B4EuGiZ2xRlDrFbuA
S5sCm4TFcpa9eOA4Xnq52oduT2t24jHEw4bWMaog0O1DfadVVvEwxqI5Jlo2
CTb9y9VhF5vhVFpwyh56f1ZIXsiL0myGFO9lhR61Cdjamf4ab9oML9t5EoqE
JyBYtOt1hn8zyF8USikkjEOha650eX8zfK6Li421HIetE5bJUy4tkBFXNqFs
OwafzQa8pUpa4eXKyz8PAkbAn3IFiK2tYPlu0tLeagQ8X6dsx023Qp7+i6tf
NUbAPtzOjZm9DXZ3sv26s28EtBv3OI9btMG2r7H4maBhYD8n9TSIsx0M3r94
mhs1BME7AVrzdh3w6tuDGpXGAfAPXNyf4NsBla9dj5thBsCTW7/5TkwH9JSN
2X79PAD2MifVUV0HPCzJD/E1GgDRQmP+lgOdQDqGfbi+0A/Mc7+Op1h0whVq
bkT68X5o1OY9oMPZBQzhtIb8kF44IjnF5cTTBWrFnxsobr3w8EA6/xeRLvgQ
ih++8YTen7dInB9U6IL3vUdsnir0gsSNK/cfPqbPFY/rRKSJPRCicv+bA64L
rhQseDKp9sCowMnCsKouqDm0J8BQqgdO756pwOG74D+b1+xNPD1Q89Omc220
C76O9cRxrnbDiqLbX2+mbiB/PtB/iT4/GMvFqIUadkOv5wB7M3c3ZBwzuV34
sBuYORxHZJi7YWmV/16fTTcMJ51hq1zpgqCiLNvjXt3w3Paiz4WeLkAypeHJ
md2QJf1sp+hzFwidHZj6vtYN5wRaPlLZu4B66qhnd2wPTO+6oj8p3Alyb5Id
q9PoPAcOFaVxd4JP79kn+dgesLGr+/J+fydwvb9+K7imB8afsOzRnemAq5RX
JzTJPSA+Svv6zbMDPhcPF1fL9wLxq13Z+5/tMHzwaXa+Si8Ycoccb/3WDkLW
S18TtXpBw7g4STO8HUr/Y3vn8agXRET476xYtsOkt+Id+Q+9wPr419NNjnaQ
vplEyx/shWDHwtno523gniYxk0ifi7oP7y54Qp8/G7ZL+4MpvaD5ueaANX0+
Nc7vqrRl6IM09NORU6gNfI7uCRIT64O9LyTa50daoWfiiVjSqz6QjCxqsDJo
BT75P/99ftsHriElSl5qrWAV5nnQM7AP5N6q7MZfaIUt1ajV+4l9wDP+5jyN
qxWEsvF1xxr7YLXXiD9opAVeOp959JmrH7Z4hMMr7FqAi5VWYXiKjsXjwirM
WqA0sYCbT6Ifkpx4dv2hz/dbzfKdOdf6wVl602/zYgv4nNKAZtt+qNkpcvHY
0wJCJfvjQl/1Q+vG32u8f5uh/lbHqpFPP/A9jOCjkZuBxdUQOxfVD73Fu3VF
e5ohuN2Cl6m2HzR8HpmMpTeDlIXY69a2ftApD3fbFdcMXesLveGD/WChrCmk
H9IMR4RfBAks9kOKv4GO75tmiH3zdhOOD0AQU92drTvNkCGWMOT2fACI3VuV
v9maQaPq0QVltwG4mKOS0MjYDKQ7wqH73g2AefaPpLqtJpDwytGIjh2A6OsS
gjLUJijoKy8pqh+AX7dcMW0dTaD/zPOwR+cAjNdqne3AN8EKg/Jz1ZEBmN+X
ubpU3QRyEk0iPb8GoLhtyjYd1wRVPgMRv/4bhOgM+fOxcU3wkDvud4nIIMyX
0ZSyIppgF9ZUy0t6EHSIhKy+z02gOkRgZFcfBG2WFyFBfk3QfG75xRmnQbC3
2T/+4HkTDIxy6FniB0HElsrMrdkEmoFPBp27ByG4VK9J9HoTVMiXm/qODsKV
9ve0m8pNkBT5+FnSr0HQuFvT3KHQBFY6xX7D3EMwOlFrVi7RBMM7rPtJAkPA
wLdv4dHpJriVaxb6V2IIfBas40+INIE0676EI9eGIKfgd3flySbYqL1fom09
BJRXnoo/DzeBrVO+osmLIShQzcokcDTBOD9T3TOPIWBx2LfNzd4ENR7YzoDQ
IfAXZ2rPZGmCwEu7yDVlQ3Blo6zTeAcPPJgMnkuswxDE4SBQuYCH4HubKWpH
hkFZ9INJKAkPDCw6p/X5hmF0a/S03RweZi3WLzjKDEOpkcCc/Awecvlu3cYY
D4NQxtvjvCN4uBr2++1JzDDEspq541rw8PD1lTnmWyOgcOoa6UMuHsQt8H2L
hiPgHW6XZoPBw5K2fl3foxHgcV6+c/0bHvyFbJNTX43AWrt7FyUNDzkdcfch
dQQmm+h9XxweXMpFNUWwI+C8K/5tZAweIL1Q7kDZCOCkHQ0UovDQ7dZydKR9
BIg/Xv6xDMPDuvBmh8vGCPx88TpGIwAPNRzvq0yYRoHDpKEm/R0egjY5c1U5
RqEgR+fyXx888HeJf+QUGQWHM0+DXnrg4bq7sSpWl27fZbnT64iHsO6yUkLW
KDRYqrxrvYuHB5Vqma2FozBtutx/TR8PwlldEYVVo9DOy5iVroOHUg+S09ve
UfA023G+dAMP42LHz/IyjAF7OC/GTR4PWYfT/2PcPwZSoZVTxhfw4Lgjxbpw
dAzaBWC/qDQeGHs15krPjMH8fk5GHzG6fl6uyfp3x8Dm3AkLBS48LFszhlw2
H4PzTd93OR3CQ6V+iKeA3RjsZB3TCGfDg4541v1fPmPQL9ZUFbqbrl/fwNEg
LN3+8Wm+V9RGqD1z6WMt0zhg3z3D5FY1guk5/dTfB8chq6gVtZU2wpqMY8XJ
E+Mw+V2de6CgESQUMItvJMdBQ5RlKS2jEb5onLxz/u44XIzjPHT6UyNYPmHm
Sc8aB/H1i/8KDRqBwVroQk/BOPB8kJK6q90IsbbXtHdVjkPvPWH/qeuN0P7i
jbdJ1ziYDzIqYC41wiUf2uzRjXHQvy8owsHdCCyJA7n+NyfAcz7GmL23AVJS
lhtwBhOg4FJff7ylAa5kHJqcNpuAaaUXp/fXNIBD7k1OcJ4A936B74l5DTD0
A7n8jZ8AHxahweAPDZA9mAlWixNw/YCx3Vu5BhDZVzGfpjEJfn8DzZt96uFP
o6RNitYkxFYU0Dpd6qHKP3khUW8SGAecmJBNPRgyvl+MfTAJ46KGx/X06uHd
P93lEIdJuOgtrWDIXw+zK3P/3KMnQf501eZCSR2kTB8+ZjA/CaXa8y6i/bXA
W/FMjSlgCubNtDrvD1bDkeMbG/afpoA7j9i6Ul8NbC4B+YNhU7D7wlc978Jq
2DqX/l/O1ylwyJS9ZxVcDcOp47/1CqeA3f3syWCVaogM1E1IGp8CfKbcm7eY
n8B+/9LfyxemgVmlPuopqQoYNnZ/c5qeBstSt+tze8vh4czQktfcNOSr1PMJ
zpTBz9b8q58o0zBfJbbsiMrAO9G0N3N1GlzZkjxuupbBP9Xyf6OsM7DFpGS0
e6EUdj47GqrLzkB8jtMif18JbAlP7uHxmwFpmeIc3uoieHCwRPv0hxkgxhM5
eZKLoGLtU8zF4BngGTFkU31bBO4tl8/pRs+ANgebxWHlIth0iDL2z5kBbskx
gZp6HGz80P7+q3cGZuvro371FcJf3UrTeuFZ4DvJvjPKWwA+Q3ljV8/MQlBy
eobg7gJgf5TyoFRqFipC2vHfid9BxOG9MebyLNwS4juamf8d9IN1DUN1Z2H7
mW0ws8Z3yG8maJl4zIIyp+tNkXf5cEVvsLXfZxbIBkHks/b50DjUfFMnYBbW
3E0e+N/Nh/H5vBsq4bPwU1+r6MyZfGBjdlM7/Y3un8DzHyYnD6xV2BVXemfp
/QQPG/trLCw3M1TYD9P5Dh1STDXAgrfeksL8xCz4fZown5PBQvSjQbkR8iz8
9ak84EXLhQavlPPVDAQQuFbxb9MhF/SYIwoV9hLgou8jgeN3cmEs+L0Mbj8B
bCZKwkkXcmH5q61UJhcBfvu+4czYzgGBiosSnyQIgH9e2hcVkQNYldOYvTIE
iFYo+H3PIwfkW/4T975EgN3n1fZ3P8kBnWEGsRfKBGiV5NtPkMuB0UdLGRR1
+nl/sDe8BXPgKYkgbKlFAE+P+xe72XPAc61Z8J4RAWLvnXSPncMAq3dlSrcJ
AbJ2X13m6MNABHP+qVuPCWDE6WJ0tRYDOVwRJ8GeADJbqj2VyRiQS3ifUPqC
AGH1we0HwzBQJ+zGK/OaAJMFGl7cvhgYOW/2n7AfAd5yCL7WeooBy/ULB999
IEDXubvphcYY+FXJtocQTABBmQz1Jm0MuPnMrKp+IUAaX6hwsAoGGNXLyWkx
BHD6E729Sx4DwftDxxkTCaDuLv9e6BwGjnc+7bZII4Cfy6zTijAG0r5cbaj7
RoCAyS1rZ14MSN47Wi6UR4CtiG725KMYKOWl5PrhCEAWTdhyY8eA8nRN8mwZ
AXDL7o27WDDQmhEToYoIUOZjt3V+NwbuPnP4kFZHAEWeW6eObGfD9Dl1D8Zm
AtTkdOAS1rLBdpnX0aKDANyz/pK9y9mwWrpsUddLAFpM307h72zw9mgxEhom
wC3XP8GKtGxgVU655TdBABfSrJIDNRsi9r5Wmp0lgOQt3+c3KdnA33L7vCqZ
AOftHO3bF7Ih+7OIWBqNAN0aJz+u0fF5/e3/GFcIsHwIz9xEt6/i7j1osUHn
7+MoDIvZoDGWvafuHwHKnReOG//Khp7kt38FmYjAVttx5PhSNphaGi34shKB
2dnn6pvVbJgXPzcxc5AIhlIypV4b2eBEY+pROUoEpt1DccL/smGrcLQh9QQR
8h8qbj1jxIC/a2H5Hn4iSAYszeqyYuCQ4kfsY2EicPOaGfdzYCB+16OUWnEi
1PEW3Vs7hgGRhkuRglJEsOffGSw+iYH8j+wffS8QIWuCu5FDDAMKtwkeMwpE
sOOq3M8kjYHbg2FPUtWIsE+w63CpKgaG463v7blJhAC3/b7OtzFg8UhJ67EO
EbxkbJMa6PXjSqaeF7xPBD/jyzSRFxjYnVcn5vuQCEnzfeEy3hj49CKOZ+YJ
EaTer7e2f8JAypYGY6oDEV784IoszMZAC3taj4A/EUyZL5xjXsCAYY9bo08g
EZ4r9Os7bGFgMkr3x3QIEQYc0z+Z0u/LMv+/lJQ4InyuubBaK5sDfLLGTgLf
iWCmmXvvnV8OZP2VsvQpJoLKMYZ/pJgckKlgNp7+QYQ0n0vvGvJyQF2t6FpK
AxHMpcv4JkdzwOEuxyGBUSK4nuT5aH4lF8yzzsUKTxHB8hMFUunvg8G6tuBp
IhFynFO84VkuKMR9uiD1iwjqKwwqfnG5wDTBanyVcQ6MD8dLHWLAQqwlU6qx
xBwwxDx3rxzAQt3LjfPh7nOwInr4mJFxPhQ3HK+KfDsHzl+XBHes8+HbMXn1
WP85OL8YFYF7nQ/BZa/uJYfMgWWkvr92TD7c2172yEubA1Xv6mfHh/Jh0Y/W
0NI6Bxd8eTJ5jb/D8S+zRoy881C5Pfn9qnEByPme5HIWmIeKU+e/SzwrAEMn
454Z0XngO/+hmdG9AMJ0OrVqZeZh64ltkN3XAmA7UKHsozEPD11LudfHC2Dn
XbgEg/M8aOwyzi+5Xwg8L9tJz1/PAz7igUucTSFcttiXOeE5D+wCvn56rwvB
VdnrFPowD1yc8pPnIwvh944Nl2fiPLx9NPk9ta0QZl8pM2w1z0PeE+PNv9I4
2P3Uo/JZ5zww3G0LfaaIA37DUreRvnngEPvVEa+BA5Pzkivlk/Nw5VKfiKwp
DtwFrQvEifOwekCmNMIKB7Gcac9jF+YhKSF1Ot8JB/2/jpNfr85Did6eI2z+
OFie0M8kbc5D9Lwe5sZnHBzu+GxxbxcJgl9lbylE40Anl2lCjo0Ef5bNOc9/
w4F9vFJ81iESSOb+fiv/HQdBgW/ucR8jgcP4kd2LpThosv7d8/cUCZTO9BWY
NeJg3kgi9KkoCSb3vsk50Y6DvRpPtQckSIDk1jS8e3EgdCllv7oMCcIbyC6f
h3GgLDKGL75EAoFA0dPqkzh4eJTbX0SRBJHMPYYZBBx4Md5RiVQmQfxlt00M
GQcV0/jKl1okkPX4jM1YwsFw1543BD0SXCTvaon5i4O1n1flDIxIUBWf7XZu
Ewdc+a9X6kxIkGUa2Wm5g4MLibiC849JkDfcMXt1VxH9+5n2PM2KBNa+rr3F
e4rAyUP87BF7ElS7TFW3MRVBiO0Tsu8LEhwzSev1YS4C7P2kzCVXEvReq5Dv
YymC8wOCG3s9SbCYhDEr2VcEP/Qyb53wJcHNZpuzyqxFoNwunng2gAQ4KSvG
//+eHX8D+1vpEwnarRw1b9Px7XppVf0wEmyl1Vl30vf3KxVFPo0iQb2IQNdf
uv8HFXIkt3i6nuljS3X082cuVVwOTibBSQztxqW9RUAvueDkDBKoj+/V0WYs
gl+StZM4DAmKr8e837e7CF5lX5fF55PgvPuKi9U/HOwIN78bKSKBO/7ON5st
HPglaw0ulpPA/NIc9uA6Dvbzdonv/knX/8siz90VHIRF63scrSdBkGH9J6Xf
OOA+Mtgh1kwCsvNyUgsFBwmf7wtc6aBjMZXNv3M4EN4/4Xy7lwStUY3X66dx
kONv3mg+RAKixV8Z2TEcyOwmHncZp+shXup8bQAHZR7Wth9m6PEPPW+gdeJA
aZ1SFT9PArbtI/OKzTjQ+rX0uPYPCU6EXSgp/YGD3mevivv/kuDW+K5dhEIc
GM9tsJC36PXJ5HEnH4MDq/FdWA5mMtjyGQQKxuLAp/Xgmsl/ZOB8IjL7mn4f
WDTCNR35yWA8qDEf9xQHn2u5vvoJk8Fw2M7vwQMcxJfzKmPOkYHbLqqmTBUH
JVlngtZUyODUU+Iwx4GDq0J54/tvkKHrobTgrj04qEuUkT6pTQbmHzVh5cuF
0B0p369mRAbhQ2fLDwwUAsVPnT/clgwHfV1PXY0qhBcMLU4ZjmRQ2FzzM/Uv
hI032vVlLmTQrLB8c+plITC/MLCZ9CaD3Hr2oVDdQhB49Bh3NoIMGV942ET2
FsJdRS8NfCUZqsxfvS81LYDRFeLBmhoyZJ8u56zUKABzrPbAj0YyBLBmv3KS
KQDbk3yWeV1ksJHqe7TESO9fd1X6RhLIUPB1W9Mq4ztg6zeqLA4sgBzhoD3r
aD4wartc2GWyAPWFcbvYL2Lh497xrY1HC+A8xqbTtZELB5Fa3bLlAhDtxdM1
US6ckDp6Z95hAYKzbkfbqefCOU6cQ4ffAnBLCe0M6+eAcf9vzNecBXiWan/x
+9Ns+G5qJ6CwtQA/DzjvKvyYAeZFv3iTdlGghbdnbEIjA46wOR/fy0yBP+ax
K4bMGeBS5sbRfYgC5bWz+xXfpYPCkYAda1EKdL4gWca9S4OappThGD0KKE38
eKSTkAIv+EX6/92lQKACy5CIRQoIvfrW9cSEAm0qFFKeeAp9vs/HS1tR4MHd
2JjlsmS46VVZ3ORBAZxIROSe6STY6lMqOOdLAb+yprPi2CTAStTlRgRQwPba
0SG5N0nAMdyc9iicAvlntLfFuJOg9/xg2HoWBYrORr88UZEA7wLvB5thKZDT
yPXJ2CIBLk6Pf6gvpED3SmxZD3sCRH8meIdW0fmbbzr8//+gN4hW7n9rKcB4
L1nrKM9X2LhCeWXSRIH58enx3x3xYEL+Y3+6lwIJIdzvT0I8sF9zsfk8RIEX
ebH/Ajbi4GfU+pOVcQqoDnovBJfGgYAag0k1iQKa9rPXE6/EQU+cr5EojQJX
xl9N/9wdB35/9up/WqYAZxxVNbclFoiJB27e+0ePp+xGCadlLESthlxHjFSQ
uEFxSr4UCxpaR5SFWakwPzKqLsIWC+upUYqBB6ngvjHuhZmJgeyNE/K/j1Dh
EN+K6OWqGLivm3D+7gkqMC9jD0zExQBb1impypNUWJq0qUhyj4GqnbQzgsJU
WE/46PbhYQw8NxAT/SBOBae22Osp6jHAn4MRoJ2jAkXTrndJOga6d5/jM7hA
BS5ZkY/+J2PA917B8R8KVLisf9nDnCMGzudfOHpKic5n/cnzz4wxQNhbxvFe
jQp/hvxbjmxFQ6TJFTaqJhUmWdxT1lejQR2HmO/oUCGwfD5TdSUa1lhV9pQZ
UGFf4aLkFn2dc5H9S40xFbKNpU4UbEbDma5hoVYzKmz+UybY7okBVVxGUZ8F
ne8g4acUewyYRjldn7CmgsZy9sddPDHwyu3qwLw9FYi6L5qnJWIgxITV6s8L
KpT2TOYNKNH1Uupf23SlQvK/yoMzd2OgVjDlA5MnFcJs7j9idoqB0b32Jw76
UqELN/3j+ucYWCXJY7gDqHB8M6EnJS8GDrYxXflfC+f+FWMegHHCUckcRXaJ
hmoL5ZRu69rToYupLEWnmBSpFhGKoRiHJKKpoVa3qYbpMs30vmlXV+myhVZR
7ZCUW7Vsbd73qwy1uuz2g7/geT7P8/vHSMQgTLp2eJEqDctutzWbX2cQtX1I
p/TfNGxMytxtlzL5l/W9xDDjdPAFB1gHyWS++kC8y7Z0JKyfOscrj0Gtei3r
9Vs65NzHUr6SgUOdxPxifzrqNNJWBd9moNJ1Ff1plAF1o9X2U5UM5qxSLS3J
ygCncKz3fA2D/PGijU7dGTBLbDx+tYHBNfedu3tMJdjlvedG1hMGneKQC2vu
SnB8tYWZXDW5v7W/eZ9WJkQGI2XFHQwK3QT/ifmZqO1OfNHQwyBBVTFcMjML
P4TVLBr4wkBe2/U5NDobDp5XKfUog9Ki1uecvmz42vo4TExhIf7L9oT1Mimu
fGUDdGexmGf0yFiYLMXHWEPZj1wWqQOCMdMzN1GVLVwe48ri9NcmD/E5Gdqj
eZUiDxY+wqCtPvdlIMH67imeLJQpRll92jkwsqAOKfgsBgc5N+1Tc3Cp/GVR
6xEW/MT3K0crcuHVts5+cSoLc61ed3sdObTXzaqrymRh8eWMzeIAOepyOt39
ZCz2zbs8WFgsh9WpyL0ZNIuR0eweN+8CcLhl8Qvus5hoHI16e12BB5djv6t4
xCKq//PE8TcKCIe8b/m2sqgX7HubuUKJD/fVZSldLJ7GO1M6VUr8cci6V3+I
he/OXFH5g0Kca58aVjLMIido/U9e0yisdmwb2THOIvljpi93AwVSIL2g1iAY
su0TMxEU8uYe5SRrEkg2CMPj8in4C5FmwyGYrXFJq/o5Bf2/OSaquQR+WiEf
dkyn0bztNR2+gKAnSNtu/koaMZXUGj0uwaMm+6hGTxrrTYQNxSYE0VKq0Sqc
xieRx1bPFQQChatEL4GGYtig86MlQbm7usMyl0bg3oEgsR1B/Oit+s1lNBY2
VRLLdQTSx/qjBg002myvRLU4ErCPJUdDmmhczto544gLQf5wxZ7eZhqOmsvF
HA+CmATK3aORxsixkYW0J0Ge9B/G9x6Noq6HuVt8CJIEEatrlDRCnFOsGD8C
i7LFPJskGoZFIXfjAwn406sDf46g0f69vYvFfgKzZ2vfzfKgIYqe0dZ0eJLf
M9a5ZRENpw9P+aERBPOGzs7wf0dhzDvnvXYkgZWyVBmYS+FOTcQxxVmC5oAl
dy/6UQhdvmmMF0NwPUimOKxNwShJ71J/HMH8tuSX1XQhOse6deMSCZ4Ezr9h
wCvEtZBiybJfCGx/3fXF+IUSm1vPmTWmT/aPfMKc8FeiTMbdMDNvktdmu9M+
VwXCZpOHeUqCE03ePEtFAUxPVnu5FBNoWI8b5k0pwKu3olfvSgniPl3YtMtN
jmS33fsvVhHsaNdr6YrNh8cdi08mvxMYi0wPlpfkYZrhmLDhIYHz0Y5Np9tz
ET6YnjRNRZBSGv5MNSDDCv5BQ1kHgX+82CbizS10N6wp2PiaILUlUlOn/ia2
pXZUn+8nOHXMLtnxdhY0NeS8pYRgSWDwOFdXgprQk09r1QQT5Zt59VvSIHjm
ErDnKwFP36lqTvANfPO5oPybz+V/zLOfOw==
        "]]},
      Annotation[#, "Charting`Private`Tag$33341#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,
  ImageSize->Large,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{-1, 1}, {0., 0.9999994076421085}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.764563697032604*^9, 3.76456370334973*^9}},
 CellLabel->"Out[94]=",ExpressionUUID->"042c00a4-5212-4230-8faf-9b5b609b952f"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1344, 796},
WindowMargins->{{141, Automatic}, {Automatic, 0}},
TaggingRules->{
 "SlideShow" -> True, "SlideshowSettings" -> {"Toolbar" -> True}},
FrontEndVersion->"12.0 for Mac OS X x86 (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "SlideShowHeader"->{
  Cell[580, 22, 4399, 98, 75, "SlideShowNavigationBar",ExpressionUUID->"c443756b-bc88-4e3b-a9dd-ebc17ad287f3",
   CellTags->"SlideShowHeader"],
  Cell[5697, 149, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"65af3d05-da9c-4201-bc7e-fc9881fad03b",
   CellTags->"SlideShowHeader"],
  Cell[64767, 1276, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"ced35eef-05a9-4c5b-beeb-93fbec39c1d7",
   CellTags->"SlideShowHeader"],
  Cell[88964, 1743, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"1f165dea-40d2-49ed-a901-a80aa8d15f1f",
   CellTags->"SlideShowHeader"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"SlideShowHeader", 126113, 2432}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 4399, 98, 75, "SlideShowNavigationBar",ExpressionUUID->"c443756b-bc88-4e3b-a9dd-ebc17ad287f3",
 CellTags->"SlideShowHeader"],
Cell[4982, 122, 240, 5, 30, "Input",ExpressionUUID->"825030e4-751a-4878-b036-d7b8185725f5"],
Cell[CellGroupData[{
Cell[5247, 131, 259, 7, 160, "Title",ExpressionUUID->"db09f6dc-cd51-4f60-b989-d0bd041c5d56"],
Cell[5509, 140, 139, 3, 67, "Section",ExpressionUUID->"12ad6237-6c23-4157-9523-c690d4971f36"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[5697, 149, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"65af3d05-da9c-4201-bc7e-fc9881fad03b",
 CellTags->"SlideShowHeader"],
Cell[CellGroupData[{
Cell[5844, 155, 157, 3, 67, "Section",ExpressionUUID->"9f6092a6-ba30-46a9-ae7a-28eaf8d86d3b"],
Cell[6004, 160, 292, 6, 35, "Text",ExpressionUUID->"a6fb594d-3a5b-45a9-9ab7-4c59081a4f4d"],
Cell[CellGroupData[{
Cell[6321, 170, 982, 27, 73, "Input",ExpressionUUID->"e2944af6-5975-43fe-82c5-734786dd8e16"],
Cell[7306, 199, 12107, 223, 367, "Output",ExpressionUUID->"7ad71d8b-f25f-4735-89df-6af91bbc5bd5"]
}, Open  ]],
Cell[19428, 425, 417, 11, 30, "Input",ExpressionUUID->"c750fc2a-6f2f-496a-b487-242d7e074886"],
Cell[CellGroupData[{
Cell[19870, 440, 924, 24, 52, "Input",ExpressionUUID->"a2082dcf-a84f-4748-8969-810cc7ad8e0b"],
Cell[20797, 466, 21174, 374, 367, "Output",ExpressionUUID->"a185c92e-4c1f-4a66-95ba-d51bd36ac4c9"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42008, 845, 541, 14, 30, "Input",ExpressionUUID->"a3adec70-cba8-4948-99ec-b09fadfceb22"],
Cell[42552, 861, 1011, 24, 56, "Output",ExpressionUUID->"90025396-220b-4f1d-95ec-175998ea7f79"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43600, 890, 793, 20, 30, "Input",ExpressionUUID->"2f94ef69-4336-4f58-9a12-5021d868d4fd"],
Cell[44396, 912, 20310, 357, 360, "Output",ExpressionUUID->"66c96ffe-6c76-4900-8672-080a253bf775"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[64767, 1276, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"ced35eef-05a9-4c5b-beeb-93fbec39c1d7",
 CellTags->"SlideShowHeader"],
Cell[CellGroupData[{
Cell[64914, 1282, 218, 4, 67, "Section",ExpressionUUID->"0827e0f6-3063-41ef-99a7-bf91ae511d9d"],
Cell[65135, 1288, 404, 8, 58, "Text",ExpressionUUID->"bf486cab-c44e-4920-9d6c-04f7f907337f"],
Cell[CellGroupData[{
Cell[65564, 1300, 1802, 54, 136, "Input",ExpressionUUID->"34c62287-5f88-4042-a6ae-06271f35e9c8"],
Cell[67369, 1356, 21534, 380, 374, "Output",ExpressionUUID->"e1d77504-d688-43ce-9669-7d652054b834"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[88964, 1743, 122, 2, 74, "SlideShowNavigationBar",ExpressionUUID->"1f165dea-40d2-49ed-a901-a80aa8d15f1f",
 CellTags->"SlideShowHeader"],
Cell[CellGroupData[{
Cell[89111, 1749, 273, 4, 67, "Section",ExpressionUUID->"1f0952eb-56fc-424e-b45b-c771242e2c61"],
Cell[89387, 1755, 288, 6, 35, "Text",ExpressionUUID->"09eecaae-e115-4c10-991c-404b57942b51"],
Cell[CellGroupData[{
Cell[89700, 1765, 2018, 64, 157, "Input",ExpressionUUID->"25942195-732a-4c7a-83bb-aafce7829194"],
Cell[91721, 1831, 33974, 584, 374, "Output",ExpressionUUID->"042c00a4-5212-4230-8faf-9b5b609b952f"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)