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

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

(* CreatedBy='Mathematica 13.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     46202,        920]
NotebookOptionsPosition[     44928,        893]
NotebookOutlinePosition[     45326,        909]
CellTagsIndexPosition[     45283,        906]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "Input", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"kin", " ", "=", " ", "20."}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"mass", "=", " ", "936."}], ";"}]}]], "Input",
 CellChangeTimes->{{3.8698997367340517`*^9, 3.869899843251645*^9}, {
   3.869900162194075*^9, 3.86990022514108*^9}, {3.869900263296324*^9, 
   3.869900278820773*^9}, {3.869900716799706*^9, 3.869900721797647*^9}, 
   3.869902565634445*^9, {3.869902624866561*^9, 3.869902626057645*^9}, {
   3.869904698619644*^9, 3.8699047019298077`*^9}, 3.870142822697193*^9, {
   3.870144764448131*^9, 3.87014476799769*^9}, {3.870144888379427*^9, 
   3.870144970243932*^9}, {3.870145792817855*^9, 3.8701457939369917`*^9}, {
   3.8713910972634153`*^9, 3.871391103619568*^9}, {3.8714518138981028`*^9, 
   3.871451832796793*^9}, 3.8714523038600388`*^9, {3.87145238207703*^9, 
   3.871452385670923*^9}, 3.8720846338697147`*^9, {3.8735416074490557`*^9, 
   3.873541607808901*^9}},
 CellLabel->"In[7]:=",ExpressionUUID->"4a74809f-4d6f-4bf5-84df-1c81e845ccce"],

Cell[BoxData[
 RowBox[{
  RowBox[{"sps", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "0", ",", "0"}], "]"}]}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "0", ",", "0"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "0", ",", "1"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "1", ",", "0"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Dashing", "[", "None", "]"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"0", ",", "1", ",", "1"}], "]"}]}], "}"}]}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{
  3.871701793434409*^9, {3.872065693342557*^9, 3.872065697733029*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"4abccf8a-d47b-424b-acb1-36e9461ec53e"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kIn", "[", "e_", "]"}], " ", ":=", " ", 
  FractionBox[
   SqrtBox[
    RowBox[{"2", "*", "mass", "*", "e"}]], "197.327"]}]], "Input",
 CellChangeTimes->{
  3.870144988515087*^9, {3.870145025281687*^9, 3.870145047838258*^9}, {
   3.871451965683363*^9, 3.87145196928242*^9}, {3.872063208516612*^9, 
   3.8720632133894377`*^9}, {3.87354165429843*^9, 3.873541682028607*^9}, {
   3.8735428967851543`*^9, 3.873542921291761*^9}, {3.873729554724765*^9, 
   3.873729555224043*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"ddda9570-52bb-4f7f-8d3d-86bdfd8586a4"],

Cell[BoxData[
 RowBox[{
  RowBox[{"kOut", "[", 
   RowBox[{"e_", ",", "U0_"}], "]"}], " ", ":=", " ", 
  FractionBox[
   SqrtBox[
    RowBox[{"2", " ", "mass", " ", 
     RowBox[{"(", 
      RowBox[{"e", "-", " ", "U0"}], ")"}]}]], "197.327"]}]], "Input",
 CellChangeTimes->{{3.870145798453801*^9, 3.8701458265877657`*^9}, {
  3.871451836727564*^9, 3.871451840179946*^9}, {3.8714519734346457`*^9, 
  3.8714519739887877`*^9}, {3.872063215936698*^9, 3.8720632175080347`*^9}, {
  3.873541677338694*^9, 3.8735416793910713`*^9}, {3.873542899115011*^9, 
  3.8735429242501287`*^9}},
 CellLabel->"In[10]:=",ExpressionUUID->"86a69de3-cbb6-4ee5-8cee-9244eebdf060"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Delta]", "[", 
   RowBox[{"e_", ",", "r_", ",", "u0_"}], "]"}], ":=", " ", 
  RowBox[{"ArcTan", "[", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"kIn", "[", "e", "]"}]}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}]}], " ", "+", 
     RowBox[{
      RowBox[{"kOut", "[", 
       RowBox[{"e", ",", "u0"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}], " "}]}], 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"kOut", "[", 
        RowBox[{"e", ",", "u0"}], "]"}]}], " ", 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}], 
      RowBox[{"Cos", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}]}], " ", "-", " ", 
     RowBox[{
      RowBox[{"kIn", "[", "e", "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kIn", "[", "e", "]"}], "r"}], "]"}], 
      RowBox[{"Sin", "[", 
       RowBox[{
        RowBox[{"kOut", "[", 
         RowBox[{"e", ",", "u0"}], "]"}], "r"}], "]"}]}]}]], "]"}]}]], "Input",\

 CellChangeTimes->{
  3.8735413196128597`*^9, {3.87354136431686*^9, 3.873541446430656*^9}, {
   3.8735415717184877`*^9, 3.873541576826765*^9}, {3.8735424454621572`*^9, 
   3.873542574461391*^9}, {3.873729533939823*^9, 3.8737297084349613`*^9}, {
   3.873729740827138*^9, 3.873729822160321*^9}, 
   3.873730350881184*^9},ExpressionUUID->"f9b93666-a5b0-48b1-94d6-\
6d1a85496d1b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "2"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "5"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "6.5"}]}], "]"}], ",", 
     RowBox[{"\[Delta]", "[", 
      RowBox[{"e", ",", "3", ",", 
       RowBox[{"-", "6.8"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"e", ",", "0", ",", "5"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", "sps"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.873729834304262*^9, 3.873729962816473*^9}, {
  3.87372999361537*^9, 3.873730015311598*^9}, {3.8737300586286993`*^9, 
  3.8737300653434553`*^9}},
 CellLabel->"In[15]:=",ExpressionUUID->"aa0fae4a-3cd8-46d8-ac13-32fe6f27924f"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV13k4lF8bB3Bbxcw8Y5iQkiV+UiRZKpJzi8guSyJLRpYoRUhosyQp2dci
S0mL7HsyZI8wspR6lRZEsrSI8p75a67PNcvzPPe5z32+I8U4ZeHKxcHBsZ+T
g4P9evsCY/WOY8Fag7MyL+SOrEPPWZLGNpK2yC7c5v67g5FopaU+TFLyOFI/
Mtb3OTgWUeMrikhDvujRffpiumIiut1wQ2KV5Dk0Y3zo2nejFKQwcyxmxf0i
ijORPlHNl47eUjP/LQ+GorZfYrrxYbfQib3+J39JRKHGQ553n5y5g0r7R/MW
dK6jilqvnmqfbPTHy3hk1v0memPVJNTsnIN0u+L++zkYj17fKeBMkchDrCWP
xO8SaajFV44+ppqPAlrFlWaK0xFlo+ND5kQ+Wp/Q3zmtcwsZePq0q6XdR4yt
wP3VPQtFn395O2GiAAU+IlfPDeYgC+EPRsdsHyHD+M+r5+7kIsvTVhseTj1C
YoFMq9njeUhMQdN+LOQxatA5+31m6S4yG4iP6k8sRGvefNg8LVGA7sUZ/7W0
LEJDDU/9p8YLkA2/Iak5vgg9uJfa9LX4AbrqJGyv3FuETH1NnSZ1HqFvVI+i
bUbFKJmvJvmL+xM0yRyWLlEtQYK32rhmB0tQVJLJ2siPpWiJWlFwWrMUKZKl
awtEytDHy7nm3++Uoj2Cv39oGJahCvcLmTPHy9DOz6feVT8uQ3YqqhrflsrR
WCmzLsG7HGV3ZPl8lahGp6H04JXXFWghtXzr/65XI4mv8mcH+SqRnnvnWN9i
Nco5JJJYubsSTXH/sq5h1aCEnrNm25Mq0e69phpRkXVocreTq7ZBFep5ssy1
eeYZEkmotOzKqkYXOxQMdzg1oeWopqYNgnXofnGl7MOkJlQaHWbNvasO9aXu
4/7vRROim+mOPz5Sh2TcD9etU3+OzPctdprl1aEO7nBFDnoz+k66YvM/5adI
eO8bwZctLUhCL0j8w/56VPgk6s0JxQ5kOtcpvV6mAU0KTl43dO5Akn0bOnyV
G9DmAEMtucQO5ESKEayBBnRHk5w9ttiBqLrCekIODSi+/bqbXXMn4jKKU+RO
akABYzdn9e270A3VnyofOJkoi9wZJrC/B2UkcMeV9jCRO+pNm/JlIamwg0YK
Bk2II+dnzZtoFkoK03u6yqoJpfKIjXTksdDmmdnsV45NqK3dTbxggIXucT5f
Y+vXhOSslnNcNfrRyYSz675nNqFxj82P3nG+QuLBU3on55qQZ/z5+p7YAfTe
XCOGlPAceX/aMlZWOIwu3HTXjnnajEYdh3mZ7cMo6oaC1Y7WZmQ5dFWx6+Mw
agnI6X7R04x2d34592n9a7SO7Lp5eqwZcRffpQlHvkYRC0kZ/XwtKDVESuus
0xvEyXm8T9+qBTHpomnqtLdIo7UvjedjCyJJhW/+3+5RpOR3UbfueyuKoXVu
zT46hv6Yt5MqHrcjxhurZ9wvPqM31ZO648adyKc5nHFq8yTqEOvIjS98gYpP
UyzXGUyj/TY+11b1dqHsYK9JJZvviPBVWpn/0o1Cs5YPLkzNIvLBbdNBL18i
DtYo1AXNI6liqkekSQ+S9XMLtNr6A9nrm1AiinrQxx8ph8hPfyJTrdHHHty9
yPCT5/l4n9/ons7GDjGrXmT6MVZeXPwPsn8/+FskoxdNNjEutT1eQvpZF6yK
B3qRosDARh2nv0hFKFtomNaHSOiH88j8PySZUd2Uod2H9Kftfm2y5oBcbhXH
5dN9aKjxGKedAicU+1yoDEzvQ+VXnztK/uCEDb8UZgSe9iG3Vf6j3d1ckFJy
dFvwuz7k4htcTknnhvtCmjNSf/rQdnf5bmNnHuBnVEdtE2Shl4IKZH/1VfCu
W3tRZAsL1RaGZihwr4YMvx/k6j0sFMy51fzF29XQu2gksMmEhUgjQQe8y9bA
SuzG+dtHWGjn2QLiZuUaiHwfZvMS21lwqO9JzRo4fSEzncOehaoOqDnMNqyB
qEShpwxst/IZH7+uNaAKhN9/DizEjDmWEfR5DXxIkivJc2ShQG3Tb1fW8QL3
eFrhNWcW+nRPKjErhBfa27LSf7qx0FT4+0XuS7zAwbA/Je3OQnOMbCePMF5g
Pr8uYY7NISElv+MaL9hu9RMswF6fItnYmMILlj7ROrYeLGQeKfH9Uwkv+Ixo
jZQcZ6E6943GChO8oGHKo2x0goUS5UR5qq35wLbsS5mgDwvZO8rmjx3mg02r
xaxVsWUSVQyp9nygwFk1aI1dtmIS68Lgg1ZR/fBU7P6BUDHqKT5INs8JFsP7
kB4xpepylQ/OzOlPip9hofj3z1yJWj54XZv5j/BnITuRLt7d9XzgJbplSR57
k8nrhwwmH6RbiLwxwC6pWpitbOWDHCWt/RHYfTe3XGD084HI8vWmP9gCWgnJ
ldN8sCrmvfdoAAvFpru1OkuS4Pr5c+XJgSx0NiD0TIg0CVrN1t4vwna0yJRI
kSWBhaDopQ5sBdLA2RcKJKgzCRn9i91+br/cLnUSBLufcHE5x0JctjJRFAsS
OClzBW4JYiF/kQ+GlWHYuXyUnGBcn/m/P3uvkKBUWLO6Clv3pWjuVBQJPM+V
Gr/EFow8uCQVS4KSkWjZZezCnw0Pr98iwUmLzHtWISw0/iqLwignQVK/2DuO
87geiQ7dlC8kUBwuOKx3gYXESkkbvCZJoK2iFnYE+3+9le7t0yR4pHgy9jT2
MX5BzisLJJj6xW2Zju19rUVlhZMMHM9Kzk1hh17anjYrRobFo3usoi+yUMEJ
zmMDFmRQOpgRUHqJhU5EFxapHiLDLd+g1c3Yig+O/I23JQPZlzdoALvkS3my
2VEy2GSaUxax6xie7W0nybBbyvud5mUW6jnMUqy9SoYLw1XoGfZv3buLWfVk
GFN/q/AklIVCarPCnjPJsCpa6GAdNodyOjHxnAzbQi84tGOvkbgppdxJBubf
P4pj2MKLgQbPB8lQ+9egXCiMhVQeG6eNfydD/81UqwDsCml9aWKBDLsiq/9d
xtZI13684xcZnlaciLuBvS9yJzPoLxmc32efzsM2d5acoJAo0DI/I9yL7S20
sHuHNAV8gm15ZcPx/on+1mgtS4Hfut5aStj+XBPGQVsoYPTPz0EDO2Tm7dGm
7RRI/nbqkCl2dHvrVWtNCvyQZqr6Y98/nzF0zpoCGd9/5NdjKywkMTIPU8DD
Y+pUK3aRZ+xU4xEKDFUpS/VgV9pEcFAYFAiRaEHvsVt2nJLL9KaA+f4HJ7ki
WGjs077AxkgKzGr5xQH2vqzunK3RFJDlaY45gH3nsF1XfAwFprmGfcyxHTt9
NrkkUWBnpfuEE/brojsvuHMp8OnpZOZ5bHVPhV9e9ygwukmyMxw7VbpKqr+A
AgdSLd5GYx9KfhmQV0QBtTOdpWnYfcH/JPfXU8Dt9qstZdg71K4ZPWZSwNWw
7koNduw3oQChZgqkHhhub8A2dd7W+amTAjmP+ri6sDv07P2vDOP61EfBR+xG
gZr21gX8vO93zhJX8LncuX9h+28K3AlilQliXwrvFU9dosAxeREnEWytX+Nn
PLgIOFH+xF0Su2ZERJyPRsCTDeMGO7BFk3MP+NAJaFa6pKOGHWi2/cywMAEz
n39Kq2PvbNRrK9hIQFpPRhpgl+T7+xrKE7BkFX7YFFvAmfN2iSIBl/f2RR3E
Pr3+Rut6ZQLGmGo5VtiKN/LEvu4m4NPs4at22Df0lPQtNQmo0JY95IA9tVLr
U4sI+OMvQz6K/cCX1RKtR4BIXfdWV2w+BafZeQMCqo5eTXLH9vg0ucHehID9
/DFfjmPLHubyUbDCv5eSqu2NHSEQk5FoQwDn/AOD09gfO0Rblu0I4Doio+6L
nbN3x4YuZwJ+7XHv9sduNCWFLR0jwFmDZnsW+4PT2OQWDwLWNTh0B2Jz+dRZ
HPYiIEjbWjEIWzo0qeaKNwH0HcuBwdi6Cd6byn0IGFmxKQzBPpanf23Mj4AD
G727z2OHl0vOCQQSwP1LZ/gCdl7Loi0EE6C5+K77IvbzwT6m9wUCQmJ0n1xi
3+/4wy23LxNwVSIo6DI2z5/wuM5wAlj0CKVQbBmy4+JiJAHfejx72dYV2+Us
F00AUavkEMa+/jZa+6EYAkJVR1hsh2tNKEXEEbA301ctHPuuWWNqaSIBV2wW
L7PdcjSD40MKATz1/tVsf/bx86BlENC7bfoN26vDTHq0Mgk4yuk6wbZsouzu
k9kExMR9eM+23t2VrIw8AmrVvJrZdqsYWtORT0CqDk8i21dai0/9fkCADn+p
Cdv3hq4NyhYSYPg1ZI59fy0TLsi6mIBHcu6hbH/+o5kfVkbAqPDZv+znW00R
5i+pJOA2d4kL25s3zgSM1hCwbc/GcnZ9Dii2vaPWE6DLWz/Lrp8HytbbyySA
oyxNlO2r5kGFXs/xekeUK7Drfd/ZUji9Fa/fc2F59nq0+SpcaOsgYKikcS17
vcbDVn3+2UWAUFL9BHt95e5Vllv2EzC4q9b6HPt6lbEbQwcJMAhomWL3i0fb
8Yii1wTsydt5KgC7YHKDNfGegLPVdjvPYHcsLdTt+Yj7S1A02Ad7ktIt4/kF
94PWwYensLduv7TQMk1A9v90ur3Y/X/mY8KlPwS8qmgSd2F/P/zpUuFfAhjr
Pzxj74/JpGSXtxxUCOnKNXbElq86oKKxhgrMpv+22GI/Wn7UN0+nQsWqsm8m
2C+IKxqbRKhwZtX5aUP2fhN3yjFfTwWn04JD+tgK2gK+jyWpUPq74rQ29uMI
fwH3bVSoTvEvUmF/XtMhNEKJCgO+mZJK7M8v6M7nqlChPjrsvAL7fhlrX42q
U6HxRtJfGXZ9UWmqnR4V7GUcRYTY9fyZznfOgAr3Db5N09jrVxgalGJMhfSK
o8UUdj+IWdj3W1DBxWLlHzd7fyx+lzB1ooJVaX75LJ53H4uGYk8wqNAE3fNT
2NIeDZzRrlT42rYsPo6dM3BzrNWLCpxR4ZrvsLNKFfO1A6lw9ykzsR07/eQJ
RbU4KqjvW11/iz3PZayyLBOpML1zpDsZW3RkD803hQpj0Txdsex5bkieK7xN
hV9Mxzj2fE/a/KBc7gEVyDw1Gl7sef3+i6ZYExUkWilr1bAjrV2MuH/g5335
maMWnz85ZiP2rr+pMHQ3SKAEu87A2rt1Cdcv6BmlgH0+7tWPi+bih7GXvk3J
2Pay8oN0Gj8YJogr+2Ar/5pjyMjzg5t2oOgm7HepoUH7nflhoPvALj983i7G
L0XnH+MHKdXh+x7Ya2/43ebz4Ietr16scsA2uOzW0OXND/ez5y/vxy47brj6
UAg/vJscDGOf/9c0BOPdUvjBjL+mpBjnA7WR7ILILn7geR/x4TXOGzRFj99P
eviBEi5f14U9dVFRf4iF31/wCW/AzpOu/Sj3mh8mrOZG72Kv9eqX6PjMD8mx
uztOYf8J4TZazUED7qp4Zw7sYfJO5gFRGhzd969hLc47Lapp4yPraTA+NxOw
CrvMYZnfR4wGWTx8Yj9xnop90uiYJkED0x6myiD2AUuz5Yn/aBDkyUxLxa5K
99gVrUyD78s3ddZhp27JeNxlRANJyuRWMs5nC/N6d8NMaKCpQ7H9jfOfWf3c
LXUzGmxYt3LmE/ZqS8Prdy1oUDzhefoZtn/Ioud5WxqsLI+H+WJbvLSR2+ZO
g9WXx0Vf4TxJ8V+bez2UBrs+zctcx3nzYsONVKNqGig+ON7RiPPt5NuGxMUa
GvA8dWE+xD60NB+bX0eDtjEiP5Gdh3faRXE30KBFP1zLDXvooWxQbQsNAjYE
ZPJiK6c8OyLfT4Oqw2ayhmfx/Ds5K06eocFno2y/BpzHTTZY3+2QEQC50eW/
QTj/a3u3HhGIEoDkKDeHMi/c72IMmRc/BEBhm92Cpws+H5ryk02PCMKoYtfU
sB2u57bGrX31gnCO5pBUdRDP53me3vatdOBsthQn67MQo3hOMuomHeRd8166
auL8Gmq5bBhHh7IKoUFVbH/LskFKAh1u2Sj2cWNf/eEfE5tMh/FgRm42/j9X
qL64lHKbDk6dSiUjGrg/mSsD9x7S4b+MneHm6vj6fZQbz1vpeD330bbvxPsx
96THlXY6PJDtj19WY6Fcv26dA5102N7ly9eBXSEc+6ezmw4q8R/7jmGP2NE9
WK/oIBA2O5CmiufHmKjOhzE6TF7yvMepwkJKZUHieZ/o4L9G36tbGefXiDeL
rl/o0HeLIZOBbbz5dtHEJB2W9VzOqGL7ekmKz87SQa1oNNxtB86jmpcXS+bp
oNNBnlLGjiA+9Pv9oMPPsBj9FSXcX09yo3//poMtf/6bFOycSzzuNX/ocG1p
39pj2A8Puu4LWaYDo26fthJ22aaWjVr/6KDkUMBY3s5C9fOyiysrdDg8fi6w
Dfv/5P+Gdw==
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#1"]& ], 
    TagBox[
     {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV1mc4lm0YB3AjGmh4WiQrvSRpkKb+t0QZLVmpVDLKSkplNFWiSEJKIiGU
qFAkEdmK7IyMZ/MMUUjlvfp0H78P13Uc531c5/88VRyOWziJiYiIrBEVEfn3
jTvvILnS0X+T7dM6xSfMPShtUDa3Ud6LLeK10nK9HzBRVhiorHwM83zGLyw4
XYXpEblZ01q98cD2iJLRjDrEFYUqSSj7YvnYYFX14S/QEjiGTbhcwOJYb96S
Q43onP7w7++Wy4hdEWUqqtoMtZSper/Gr4Lbu87zak4L3PV9PEaUguEw+UXg
LLU2vGrsTho2vAm6wkyu3qmv+OVm3jHocguh32csqkpvR4JoR/CI0W3Iro85
41vdgS21txf/bIlAuGKdh3pbJ9gxW4uHj0ViqeP9xWrNXbjp+Gf/0HgUWF+t
9u4p/YaG8aORQqV72Fglc1v9XTdOlyuuELy4j7JXsSbGqj2Qv9NYzTN8gPXJ
6p/9LvfAQZMS73eJh+v+9pSS1b0QG/f7PByQgI+VzMSg4F40dlstHDZ6hJOG
9C+lbb04+0wq73tLIq7Ez/mb4tUH0wim5PeEx1jpuG4w73UfFM4WWw4eS8LB
FToVi373ocjwjFAwnowz5zvcEs/RMbm9V52nlIZDE66z57kw0Fr0zmeAnYaT
1Udiah4wkJ4SU9L/Ih3HR+3qmusY2OG94yDX8Bn6vyi/1NdlInpqfjTLJRMx
bSqZllwmss+Yep+cn4Udc8dfv53NwgPuGrmhgCzYP1k7eUKfBc86WZfvRi8g
WXbRwCyUBdkHFWKDLS/xRDFx0hJlNsan56Z5bXyFeXXGk55sZoN+6fEuYcIr
1EmdcbJwZCPX5fxDwbFsNNeWSe5KZsNOR3c9fzwHIva0jFRFDgxTVHo8DuXi
i0VdhfUGDrTkZlznlebi5G7WXWcbDv7+5TQNhL6GZIt/1/tQDh5VxZ/oV8rD
g9Jhm44hDoZjcjS/3cyDZtVIjcp0Loxdqvu+jOWh1PjwyRZ1LgbER6zyG/LR
XXTuZvJeLtbq71gfHFSA2I+SCQavubgh5TgUMFyA5RsP0c/XctHV5vvM6/A7
VC+1n2nex8WV08mKthsKsStdS/SxTD/qMn+LqQve43L6+wvRB/ohi6n2MtpF
8Je1W//Eqx8b/otcH+RYhEA1tUj7wH6EDacN+dUXwXjJeFjVk37o3m5ycnha
jHyt9q7GgX5cqNIyXXmwBPfWJl05dGwAqS9e//c0qgQr1J78rj0zgC8xm8UX
15Qg3+5a2MTVAai52BbMX1cK2Wl5btmPBlAlfkVbhPYRHWLBL683D2Cufrvs
57IyBHOl9b7q8UCpOQu2/SlD4OybtPmbeTgmNVj9Qacc2qKJMUu381DQJnk1
N6EcMp17RpsdeHA4vWokzq8Cl0PHj1jd5OF5ZnC7u3YVMkXbCiNaeeDKcm+a
Hq4CK8zZld/Ng/pp000akVWg69v0aXN4SNgo9ahvrApNP7ukD43xEFF509nu
YzUs9HPtbs3n43TfrcGt+2uhNeeplL8FHy+NhYmLb9XCvqasx3IvH4K0XZbi
H2pRIbshV+4QH0e9Zr0uVP+EzPqABCcPPuz+RASsHvqEisa102YE8REvVR04
y6gOl/K/5LBy+dik77//on0d/rY8PxtZwEeH59LVgjN1+P5rW/KKD3zM/3KD
WZNeh//EP2isreXjdoy5yfWZ9WjuuTw3pJePwMWfpk+01+PUrymLjk4TQNXm
PMvjRz0iWynVxTMEKLquXdQx/Qs+V0ZXNdIE+N0fdiLf4As8bKzLpRcK4PNy
Z6PPky+w0a9YQNcWwAX19wa8GxDqME93zy4BRBJ/5rffaMDG0AiF9ZYCxExS
6KhKasCsfYU35toKUFHprJjW3ADtr8rnsg4KoGH5O9FpfSOm9chMdfcUgH1U
/VmXaBPWJjUtXhAiwMVq89pa+SaoBbgfPBoqwHxtb36BThMsFOcyn4ULsG2o
YEWsUxO21mfpK9wVIO2cRY5NZRO2Scq/jn4sgGvEucK68GaIVmkUTH8rgNiP
xG/v05px7+rz9PZ3Aty3qRDJ/NCMIkl9s0dFAlQp0AxDh5uR1TkknFsmwNIn
qeUmti1gTCgyMuoEGHjbUPdBqRXP+/XM79MF8GQs6ct+3gb/LTv5AnEhuu3b
phRXtsHYplhsioQQe1qva9fS26DQbfVLQVKItdUsX4b8V2yJDH+lP0UI8RfJ
M+cGfcWNBW1a+6SFiAlQ2XTmYDviH1q2Ss8WYtrPuiOBfu24NHDo7U/ic8cv
Bt+Kasd1lb2jnXOEcHD41vikuh06JWfvJ80TQmtbnFvrmg4svPdGh7ZAiGKa
3L11MzvhNT9Hzl9FCN3QikKjpZ0w8r8crK8qRIrkWfpu407YplS9+kN8c6x5
uWtAJ5oCXLafURPC+ltk2X12J1RundbYpS4EN33m8HhRFw7rKbEeaAmxX61I
bkpHF9Isi5z0lwnxOe44Zo90Ya7XSn47cXb4pxCtZd/gvMpem7ZciAunb6oc
uPcNI04iTNeVpD6VK+rf1nbD0tz964fVQlRsvL38vFE3zludNtDXE+Ka7cM1
Cy26Ib5ZVyuXWCz8zdZ9bt3onMyJebxGiF9/B1xa4rpxVn6bqsc6Ifo7rFI/
i/UgfM2lU+82CpE+6pDlOYPMtU4rbxV9IY7O9nojo9CDOzSO8hViullIhdnq
HqQeuXvfaJMQHfmF7HKXHmj5vDfOgxA1MepLimp64L2Jeuy0WYiQbN2VB9t6
8PPO98Yc4m11Buv+Mnpwzaeyc5KhEKWT95tsnCDnt2efiycuOB1+7M3KXihJ
iuqUbxHCLyLuhM2mXtgKLtpNNyLv4Xm670/TXmD2RldL4mxGabCuYy/2iNwy
7STOsBxLy4ruhbT7R+NeYyHcvSRf7nzcixXrDrGVtwqheZOWz8/sRYbOzMv2
xCkly6q0Knsxo3dLcjPxw1UO3NTxXhwbUfyWv02IsJnVmo8O9WFPXvQrS1Mh
8tYszNnh2Yd9SeFRQcR0++P47d8HH3elI3nEGzJoljZ3+6A3bFEpZyYE22T/
OZlPfRh1sBuvJZ594vnkt+19eLH/6MlRYsSIRBzl9GHGtrUsFXMhoplJKSWT
6GjKbq08QWwYyPvsu4GOsy9WZ0lsF8IzHXbqJnTM2x22QYP4fv1teqM1HSzP
1HoTYqGS3thybzpeHty+OIRYYev1wM4LdBiUHmGmEm/1/CpzI5QOXq9FXhnx
w4LzqsxUOvxnJ8dOEFf21T27k0vHrxHv5/I7hBietmiNQSkds3YI2nSJzfaW
mz34RodmXNwlZ+LRHzN8LOYyMPAq0eA98aKFDhMTixjI/1C/tYF4x5bs4IyV
DBTkPDrJJE6OsH04ZTsDLzg2m6ftFKIuL10j146Bpct3suWJf3f/fnnkKAPv
fKbmaxJbLH9UXhjIQKl2CXMr8Xnrod3utxloNH2zxYo47ZxRh1w8A/M+7Ws8
TNyYdNe5/BkDbebKdz2IRWo4wlP5DAgVfcLOEi8d2uCvWsFA1tGyt5eJreXD
JOqaGOhadE75JvFlg+5b5/oY2K0x8S6SOOPoKvmlg6Qe3cyoB8Stt64ktf5l
QHVkIPUx8aTXzdrXpJlYO2fOWBrx8i6NPB15Jr7qXgrKJLaT8DfsUWeCNum6
VTbxNa3a2rDVTOw2OOf4hvjFHiXbjYZMuJa8yHlL3OF3opeziwnPSydNC4mn
JJa437Vn4sx9GbUiYp3KOSNb3JlwOlBhWExsL3S59N2XCb9dnPR/DpmXL5UQ
RPa4ta/s/jlnk3T09igmTqg4W/473+1krzyeyETIEdz7d790aFZ6ahYTcc98
NQuI12SLrbYuZMJEcvvUPGKHdsv34jVMSHTy9XKIp8X2eJxoY2LR7xXVWcQv
7TwVvjGZuFcszH36r1758SrzYSYml5uNJxOLfw3yzRdl4d3uA/fiiZ/dm62h
MYOF27abImKILfc+ao5SYGEgWIIVTpzUmq9zYg0L7K6o1AvE5jFbe7u2sFC2
5NK808TDNo3h5hYszB2LZroRG7XweOoeLOi4eab9ew/8aL8HUX4saLdQ8SbE
0daTzcSvs1DjeHBsIzGzSTm1K5EFt0M3ulSIw6IyrM2zWEhgxbjRiPWs1kvk
v2PBvfb7oUnEQY17HKJaWTDJ9blOJ+9Xo+Gagvl0NlKGMxlxxHURtOq8BWyI
3/6vLoT4rEWCr/oSNp5L3NA8Q1xZn9cstoWNlRuqVu0gdq0bCM/zZSN2v0/6
D9JvsuG+UA9io8L27bIu4vydkvzISDYWDMlp/OtPqc9KZl6ZbKTKRV2NJH5W
ayGhzmBj0/Xd45rElqHfXkV+Z0OGr7dzBvFvc3cHMREOlBsGlYdIfmyvufq+
U56DxBfqBnnE/Ko3vpG7ONjuPRaCf3kTYqQhZs9BwvvxPiXiTaZfmo+7cWB6
W7N4guTTrcp+HbNrHJxIVpAqJF5RocgXLeDAOvnYd13iEx+vOBwne3VoSPFL
aZJ/inRjyZ86XATppJmyTYSoFpuaHkBxMdOj+2IJ8WIqdDCY7N1WSd5VZ4nb
8qMuJt3ggnfcsaeL5K1BVsrDNgG577Vz1D2Sx7NiK74avulHU+Ok+BoyDwrz
Qs5Xlfbj7qUQ0Thit1Zz1d31/ch9t1fWg/jjnPpj9tx+DH6Po6SJfcPbRs8q
DCC2gNZhTOZLz1Xu3OcXB3ApZebxDAPSf17Se+ab8NCh4kffTOZZzX2ab4Y1
D625IzISxMxS+fjNjjzsZ2Uwy8j8k5db0u9+gQfB6/gnJsRXio0Ci3N4mND3
PGVG5qe17IVXrqp8TBn+GketJ/P21SDt3Tgf+R65msNkXj/yd+OusBHiRoB3
UxTZB1YUTR0NOiDEaosbm82JiyalSnw7IkTC2yvnxYi7Q+nKoV5ClHvJO3ou
Jf834YANJ1gIFu3jgKEmyf/SnaWPCkjdd4M1mWTfiJDRfSirOogLO06PTVlE
8iz+9+7hgUH0J8lyG+aTPGzopgr8hmCifOviFrJ//XfK+ayl5g9Iq/TIrGwU
gP7jrrXUu58wc9MQq04UwJThei7ixCjauktm7ncTYAc9fKmi4i/09U7PjNER
gFvicLEiYxz6Ti5dH//yoT2reaHhwT8oTzfdyirlYxp+HO4Y+osVL2UiIsL4
2MqzG1G1EqFG+kd3he7ho/WDo6idlii1fsE62QdyfORcL7VX/iFKBSoZhaf3
8eAs4dP96ZMYtdb3cpdjGg9HvP1zpO+LU14JN00MTvKw3GXpJ/PDk6gNsZVt
fpt4+CyrJeWzToIal2lyjZXg4e3zy7Fa4pJUoJn19oG6AfiLau6q6ZSkRP/g
lXHMAKZ1+G3zzJ5M+VrLFKsdHgAjRSUyPmAKxTi4nb9EewCRGnKT8qymUhXK
e9sOjvQj/L5z+WHladS+lsauV+Qd2kUe+CTNmkZ9KNVb5XenH6NbksfiC6Uo
34jWO7ft+9HH2Hz2Q5A09Xi6gdLzZf1I1F+5oPawDOXbIh+tM8ZFxlWfWS7L
plMP+AW2mRWkb6yOmIn/mE4xtBZ3R9zjYnXHo7Sg2hlU6D7T3b+cuYhZEptR
azaTcooN2rh+LRcXikJjzPJmUmYHXPVoklxsX2CVXKU2i7rjEdVm28SBgWf5
vlnBs6hU9a5mm1QO4hUc1Gp+zKKsVz1eOeUsB2olT6J37JOl6hdpXjhsxkH4
sg+aXwplqRkdElmSyhw4D02qr9SkUee3Tkk8OciGw4vvysG3aNQnMacGtzI2
3C/v+W16m0axpGIePvvIhs+e7BbpOzSKkz90glfKxvUfPmHh0TTq85pFK4+X
kFxcNzZ+N45GFXCsZp0sYmOseKI55SmNqpqRZBuQz0b4F+nQ0nIa9dro5UTk
czbuP/Y4eq2SRh1gzvvcnMHG41OfDLdV06iTUXdT5hPnzg3/Vf2JRr0z6/N4
8JSNDjva0YYmGhW4J/5AUiob6n1yhr19NMo6ybA4J5GNFdl+ikkMGpV8nCY7
+oiN9Vfbx5xYNErxzw/X9cTm6nFZHC6NOndjUPd9PBvebsqKg4M0qjcrZ0nl
AzYCNl4aezlEo3ZWTH8rRXxVprfx1A8a1T5+2nJHLBsxmY9vjI7SqH107/iG
e2wkXpzkkv+LRrVkSuyeS/x0t9PmgN806lRp8rS9MWxkq5Yt3PSXRl3IsKiJ
vctG4dB/YxMTNMrgp1R0VzQb/wOdwtfF
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#2"]& ], 
    TagBox[
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVUHc81f0ftSIj6/ajYfVkhaR6SlG9v6EIqYyECjcjiRbJSIOMn3AJScqI
IqPrScN95GNUZO+RkXVd8w73Xr7XyK/fH+d1/jnjdc424jVrdwE+Pr4ff/B/
fhFKFN7tFnzkbZe2BY18Fn3tULG0V3GA49UFYYqJ39Da94owFRUvCDBMknYz
a0aSiR/JYr03IVa9cZ+vVzt6URmrvE4lEGacrQ6EKXQhHYZb3JrnPXDZ1rhp
5XoPGpR8+Xul5yFsWgoLjn7Qh1Rfi+5fWn4EIhajt3Uu9qOrh/19FpWj4cve
Yx1XZAbR+87hHI7xY7jq3iCakz2ElrwtB1ie8XCE1TeVKT6MTJoS1BZ6EiGa
Qq1RvjCCJlNNqzheSRCXuVi7f3EEPXZbPc9eTgbruZ5VwehR1LF8OYmp/Awc
P9llbn06hm7XKukxStLgn/BbJVUbx9GWJ50Nc8bpQDts5zL6eBwRtTDBGc8M
sGnyl3bxpSKB5aAWTkgmNHLc99zvpaLOYTtFzrEscPNpqco9MoHuFIqXzfdk
g/Fw19+b+WjIPHFCeD7zFZwYjD8w4kRDCneqbFleOfA0VWQwp5SGKo0DmIzl
XFjIyGjouTCJRPpHNeaU86Gz2cBF9ckU6q384j87mQ8zOy9PqPZPobevU2tm
St5CS9RGR8Nt08jqppXztHEhiN6/LbM+bxqliFJSaJ7v4IJfwBun7BlUGmB+
89YmMlywunTcdnAGpU/rb2aHkGFJMP7oktws8m2V9Zw/VgKa8jbL4RGzSDa9
ToDV8w9YxHsFT9nPoWXJj/nXD72HcINDD95Ez6HxB69OMzPfQ9rBMwK+lDn0
0TP0JcOrFII1WpT1N9OR496/DejLH+DTIeJwUiMdZdVn3JhRLoM57lrf1vVM
xEn9oPXrcRn4kvxeBygw0XHPhrF2XhkEHD/RHqLHRLOCi3aUDgo4NduXxtoz
0YHDVgbRkeWg+OVNuEQ2E7W+WxHQYCDYPKU2Y6bFQrIgenGDbiWEBuS8xQxY
yFA9ySDSrRIi6kqU+c1ZKI6Tzw5qqwTLCng1cpmF/k7ocicWVEFo17V/TXJY
6F69jvlu5xoISIm2jZOdR3kln9QLkmtAZbhGJUN5HrWnGgmqNdZAdqb3aLTO
PFL1PFe+6eBXuLRzLJB5bB7VC4br8hG+waOGJwT5gHkkd7hftuX7d/igJmfh
2jGPit9F91/VrYdnI2TK1wA2mpadfmzuWg8iU6+ycu+zkcZt8yOaSfUgpVNo
7xfNRpmHxLPGePXQUpxhMZDGRok/Hns4fmuA929EygbL2ej2WDzL9HwT7LT9
emZwlY0yxBvCZI61gpHf8j5Tfw46cjj4/P2LrZBzzf6LdQgHDfhq72MEtILs
yEimdRgHbWqPmWh82wo3fzrI6CVwUEKq5Yko6Tb4KTwW7l7IQWFqzZJr/W0Q
Mmuk2DvMQZ7Q9mz2ZgcUbGk7WWDERXzZC5T+mA74xTW6RjTjolQhhYH6nA4o
Y6GfMlZcVPfDQym/+4/eeWbEwYGLNG1Xst0NOiE6/vYQ6RoXTV7WKBzi7wKp
JrZ4QRoXXUm8W9FK6oZJvfQy4VkuEuBm/0L53ZBEl6+eYnJRmn0d37vqbmjj
u6tRy+WiegWCcSynGw7W6R2+scZF2m/yak+c64ELt8TwMNkFNPtvR2u1ci84
ikYbF+svIF/qjrHS4j7oG3RzuhS6gIYv9q2v+tEHvLF9fUJhC8imN0q3abwP
Jq8ePJkVuYAONNACqVt+gjb2wKSBtIAES3Kl5SJ/AjVrLnM8ewGlhmw7EuDc
D6sNx3bs/76Aqgibnx2UHoTZPGlWiugiEtsWrvHrwDBMPd9dIPVoEdUdStgV
emwYdCcFPG5HLaKIcy/1Fa2HwXzgif/PmEUkQPps6uQ9DCTnxzlpiYto6fes
Z8+LYcjcTqYsZyyimQG7vBaBEShxsfYxoSyixlSNHZWNIzCYoq8tPreI4qQb
tLJcxmBjHc14pyWOyvQVP1j5joGHQ7+a8SkcjV+8BivBY/Bs9dvoWWscGRYR
bO2fjoExv/VSwDkcTZ44f3dD8xjMBwXgmW44Mg6bawk0HIfqbcQBUgiOcK6U
v7UcFdLv51qvvMXRdkXi2tp2Kmi6pZ3sKMKRlUlpdNFuKhBPL93NI+MoN/Hc
y/UnqUDWbok6+RFH1ruyaivC/viNwv77sApHRZf3bNFmUUGUK2uZ1oOj3vjw
nN7ff/IkSBS7nzgS+tStGyExAfZDGWFSgzhyXBdsPKIxAT7s7bqhozhan11z
9enFCThLjFIxm8MRsd8WCTZOwPSRcguSAA+JPR/xudE3AVrUu/Lq63joH0df
hV8TE2DjTb5KEeEhwZ+RgRR+Ggyfe0cZkuChwmcbNTWlaHDqt9wnXykesnXI
6k5WoMGslLrJqgwP5fRS9t7Qp4FguaW8nDwPWaaajg6Z0KDYneqSuZmHOPad
JEtrGljlaxjtUOChdHlXoDjTQE6psZysxEPHeubmNHxo0G8Y1b5/Gw/RU4LS
k4NokLIYuN1QhQcpZ0UsBKNo8OROUPsnRR4ckUviXU+mwS4nu5F9W3kw0aWS
N5RNA5KTq9WHTTyISy46a0n+k6eHy++X48F+O4N1lC80OBpx+9BnAg8iO22I
yb00aMLfxldK8kCzI0LBUnISarsU9owI8eBK6yypLHAStGOei7qycJAlBYJG
5CR8jc7+ukDHgXJKmJ6UNAnlur1VsbM4iLcoW1x/Nwn6zashiIZDYZP1Og3q
JJi4CDTvHsKBXv85MOn0FITZumSdbsDhxrdw4jWNafiPm2mZz2sclMaPCy/s
nYbMyTuyp3NwaBAQfRuCTcPDYRfnvdk4qGGxrGiHaaiICG1deYFDHyX5fk7M
NGjf9eKmJeNwlPz6ZR9jGtLznYIEH+Eg87zup/HnGZjPJB+854ZDyXUJm00n
5uB4k3UXWw2HrGDvaT17JiwIRtoY5S3Cw4yVM5xZFhBjXwfqHl4Evo5hrDyI
DcwXWTNbBhZA3c/jjq0WFwgT52p0ghZgnPv0rPiXBeD5pcm8V10Ac+qVu4k3
cOh0+WRIauKC1ThJW0lpCcZ99lwfDOHCdA3xfl3RMgQG3Ch138cFXZluRWPn
VZjh6L+MmeGAGHBdB9i/oVe8eSM5nwOmc46Lf9nxYSvHhcXt3DnQW+3G76jD
jx0dVxpcVOXAh6ivF1W4/Jhj1Plp9Rk2eKzzH25uFsD4/2MqTy5mw6WbwR8k
0gQxoaGuavU7bNjlqd1s6SqESU1G/HXOhA0tsjri/gfXYVHrLTMGNrDh3+KH
z3UEhTFOh5PLo4F5CObXOt04KIxdXQkvSMibB7GBIDPfUhGsfrtNBH/APFBf
b0vKCFmP3dqT6798Yh6SNDcLldmJYhem/xJ22zIPpDSPWlcVMcxdUbhIksEC
x6QLzRI0Mexw2mpgzDcW4Ca5vIwKcazQ7MzjiGcsGKMa3amOlMCeE7faX7zB
guzDu7c2uW7Afie9E6o7zoKiR/4ynjslMa/8OP4BZRZE2l2yEORKYm20X8l+
S0zYN5CVH9kkhXl8eUkmtzMhdcfzoiYLaayYVrH9ZDET7lXGplqUSWMHNsj4
Z0Qy4eRWu9x6VRnsVGpfqCuRCUd9a51komWwJpkuYxOMCRkKRNVGrgy2ls+u
oW9lgmrNmxQrJ1ns1Xm22T4uA0g7q7XaK2Qxib1cb3opAzzYQm0/tAhYYJBV
iMN1BhBL5lWi4wlYy3mpcI/dDLj60GbFPIGAnSrfVHRKjwH+NqU9Ek8IWNqV
xKkDuxgQxfWPI6UQsHTZhnTxnQwoPshbfvqCgKnZ8M2WaDKAV7XW/bqAgMmp
sbUEVP70t0vEfq0lYDHPqXIlkgxIe+VzOeIHAfPiuOWlb2DAK79mY7MGAuaw
hltGSTDgoxxpqaGZgE1N2rQ4izFgwJFwuaOLgD0jRnyWEmaAxthm49ExAsYQ
8d5ya5UOeqVBSjlUAha3n9rvvEIHg0f9PHcaAaO6mJItl+lgqfGCPDVNwHhP
Pieo8+hw01tFicUiYHzi/Pr9HDqEHHrA+4dNwCzyfeLr2HR4tGG0049LwE7r
FC59mKdD6rtXMThOwIJGn4gmMOmQfV/Ik7L0559U/aJQBh0KzrgbhawQMMMD
2Zeu0ulQ+td3xSO//+xZ16bjOEeHCrY6b22NgHkbV603m6XD/wD6wp1J
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#3"]& ], 
    TagBox[
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwV0mk4lH0bBnAhWSKMhDQmRFEq7cR1l5RsLUhpkT1pkRJCpOxkiUKWUApZ
k0rlT6lkKYzd9Bh71rln4x7E6/1wHefx+3KeX6719tdPOvHz8fHVLt3/M+2u
vdB2R199yumiNKbuCVRDpZhaU85AvZn1pmr3KrT4vfI+heIKv93qKyM2/kAS
8eXFop0eEHAsyOcl1oDSqqKVllN8YMjYRO8NNKHNDMeHiy4BEFirI7Ig14L+
SKQvzHcEQfBZ8WfyVVSkmiOye3YuGFpOxg4MY23oip7n1RmlcHDn/hPakt2O
3rTSn3MMoiAq76yQzUQHmnUzpTFdYmDmngjDmdKFDjXGbZjuiAfRV2/8gy17
0N+kI9Uc1wSoz+tGLadpKMrx3zn2XCKcmY9S/nH8D6LOXUrAlZLh9TJO4H1K
L7r9g7yNUZICVt78Nwxme5HCo9b6SYNUcLyt2tfuREf2GpjAuEsGRHrrb5fY
3Ie8X4t9YHVkAYwmONTF9yPj+GEh1rNsGAw2iagb6keK3tWWTNfnYNp+9bv7
7gFUZeCFM+ZewKjNrRO7mwfQip5+9UmlXCDtOFT6iTGIOqs+e078zYW5Bf6e
dTuGUF5O0tfxkjzwXJVwgHNrCJl7mNuOGbyGMDez5MOcIfRYpOLxiEsRTH4P
99McGEZlXsYeN+WKgSu2wW9s7QhKHdsjz/YrhnUFBSKxFiPoWpO0C8uwBA5H
eZu/RyNIOrWWn9lRClYlUfOBUX/RnER5rvv+NxBimJjh9/kvGryXfRx/9gbK
/erink/+ReUud9MZrmUw7rjiZILJKLLZsVNnau4tyKZWeKUvjKLMuowb40of
wLNMxb7fcBxxkt5q9EZ9ANaY1hP36+PosEv9QAvvA2QraJMdksbRhMCMVQW1
Ar6tvqIQ8ncc7dUz1wkP/QS7Nw5abg+aQE1F8/zqDARXqUNlpPRJJA0iF8S1
qgBTigl6jCaRrlqCTqhjFbjjf1NL6JPoISeXfae5CpaR80Uq1k+hnXFtTvb5
1XB/b1br8YwpFFC32Xi77Vf4k/xxjVEEA70qeaeWn/gVZkISgg1yGKgl6aDA
hoavkFJUS+ZWM5Cqy+lPcvtq4DcyCAgkGKhO4IEWH+kbKF/qrdLajiNZvR7p
39+/w92jNbhSIo4Ki8J7rmjVgbfy1vutRkw0Jj0WZWxXB8L8fuJJJ5lI/bax
/saEOlA5ILji6DkmerZfLHOAVwdONPzEretMFP8zytnmWz00m2oemklgotsD
Mcwj5xrh1BvBoH4aE2WI1d+XMmyCLP9CarkdC+nr+Z4LvNAEO4R7ah+5shDt
muYuhlcTPODqeTneYCG5lsjhhrwmOD8zE9cXyEJxSaZHwySb4YV8QJp5Bgvd
3/BLYrGnGXTW1KtIdbOQCzQnT3hQYVeXdbSQERvxZU1X9ERSYWOqefJtczZK
ElSk1T2ngkxB9fH/LNmo9qczObedCiW5vqsT7dhoo+V8lpNOKzQT0Vuy7rDR
30vqr/9b1gaqmluumOSz0eV4/8qm2HbwrqZ0W6zgIH5uVi/KbYfYrxaNtJUc
lGJdy1f0pR2M2K+ibaU5qE6RZBDNaYckM/nXJ9dxkObLVz+Onu4Avoenqbg2
B018pDZ9UeoEA58g+przHHRtaNNAWWEX6PbclrYv4CD6hS7h6p9dQPWp3PGi
lIMsOsO0Gge7ILxN913/Ow7aWz/iM6TQDZQkaxmzLxwkUPJCUja0G8KGS9zo
7RyU5Lde38u2B2yLvAzEFjmomiSfvE/yDyh5p7W+N+Ei0fUP1Hv30sH2DznZ
upOLavfHbb1rSIfFheGgPTQuCjmdvmfdSTpYSUQ4yNC5iD/2/ZGzbnSYln2j
+G2Ei2YXJlw60uigTbZn889w0TjN6tVv/j6IGszdxZGZRg1J6puqGvrA+XoT
ud58Gj2UrNfIvDgAt5rpkeSP08i+xxIJNAzD24CYCL+IGXTj2wP76+pjoM+9
oL1Jh0DkwcNC0zvG4FSMh/LD/QSq5xfJ88PGwGFt6TRTn0AbsGhm+Jkx2BNI
Pl5uQKCuisTA55FjUOO8c2i7GYEOFOekdzHGIFbYNZRrSyCpp7XdBu/HYZ1k
WtvbBwSq/BBxt65mHN4Xap/ghRDIrdNU+UTzOHzOsW3WDSfQt9XNrhfGxsH+
ydaBz9EE8ontIrwVJ4D+eatz+WMC9QWPyRYGTsBzdk+JzysClbivtJA7OgmO
H6Ltj9UTqCGF5FNwahIKA+x9shsJNFyjkHHQcRJGdxaGTf8mkIL8pvErAZNg
uDkzIbmVQA+qDe9Xv52ES6+c/Kl/CHRKOuDNZeUpyHB3cV2DE8h9f0jX4tYp
KOW+8z3LIlCEc/Rigt4UrL3+MTGdQyBU8dQEnZ6CTWv9BtfzCLTR4f0AKWYK
EhUEhdbx89DBaCSSmzoF/+Zlrp8R5KFz775v1c+bgjXbZ4cThHgoTqzN99K3
KcgW/SgkIsZDs2+YpM9zU0APUSMPk3hI5j9i30kRBjy6yVq9VpaHtIT5Lo7I
MoBY46JyTI6HjLRXhPipMmCw8cyRIAUecjgn8VpKmwHHfqUHlSny0N2Q1S05
wIBdO0Tah8g8lFSsSOiaMaCCd/mg7HoeKu1WITfbMGDDKXZvMYUHDYKah5wv
MaBXfzsEkXkwrKV9ec6TASGXhO9bKvKA78y+2Nj7DFi1Oz9PXYEHuwqO0CrS
GSD/UjmraTUP0q0cb61qY8Bw3qe1W8R58D7ALeV5PwPiNRpdhMR4QM31qNqH
MyDS3z6OLswDwQBMqHQ5DtfaxO48EeTBXqr5SRoFBxVJlahV8wRk+rqNbbPG
QU1nfOXGcQK2VYkQoedxqNk1E6Q9SkCV4KvlvQ44uJ2jFumNEECPHqREu+PQ
pFiseWqAAPKz89aj4ThkeoinxvQQkFJzrCbzEw5iDYr9Kg0ExIvvTJdWZkJf
5CfdXwUEBGXMn+BMMEHXWoIufpMAPiod+3SHDRWxv70jhAhQu+XsbanBhcyY
C5wbyTMwyH1ySuzzNGQ5/XczXGcGjIcu+8ffIOBW4funst3TYD4Yq0kmz4Ka
zMdNovemYeyrfWBtwRyU/rwY4qQ9DVpS7esMbP+B0cJNi9N9XBAFrh2NvQDd
FwMGIZkLRyZtZpSt+DBXlztmyeZc6PziuMxm8zLMzusC8hbmwtuwmgsU7jLM
2DG4OPUnB5yXe9J//eLHNjz0i18WxgEHD9+3K1MEsAi1F9tYphzY6qL5y9RO
ELvsn5WMy3Dgt/RmMc99y7F2s5v+qt1s+FgY9HSzgBC2bSLffO45G3yXaRxv
+COEaZI93eWuskGUdsfoWtkK7Mg1xQe0PWwYylmfkOEnjFGWMbbFr2BDwkZ5
wQ9WItjVVazbSa0siE1x/mFHEcX68qt/6uWwwCbh/K+VI6IY+wf98IAPC4hD
L3gZlWJYYKXRoVFjFgwMHfT+EroSe3zxSIcChQVZetvXNtqJY1D1lVXNYkJB
sKeUyxYJTO+v7vmzP5kQauVgIsCVwCw2Ss49zWTCLlpmbmjjKqwlcTHT0YsJ
SZueFjSaSGKvH4TJNRxnQkBVdJLJB0nM+syrY7/VmWC21upFnaoUZsi22mHH
x4QD136clQqXwnxD/W8b03DIULRXbeBKYbuV1jj1v8VB9evLx+ZnpbFDrrtj
CuNxiN3yRaOlUhqzP1pc+WXp75zZgs0/NUjY9d02Y09McbAvYVHCY0hY8IE7
zm6aOFwJspg3jiNhVLt7gZ0aOHhalHWsfETCcry1XxguOYzr+TD2MQlTmYyU
oGzCoXAfb+5JGgmbqZq82KaGA696sT0nn4R1/pvnP6i8tN+yMrrmBwkb/NK7
TUUeh5Tsq5dCfpIwyzleapwcDtm3fhkY1ZMwkzURJL4ll8vGztb/ImFdj5ZT
/sjiQLMhXaK2kbDayf3lT2RwUB+QN+gfIGFpHKN2CUkctpXdIT8fImHZJ/hH
/VfhoBPcw3MaIWG9GnJCkxI4mKqnFY+OkbBCIQvbenEcPNwoZCZzqf9yV3SY
GA5+++/xStkkTKnQmD0jikOweH/rLS4Ju9cn5uCy5KSi7EiCIGF63V52h0Vw
yAoUdKmYJWFu95Twt8I45J9wOug3T8KcaZLhG5Zcpvx9nf4CCVM8b7w5cQUO
lWw13uIiCaN51nUJLvl/vLiFRw==
       "]]},
     Annotation[#, "Charting`Private`Tag$11606#4"]& ]}, {}},
  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, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  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.8737300176193027`*^9, {3.8737300596459723`*^9, 3.8737300665241117`*^9}},
 CellLabel->"Out[15]=",ExpressionUUID->"3524f72a-730c-44fb-9dc6-30e13e9dc263"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"CS", "[", 
   RowBox[{"e_", ",", "u0_"}], "]"}], ":=", 
  RowBox[{
   FractionBox[
    RowBox[{"4", " ", "\[Pi]", " "}], 
    RowBox[{"2", "mass", " ", "e", " ", "100"}]], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"Sin", "[", 
      RowBox[{"\[Delta]", "[", 
       RowBox[{"e", ",", "3", ",", "u0"}], "]"}], "]"}], ")"}], 
    "2"]}]}]], "Input",
 CellChangeTimes->{{3.8737301794533577`*^9, 3.873730202629624*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"1d64fb32-c3b8-4ea4-8f06-4b01e68febc7"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "2"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "5"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "6.5"}]}], "]"}], ",", 
     RowBox[{"CS", "[", 
      RowBox[{"ee", ",", 
       RowBox[{"-", "6.8"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"ee", ",", "0", ",", "1"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
   RowBox[{"PlotStyle", "->", "sps"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.8737302193640547`*^9, 3.873730285086452*^9}},
 CellLabel->"In[38]:=",ExpressionUUID->"5bc7a9f5-c5f5-4792-852c-3f1f56196559"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVjmtUjAkYgLtRLjuJsCNlsxVFrok2vW9klcas7aZtNXQqoSVKy3K2Zag2
umGTUqhcYkeNdFmrccuYlPTNN9/caJrJNKX7l3NkpDPNtj+e85zn3+MYfTB4
t5mJiQl7gv8duLuHfNobBhLfnPP3SAVY+M07pGF6oAMrY/MbsQJWXysormZu
Qqf61ZZmE816rbh/hRmKpGMRL6RVAeGx+kcXmLHI33vhpL5ZAe3L//TLZibj
S7vtP28QKqBk9rb0s8xU9J1FitrrFMC1ql50hpmHhjDVj17FCvgpopQ//OEG
Si/XGWxiFDDNyI30d+Cj2UiE1eFlCnAOGup5nlmDIR7PyjPH5HC01QsLPv6L
hwynfWxfyGFujY3a1PQJMtidZ4cvyKFO69dCrG/Aqfyk3HCOHBrir9hFDAix
ZY6gPNlVDiMuPtYVBxqRyTE33B+Vgdkob4adczPa7lHqC0Uy8DRG1YVUtuBi
97DAjEsyaEn1aBW5ERi8JYD6vE8GC871j94ZEiO54Ga7u5cM+vUatWyFBIcY
Dymf6TK40VffzMih8PPxF1xXlRQY9uW0TifFf1ZzmvIrpJC1rzauxEmO9oVt
PjWnpPDJl6fbe0KBa+FIVmOoFJJXCWvmPFfijry4r391k4KmUphJOr7FIHZK
m7uBgrF4aZHr8TZccuBdNJukYMkGy/G4ehXOi7R27bhNgX+Z456quWrMlKZt
zE+h4Kj9roXvIjW49bfKi7wwCpYOhE1TsDvwi7NJvnEZBbOsh62EXzoweGFZ
fPYkCiyeFPG/LX6HnkkB2nCNBH7xqH8zhaVF0fCazoQ6CfC2RJc8/qRF5fwG
gTBHAiP982Z053filFy7pjXxEqgN8TfP26TDMpFlgZWfBE4EpF2r6dbhzDmz
E7ocJODYHqi8ntaFEeyRh+oxEr5bm5rHXNWNd/nvS9lyEiIdiFFTqhs59duF
t6tJGGnO1h7ivsfSKdNHz+SQUGCxMznYqQd7lIpC//0kzFLXlpkJe3BsXaPP
DBYJk1K4oV6JvdimmRuXuIgEkzs7+dqZfehTlP7HX5NJSLrMz57Z0Id+U4v2
cTrF8GhJhpVNQj+uDPIeWvVMDDaBSUtPMQbweIpDXWWpGFTuYuqYYADDztu9
6jwhhntX0Um3cxAz7rocrIoSA4fvaWM/PoiXJlevP7hRDCfTSrmqv4ewyoX1
4JOjGCp67IsYW2m0kr9UjpuLIeNKWpT6BxpVXVueHpvo2OBB58ogGm0DbSo/
molhvkBwjx1O4+7YxOI+UzGczd0hyoqmcTxGnlFhJGCvZ+GHqcdoTKsnuFF6
Aham2gZY3qJxba6XhOgmYHzd718pbtNYxvBY2dhFwJvBTsktHo2SspLixzoC
zoVXczZX0ZgjaCmo0BJgdAs+nCagkV8Y4ZKlJkBF5l41p2hM5mRfZMkJeJCu
j6FkNPZ7ba73kxGQ573L9bqSxkBDzKC3lADWzWU1G9U0ZtLNR5ZKCHj42+sm
bi+NbhnLHRitBOS7r8kNGqAxoeTzrcmvCUjUFoc60jQmWX+z3viKgMVb92ue
fqSRyUvPHG4iwMJUeuOcnsaTHsLve18S0FHrHR/1ZeK3NXK6tpEAQfz15SsM
NJ4u36Z5KyLg0oJpI0Yjjb5DRY+pFwT8B4Difmg=
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#1"]& ], 
    TagBox[
     {RGBColor[1, 0, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVlGk8lAsDR5HiITtZKlK4CpPEkO0/4RpZcm0pISERkcstJes1bVLShFCI
VCRljRnRZKmunVL2LGXt0SJR8fZ+OL/z9Xw6ip5B9od4uLi40n7zf1seGm+v
nXAylm/iZCw024LXVO7YoKw2rmTrX/fNscX2zNSMElkzyMw9p1qctIVVc3fx
DVlHTAr5+dur2MLZe746SdYbGRF6dG/GbvRvPWuaIBsK3q8zhPIeGzAzZWqV
9/+DC9wy3EPbbGApkm/4JOM4POVn078I2aD8Y5MOqXAS5ZEGZQqN1rhUKK5q
pxIJ7rHaRD8ja2RJ2Z65IBsH1rq4zw16VghvTqnbUxWHkhPshWhZK+xhDPFs
2s+AkMA1nthFS6yeC45iZZwB5/p+F6dqS4R1McOmFc7j88+Mjj/MLWF3tefI
bpVL+O7tLf3Eaxdi+Ev+OC/LxHdDu8aQF3RsoV85V3yMCcXcC23RRXR0MIIm
ehuZiJrYWFt/jQ6lFeoFlOPXUE35zvjqSUfDr1xKV0cyoqaYlR3L5lj9JVlH
4eJ11IXbaq8wNUdq/ynTsqWbKP3FP4opM+xcty93wCkTa84OlSp0m2HSRXcl
f2Em2mZrWMrPzCA/fdKrUCMLwTcYrFvpZgjyoypz9WVhr1K4gNNuM4h7F929
rXcLzv/WLZpWmWLvvuyi2U+54JX293l7zwRkr6PALbnbCJqJe/MmzQRn3PkP
OZjeBo98qcRsvAlKvILkypm38WjmH95DgSYQDjRihOvm4dSNhlZPbRPUx751
WRVxB34+Z7Si63ZC677YqnX8+XgdrBH5+QsNLyj1Hs2a+ZA50Ttx/gMNHo/C
WJH78jEfOmmwvZeGhPLBY+/y83HQWX6hikPDeG1hb55NAfjl/SrTkmi4+WrX
o21X74OP2lsarkOD4HKMK12+CJ41li2lQkCqn7O7unkROoOzHbBgDKUudQ+x
wCLQjtl5j44aw/het1dvdRGM76q6RLCNUXkvuayC8hC6fgekxgKMMeMwzbZq
eYg5jaNOBzuM4HQ35b8QoWJEpJ3d0/7AEBlusdY+isXYE0gtsL9hiGHxo817
dYoxULGwNB5viKAIk1Yjt2Kw/M2EHP0Mcd5+poOvsBjd5yRYCsqGqP5p0pNu
XYLj8/4qgVkGULb7OP4svhQLrsELv3L1IRdF1z+cXQpfB6VhJlMfooVZ8QIV
pcg+x7qLOH0s8jtQ7IZLUWa2iDYvfbTVVoQM6JVhTfg29ttN+jilGbP0fawM
7xg6Ksfu7ECLiKQkhVaBmpeirk0sPTwzOnqo3akCifLyDpmFenjs31Ae6l+B
CXqKcWymHnIaw/axkitQFR9/9sS/egiL6r9pMVMBrguPl9dZ62EjmbfZK+0x
nFe06xoN6eJEyw6kfq0EX0wtXU1CF8KUsrEX/FXoDnHqyeTTRW6C5sUf66pw
xUKSo/CTinbrP966/1mFtZNXd+8co0LtP8kQ5eQq/JlBv0I8pmKwkcwr1mXB
4HxSe9YBKuhPbws3n2Ljpv/z6ehyHUiXig1wc9dAPGdzjkS0Nlwv/nAcXVOD
aJs7aarHtXHLe+y/BvUaCHwOnTcP0Ia6VGVl/N4afNRuLsjcqw2T4x7XpB7W
YJNjdGLNNm1QM5zzmFtrEbe3WfHB2HYcdM9zuztQi4klq65PDttRPmza1GrI
QfHBR87FdC28toj5lEDn4Gek+J1JaOHbgydrrO05uG8wGbhZTwu64TsOPj/M
gYZBG71KVQuPJbZ+q03iQK1TnqYqoIVKM7kNxeMc6P3jG0tp2QZ23mwIk/kM
ggkSiQ9dtoFz5MbafdN1SLJrPCsWp4mSCX9/9291iJ7lerQ+UhO5vvosL656
iMS4E5QwTTB8ul0CJeshFTmjsz9AExaeYmn/Gtaj/u8LvB8cNdG6jyHz4GI9
5man2LUqmuizOCq1QqMBjeYxIzVpWzGnYiRSeLQRF4tN1n6xpIBnoUB0rfJL
qOwaP538XA3UZY9yhwdNWM0zpvK1cTOa4rRbGra04uvp9Ub3xVWhkDi1cO9j
G+Kf2jwtDFHB1PzgwCvNDjx8mmV9Ul4ZuZOsl8KXOlHSGwva600QXn+HHB3t
QgOtc4dXxUZc9CvzyVJ6DaJmdii9QRHfaAWjvlHdGFwVNR/QswGhWnWla569
wdirncFG4hsw+KAuvl2xB55smuN0qjx+HOlK33yqFyMVmwTcqOuhtpNvyYfV
B+KDnYbMqnWg31I8/Eh6ACbKFtmVr+RwYv2Bje9cB1Fs7yPZ3yIL9WknwW6b
IaxN5D/NOygDCZFZ/rrFIcSuHGJxpqXBW5NetCnjHbh+BrT1yUrDX5v1lrAa
xnZHwcPbzdagYJdn1pNvw1jZ7NUvHS2FuSk50ffJI3ixQSIFXZIoc6CvYJqN
oj2ekjuiJIkoC0Zm6ftRNF36kp17QQKK/ZZvchhj2N9paXb5hzj0deOYslrv
MSJ16mWBrzhc5VsXuDvfg3t0ern/vRjmXiYMH4v5gFivxLSgg2JI5XUPtVca
B/gMOGpTopAYKLvFUzeOGnWLvsrTolgZEeO4I3gCTrBbmycgCq577kXD4pMo
/i6a6mopgr/TihLEOZNwv+rsmBImjGq1c/xigVNgn9Vypf4eqpjl3+qxwtOQ
ExET5F4phD6Nts6T7GkUZsR/2WK7Gg9vQmnUfQas50P9V+8Lwq2IKrZ+aQYj
rTP+URKCiGZkx/Tlf4Rq57Uc0VABFI6vTxe2JjFvOX5TfYzAuRsMj4HdJGiv
lL5XDRPwtp9RfmBHwmv7Z9auIQLr2OyHNs4kmtQ6Y317CVy4vL/hoicJLdsq
gYJ2Ar7U658ETpLo+evGlOkTAhvjJC348kiosS7p1qUQWNI7LdR9l8TTacWi
A9cIvJ0Z6cgrIPF+41zCjyQCic4lbuaPSFh2fYnWvkRgeYt9CINNorYsYkdB
HIG+9ss3V3SSoMgZDBQEE3h8Zt6r8xUJo8XyapsgAkyDA5tz3pCwv93NMxtA
wOo2pdRkgIR5ZXa1ji+BqrDmFzETJKIMr36udyeQrKFz2W76d29G0M8jrgSC
hzMcFUkSAW4q/qIuBFStAwZrv5K47/fkuKsTAV7urtzEeRKOUUclVjgQGCoz
OOKxSML1fJJk/l8E2Edytmr+InE+RvKU3W4CKQqCc8vLJASP9JguWBH4H0gX
d+s=
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#2"]& ], 
    TagBox[
     {RGBColor[0, 0, 1], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVVGc81QsDVtybVPib2ckoo5Pr5pD12MSJDkqEXKsiVEooSaEUik5GWVmJ
7D2SZJRQyLVnchH+dMOVyuv98Pyez88Ud/Qyd9nMxMT0aAP/Z2OXqY666aOa
NzctRvs9twSLruC5EYEDWH/w7oxXtCX+TI5LKBbQw9NaY4FMX0uYtPUUJQpY
Ymv6U51OfUtYOa+8iBZwxkz5eOD2TxYY2n9LN0LgImSuUX1pey3ASN5ZJ3Xi
Esyt0gymCAsYc2Sr1yb4YPWUVkf5mjnK5luVSDE/nLU1sCTfmyMyl2svXfoa
PPqY6Ln+5kjhNQu9IxCM9itKeb29dFxpi204VhUMEmb8To10HAsZ3SxxIgRH
9ph68RfRsX3pfGB1QihqVaw71+7S4fuR4TsrFoZX063yGdp00B/0u5lKR6JC
/6w6R8kRBLEW7wkTYGCKnvgzytkMsoZRt4vOMeCxrf/ngK4ZOkO8pgeaGVhk
i6WdkTCDJLN8DsXnIQo+GiQafzJF0890ysfOGHzqiqZyuZpi+78xSmLh8Wiw
lY7nvHwYcUP+uqW/kvAjtOykYr0JtIWt04ePJuOWfal4RJYJZmyUf2PNTUZ8
nE0G3z0TiM76OeXuS8FEf/1Ssa0JvM5QpZgGUyAgbaPG9t0YXM75WRkqqUil
xfbXqhvjuPWT/IXFdIzK8jKcRoxADliypQpmwHOreYtQqxFC7VldLHQzQO8d
8V+sMEKxk5dgGSMDgosSmt+ijcDuqRFyRTkTa4pzYkeMjNB4o8/m94CnqDni
aSdbYQjF58TvwqzZKEx/aNKdbYC3lEaHNoVs9DIiDrx5ZACHQt/qa9bZuJK0
NbjjjgEiykbOjWVnQ0yBJUHa3QBTdbkDmYdzwM+mU9Qtb4Ck7kOFfzx4jqFA
R2HxEn1sWw+yNRTNR+HZFXV6hx7izljZyxvk45vDtCFzgx4kP8o7EJ75OMt8
k2gq04Pmsx6ngRf5yO8dcAxI0EPls5jSckoB2i77J+Wf1sOcxWyNSXsBnOku
F1hZ9HA0K/ad944ilPY5P3mrr4sEuxs0V/EiXJLLnjyhqotxLo+240pFqBKs
Clqj6MIrQOe9hl0Rkj/TbB136iLMfK5zS24RFM/5SLZ+0cGLHzr9j2nF2CM0
VvyDoQMp+vzU67sl4C2MN9Je0IZgoKHqqSclMDlwQ/zaZ21w5qbcZSsvgYh1
aE5zvza+s1pQ6OMloFQG+l1v1MaHunLvYZXSjR5Ox4Q+1oa/QtCv/z6Xgv12
8qHOQ9po5+DhoWiVwysrt2wwTwuvNTxcOo6WI0rT+/utDC1UuDeVXXQvR9Lz
qQCNBC2kNftaV8eUozoxnak5TAu+gUNJRnPl6JfPkDZ30cJuMlPG6VEFXtVH
fFUV1cLl9oOI+1aJElubgP59ADul9PNb1iq0cj89+VgUSI9QCF8TrgKaV8bd
OIAO2p4+e/0qOHTW+SotakLuHY+3VEwV3ItPJfmUamKkmcwsUq7Gkfir6rFa
mjB8lcHe5l8D/xdrD1odNMBfQgxv2vQScb5LNw81qsE2fM1ygu8lolZ4k8Qq
1JDq/Pldk/xLhPBbLWzKUYM8b2Xl3eMvYcnntTpzXw06Pg4PeQteImJhfwvN
Tg3UBKtMxv66jZ2uy2r/p4q/7DPtsobrIHZLZ/E6VRVl47qt79XrURp5C/QO
FfxtFLQYYViPTtG3Xl+bVLCcV8tHM6+H7rzszoQaFShfOfjXm1P14A349ylr
lgoquPcv10XXY8W4NkM/UAWVeoK7iqbq8VNo1sFgvwpqMhe8GYzX0KeQqrse
KqPeLVHIerYBxhG7dpj5UFE87e5uv9wAau6VqeOeVKSfVq12YmqEqcsW/dOu
VIS49th48jRi6ecXt9hjVBg5Eo9uqjci7J5mK1WZivfWITvzwhvhNlLHeWlV
CYNGHrzM+5rgyb4Y2BKkhCVpDY5cj2bY6ETuako5gM2rOZxCUi0ITQ74dWtR
EdR1hzKLvFa8KVz2S477A63BB9qbZN+jZS3y6th5BYjd/7L6bP4D5KKuu5wp
puDLyshwt0InFqZY57ZMyiN9prqFPbILMRJiP1Jl5MAu8pScmPgIUUfrB2/S
ZRB+ptQ1RfJv9E+zjhH2e7GslTNxOrAHMtPjvjSRPbio2FDC97oXchShJIUd
0hjJa7jbId6PuK8d7VJMUlhz+/hYxn8AM6B8/cAuCTntLb9cqwchnDi58bIS
MEwVP1XIP4zg0pqTMud247LIyd1jtiMoEilOnqsRh/zs0W09h0fx36gms5mA
OLg5Flgbvo/iH5mdz1e8d4Hl5eN8iYQxqOx74jzfLAb3A9V9W03GQWWrjGRN
F0XOIceU2uVxqFQc+3tznAiWvghyTsZ8wp9P+ZM5c4VRamHIzNCbwNQ731S1
F0IINApJLpmcwI+SRA+lGUGIDxn3poV8hh1PjIqrtCBUlYMZAoqT8EmzP2Fh
JwBb0ferm7omQXmWqRGbsxNLLRHj54L+wcC7prs93/gRx2J/0VxyCon9zvJd
5vzgHi5N3dwwhTb2HelXq/nwW0CQ5cHz00i9byggIcYHpmf2+eNcMxBqEuVk
Y/DiwqP8CK76GVit+2iqsvDihdxtVsLzCwaPs7RGBfCAML4gf4N9Ftw3Li2e
+o0Hg/s+dPnVzEK2e7Xw+W1uFCRBcsJ+Dqs/bFR3CXDDLp9KiPyaQ2/E8R7e
Ai5cD3kSNJg9j5kSGl+ODhdyp0Qes9NIcA/1N8WMEridGOIwbEqizBSHp4cJ
OJvPSeXRSSwvBwurDxEQrqkpOGxFgmOP28PxPgJ37p1oCnckEcYcP6/YReA0
NX6RzY/EhYlgylgjgd3BPEZbMsmN/CpcvHMI/FK5uqMni0T4tR8T7c8I9M19
6szMIXEoeiRbNovAfatiO4NCEjxqamyf0gmsy5p7h9SQEK6w4bVKIjDYcS+J
uYvEsQ/Dh2lRBCpCV5y6uknIK3dQnt8jwFA7KZPWS6KJx/jo9kgCJhmUEp1h
EqFzYRfa7hCo8m17GzRNQuuU/ax5MIGYfUr36LMk1nJTVkpuEDg/nmApTm7o
vTBiyB9EYC/t7EjdNxL9/+woHwwgwLLpY/r9FRLs9gPduEpgtFTNzeH7hh8j
DgfT/AnUuKXtV/hJQto2Yux3PwKxYtuW1tdJ6LeatbtdJvA/+v6BYA==
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#3"]& ], 
    TagBox[
     {RGBColor[0, 1, 0], AbsoluteThickness[1.6], Opacity[1.], Dashing[None], 
      LineBox[CompressedData["
1:eJwVxXk8lAkDB3BJmyg8KM241jHkmqgwUX7OJiG5thSTig4dypGjYi2jFKt3
m1CInEmMY6bCFJVYFIlCNIQswsOWRMX7vn98P1+1gwGufqIiIiIh//P/d/iN
tNWOeliMFM1eWfmOCTEb6uk+yiYEePZfjmlkYmNmanoFxRYVPwzUiEomHF52
lmdQ3HHR9fm+gFQmdvvOPvqL4guTSynUjXuYeL/+ok0iJRhmSstel/VsQ9Zq
57jLlFjU3A+2CJq2w7mXKXW/VcWCTX3gfWLQDr+x+0U19rFxfnQdLfSNHVbO
nImqTo9D5L0R7SeVdgjr4ISNq8ZjMn/+pEaMHVyuvfPfqfUnJJ9uiV9GsUO0
eIV2PIUDPdIoLGiXLXSZ/7lUfpoDjXNvdXxtbPGaHTDa08DByjj+fj8TW2gu
1S+in72OY5kBPalKtqj/mUvveJ0M/9CxqtpRG6z8nGysmnAD9owB5/k4G6S+
j7DhL9zC7Ys3BO4vrGGl5Jkr9MiE8ahTn8oTa4ztNV0mXpyJ+oVS2le+NVTG
ww8VG2QhvP/yk+ZMawQcM6GJ9GZB0HjokUywNWR9uXfyGNlouZt2b1LFGns8
b3OnpnMRbzR05vgFK5A97hLZ1DywxHq6PIKtEMcS93OzycPa7q8x9setUHEo
gHqfk4e9thqnHD2tIHVqK/ucaT68pYxW1Rhb4fkf3Xt/uVAAe8niAocpS2y4
R/yiJH4Xmm6i9MLjlpBcjPZiqnBRsIZfOjVtgdRju1n627jImnCUnhq0gGaH
vg9xiguua1L39zcWsCjsPNTziAuepUkAo8oClYXJ/Af0Uvxs3tVlE2OBCbdx
gUNLKdq91f++rmABjzspzUGryuFCrWgSbN8KmsvkyLMrPLhCtkOxyRzUKKbZ
kds8KGnJki2PzCFTnHVF4gEP6XbzpQll5pgXd6O7DPAw99XrX+0b5nhV+yBI
yOCjSCdhsueoOSIMoxe+feRD0VfWiC5hjhZpeXm65QP0rI6N1HA3Q2jLZqR+
qYSQdV2rVXQzpOj8j43iVTjCqzwaOsdAbqJhwnelKqw4cfGc1hQDbY7a3Sy7
KjwvOziT+p4BvWb5IFpyFR4HOp5oechAXwOZX25aDTO2uGLIaQaYT/KkXkYI
YNZ22iB6yBQKPEK4ZEkN/h5krTAUmuD+gM2L1i1Pkcbuymoe24Sn/hmKnuN1
8FE9udfeYSNmtLZKF59sQODya684P4wgOlcko0hrAr2ezFavNYTJos99t5IX
UPDLClIOX48XsZta6nVbMW3OLRflGUD16qe5wslXuLv67PDsDz18mu0TvjF8
DecPIyrGibrIHatukvqzHfu4D7dH03QgpVxADg114MzHgEQMayPhGP9wluZb
sOIZ7NK3WvhqWTR0NKoT0dWbdhGNNARvqOOtedaFap+0pPkPmugrqbvSpvYO
6ZXx1a7fNfDdvyNNJ6IHpWmyy67qakDPavnC4epeHLj0zKQ/RB3MbLUjZQpC
uFPHF0pr1BCqvF/9g1cflDLOH6CrqUF/3EOy06kfbZ6PH6ck/Qo56Snxuvl+
aOkVFBvPqkKsJo2rkf4Bvdeb5xY/q+D4puruFQ4DCD284dboJ2UU2R/Mevx1
ABrFmXIbvylh5hNVZjh5EA2Tv4qnqCmB78ZcyrEdgtPNiJG/7BQRtZ2dyRse
wrQwUcUohgq19zu6ctgfQXDK5de1UGBmGsuhbBiGyBZtAUeOAi+V1rkl7cOw
VinQVzi1FjNNiQOno//B+bkIHuuZAlLFWMGumiNw+/LySbyeAuSE/GzRuhHs
c/dy3ZmxBssuRLtvPjOKe/qtOf+uWAORQhZ3QHYMZ79F6kXGrUbgTW6i7NMx
iJwMOWu2II9HepfEiVOfsPOfOV3nSHkQOwL1/5Aax5rIkqlBSXn0GrxqDxeM
43bmT+d3HDmU3oLmEGsCnhz7Wh1dOXhzTQjlhQmsPXpAPOq5LH5n347uvTsJ
S4Yp7aSXLIpHlNOkHEmsMnAekl4gcCmD7SPcSSJ5l45T4A8Cvq4TtBIXElHq
w9865gkoCQSlTrtJUN0b1qbNErictK8+4SCJnmy/MN1pAkdNbkxLhJN4GyJT
6DFEQD1WfvvyfBJh71QOvG8isMA4v6rzDom+zsDl2xsJdE8Mvs4vIrHalTtZ
3kDg6u4K721lJIT1OVHxdQQWdV2D2AISrRJNhZsfE+htS7q1tJ2ExNtkIr+M
wMO42UPtb0jETG2ukS8lwDHfr5PTReLHcOu9mBICDnl0nrWQRG5hCvNgEYGq
sJeN0aMklrImw9XzCCQbGCe5jJMwkP6ddi2HwJmBdHc1kgTSr8iIZRNY53ii
r/YLiWrNjtzhWwTElnTkXp0l4Ub57LQng0A/39zfZ56Efq+zSVMaAYF/znrD
nyR4IT3eW24SSFGVnFlcJCH1Oa6uJJXAfwHuv2rW
       "]]},
     Annotation[#, "Charting`Private`Tag$17084#4"]& ]}, {}},
  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, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  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.8737302421940727`*^9, 3.8737302865099287`*^9}},
 CellLabel->"Out[38]=",ExpressionUUID->"faad53b7-2ecc-458c-bd43-c2e2528af9c2"]
}, Open  ]]
},
WindowSize->{1389.75, 768.75},
WindowMargins->{{0, Automatic}, {0, Automatic}},
FrontEndVersion->"13.0 for Linux x86 (64-bit) (December 2, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"a31d9556-d2b9-4444-bce5-a086ed60159c"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 1110, 20, 71, "Input",ExpressionUUID->"4a74809f-4d6f-4bf5-84df-1c81e845ccce"],
Cell[1671, 42, 1195, 33, 51, "Input",ExpressionUUID->"4abccf8a-d47b-424b-acb1-36e9461ec53e"],
Cell[2869, 77, 586, 12, 50, "Input",ExpressionUUID->"ddda9570-52bb-4f7f-8d3d-86bdfd8586a4"],
Cell[3458, 91, 654, 14, 50, "Input",ExpressionUUID->"86a69de3-cbb6-4ee5-8cee-9244eebdf060"],
Cell[4115, 107, 1856, 55, 49, "Input",ExpressionUUID->"f9b93666-a5b0-48b1-94d6-6d1a85496d1b"],
Cell[CellGroupData[{
Cell[5996, 166, 950, 24, 29, "Input",ExpressionUUID->"aa0fae4a-3cd8-46d8-ac13-32fe6f27924f"],
Cell[6949, 192, 23640, 414, 230, "Output",ExpressionUUID->"3524f72a-730c-44fb-9dc6-30e13e9dc263"]
}, Open  ]],
Cell[30604, 609, 543, 15, 45, "Input",ExpressionUUID->"1d64fb32-c3b8-4ea4-8f06-4b01e68febc7"],
Cell[CellGroupData[{
Cell[31172, 628, 787, 22, 29, "Input",ExpressionUUID->"5bc7a9f5-c5f5-4792-852c-3f1f56196559"],
Cell[31962, 652, 12950, 238, 227, "Output",ExpressionUUID->"faad53b7-2ecc-458c-bd43-c2e2528af9c2"]
}, Open  ]]
}
]
*)