(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[    109212,       2301]
NotebookOptionsPosition[    105123,       2165]
NotebookOutlinePosition[    105489,       2181]
CellTagsIndexPosition[    105446,       2178]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"(*", " ", "\[IndentingNewLine]", 
  RowBox[{"Zabawa", " ", "z", " ", 
   RowBox[{"tzw", ".", " ", "problemem"}], " ", "Keplera", " ", "dotyczacym", 
   " ", "ruch", " ", "w", " ", "polu", " ", "sily", " ", "centralnej"}], 
  "\[IndentingNewLine]", "*)"}]], "Input",
 CellChangeTimes->{{3.498647034065316*^9, 3.498647073544476*^9}}],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{{3.498647079755879*^9, 3.498647091854813*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Zaczynamy", " ", "od", " ", "parametrow", " ", "takich"}], ",", 
   " ", 
   RowBox[{
   "jak", " ", "dla", " ", "ruchu", " ", "Ziemi", " ", "wokol", " ", "Slonca",
     " ", "po", " ", "orbicie", " ", "kolowej", " ", "po", " ", "to"}], ",", 
   " ", 
   RowBox[{
    RowBox[{"by", " ", "okreslic", " ", "iloczyn", " ", "kappa"}], "=", 
    RowBox[{
     RowBox[{
      RowBox[{
      "G", "*", "M", "*", "m", " ", "w", " ", "rozsadnych", " ", 
       "jednostkach", "\[IndentingNewLine]", "G"}], " ", "-", " ", 
      RowBox[{"stala", " ", "grawitacji", "\[IndentingNewLine]", "M"}], " ", 
      "-", " ", 
      RowBox[{"masa", " ", "Slonca", "\[IndentingNewLine]", "m"}], " ", "-", 
      " ", 
      RowBox[{
      "masa", " ", "Ziemi", " ", "\[IndentingNewLine]", "oraz", " ", 
       "warunku", " ", "ze", " ", "sila", " ", "dosrodkowa", " ", "rowna", 
       " ", "sie", " ", "sile", " ", "grawitacyjnej", "\[IndentingNewLine]", 
       "    ", "Fg"}]}], "=", " ", 
     RowBox[{
      RowBox[{"G", "*", "M", "*", 
       RowBox[{"m", "/", 
        RowBox[{"r", "^", "2"}]}]}], " ", "=", " ", 
      RowBox[{
       RowBox[{
        RowBox[{"kappa", "/", 
         RowBox[{"r", "^", "2"}]}], "\[IndentingNewLine]", "    ", "Fd"}], 
       " ", "=", " ", 
       RowBox[{
        RowBox[{"m", "*", 
         RowBox[{"omega", "^", "2"}], "*", "r", "\[IndentingNewLine]", 
         "czyli", "\[IndentingNewLine]", 
         RowBox[{"kappa", "/", 
          RowBox[{"r", "^", "2"}]}]}], "=", " ", 
        RowBox[{
         RowBox[{"m", "*", 
          RowBox[{"omega", "^", "2"}], "*", "r", "\[IndentingNewLine]", 
          "      ", "wiec", "\[IndentingNewLine]", "      ", "kappa"}], "=", 
         " ", 
         RowBox[{
          RowBox[{"m", "*", 
           RowBox[{"omega", "^", "2"}], "*", 
           RowBox[{"r", "^", "3"}], " ", "\[IndentingNewLine]", "    ", 
           "gdzie", " ", "omega"}], " ", "=", " ", 
          RowBox[{
           RowBox[{
            RowBox[{"4", "*", 
             RowBox[{"pi", "/", "T"}]}], "  ", "-", " ", 
            RowBox[{
            "predkosc", " ", "katowa", "\[IndentingNewLine]", "             ",
              "r"}]}], " ", "=", "odleglosc"}]}]}]}]}]}]}]}]}], "   ", 
  "\[IndentingNewLine]", " ", "*)"}]], "Input",
 CellChangeTimes->{{3.498647274580111*^9, 3.498647292747118*^9}, {
  3.498647330491349*^9, 3.498647362041929*^9}, {3.498647546250313*^9, 
  3.49864766995125*^9}, {3.4989776880169353`*^9, 3.4989779075059357`*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.498647373117508*^9, 3.498647386406022*^9}}],

Cell[BoxData[
 RowBox[{"\n", 
  RowBox[{
   RowBox[{
    RowBox[{"dist", "=", "149.597870691"}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"w", " ", "milionach", " ", "kilometrow"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"omega", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.", "*", "Pi"}], ")"}], "/", "365.0"}]}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"w", " ", "latach"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"kappa", "=", 
     RowBox[{"omega", "*", "omega", "*", "dist", "*", "dist", "*", "dist"}]}],
     ";"}], "\n", 
   RowBox[{
    RowBox[{"m", "=", "1"}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"w", " ", "masach", " ", "Ziemi"}], " ", "*)"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.498647418772759*^9, 3.498647498230155*^9}, {
  3.498647640598517*^9, 3.498647643883833*^9}, {3.49865202825113*^9, 
  3.498652046129769*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.498649068888134*^9, 3.498649070379031*^9}, {
  3.498649119000066*^9, 3.498649119205114*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{
     "parametry", " ", "krzywej", " ", "stozkowej", " ", "wyrazone", " ", 
      "przez", " ", "energie", " ", "calkowita", "        ", "Etot"}], "-", 
     RowBox[{
     "energia", " ", "ca\[LSlash]kowita", " ", "oraz", " ", "moment", " ", 
      "pedu", " ", "Ltot", "  ", "\[IndentingNewLine]", " ", "rozwiazaniem", 
      " ", "r\[OAcute]wna\:0144", " ", "ruchu", " ", "jest", " ", "krzywa", 
      "\[IndentingNewLine]", "   ", "r", 
      RowBox[{"(", "phi", ")"}]}]}], " ", "=", " ", 
    RowBox[{"p", "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"1", "+", 
        RowBox[{"e", "*", "cos", 
         RowBox[{"(", "phi", ")"}], "   ", "gdzie", "\[IndentingNewLine]", 
         "    ", "p"}]}], " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"L", "^", "2"}], "/", 
          RowBox[{"(", 
           RowBox[{"m", "*", "kappa"}], ")"}]}], "\[IndentingNewLine]", 
         "     ", "e"}], "  ", "=", " ", 
        RowBox[{"sqrt", " ", 
         RowBox[{"(", 
          RowBox[{"1", " ", "+", " ", 
           RowBox[{"2", "*", "E", "*", 
            RowBox[{
             RowBox[{"L", "^", "2"}], "/", 
             RowBox[{"(", 
              RowBox[{"m", "*", 
               RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], ")"}]}]}]}], 
      "\[IndentingNewLine]"}]}]}], "*)"}], " ", "\[IndentingNewLine]", 
  "       ", "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.498649134278972*^9, 3.498649194907016*^9}, {
  3.4989779413149357`*^9, 3.4989779526009355`*^9}, {3.4989780329679356`*^9, 
  3.4989780892319355`*^9}, {3.498978618315935*^9, 3.4989786838079357`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Ltot", "=", 
   RowBox[{"m", "*", "dist", "*", "omega", "*", "dist"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.498652110081*^9, 3.498652170699556*^9}, 
   3.498744087668823*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "kappa"}], "/", "dist"}], " ", "+", 
    RowBox[{"m", "*", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"dist", "*", "omega"}], ")"}], "^", "2"}], "/", "2"}]}]}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"p", "=", 
   RowBox[{
    RowBox[{"Ltot", "^", "2"}], "/", 
    RowBox[{"(", 
     RowBox[{"m", "*", "kappa"}], ")"}]}]}], " ", ";"}]], "Input",
 CellChangeTimes->{{3.498649201277267*^9, 3.498649248023225*^9}, {
  3.49864953412665*^9, 3.498649535304441*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"e", "=", 
   RowBox[{"Sqrt", "[", 
    RowBox[{"1", "+", 
     RowBox[{"2", "*", "Etot", "*", 
      RowBox[{
       RowBox[{"Ltot", "^", "2"}], "/", 
       RowBox[{"(", 
        RowBox[{"m", "*", 
         RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
  3.498649333396809*^9, 3.49864934323158*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"rho", "[", "fi_", "]"}], ":=", 
  RowBox[{"p", "/", 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     RowBox[{"e", "*", 
      RowBox[{"Cos", "[", "fi", "]"}]}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.498649594720883*^9, 3.49864962916204*^9}}],

Cell[CellGroupData[{

Cell[BoxData["Ltot"], "Input",
 CellChangeTimes->{{3.498744118332438*^9, 3.498744151224794*^9}}],

Cell[BoxData["385.2457248300915`"], "Output",
 CellChangeTimes->{3.498744165915021*^9, 3.498744986179862*^9, 
  3.498746068598442*^9, 3.498746342998448*^9, 3.498746673910418*^9, 
  3.4989781272799354`*^9, 3.4989790019789352`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["Etot"], "Input",
 CellChangeTimes->{{3.498744131577819*^9, 3.498744148176283*^9}}],

Cell[BoxData[
 RowBox[{"-", "3.315849695761892`"}]], "Output",
 CellChangeTimes->{3.498744166006375*^9, 3.498744986423355*^9, 
  3.49874606868533*^9, 3.498746343099041*^9, 3.498746674095195*^9, 
  3.498978127294935*^9, 3.4989790019959354`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["p"], "Input",
 CellChangeTimes->{{3.498744137317296*^9, 3.498744143569849*^9}}],

Cell[BoxData["149.59787069099997`"], "Output",
 CellChangeTimes->{3.498744166158971*^9, 3.498744986620608*^9, 
  3.498746068959605*^9, 3.498746343456615*^9, 3.498746674438343*^9, 
  3.4989781273139353`*^9, 3.4989790020129356`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData["e"], "Input",
 CellChangeTimes->{3.498652277182697*^9}],

Cell[BoxData["1.0536712127723509`*^-8"], "Output",
 CellChangeTimes->{3.498652278667227*^9, 3.498744166313408*^9, 
  3.498744986787569*^9, 3.498746069196193*^9, 3.498746343660403*^9, 
  3.498746674723174*^9, 3.4989781273309355`*^9, 3.4989790020299354`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
   "Dla", " ", "znikajacego", " ", "mimosrodu", " ", "e", " ", "mamy", " ", 
    RowBox[{"okrag", ".", " ", "Przy"}], " ", "ustalonym", " ", "Ltot"}], ",",
    " ", 
   RowBox[{
    RowBox[{
    "Etot", " ", "dla", " ", "okregu", " ", "jest", " ", "minimalna", " ", 
     "dopuszczalna", " ", 
     RowBox[{"wartoscia", ".", " ", "Inne"}], " ", "wartosci", " ", "Etot", 
     " ", "moga", " ", "byc", " ", "tylko", " ", 
     RowBox[{"wieksze", " ", "!"}], "   ", "\[IndentingNewLine]", "Ta", " ", 
     "minimalna", " ", "wartosc", " ", "energii", " ", "moze", " ", "byc", 
     " ", "znaleziona", " ", "ze", " ", "wzoru", "\[IndentingNewLine]", "e"}],
     "=", 
    RowBox[{
     RowBox[{"0", " ", "\[Rule]", " ", 
      RowBox[{"1", "+", 
       RowBox[{"2", "*", "E_tot", "*", 
        RowBox[{
         RowBox[{"L_tot", "^", "2"}], "/", 
         RowBox[{"(", 
          RowBox[{"m", "*", 
           RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}]}], "=", 
     RowBox[{
      RowBox[{"0", " ", "\[Rule]", " ", "E_min"}], " ", "=", " ", 
      RowBox[{"m", "*", 
       RowBox[{
        RowBox[{"kappa", "^", "2"}], "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "L_tot"}], ")"}]}], "\[IndentingNewLine]", 
       RowBox[{"(", 
        RowBox[{
        "pamietaja\[CAcute]", " ", "iz", " ", "energia", " ", "jest", " ", 
         RowBox[{"ujemna", " ", "!"}]}], ")"}]}]}]}]}]}], 
  "\[IndentingNewLine]", "*)"}]], "Input",
 CellChangeTimes->{{3.498744203848306*^9, 3.498744251553463*^9}, {
  3.498744282228694*^9, 3.498744400653491*^9}, {3.498744433517811*^9, 
  3.498744438210105*^9}, {3.498744500616396*^9, 3.49874452544253*^9}, {
  3.4989787859419355`*^9, 3.4989789293929358`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", "0", ",", 
     RowBox[{"2", "*", "Pi"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw1m3c41f/7x2WUIpvsPbJKNhl39t7rmOe87RYtDQ0UUillJBEZZaSyUhK3
QkXSlpRRUpTMltX39bmu388/rsd1OO/X+/V+3ff9eF7XOTJUlFsoMxMT024W
Jqb/fj8YimiSqaGA6f9+BPI/75UifJEvbqXGrhemj0P8i4UJi99t7VDf9dFU
7f1T89WEn9q68zdf+2bK02HqtVhNQTNTb7f8rmnTTq6avq+Ef0ZkZH/m/m2q
d0X9VTdh6fG41oJr86Z1V+tMqwmbH5LNKZdnAmUrN8k0wmyuFsNMu5hhppsj
OpIw+8TcTbcWVpi+/k3bhLDB1tzofO4VkBSxsJWLsJFwXPpIwEo4lmcr/K6K
AveCpRSlaxxgsGpQ9zLhMM2drYy51RCW1PUomPCs80WanTwP+GeOWfXdpICu
mqm3YycfdNA3aqURbvZO5Tp8gR8mrimkA2EN+SjpJBSAqVDB6+k3KHD83fg6
iWsNRI1YROoSZus+N3pYRxhGuKjsV9cp+MKhKLLDXwRCKtrcWAgXF4ZfMysX
g6jVpv3plRSstDoZJ/NcHKgZjn4pwpms3q5/f0tAsBvfRpVrFPD1ToU8lpUG
/zj17ezlFHSd4c3XjJaForDz+bvLKNBsCONOYJGD0AJPq/elFPB2KQw8yZKD
E3tDZvOuUuB2IW61xz15yIuR6hkrpoBZ2aJBfpUSvFzcW69JWN35oKF9nhJ0
hmzliimioM3u4PwWjbXAmzEkPXGZAuWRizcveimDbabNh/JL5Pl3yi/eKFQF
4UMP9j7No+CqesCFYh01aN8fcWI8l4LEVwI+5x6pgWsLp4jkRfK8mBcqfX6o
g/sNpo825ylY9y46sl1cA9z3Wvi0naag/vjWQYNwTXAwsmlMTqWghfEhPe+W
JlwY6nhucYoCjlrO8Xk2LWhpa+a7kUKBy/1endJiLWiuHfplfoyCJ6u+lTQO
aUObX1dB1D4KLmnZq7b46kFli5/9hxgKdtNZc1rL9OBp9XCl5R4K3tSvsmj5
owf2gg25y3eS881dOFGSpQ+ZB86ts91CgcqdQTXRFwZwnj6a9tyfAn7Vs4k5
1kZQM83/c8yXvP+a/UJKWUawK1JP5J8PBXu3qx2oGDaC8+Z5RWs8KchOkqVf
jDeGH+FXnq10pOC8D5em4F0TKGO7yjWzkYJvCivK12UB3L6THfdBkJy//YMl
z1nM4RhzXXIsPwU73fZHTyiZg+DV15yCvBQI3VbYze5gDn6CJ9INOSmI2HKb
WS3DHGY+PypyZibn6wvLJn15C0hbyvIrHGeA2UudiABzS3iuNXT/BzLA0Xrw
1cVD1iDNLfF48R4DVtXtiuS+bA0pq29Q7HcZIHY5ReRQqzV48nS48dQxgJv6
UmrFYQM/VDs2fitlgM3jsNSz2TZg47jU13WGAazizg1LVbbw+dPqL/a+DFAX
ee//d8Ae1OcSNw96MYBnYWHX1DIHOHKWoz3anQFpCWmvPso5gNelk+yJDgxI
kUpyrQp3AFZG5YtYEwa85rP0XvjhAI8jVAyKZBiwd7/VQfMFR2Cx2h09MEwH
i4Ij366zu8DBKtvLk4N0iLr8oFFPxAWGjHa0L76nw+1W/epGZRdIfxU3xfGa
Dg97Bv822LnANbmgz79b6aChkZh6/JQLTLSoSnAX0+FxvZRHCLcrPMpZergp
iA57MsKXsfC5weGBu2t2+dKhZzRJK0TODaqZ73Tke9Lh90M2G9R2g/Bnn4O+
2dMheCyVNdLbDUTLY7bZ69OB4wGH+cFcNxg8IO6owUOHhDVVGSUK7lDbIG+y
bn8QNAdueD6k5QHDThkKUduDQOp+xkSmhQec0z2//mpwEETbpXtZeHrAj/VH
6xccg4Cv6uy3MzEeUNVQ+91YNgi0lptcH7rtAVvVKgrFHgdCuMqK+C1GnrD3
S2KUL08gSCdcWfnHxAvOR3yQcmELhPdXLg32OXrBTy/baJO5AJCenjJq8PeC
kFpYYB0OgHtrufeGHfACf2Wz1Yb1AfCAV/vkwTovoC2JHbbyD4DK4zLsuSre
cP8Ar+NCgT8UnBWyPMPjA1PK3uEcmf6gOxj1tFHSB0QL1x7kT/GHZ9unNg+r
+cBzO63LXDv84adiUreMrQ8Mu15RrQJ/WBr8+dkj3ge4bbI+lPX7wWDmfYsL
kz5QfWL1t8OCftDsJPQvoYMG7Zd5cwtW+EHH0Z17OHppsHxXfVzDX1/QuNwY
c+oLDdyvz86++eAL+d1BK/ex+kLk95Vz1SW+8Oi86dyssS+Y1KWxzmv5QsaB
PVsTbvhCffvfxS92NBDZvijUfcoP/rov8Qkb0UDm0bI5eo4fsPXqD4I6DQxF
XEfGrvrBaZudtQe5abBfXqtz5L4fSFgFNOa99IFthXq0w3/84K3OGoEUXx+4
s22PkjLlD6yKZy/JBXuD1JYsXvf1ASCwXuBQo7s3nIy7Y1i8MQB6JH/3O1l4
Q/50l9m4dQAkNK8woSt4w/XPsh/D6AFwvf/zkM4XL/DoCV76nRYA+mrzJ25G
esFU0Mzd05MB4H/+Ai/zZk9g0a022FsWCF4u1k9f0jzhYObI03u1gSD5g+Vi
vq0nnOdd/26+ORBMLlfFKCh7QmHUIavAN4HQ4bSjbuyLB2idEpdqZA6C7/s3
2wWFeIDX6yLpXloQeLiZ/+ShucOLkHvVLUx0iLjBa+9i6w7FI4vT8Wx06CqC
XycM3OHEx143k1V0EPrK6j8t4g5Lb7u3VvDT4f4jy2cp79zgh0/vqJciHT4O
WNzu93WDL/NXOaft6CBzR7N4o48rJAeb2m44S4dzJZzS49ausG7bd69dmeT/
WY9sy9FzhaAJOqPqAh1CQuZOfRVyBb2NSREyhXQoHpZ2dH/tAjOxm+SeV9NB
1UzYIczVBf7dKo/c/pIOFzhfpBqlOoF7vlv4BV4G5LYt7eyKcIJDXU2jZYIM
6CsKqvK0cIKvO3vP1IkwoGsnBDvMO0LOVJxZE+kz4bQLTS8jHIHxq04mcwMD
VngxdrubOYDkgkaXtwsDqjT9/OkTtlAynuKofIIBtc//3ph8bAuLIp8UclMZ
UJLlXrm/2BZ6Gzv+rTzLgNMuH7r30Gxh3fOg2jfnSd+0NRlY88AGZp+VsWmX
MKCXLSGRM9MaeDRv7llsIn2u/7MNi7YljHp6mbT9YEBAzLZXNjyWkKwj3Rg2
xQC9NxOKKd8tQM92rp1llgHSA74Wc8UWcCLbJl7jLwOsTexp+YIW4BC/PcaD
hYJfN27XLPw0A1bppJ42IQpOatF6p5+bQejd+XdSIhS0V/h7DFeagbZoKWOP
GAVBB1y23wk1A59Lnr/4pCn4G5yspv56E5xMFFQXVqZgSrC+3KwGII2J7mFn
SEFgazdt+jTAOGVq429EQe79vvG8zQDPRrNvRZpQQA28+fpVBiArpJpprxmZ
ez+KHC9Jm0LSYKKRnx0Fc8vz3R+/NYLWoMwRWxrxjFdR7ZVnjCDQ8YA1vx8F
9gqu8ietjMh+bhR7S+auidazUY3ajbBCIUPEg06BSP738ypnDEFiR06oYDgF
Pl37w0fM9WEyqrK/cRcFz/It/F7N68Hgwa5bOmTObw1dOt9Yowdin2bWlhMP
6Iv/rXFQTg/g4VRv4n4KvI5437/DrAv35k1C2Y+QOd72pT+2QQcaTqxhDYwj
npQ816W/UweCnzsa3oinYJKWElFEPKNcKX2LBfGOuAuNO4xRC7rM/r41IF7y
zvMR69u9WtB+RENl6wkK1DgNZret14I+LjpPzkmyXmcX9VOXNGHX9pScr8Rz
qKv3v6Uf3ABZxpGXXM5RkLd/NLo5ah087p9WccihIPZnkmd0vzpE5Flc0yMe
lRrKIy7qqA62CqebJIlnUdHRvT4qapAcpJ45QDzs3c34gaRhZVjr2cuqRryN
rnMmj9ldGaa7JdlnCa9oanXb17IW7qm35NQXkv0/Tqm7XVKCWL3gbRuIBx6+
+mz7A28F4HWwejBxhQKtK7FVfO3yUGKaaZRBPNITR5gCtOVBYGTUXJt4pgT3
kMsnHjkQfsG5EEY8dGBxOEvqiCywLTt37S/h+i3ndnuPy8DpuhNjycRb45gG
4hseS4OdigvX+QoKkm3mf33WkwYrp/cza4jnBibESIg5SEFcs29YBmG5Vk+q
iFUCcO721BHixb18TcwfBsXgVxF91QThwEtHlrjviUKfz4kbNOLV/FaaksF7
hOGu079LMsTDLRQD52v018AZuFqYRPj21hGPpXlBSBSMjBghHBJZdNSsWQDa
J5Qn//N6NqWowbgE4vmRzwyyCKe5KbTeseQD0fMWliOEOWWLksfZeWGUOZpv
A8kNP3Q/yZ6WXA23pHKmbhGenl/O31e8Cu6b7WKZJiyS3JQircoOhcdvv1Ei
uSSn+pl+UBUbDPf4bvEhbGG6/MV5PRa4e8ah5SjhpCLlFZ33mGD58fi3ZYSH
qlYLzCQvmDrxTtzsIHza/Y2wTcZv0yezJ+1GCE/9s43LKpg2jdjtUrJA+Gsn
l38z7ZvpfTvdlv9yWLx3AIfdziFT1Qr9SyKE1/iwe5uqPjWt3ORhJE2YW+/v
t+DP1c3sWskX/8t5SSPmTi7H3zbXiT+/K0G4he9jUIHqSPOogFaeAGE9GRON
lJWTzSfLb5guJxx5uOem8eJss9Fnu6Jpcn3nyxJO3yf/NntYrnrUS1ikamdx
+vBS84cHs5V3CT+Oklgr8n0Z5rzg88wmfDvsR8ynUFb86bf9XhRhjV+DnKUD
y5HzAd+XTYRndqjqh9FWYsknzhf/5TrleTZriZccuI5vy+E3ZH8Dq8qTnzhw
4W0x4+/ZhN2mVYoNgnlwwjJDzItwk1zRFvYCXhRzPsrJRdgnLXDTi/d8KLpd
shnJ87xtcuRJlogAerjFam0nvK5qba2nlyDelauKFCI8v3QknytdCNsT3jLu
kPNyj/21xDZOEfzVFpo9Qc7bjkaFf7y2ouhsydqTQHhia0fcjUQxzD/O9paH
cHB757KXSxJYfa9KRoKcV/9za2Vp+2SQ9ulUzy9SD+aB3/RThGTRy4BfnU54
lKPGqLZWFs85KVi1kvoxP9lDm5uUw6eSeQ8Oknr7fXsvj1mkIt7azDR1vYQC
73O+dRYrlHBH0r/wH6R+S233DUOJEioE8havJfwSpTQUhtbi2UKHramk3n1a
y4qzfVTR0O6MuXA+BQ87fdhefVPFEyEdLcYk142925TBfkQNl9Z/rwkk/SSk
buaaf4k6hg02zZ4j/ef9n3N916bW49fFvh+tJMfdvD004GWmgY4b4vjasyhI
uBZEe3lSA8cWP1TczyQ55el9UZTcgPfHWuF6OgVPm98OhlloYsv7XXNOZ8j8
eN5xvi1VE6UGeV6qkhwowShrl+rRRGP/eDlW0h8Vbx61aI/UQq47S2pXST+N
nXHr7D2tjaHPs8ZuJVLwKmlvIMdbXdTfttFz8QDJUZe+aYjI6GG+o2f3WdLf
s5d0A2Q26+GpVrt2GZILx2/Mx0ov6GGFxAY1bTIfctrlB3qkDNAT2T7KRlGg
HbS8VCxsIypJGNutpcjrbQq/D4yb4uvM+dHXFhTEzPPVD5J2dORE398bZL4V
maS6/bMAFFyp9+QYUKDQm3JwXTpgxfaBBkkyH6+rFlLu6zehSkrLc34tCjrj
eioqw8xwU4LAfhYybxkVb67Ep5qhmUbxUooERbz0rZNLrRlyHvdx5SDz+VaY
3/4Py8yx+vVqxhLJiXn8X5lu5ZrjxfsCJlc4KCgefyjn9dwCd9m0Zi3+ZEAi
K2tD928LzFDvsN4ww4CsyRf25pKWGN4s9iJokgFl2St0hbdYIlvqwJ7KMZLT
dphN7mW1wtq/33g4BhjQs5QWwKZrjc6i/1yftTPgb/36ExwZttiy21U75RwD
FH//0fx32xZpsX/fK5McaLgzoenbB1ssGKGXt50kPlO/5HFDyQ5nm9fX/zjG
gDvKa6mBBjssiQytWb6XAecmQfjdgD0+aXbndia5cU16rMGMnBOe3msQYiPB
gHd8n7InrJ3w1OqRM87E1w5ess/+tMUJIWyMyY343Ns2P+fqGid0fR1ab7Oa
Aem3Lh0fNXPGikV65uICHXYRN+5b74JGh/fxJ/XRIeyV3xczW1cMvHHRQZT4
Zurh1NljNFd0H4xe05dGBwcjt9f3I11xq7/kyPlTdHBT/7lM84QrdspIpLMd
o4OtXGzOhw5XpI1ctKnaSQel/FoqzN4N1x1w5rR3Jiyi5iNn647eFa8OnSG+
LHfvK+2tjzv6SFS9eEF8uodPv/Z4hDvybVqTK/8nCLbR9gb3JbvjZFnufo7P
QfCKn9/cpt0dHwin6hxqCgIHKddmIzMPjBF/1jAdFQS+DT6VMbqeeGppbZlu
VyAcHhuxfmPliQ/O1SWOtgSCd7r4h/XennjQUeR6+q1A+Gm5NPxsryc6uhtP
dOUHQubYgdHHtz1Rb4pKfbcjEGyPsN5PN/DCcPtjAt4CgXCjb91MmK43Zueb
l866BYA6/5XVi5beGPn0Xeh5qwDgjJU6d8rTG+MMC95tMAyALbWHHp7f7Y0b
PfC3o0wAsBTyRQZWe2OXimyw3IQ/RNUq1Zeq+aDAoR3iViQXrpJmofOL01B3
WCrhWL0feH/d4P5SmYaG/pt3rKzwA0VPa75UPRq2FiafTr7kBxsrL42Mu9Gw
/vbZwKhEP9jhfDklJIWGn8qLro27+cGTCeFO2Z80HDLUarz63RdmTkpKPXno
i6dBpk9JzBdsUyckm1/5Yq72pktvV/uCeFp4fMWQL7Zz8f5MYPKF1odnGNHz
vqjn/HVz62canPwyrVO+zg+5rYdfzt2gwaL93La+dD+Mir95VN+cBnaN6psl
af64TLBtkivUB6L8VDYahvnjMQ6/xQ5PH7A83/DScZc/sq24vnTEygfm3Xc7
+Kf6Y/rCpqlXSj4AJV+mBVv8UWx1Bp/qmDe0CyU++6MYgLYfODVObfMGm8RH
WsHfAzBG/9aZTVFe8PBK5rPoPwEo3GPyQz7IC95+S8E9rIHI4cXycpmzF7x5
/u9kqHgg1k5Pbipd5wWPaj2j39sHIvzJ0y8d94Q1e/Ki5csC0UXYtmwfyY2e
XkGf2QKCMMdt8++DQR6wZ7F2W2t4EL45vJ171MkDEtf1aB7YGYSiCh8POpt4
gF/FmfDO5CAM0izYzCHhAeItzUqs1UEodiKFH965A+tbN4MwZjpOfWq8WeDu
DnJlJ4VfuNLRP+GPIN3MDWp1ff0LveiYwsQ8I7XBDZI1FHm3+9HRyL6k8Z2U
G5w/EvDodwgdnywyKW9adIXmVFmbV3vpeCHDdfx5vSus0l0WcTOXjlarOYNX
qLlCAPdxi+8f6bibucF0mMcFZrVdMgdG6MjIvvnGmMkF9JciOp+O0VEzY0pQ
acAZ/FTufiqYpuPDGwH1YbnOcCObXs7PzMDgJ3slTwg5w4/AH02a0gyUKQ3P
jVjpBHyywzUuNAb+YrWtuDdqDxo8nl+b/BnYreHzL/mhPTw02Pt5LZ2BnO1J
GfYl9rDN74juZBgDmx4eK2qi24M0e5m4w27yfiYbjD177GBsgM9S/zQDjQWn
K+0f2IJh2/yF900MNKlKtM/OsoZ3j+0uM91n4Lj1nt2l0daw8eOvLKk2Bk66
H/9eZWcN/4722Dp1MpAmv+vnzSUr8FaUr9vew8DHFcY59DArYA+bH2r4wcA3
Y221Z0hu7FzDCE+dYiAHk99XAy5LENE62+U7S+7Xxk2k/4sFVP2TNBj5y0B3
I518gYsWkOD3JaaZlcKpsc5tG5dZwHn69aEDIhR6D1yu2f/EDOL5poXmxShs
9fz68VGJGdje/yQdI0lhw6ekZr4jZvAyfvGGvxyF80zLKrM2mMF5toyr39Qo
VOe+fpqetQnEhXeH8phSyOT1Vmp5EICng2xNdBCF25vCqp/0boSr76/PqzAo
FPBmPbY3bSOkV0VXDVIUqt14sF/MeiM8LxyiTMIojD9eKepcawg+jW+e399K
ocdw7Ybw0wagadWYNLGPwtwdQrc2b9KDFzGrd9IPUNgXOOh387cuFIg4tT6J
pfD9jTy7yUpdOGmknXv+MIWM634OQaK6MHd/ivvXUQrbxPg/zU9rw/LKj/de
p1L4e6PRLtMiTbg3qz65cJrC2GlvS02aJnA2WB+TSqNw+FZxkjS3JkRVtD/z
OkehXN7JgokDG2AZp9tMahaFDl8/6fi6aUCmcKRa2nkKtzK/TlRh14C8ELW2
1GwKT1tzV43HrofsfY0CsTkU2hSJFqxwXwfDRiEKkpcodEosEfnIpAY2n3Q9
5wl79jJPWm1VBRO9Dx9f5pO/Nz46U9yjAsXHngrtvUzhjal+QcfrymCfebUv
sZjCF2Z/X7r6KUFYzK5Z4xIKuVV8pY61K8Jbo2ORU4RNu8aab25QBEGWB+F2
Vym8OHQv4+8KBTgwJHaTeCbKbo/IcayVhTpjnm6PcgrbX40lm8vKgm+62d5B
wjl3dKR00mRgTdF7aqyCQjPb/jccW6XhUivrlYhrFK5U1j3wsUYKDl3/tuMj
4RXi96pKzkhCYmRC0aNKCv/9O7pfyFocItXn2XWuU+hak5bZKisG6/Nbf+US
vq1oob11SQSCNb9sDrpB1kOz4y2rWwNVL3iniHcjt2VMn8lZISiT1PjHfZNC
GeV50+6tgvBh36/q64TZ+55NDMrxQ8e5fT6/Ca9XZntI/eOFF1GFiRuryHl1
orkMvuOBny3BhrGEeeKPFtVs5YZPDpUxJLeh8Itw2Xy51cA6Em04TjhTy6Hs
6LtVsHP35USparLf05KtjLPssGpK08eR8KPP3UsGNsshMEy0Ooawj53LZ85/
LNA26Zh+kbBDyTmWvrplQGtqnSe5B6+BYKdt75LpT66933oIX3xKW/fGas50
k3YgNUH4O0+/jX/tT9MrZ3fTmWsoHCk+rflBZso0MLfuMw9hnc7q5bWSo6Zh
fRLTooSfFNU6OM8MmJ76Up5IchzWv+/sDnDrNI3c4ldAchyWfRVhe3m0ujlI
br2ZJOGFe9zOa769aS7+oLhFiPDGsncaOTnDzRsOWIitIvxl/u6Toubx5udj
J9z+kvVEVl9RcN850yyitcQz/N/6+WhN3ax/msMjirxIDkW198nZVucXmoeq
4uSuEU7PknPqy2DCAzZ5e1MIi1BvfLj7mfFjzDJ3inCNbMSgiRIbZly4VadL
eLD5s09k9ArslHyYv5ywqVeA7uk7K/FQkgnvC7L/+bT+79eYOTGYUl51gXDs
rf2TrfZcqCV09pwI4YtlJd4T23lRdLOa7EvyvM3mPMui0vlQVmWfwXHCpnP3
Ur7V8yNzWeS/YXJ+rEvG+l8xCeHaHbrXTxH+xtT1cJPCGlx+e9nn9YRfnS5U
KrMVRq2s8bEIch5v3YvWDjsnimceGDcskfMbceWs591bYlh78RtvGmHmnb9e
r+wTR3WDTt8Sct47jMbfn5WTwlW773vmkHq5ZOy4Xu6dDGo+61JlIvyhbuh6
+j5ZpG1f3kIn9ZZ16dT8kpAcPgytVeYrpdDO5Iltm7s8yjTa/NUh9Sqe1dPL
0aWIrrZ/u3aQ+q46PjfkuEUJR6LCg0uLyPn9Xe6WsnItPm2nP2crJP2S4+zd
b5bKSP/tjfGkf7R3fX1Ca1LFcx6xl9PzKFx9a8iQ8lfDIUV7y4Jc0v+HPCup
OTWsL2B5U0b60cD57nE73XV48Ua13GHSz4RDtZneeWngvbe1NaqkP2Y21j6k
d2niR/3Zj6MnyfXqp26LimthGr7tLTxBYaCjiW33Zi0snmCWZjtOoUvJlh9K
7NoYUmorfDKBQttV2gMZZjqo2DbAkxdD4cN7bUx2t/QwWb3r64o9pD5r9Y8V
serjwvFnp7ftItdv47H85aaPMalfrZSjST1ztSomT+ij5VO/3VQkheezo5Z5
rDXETqOVG977kvliY/kqKdsIQ3p2fhGgUbjP4YS47pARPso5JWnrTebF42KX
fmVjLIrP4Cl2J/UWLKUjctcYG9tP88g7UAiXphd0P5ig2kLaVwcjUn+58oeX
pQHO2RgMGRtS+OY53Uu8FtBgVJpS0afQ4uxBIa23gBu8by/91KKwNH1KxUN6
E2plnKryVqVw3crgyuCbmzCnKHGPJZmvp+8G0m92m2EZI3vDhBCFc0K+6Ykz
Zig7dfZwhgCFxtdVZLzWmOPgsZSu59wUPo1crzIaaI5jzr+d1ywn+9totvB6
3Bxjo+9ve0rm/XtJ5z/3V1niH509h6aIP1jeuutUCdb4ODIn/Bnxi4vvGVe7
w6xRcOmSdnkzA+cMFMy/nbJGWZOvB5zuMLBZTlaYq9ca76ptfEG7xsDtB/XZ
lu2wwYpTNYcPnmVg+I/j/ufybbFmhVT1H28GOW/7p4V+26N0rc0bHQ/iRy1q
GQ7iDpgrp7K4zYWBP4tGJA5scsCe4Gjvbhuyvo6CpsaTDmjGGPD3MmDgPbNV
qwakHDHUp+BKjAgDazsj9p62ckI7mrQf+xs6pqs9dL1q5YK3rQpPNT6n49LU
EhPl44KyzVb8W7vo+JRariS82QXzL9yXv9dKx4sp7x/EpLqgt2fAAaUaOmq9
8apfeOmCuesEWCTP0FHZ3txUie6Kh++Fhl40p6Oh48aRV7vd8GfnVEmjCR3X
GTxi/pnohpW7toT0GtAx40/kPPd5N0xuyKSzadCxJH9ym+4dN1Q04j6oJE5H
nrbYlxaLbji2aeSbypcgDLvfExd8zB2dxAJeWJgEYfpxjrrk0x6ofOiEu5J2
EE4eX5QfyPPADFMdB2aVICxvacrZUOmBqQXPKi8LBqGB9vDxh50eWJ8gvzfp
eyBGa+9mr2f3xIEnfi5j2YH4aXhH4OV4T3xcO/NNjeSDYQ+vnylRXhj4ZruJ
3FAAunhZcW487IXyb0x9+d4EoFPy6rqRU174pC/jz7vmALz973elerkXqq2z
EhNKD0CVjpB3+4e98KyavY+ufgDeTdHhYPfxJv3hx3LNQ/5o+hjiPhn6oPGM
60uLaH/sPjaqzmXng028Vaecg/2x3lbu1waaD7LvFXvtYOuPwiaibSF7fVB3
i7b8FwF/bJbX/RxT44NBf5P+5ZX74ZmfZXtsVWiY8Ytr0/pnvrjntumhE3y+
GPdsz1+r+7449PGa8zNpXzyYNC7sXeuLz68y13Ot98VVVz/n0rN9seenkfYe
e19c8c1+rQid5LV2nX1fj/piZvmJp3I/aGjFcfhqyIwv3rHWf3WXlYZeo88l
s5/4Yc81tcnhXz54PuR16PJeP1Rs/Nu/YtQHDwvxP9722Q+dC3vMNnb54JnA
Pz1KS35YZ8J5zCfTB5M4R1ts1pNc535xNkreB3cuzkttTPNHXfZLkqom3qg+
80OtyyEADWDOO3O9N3brn1eh+QTg7pevXv2V9kbbHUWfPwQH4C+N2KtVLN64
gnH0+PMDASj2SUjowyMvVK9cockoDcDC7l22im5e6GZvlS/FHIjpAXvfXqR7
4i3JiotsNwNxT7XctUJXT9yigNu2NwRif6fu9mIzTyx6kjjV3RqIXzpfdWbL
e6KityXrod5A3JXL4Wr51QNlAr/wHmEOwgllpdKb2z2Q4cx6Xdo9CHn5dnwf
3euOXxoybR7/CMITDz6d9Ixwx37mM+Xmf4KwT4z/dKOPO8a9dDfoY6Iji+vh
xwcN3DHXWDLLbBWps47Lx9rn3FCBFY5+InXg7X/EJeGgGy41DGlv20RyXrr6
hbqDrvhIeml7djIdkw3brXS2uOIijG7XPUVy2qYVAjdpriixcuxJdxodR81l
Zc7rkrqMWXFp9gIdex+FnVGddMGHL39FL79GR4vNH1fvpFxQfv8atO2mo01b
1WONZCdUqngYa8XHQNu8a+YcwU44mK3Y6yjEwPVzPzZ+MHFC9tLpbS6iDLwr
/6w7/JcjCsx6GljLkr50csNl4VBHzPm8Tn9+AwOrWGJfToEDum5tGYolfWtz
T492xzTpc2eUjr9NYeCebbuNFZ7ZYnkhXXdVKgM3fCuf319pi9cvXMnTS2Mg
b111A0+ELX723pSRmMXATDFJjdn3Nth8zJLzRyEDNRO6QnLarHG079yRyQYG
VscdKc4/Z4lKNXc23/rCwKOt/n2boixRMfvp04AxBpY3L7ztt7dEGzbN10zj
ZP3drPYr2SwxPtG93WCagezU/Lz0Xgs0l7Nau3eBgWfZG+WCfc3xJP0cXYqH
wvLl4BslsQl3O286LqNNoUpKKG3PL8BUjVjxSF0KY/q+HIzpBlSzjN5fof+f
p94M2BwPiHO/TkgbU8h5sfawVocp6k7cmuu0pDChen/+Jm9j1MgZP6rkRebW
UPLVRnZj9As7SBPzIX9P/x64vsEIH4xL/VhJ5u4OrnqdBXEjXC8R7dYTQKGE
QNdfnY+G2Nwr+EczlMIZOfXzxlv0cZlktPn4TgpPSq5TMBHXx7W7NnSe2k3h
nvKjnXpdesi17uNKJeIF5gl1d4XW6+EQZVxuv59Cq51BxYnTOjgWS9/seYR4
rIZslGuRDpaKCud1xFF4ePLgijUeOhhSYHzUkHjGqvGm+TN12rhm57qIVYkU
CjVpb9+4TwsHb8albCae8rRrs8zbtVqovCXGoY54DH4c947q1cSSyjvNc6co
fHDsVVmqoSbyHq6r23OGeMQ1u8BTCxpYuvnma9UMsj/9WRP2/9SQ/kIuN5t4
VEMkz879m9XQf/URmcPEs6Z4mN/nv1bF791vmYKIh3HaLM71VKhg697SJ3wF
JOflHYn5670WF9U+nDMjHidzuddy+r4Sngmy/chMPG8tm67iZ3UllAvhP9FE
+GL12+t3WBTRIXf5R6UrxIt+WU/L3JBDk9lz1+8Tb3xZruY/ISKHnc1cjQHE
KyN6FphuH5PFEo5z0rOEr27amrbRVwaHBUV/8JHcFng9o+1XqzTO6/R4XSS8
q7+ho3K9NL6+wCkkRTx23Mm4VmqVJG5pCjgvTDx3/kJfQMyAOCpt2Eg7Tfjv
vrvP22rFsNFeYMc/wgq7csxcgkRw2TXbwpfEm9u5Ol+d1BbGpdnou7rEq+dE
CvxaVq7B3d+VVTIJh+VU7BCrFcC1TTMrLIinGxRJZUEKPz4dfxuTQfi0RE0v
PZAPn8T8MRwknMTCnnxQixfLtpl6KBHvj9vd8jWdnQddeHPvRxLWsF6R7nR3
NXJ0LiZeJdyY/OXLtlAOPLPO89Ig4VWjUaFJ3Cvx5990DkGSM3oCfxRk31mO
TwQruywIQ1T1+pJgVuTJzBiOIqzMsfrqtdXMuKHO1CWL8NquAv2gqH/Nov3F
PLcJ76vItw14MNe8rvCuwmvCCh7NFR5rfjXbNB08M07YVbA8xWrLVPNjkY+O
y0iuUtFcIx5xeLT5ePsI/b8cp20q89BIcbD5I1Pyo/9yXOpjnas6+p3Nt9bj
sf9y3Oxp0JZwrTUNGk7L/C/H3TvDplza2mPa8u7PnAThfW7tJ+INP5sa//1T
K0h4mP13Tdm5H6aNLukP2AlfsguQCP05Y1ok3qn4m6wnwfZ9SPTlP6a994v7
BwnHpHYl7nFcNH34UuFb23+59uWD8jsLTOCxytPhCuGJQ1UetW4s8NHekDnh
v/tZGaBaXsoGzqX9nDTCEj92X8hZWgGTNIetqoRXLZSuTPRYBefaTkv/Jfv9
Or58JKKcE3YfqFa+T1ioY8pmdQUPiD5/bWhFOOTf2WtrnvHCwTMnSpkJVw/f
eCzxkw+sFVsuNJDnXVHEt0bcVBC4peaWSxB2Nqhq4A8RAhl32aMP/zs/ovFL
bClr4A6L8NFthF2n9zj2vhABK68HXNfI+Yv+vdR557coMIVfzbEgHJwU/zxd
XBwePa8rf0vO68klex+tMEloiWzwniHnXfRngdupgzLg9SL26zTJaYUbvr/u
55KFTlprdjjhiDdFLOqXZaEhV/5pD6mnAxkfn99tkwPLPfWPykn9HThQROVz
KcKtJI1lYqQ+2RJcqM4CRYjcZ9C/heQ4S+XcmSlNJTjtc2rTbVLPVz80m6zz
Xgv8EfcumpP6f2JKuxyUrwKQYF0kT/qDOfL5+G9QhWjp9HrXfHI/iZPJbg9U
wSA/3/0A6SdfRf+0q3xRg0qxTrkm0n/weOV40rr1wFDxVuG5QKFj8my2Pa6H
0pDUBJFscn7eFruxKmuAmyiESJ6nkNqZU17XoAHNMqvbRDPJ86w5kf/lwwZw
bJRm6U8j6//y6UyZoibIb+O6/oD0P91rM/LhUZpgfH/fcPFpCpX6GL9fMWmB
Z0vPLI30S9GmIY/9ctrg355ScTqZ3M+KbVWvw3VBLSejWvYQycV3c+O33tKF
VQum/EWxFJb0Pf2ywKoHejJ1GyUPkH5rc0hiVaEejG2u+8y6l0K5UNOWl336
EBi1rTyd5D6P4AsneJw2Qu3gicB6isJ/heijm7sRHjw5unOYTuHzrLe+HmMb
Qf3PS0/OIJLrpfIL9yUZgeSztFMOfhSmBIb0BDUZg9eHbY9Pk9zHZKPRYScE
cPsx56KoObm/wColUWOA0mlVm59A6unBxonhYIBi1qtbOkwo7GPOTdtaBbBT
QsEklOTAf8+kx8wcNoHxWpQN2ECeh11IJHO8GRheVI/1kqTQbMBuW9yIBVTu
aDPcPcvA0pHEm+KrLKFs7bIL70huc7N6016jbgmuf7RTjSZIDoqtsX2x2xLW
m4kYzoyS+T97uqab2QqWBqVuqw4wyJzN3NMqYQ36Fav7Hj9iYKeSxKiPsy20
JDgWN2cz0KWmvNpkpy3Mz70UzMlkYOxwp4Fkpi0YjLFoR59j4JXEA4Ed72xh
o/Th/dynGJjGeYf3TagdlHCMmogeYWDH4+1Xag7Yw6GQg+GMEAYyZYkIOOQ7
wtcUseWNqgz8wbp3lWarI6SPsuQ0KTGwwjEjnHvUEZ5NWBk0yjHQV3ef6k1N
J/DxiLtUKkau7+t183ibE3yA4/JWHAw02um0+dSoM6zb36OcOUrHFdJfj9oJ
usLaJ3J3TxbSsVtJokZZyRUO1dMEpvPoOPm6Z5zFwBX0im4XuBM/dH1gsqPM
3xUeq9r6sBN/7GR5JHW7yBXCHCIi9A/R8dtLo4IBDTfYYNZsUutFx9PRkwVa
1u7gmFuTK7+cjv+G6GKiPu4Q8u+TA+cyOvps/qM5H+EOZfL5PAO/g9DzynPF
yhPuMOks5Oj6KQgflM7kdXe5w37WWNfOO0EY3ryubpW7B9yo7+QTDg1C4+Or
tQxoniDJ5e+oXBdI6lHec3OkJ6y6vkFjqCwQxyOY4jP3ewJ/gGdK2qVA/PFp
5XDvBU+w2CJX8To5EMPO3qWJvvOExJwrbHO0QBT2mxeLoXnBpauRw5ULARjK
vULkgLs3BF/I0+7dGIAlbwZX9TG8Iedui7uvRgCmRW520432hp4XW+iv5APQ
M7j4Vd8pb2BJO+xYs5rkwrlii49t3nClQbdNfMCf5KXpp726PhBrZFtSEOeP
r26k1/sJ0OCS9Wr3uCY/ZPRq32OSpYGUV9I1+Wo/fDA6WZG/ngasDGmX+yV+
OPMyWa/TjgbTl+/Ej57yw3u+ihVv42jQy/rt4xdfP+Q4Wr+24xsNlnSGLa/+
9MVjDSxTk42+QKWVK1op+mL40Qr2jx2+UPI+VHxxjS/axx483PnWF/6sh6Jr
K33RatX93cmzvrA8sWj3z+801AnZLlGk6gfRym+y19XQkHXLmpnQbD94tcbj
qJkJDTMHWu8vbfGH1E+UBpezD66TaBaf3ecPBU8f+wyY+uDm+iDfoUR/cGh+
kl2m4YPWF/t2FOX7w2LwYrQynw+qqOhUtL3whwSLN/NvXnmjstlT0SN6AfBw
l4rLOMnDnkJWrJVzAeDhP7So4O+Fv/g3nexiCwR3a8HVWo5emFM8kzfCEwjG
+1vPGZp44WJTrB+zUiDZz908BtJe+PhE6NiAG3nd6Kjx7CdP9HBX/aNRFghs
85JXwjZ7YuTji4bCLkFw7IWrjtcuD/SeE1lu6RsE9eH5jDvBHlio9ZB9c0gQ
cDJNHlrj4YEa+/LKcvcFwccpmaxWbQ982KvGX1wQBFw9rVUff5Icd/m1n8FE
EPAP/1ymSXLeTjmhbbrmdBAR1SvRj3FDRqPSm1vWdJDyHRZ0D3NDnZbb/loO
dPhdniUb6eWG/Y28TdKedMhebSR1XNcNP1gO1jwJo8NOFimec79c8d0oa3xE
Ch0ECpzX9u5xRXsr79jSTjrIrDyh5LzDBQcz96SXd9MhqZn/+DDdBa23F5aW
vqTD2meiTDEuLnhSitae3UcHvUUXq9T1LsjM3FFr840Oown3NcX6nJFeEJA6
spIBHF4BOVXPnTBsQeMDvwUD3tDdwzbfcUCHlIjGAWsGFOx8K/wkwwE97Fe2
XLFnAFDarQrRDnjXzVZX0Z0BN7YXrGtRdEA3hZCaMQYD9ljy77qQbo9Wncd5
Uw4zQHr068O3W+xQW+fHDoFaBhQ/4jvUJ2SDWsqFq8zrGRB+2ybMZcoam/r3
lG1tYIBP7mBmU6c17nKb/HkTGeC8cvt0Upw13hApr1p8wgCzrT8ePRq1wpNu
Z2evDDOgjntRQ/SuJd5Plv6pyUMBY8hyjwnNHIeatRQy+SigLVZeidMxR/Pb
TaPTAhSoG4crNPGaow7ju/plEQr2ruxQV+0ww5U7bp17JEuBcNiBoKeGZjja
z+dmpk1BOG3XDIhtwjMup5MCdCl4Fv5B0Og3oLNci9VufQrMmjn4tF4C6g6Y
RWQYUSDA8rtB4CTgYmilTL4FBTtL15xy/mmKKscil/W7U7C5Z/7P2mpjnGSJ
vlXlScHtiif3CiON8SLj2EScNwUDX/hdBWWMUbWE7zuvHwUJh5apfD5jhO62
V5iXGBSIL8vL2rBtI8b8cLBZ2k6BuXH02ClFA2Tbf4jjRDQFuoul3y5+0McJ
EVUH3p0UiHiyvy3O0Ee330cs+PZQEPYo704xsz5q+PYv7z9Agd2CckZgvy7+
OnvkWG4iBZlRafYdGdr40KLMmi2Zgqisw//+mmmjwl+H0IjjFJTNZyfJTWlh
qY77dYmTFBw/LrwzxFELEwSmtlikUSA0zbe/gE0TDw7s/n3kLAU5Ewrxl2s2
4Eu+Te/qzlFw4I6UfB5jA4peMSzmyaTgZFTbpWP3NDDicYdpwgUKbmRaWK0b
W4el3Jdi03IokLjlemjAch1+TPzFceEiBVf/7TqQfFkdu3V/rb2QR0F3nfuZ
Rh81fD3aae95mVyPNfsStivj8h5vOd1CCvIq3eLVZZXxkJuAI28RBQWsknfP
HVqLjgmCWbeKKWh3TF+w0lZCkzzv341XKZCzUtP2vyyP2vUHaNGlZH0Xzjic
XZDDgfivHJJlFDgmZziitxyqKrygR5RTYEPteMHJJYv8rqPMbBUUVL5ubVCJ
lMGAV/GzJOeB/PObHuYk9/HkxLXdvkbBZ8v6Jx2zkviJK63YpJJwv2Vif7kE
LoddL5oIe8kpaY3SxfGNi6jo9esUvLshu+zzExGsNc/UkrhBQcNE8NibBGEM
FnG7+N/nKqdTLf1a9NegRzmfp81NCpRvuYbFFwvgX9++KwWEc7bXPfb25Uer
3iirGcKJJarWa3n4cPjXLS347/tZlZ7mM208OLHpwvZkwjyrKy13/ubCBdZl
vx/99/2suJKMWYoTZa3fPWGtpkCvdtBo+9OV+LuTb9KQ8Jxt5YpBgxVYxJ9D
30KY396/ya6EFd/5MiRJjgPTEVC5wcOMX3P8lP/7HKLDZrqy+J+l5psph+Lf
Ea5lZj797Onf5psX7inO/vf+K9cIHiqZbWab4xchuQouFxcXWv350bw764CP
EOFJw+flKSXDzf5rxoclCf+IhWenSl41///35az6UhYTCuTv/Q+2/mTC
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{{3.498652363276937*^9, 3.49865238272305*^9}, 
   3.498744166653677*^9, 3.498744987001581*^9, 3.498746069551092*^9, 
   3.49874634415429*^9, 3.498746675031761*^9, 3.498978127553935*^9, 
   3.4989790020899353`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"wezmy", " ", "wartosc", " ", "ujemna", " ", "Etot"}], ",", " ", 
    RowBox[{"ale", " ", "dwukrotnie", " ", "blizsza", " ", "zera"}]}], " ", 
   "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "kappa"}], "/", "dist"}], " ", "+", 
     RowBox[{"m", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"dist", "*", "omega"}], ")"}], "^", "2"}], "/", "2"}]}]}], 
    ")"}], "/", "2"}]}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, 3.498746238227403*^9}],

Cell[BoxData[
 RowBox[{"-", "1.657924847880946`"}]], "Output",
 CellChangeTimes->{3.498746344397765*^9, 3.498746675291688*^9, 
  3.4989781275859356`*^9, 3.4989790021199355`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746253474484*^9}],

Cell[BoxData["0.7071067811865476`"], "Output",
 CellChangeTimes->{3.498746344414151*^9, 3.498746675305136*^9, 
  3.498978127598935*^9, 3.4989790021349354`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", "0", ",", 
     RowBox[{"2", "*", "Pi"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw12nc4lf//OHBUeFekZBQqIruMbPE0MjIP5yDZ65xzy5bNaRcyirJKlOyi
VDTkWSFbQ0mZkTJKC6XwvT/X9fudP5zrcd33y33fz3Xf93Udca8gO18ONja2
3+Sf/323Z56ziv5MA7b/93GN7qiPJL3/6gnB9EtN+hkXonnDSD8fHr+Xfum1
/rh2tr4/ae7qYkdF9WH97gwpigdp9bv+phmXPuo3XBczsiOtzf4zfmHllP6W
8FgBA9J9cRHFyurf9Pfd1GxRIJ0zt3LFp4kf+um81gcESBMr4+MvXZrVl9xY
07XwiQbDZXf1aJTf+tKybhIDpKd6EoN4V/3VP5dh7PCQtPCBLs3WukV9vShn
Ipt0wN/rI4bqbND+7pxbCGm6xCpRiXR2KL00pGJK+qjHW3eOSQ6olo6KnRin
QRUW0povrYK9PFdm7pCWZg9SKp/nhI9pF41YpDvO+fmlUbjhhlZNEjfpK9+P
qLmtWgOdYUXxLR9psFsybt7UfS08DuazPUF6MsjKT/UeD1Doc+dmx2iwKoKv
65Q6H1T46a4oI21m/W/v4VE+UNh93m4/6ZfKsvXR6evh6PPVqdWjNDg4uvll
0OQGOHU/K1zrAw3q3OVvM/MF4PtTnnevRmiwNvyLkb+FIHyxTRfzJz2pmC0Z
+FsQrhtK6KUP00C/PZo9yk4YZlfyjNwepMGraWnTS5wi8DTvqIo6aUfuydjy
GhEoaLzjcXuABldvqNvf8xCF13ez7Er6aZCWrLeh/74YfLkzeJL+joyv3uvI
PI1toBcwUNfdR8b3lbK8bOo2EPnE3rybNG114dvasW3wIuld8GwvDXyei6q/
PSsONe7uHJTXNGBNWyRJTUuApHXeQHA3De6sd+aZL5CCD9T2+PNdNAh9MCg+
vCQFU31Av9tJxtO7+GGbyw5oDUnY+KWdBsdnfuUXCUuD4uroTsUWGlTsM+rw
T5eBHy5vXP81kPXT3x1hdVge5M9aeX18RIO2xwILlEF5kFpyKmurp0HjnJ+/
k64CqMyZLiY+oMH0WyEV/98K0LJm+sb7uzQwsBJfVRa8E3Z2yDxrq6TBqIXd
3aZqJQh4iHWRFTSA/oXpt0NKMJmZP7atnAa3d0daf+FVBvcDQ5/pJWQ97nqf
vylAGR72smW0FNLgxfm+7sOyKnDBv27s1XkajJRfMDlzRRV+qRwd5ImlgYdl
9NYb59WhRGnwlHs0De6lKPi2P1aHvs7DPdcjacAlq/ff5Bd1sNXjoUE4+f+E
5yLkTTTApER8RDeABio7Vqk0zGmARfN4WoYbGd/F1GAFJy1IDjN1uqVHxv/b
akclUV3YL5wlkKVLAw6H5qgEM10w8TRxj9KmQddRX8uucF3IPJmSqaROg5r7
t/nDOnXh8vHZoEhFsv+f2vm/SNgDxv+IXW4iNFA1az02OqIHnMoutxJ/UWGh
Ys/D2h6A2gyh7Wt+UOEFHy37zCSA7ycb6cQZKgyvunXQi80A8hNOj8dMUsFK
SWDFegUDWPOBc4faMBVOZPlyxxwzgMCbhj0r2qkQ8m/6WJyKISS/3LbMkU8F
HsOo6xZnjaDJz2dxYx4VNMdufDYoMYLEr5uDtmdTYe/CdUOteiOwMD1WqnaO
CmVD7Y7yE0bAts2buuskFXY9mZAWMzCGMyebMyMCqMBVqk5ofjcG/NGVZqlF
hXXMx49pdiZwrrRSvVCNCq5Cq35RCBPwOc92/LsyFRSu+hbZHDWBvdcT15yQ
o8K974dO2NwygTdj4gWRIlQwMUwN9dxgCl+tTzMS/9lDW5nv3u6XphC42GBt
/8Ae+BR+iA9QzaF/JmpOSN4eQoUO8/O7WcKo/w//3B324Osl7dIXbgnOGjw3
hCXsIfmAtGFhsiV0/eU8y73JHsrSaibV7lkCa3vj5Tuc9nDbPLwzdKMVlL5c
MDUasQN+vv5XOu1WED9kn7c1ww4yj7+5tEHbBk4+WXvm4AQFdlzjvnyKlwK3
PWWuO49SoDBQ6YiICAWenRA0NRqggMX01oRqaQqYvsS+FS8o4HjgbdsAUODj
0yVpwzoKvN6x/fHeMAqc/hzfXn2CAlfP81RYvKXA4KWzV0GUAvUnPhSrXrWD
KwUR0WEatjDorN2ToUsFrfSzEquVbOHnnpwRdXMq7FyGilwZW7CZuyH+nkaF
UgNxq7JNtnA05JGOdBAVwo2W+3/P2QBfUvxoZyEVPl978/1LrA1Iae1zp3HR
gP/+q59746zhj/SeVptXNIiYaa/5TrWEmJYo61VBjpCj6Kdx+rYpuIhtoLfE
OcPzcZ0vlWUGwDfA0owYdoWfUh0ODxvVYM/N9s23UjygodRsqSxJDUL7XZQE
Mj2gMCcwNMhWDV5r+V+JyvWA8hUO9kfe74bQNfw5UOIBqG4yWz+jCiH0Gw96
0APWF7B5HuVTAS9GowzPLw9w8GDAiI8CVBSLdNc5e4KUSdCzw01CMMO25kKT
jBf4s69wnMkRhHzO+E1Nu7yAfcJjZ3yAAKh/GXnUqO4FQ7d9Nvxazw/P218x
Go294L/A02ZRBWthZ1jJf00eXhDK7nwh/uQPfZkZx1+t2V7wVap+h8o7Xtyw
cF3nNbc3PCVaaix0FdHyt86/qM/eUKnEUns5pohNi6p3/nz1hq1qEgVdKTvR
xvDjuphZb+Bl54qSG9mF93N4CqM5fOA/MZHcBybKGNhaPhMm6gPsAWktDkuq
WOXqJWxv6wMHtgu1fc3Twj1FaURZrQ+UG7ztGZnXwiTfxdBNj3wg1i5OQoCq
jbyxVS2nG31AJDVJfQ+PDmorzNn4vvCBXW+vCWUf1UXlev9ogSkf+BdVUnc5
Wh/Ln2gfttniC+v5T4dt6jTE+F9mTP7jvnDnl2VQzAoj7Jt51vM1yReuVn8/
ulLbCBP7K761nPWFqgt/GZklRlgzdGZH9GVf6PD+El583BgNtdevaL3vC8bd
DoPXTUxwxPy5hMR3X4iejl+3dcocv1vkbTjv4gc3mhwqBnbsw5clJzN1vf2g
aejCrlavfbgoYCXygekH3xkKgdLv9yG1SG2rdKQfcChW7lrZZYHRn45Qs8/6
gaCcoJ7zQyt8Pu4ys77ZD0SCnfvkBm1Ra/sMqCrSYe0RviDXJVsMf1nEHaJK
h8GqV2klWyj4inNBsFKLDqrQdeOoOwW9jz+Y3GRCh935D7gdPlBwt6brlVdu
dFjY/s0lYdoOnxcVNdWl0cHsUE2bHj8Nj3f41G2YocNlg6+fWeo0jAjYQEz+
ooO1h2Xou/00fCmR4/VogQ4ftUvY+wtoeCn76MoDnAzQ2i5UvUHJAVMeXef1
E2PALeNq8QP2jtiaYyH+YB8D9vF/jKi/vh9XOonKiF1hwIRE0YEbL/bjFecn
4+tKGDD2N3Py7ux+jD112mGpglzfFqC6Ts8ZN8//lOq6w4CnUd7XBrqd8bLu
69OarQzguqKf+Hf+AFr+WOWoNMMAl+J5xYNiLujf9uBD308GSB6nyM4auqBJ
RWBywm8GhAjZWdmmumCa4fDae+xMkDpjwFTY4Yr3Pf0ffudnwvuK5DvHnd2w
1P4zB1ODCcldIyzWETdkzVncbNBhQp14jn5GqRv+SUp4xwdM0Jk737hq3g2Z
nxpDr5gxQYataFzovDvKbn2UcmQ/ExI5fyn8IzxwU9vPl1IxTJh22MIZescT
s/qLLt69ywSO5S/Gth2eCE1rmzPuM8EsenfbnlFPnGnqesh4xASt0VN1Juu9
0FGZ5+WfJiZYused4wv0wmYe5Z67PUwQtpjxOyXnjdeuztxM+cYEJbemjgOG
3ljfPjg/9JMJg9Od1frO3lhWGrksM0+eX7j6klqSNybXp/FcXmTCCav1v/on
vbE5N1tt02oCvq7mlxCs8sF47toJNXECzvPyRo4+88Eml4EPU9sJEF3ayMRh
H3wfeTUmawcBf434ZMo2+GKrDq92tzwBbr3+UuqRvuh8eNvQIzUCtLaYRR5O
98Vn+2/c3qpJQFHcGeJdmS/KvwpIPaRNwB1zXrX6fl+8VXV8eVmPgB9edbZ6
Rn74euuXshRTAlySWNdWufphJZEUlmlOgIlWSVZ/hB/2mPwoSLMgwIONsLlX
7odOd2Vee9kQYHrMdkNDox/uhYZb2hQCChIVA3oG/VDEvrqQ054ATY/6UUV+
OrbGKHjEOBAwe+rU3xBFOo6GP14l4UTAS46Iq02mdLz3opFRv5+ANzPquXlx
dHTW3LGh24WA3OuJveyf6Hj0aXbuAS8CcPyf359lOkZb77+V6E0A37V9IUvC
DLQJZz9V6UPAcnHEpJIFA1W587c+9yNAdR27vZsPAz/I3Wx4TieA7WK9bE48
AzcwJnqaGARkylTwa1QzULJvr2AKQUBafNjW3FYG7nv+T8zdn4Cdv/KSVo8y
cOzoEkXqIAEvOIQop/8xMHPp8+Uh0kOvbgXyCzLxtYD5bEoAASlfhIcqdjFx
bVCuvlIgAeQD3yWKORM3ZpoebCad3qZUvNKbiYkc5QzbIAIa98b+bIxjoqfu
aoku0gtNePjcBSZa9V0+oR9MwKDXV7OAaiamZZ89XER68VCbCa2NiRFhs8v/
SF/R/xW+b4yJbtXPv5mFkPlLHu20XGKi05wQnCbNQ/1s4SJMoKtn1bv7pKvo
chNRKgSqZRHXh0nnri67fMWSwKH+tdf+kA7M1Wb2+RFYrGB+Y1UoATeVzuiL
HSFwzPLJnZWkJ43FRQPzCDT7tLp0jtw/3tr+d/sdAscLroW8/9/x/Cid6s8J
fKayavkW6ZLxrZk3Jgm0PXZGN/5/58ebvVd1lT9uxlRBbdJP8+nvGrf6460o
m+gJ8nom4vgsvLX9sfceByWZ9FrK5zRemj8myXYmbiX9U0q4sDnIH29ub1xX
TMaLnftiRHKSP6b0funaQvrjVUNu12vkel+DmkQy3gtjlVY66I9Hl4tLx8n8
rP3BeCE254/2MsuHIsj8VRzYuCS+/iCGJ0VJlZH5Pjys9kBF4SDeG15/vpOs
h+jbxx9Hex1Enml+7ymyXsqUU+er4g+itOTZ7DGynu65NWX9yD6I3pcddbvJ
eive6/z1YvdB/B1EiYsi61HEdHX6yinSJ+Vn1Ml6NX+UHBrFGYDlNRfjxz0J
qI2I7Tq+JwA7M9lyhd0JyDt2TVlsfwASe5SXcl0JOPJSqOxxeACOaB2S4yX7
oUN7XfHOygA8KyF7qpnsn82/KvnnnwXghXNhKWscCRj2eurbPhqApfdPchnQ
CLh/VPdElkggbpTU5ogn+zPQ6JzJw+RAXGxqqFcn+7t6Y63ynvJAFNlVcrfD
hNy+u/N8e0sgcrNuXrUxJuuj+y+LjzMIHfNSpkX0CaBLindqHglCgSBBjXZV
cn5BTej2giDUT5DvblIi6/uk2/zmhiDU9HmVUKlI2nNun9JiEBYM3ryzR5rs
F+XSkd+6wWhqZnDo4CYCXruuOVlgGYxPQ17XNwqQ+ZnNMrR2Ccaa7BI97g0E
lGvpv34QG4wtYrejbMl5GZoXxjt4Pxg/dxTf//CHCaae10pUtUOwbNCxaZac
zwLJPt/U94VgM/fUgV/dTHjyBA/rO4dgdwQfvb+NCYGFN/Z4xISgn0wrk4Hk
9l/bjw3eD8Ex1nbugHIm9LxMCkDdUFR682ZyHXn/KLQ7NMFnFYoZF66fzg5j
wtPNORW+rqEo9dSxlzuACbf7RR+IJIRilha/UIE7E9w/2XLVNYTiJO/c+mIj
JhRFKBQVGIWhqq/8jXlO8n7n8vnMkH0YajzLYt+4zIC/OqlMcZ8wtJLO1BWa
Z4B119L168fDsM7tFcfLcQac6bq9fbwpDAs2CnELNjFAzyhZoNEsHJfZdi4/
jmPAZdZLF5n94UjTEFmAMAZw1gXsSWWG47rYKqtrTAZEhp4/6JYUjrLKI/li
DgxYecBzjKcjHD1SgdKmyIDX1zgyMmwOYe/87/LlXjpEFIT8XuVxCDmS2XVb
O8jnGxD0iwo+hMa0nLsRj+lAdZH96nL2EKYGOGyIL6eD4+HIW+KvDqGBbSK/
WiwdGjfbuB1xiED3YwVuIESHslcSOS0HItF1x8WHoEU+f+2+sPYLEYkrqp23
ecn7QRcrSZ83JhLfjL4VZor5QcXpcxEm2ZF44aZ8njS7H1xsW5ZI7olEj5Mz
tb7PfCEkxLfF0jIKOSSf3h6w9oWGrkZ+pQNR6OKjWmiv7wvyXeuceIkoFOnR
Tr6xyxeOH+6iPjgVhUfW0HX4+XyBfZda36snUXg50s4+q9sH5C/dCh3TiMZQ
0a6UE+Y+EC+4X5xpEo0fm/zSzTR9wC0zeWicGo1+Rc/0v+3wAdf+f2s6Q6Kx
eY6Z92WFD5yyvVsKFdHoFiYlWPDQG2aOudwtFY3BHXlfbo7KeENfDIf3+d8x
WL++v/LXtCdo+OibinHFokNhxPv0N56gvzy0PVsgFmeS4w/woCd0ihYJ+anE
4ouzFqySc57gOhy81c0/Fr/ZRD32VPeE2Ay2PdnvYvFa+lrG30gPqCz8IKjz
ORY9T1Qyrrp7QPx1lO2ajUWRkOj7CqYeoGg4Y97IF4c/fIeujQp4AJvlzsKp
vXH44fsewTZvd/g3OiPLdj0O92RZrspqcQWtcfq3s/fi0PSZHUM60xW6L2sr
r2+OQ2HhU3nx7q4gZq4a3jsYh/TYc4y4WRdgpekYJ6yLR33hjottW1yAvZUt
eV1APB5ukY/46OcMIdvin/tExWMX46/+j53OwOE7+bHoWDyGnLA2aJ3bD2ef
1b8by4lHg2liV++J/VD/5OtB96Z43P2SGI4tdIKWc1t7RYQTcFuOb6tatwOs
ex+w0Vs8AVuFP1T8y3IAhQePKs/LJaCX2dOUHA8HEKZJhD7UTcCuT7YXHL7T
YOW38V9rPRJQUcTv3TgfDdZwKVZdvpKAI7qbuSvN7OHmpGRUankCvugwm1vg
Jd/DWSpP/W4l4PtLn7ykXttBZ75VyKvHCTh1oGw3r6cdBHYNauQOJeDg2sv2
2pEU2JZyJb9ZkIXZAdtrS8j387tsGdIGYizU0vrzd/1Ta3B28c6+up2FVXI3
3+tbWMN6P47XW5VYeGDtmmcTzlYgKzpl+NOUhQ+0t2xKjbSAWJnos+VWLPwv
uGlOYWkfpLZVCprbs7CFZ4X5heP7IOxEQKiOGwuDHaeft6ebw2aLfP2foSxc
uMB74H2xKRyv817/IZKFzcqLOC5vCloLhvG1cSzcRt3b/LTaBE61xL4XOMlC
oR16znP39sI32/Ku/CwWFj+NyHBsMwKNtxXJmEfuP/Oh9z9rI6DzWrm0XGbh
2xsViudfGEI++8p/iSUsjBVS7dN+awBKAX+uvLrDwsLjV6buPdOHSZMtS/z3
WKhu+CT/5Qs9uJJIXanzkLw+63PqVe/3wMkX+mO6T1h47bRl4IOvOnB0w78S
oSYW1lOvPPo6rw0+f99Ovn3GQl5OgvqFTRsSmr408XayMF1vd4P9Bk3o6uWy
ONbNQkf1q9YNmzXA+wqX7OALFk53mWQvSqiDVl5Ep8kbFq7Zrvj+j8puuPLi
+jHqW3L/tzFZD7RVyXf5A6Fm71jYmvIvw9FQBQaLqvs/DLCwYuPUgBRFCTwP
rbNIHCKPf0xpqrZ6J6j42PcKjLDwC9H89NMjBZi52Rh7+gMLX8Vc4R/qkIOD
11erjoyy8GivbHnhOxmYHZP/vfUjC98vrnRS/bwDluxTm43HWWh/mHsma1YS
zqXM59h8YqHlHTHDVxzbyfeNM4HwmYXsjvJSE+vE4fjpVyA0wUI/5X9997W2
gGBNNU8PaeKA1yDdfjO41TJfRkyy0OSM4sqRg4JQLF50Zpl03H8GfxRPboDh
dl/twCkyPxVb+yr1eGBwZuO7J6R759+JNAavAh71KGKRtLVOrIbX3r/6H/Hu
pNg0C4fmqVNr7T7pp+12cpYk/XfvLUOvvO4GjwfG9/hIaxVm6Pam/WzI52zi
GiXXZ3Bde85zgR2/Whw1ySP9p+SFE4/Pf6iufzJCi/T+yYvEl2w+XPmfTk49
eX47n9WH1rA24kkt/2pp0kX53gIUP2GcqzV4FEle3+mNW4seWYpi2Q3zputk
PHpVezRjNm3DY2GBLa1kvPLNuZiH1krg5oC/be1kPH12hfFYLW9Hu195z2vI
eI+VP6QsfpfCoyoHBo6MsTDBIcbr6Jg0Gv+7OqtG5kefyzLg7RtZdPvo7GI7
zMI78z6Kn+8rYs/J2Dt3B1m46f2UUW7lLiylJe1aQdaDmUbuGq5uJdTbNHXZ
uo+F8WFSowKogv7Ex6tWvSwU9i7kotep4lHJLWMqr1l4e3VD5Nnq3TjTUKRZ
/ZysnycfbhIF6mi+MtHVrIuFkToj9wSyNfDLUO9QczsLB5pPJ55L08RP5YPi
kc3k/NiYyrUQr42+lI6aa/dZuI9rMcLSSQ/9J4kt8rUsnI2vWdjprY/6bzc/
PlfDwhVDfvVj6oDrTn+c5q5kYU7b8VReAwPcevtw5Tmyf4Hdiodlb4QupblX
5Mn+Fm5ivqzqNUKfMxf2l5xjYZCSqHvNAWM8OHrc0TSRhUm1Dvt2eu/F9xuX
W+UPsfDm4DHeniBTtLhcutAfxELfGsUbsT9MkX8o51wIwUKeFXo7v4SboWbb
T/u97iwUiLkXIRttjgu3p+7kkPPMYsaCbyLBAuOGpbet5WfhxoeBw9kJNtjI
CPVYXsPCTk0tb/dFGxxNSinuXUnm/6x62sweWwTJzHsacwm4Lve0zIqHtji0
0G18rzcBq//JWBy8S8FHidY+A+cT8FnbpMHAVXucKqXF70pJwID2Tt8Xw/Yo
F9aQwziegOZ2X+pzxaiY0uVATw9NwHnqD2rFBSo6bP26675VAurUaXjWnaTh
7P5QxQvL8difm1gh7eaIGuecOIXn4nF843+ybpmOWGblI3xyOh7N8jxGmG2O
yNhS5ivdF48WQ1Vvl9WcEJ4PiB29GY/pNlxfxdbsx2CZLpPbbvG4RqKfY7zS
GV/cq3/ITY3HdrsMUZthZ6S+8QyzMI/HaIXNyon8B8i8scleVI3H5/lSzkHR
B9BsuNYvkCseWZbslouGLij4uHvpZ0UcMr9rl810uOK77gZCtyAOy9XOSFov
umKvl7b2ocw4lJh1NoxRdMOXlQOpdXFxmPlApEo7xQ0vytwdTLQk1x8yPea4
zx3ngpxqDk3EouXjiveBrh7Iyk25R98UizeM054L3ffErab6e4dXx2KiSlPG
w+eeOLr0p8DiXwxetl/na/DJEzt2Skf/HYzBHLHUlVP8Xij/Pjeu8WoM/srg
3VTv74VDue9WS8jHYEooNZgm4I1ZNy239IvEIOXLtk9zst64T07lfvLaGGxx
KZiN1PNGNf8nah1folE8KKRLmO6NOx49GnhTFY0OXTa6r+9644Lq+jMaKtH4
oEaOW8vWB/U/Nve6SkRjcXp++SsvHxwfObw3bkM0VswXzVgd8sH605wyV35E
oWYnhI/k+ODhNdN+F25F4VuDlhUnRnywZ13SdiPlKJxp7o9tJHwx87WV+125
SAwOLrL39/XDizdPtUZtjsRXwdd+N4f5oSbP0hbN1ZGokm9wZ+VRP9xSK7b2
xkQErpRKr9mR74edNhL0mNIINFp37hTnaz80TLdQqZaKwMQfpxOM99DRMq3q
u41ABH4x6QhrNadj2Ar9k19XRmDq7oAsZQc6vmWwX5EZO4STWrdrKwLpmCLK
9yjjyiEMET2ZaJ5PR9X9vUWuWw+hiNmlI3t+01Hf7snTed5DKNyV/3mJg4F/
L2f7n1sKx9dqjh1FPAxcExnzoHUgHDv+VnrnSzBwbNpI0OhiOB5+dsfawoKB
fesK3hzaFI5b6qPDCrIYqPuwb277f+G4Jm8w4nYBA4nhbWI9v8NwLdvOwBtl
DLzcXP9K+20YRvQU5lIfMFDK8rWxZFYYsgS3XD44wED5bDm/aIEwFOM+sndE
lIk0e4UJOmcYGpoYfRjazkTn1+/27Z8PxaLIUZVWOSaOh/33ybgvFIm4VCMX
TSYOvenZbnApFLXqqoNO2jGR9aZN551UKH7vexy7eJSJD4ezYhYFQ/HEAGvP
zdNMdAgQiJHkDsVXqy5yUlKZ6Psr/Wv8ZAj2f2hcScthYntVNY9HVQhmG19e
HXKDiTwl9jentEIwQzJk49rXTLxZxayzlQ/BhyM/PJP6mKhDHd1zXzQEx0P9
180NMPHv89qvucvB+Ljz2PD5cSZ+989bk9kUjE5/uu7+mWPiru5guS21wfh7
9LT7xwUmGko82HS9NBjTXqZsbFhiovEZy/S+5GA82xP4R5eTwHXHG8YC7IJx
qJTB5buRwFO/Y9O2Ggej9Wt//6dCBJ6Ru019oxaM5eNqp/hECHS8a0mx3xSM
7gZxOpHiBObwjywJPQ5CmZurpqIUCZSSIN50FQXhyvc5njZKBNonX6Gnnw5C
SXMubQFVAs/Ob5pVtQ3C35tnfII0CWy56+S5NBSIgoK3vrkbErgw2/XfzceB
2Ndh6fzRmEBT3Uvvg68GIsxXszmZEkhs3GYgTA/EzKKdgassCWzt45xfNgvE
lU/tCLAmsFmlYvOsXCC++07lYdoSeELTbZprJgDZ5QJuHqMS2HnM+r7UiwBc
3C3y6ZADgdqNqEapCcDzjW8UHJwI5LrXE/IiMgB/hi3ZjhwgsKh+k+0O5wA8
+t8p7zOuBD6PkRFK1A3Ab0m0r1LuBG4Wm8iKYg/ANXXmPZJeBFapbAhYNXYQ
U7m5vZO8Cbx89JdnQfNBbNm458WgD4HR/lODbGcOYtGphFtUOoGvJVwPNwce
xJ7t2U4RDNJBRuG5lIP4I8LQ+BSTXP/R6ViA0EEUbs+jRPoTuIVt4M7BBX+0
z2y65nCQwFVVh5kxA/6Yr3eoSCqAQG9L27IL6I/tnXV+Y6Qp3zYdfnzVH70X
klZlBBKotBjwe+GkPxZzvMpRDiIwvOKyqBHhj5lxPrpIejqzaS7Hyh8lZZdX
6QcTeFX/fvqykj+eabFed4P0rlLRX2Eb/bFwZjKYJ4Tc3m4rNz9PIG3QwciV
9IDSYd3E9wRqOrJfvkTaOyRXVa6BQJ4o44tdpJ88pgm/u0JgZIav/TfSrtG7
fmafJLBgzfl+jlACS3IOdvgRBA4pfrXhJD1J3VxpROY1YV1k0wK5v1smb+4u
FbKu/LfRhkm/if98SU6QwOOHmwRrST/4zdewm6x78R/yEizSJ1UerLQeZOKz
khVpmqSTJT4einjCRKmq8fBR8noceJw23Shm4vkt+6dYpHW4ymZ/JTGR+N3y
Hy/p0Ait9RZBTPRYl/ophYzPkUyN0Gp7JvZNTeQtk/H8byZ/ixQ5J77d5NLz
Jj33SVO0jJwzmdz/xurI+POqXQrWZSfnDHtTCRvp924q4kMfGUhl33dJm8zf
B72dqmfbGJhRfqDfj8xvM49B2vZMBj5Rnj17gcz/XEXPhxXRDBR8VqGWS9bH
R4fW8p+uDHS2soxPJ+tnhL27lF2GgeqHyik0XwJvF0yObyXn7NNLA6slyXo7
FBdfZvWDjoeLPuJHsj6ffVBp631IR0Ert1hDDwL3bc+1/HOVjl52Yh2zbgQa
dnF6yCTTUSt1zfp8st7/rKMZVe2no41Rd32nM7l+a466x7wfSulEJleQ/TQ7
OBr3ddAPgyaq/bjtyf5/0e2d1OyH911/ozOFwNPtgpumzvshz7YW9WErArev
fWFgq+aHiq8F2IJNyH5pod+wEPPDUOJTZQzZ72kDrU+pq/ywt13UPZacB69M
vk+mvfbF4pJjAu56BO7oWVApj/DFdrHf526qkXn3zb4Q+MAHH9zSf5cuQeC8
7PoLxkU+GPvjZ8CKbQQ6va7VlknxwRVF3JwBYgTqzTgFCbv7IEUz87CUMIES
q1K8L67wwXUP2cwUeQg0SVHaomzjja8vSslfm2ViVrx207CmN3K8/LU27QcT
Rbin5i5JeOMK/DUUPMPE6IDhbzDnhRvy9p4TnWCiCmf0H5l8L9SKSLgp0E/e
D3pgLuGbJzanvJnMRSbyd6yWC33vifYt4/6i9Uw82FGpFtrsiX/j2tiy7pHz
e23Fs+yLnij8xCMr8hYTO+9NRrPMPTFI9VLnp6tMvLbpE2f1NQ8UK1DtSz3B
xAOH6UdWdbvhLXrXRjNjJr6vnMjmL3VD4++3OzmAiXerqn20j7jh7dporjod
8v71RfLXC1U3bD5qLb1BlYnCKSempXNd8W1Sp+ZJcSYugRtZhy5418MrOmKR
gSNHJmImtzpjHMPE6MFNBraOM9y//t6P9VU25YGVDOw3MOrnebkfn/8o4BAr
YaCfhhb98vH9OBuaVkNcZGCbxEXNNVNO2JP00bz9JAMv8D2n/nnoiBWwr6zc
mYGGv4UqW8Np+O+s+0PZZTqWTIk8kqTQ8C5/barNHzr+5ho9n6NIwyflCedD
ftLx3jEP5+fjVPTkOeVT/omOsYrpacHOVHy315TjWzcdZcVvegSa2ONYeEPf
y8t0dPIYqfqlRMGKG20TSuTzTuIVinL2Bgq+x1bl/Rp05JGqibb9ZYtzfY7r
WMp0VHDwV1yutUV5w44rT6XIur9b99QebJGHvahLiZeOG6pupFwqs8YDl+Q0
Csi6NwPDTGKTBe4WtbD5E++Hqa3Nh+X8jXFR5LBzaa0v0lNoa1pVjTEZGIW9
1b64dah4z9l/Rli4JTthVbkvqub8Es9PNcLyxYPrXC/64oHggeLi24ZY7Lpr
w7/DvnhNofGBE7cBZn7jqls098XPPQv5E7k6WPhje+bO9z5YytNnECa+C2fY
9yhL/vXGrT0vFBUjduJxfsn8s7+8UUvedca9XRFDcrNKl754o6z6k4KdkQo4
sbX7cu+wN64/6nWK3iOLd1bvcD7R7I1mv9p9l4ok8asAa+XLs95omcVRyKkl
gHX+cpH6Mt6o5LC4s9CJD1L/OWX0Ur2Qa5GSnWK4GySfFG4frvVApsVa/ydh
BhBpYPVO558rhqwfZ1lJmoLmr7Au2R5n5O7yu9dkYAkflzKCauMc8WX8I64J
e0tYkVp+5UmQI8o/yhBa8rUESWOKWLeXIzY9vXR5MdESDnkkfZ4yc8RLqaX7
7F9YwttMjhe7BRzxwvzqeg9XK1ASqq5bccMBb699MaHvbw2c6WObModpaLo0
ITu5yxakd523DTCnoqPjvrjX6rawVyk2rXAPFX03fI95sMcWVq427nujTMWI
tqTTkRa20AhrzxhvpuIaxeu7q/1s4U+Hho/ClD1mnJUy2XHRFlKpksmKKfbY
IZoZrrmCAoab/s3/eGGHS0Ow8+QTCny6Hvgi1I2C+tQf9ytaKPCEeSX4C4WC
1wQfNrV3UcAk5r4scy8FY7YpeP59RwG+9dKfvBUoaHZ7aFn4JwUotU/i/RZs
UUptlGYrYQdeCbOply/Y4n833p5/EmcH9i97Cu9dtkFPSsK9vRL2kHf11Ex9
kiXuEw441LnDHlR0P0y/CbPE++xJP2zl7eH6+vnk7y6WeP46S8Nitz2U/6gt
l99liduS923hMbEH5a3O+2peWeDA+7ufnzHsIbhxiHedmAWqvrvuLFVhD3hq
T8K3anMUTpY5L1NlD52pfSE6ueYYcfZNs3SNPdz/dGvk1DFzvFm8UXDzA3u4
1EFTlnIwR3CuSutqswfBuVeBIX/NMF/EyyRqwh7qEvxoNiZm2Hpaflvfdiok
anxxEB4wQcHvexMspamwotC6n/2ZCZ7Xm5Z5KEcF5Za3ZdPVJnjc/h8jXZkK
WoKNws3HTdB/453atfpUOJQXJnVawQQHDi8W0/ZTYXoLfdY0bi8yxRN1lJOp
UCPnpndD1Bj/WBUmKaVSYXVA24lHnMZ4zSTSUuEsFcr4Rma6vxlhRdUThU1Z
VNBQoEjPNRphZqJ9QssVKhiq83DYBhhhUu4NasI9KqDodxfNR4Z4x+QOm/UY
FTYH2Ni+8jBAzY2SWgKfqNC/rem4vLkBlu+8w3w7QQUnJ8qH48oGuG7jYq3d
DBVErRI36HAYoPF/s7JCC1Twvhsieucq4MaaLT7DvDQ4ID+48LBTH7+zNr2N
VaMBceKj5fx7XTTN3q3fqEEDE79cjulKXTxVuYjc2jQQfVJc8iFeF2kxh4RP
6tGgiHN/zeutumghf0nZzpQGq21EY1566+B/O3X6lZxosNN/+1rJaS1cS5VX
440kj/fKJGj9kjrWrTH48TCKBjNXv+S5v1THNsqcFyOGBjYr5cSritXxY2eN
5514cruDurydjTo+b/f6KH2cBmuuOc+VFKjhFbGnp4PSaRAb5rwz3Wg3VrCN
NE+cpYGD8MtbP4R24/qVw489MmgQ6SBGd5hWxQvtxd+ML9Bgl80CVfy8Km6q
EYbxPBoY0p2XGsdV8OlCnnBKMQ3Kfj21MEhSxgc5pV6vSmjgZrkt746bMv5Y
yg8RLKOB1eew73Kqyrh/pfjwuQoaHKxrLhXuV8J7KlvW+lbToK4uYWFZUQnf
zKiwB92jwWem4rv5q4r4wbiq68h9Ghg9fPN9yVIRd4dO0s8+oAE/W3IT15wC
mrLJLhfX00Ar5MKZLWYKeJRwaK94TAM7I/zPY1oOuXP45EVbyeNvUjHeqSaD
5XbrqwXaaKDR+XLf7h5pTGri0FvTToOnQz38umHSWGgWHjHZQYNK1ZDntjd3
YK+PvW/0cxqUBu1USFWUwlL90j32L2gQJ7gcX9ghifevakjIvqSBKkVg/11/
SdyQe3p96ysanPtwJHWsbDvOOm3L6HtDA49q9512OySQrn2sN6OXBv07vThC
msQx6c+c3L63NHhUemLHWR9xZNvozFbRRwP1uVzD11e2oZ5vWcCefhokG7xY
8VZvCwaGaqu+JS3//b6YbLwY+byivzpogAbGxam+cQ9E0Vy09kv6IA0qPBr3
ymqK4OHdXmevDdNguZevLldZGAcVKye2jdBg9OXP6r9BQvgm5gAtm7R6n4uv
2w1BvFaD9NgPNGCsOfNMRl4ApRIGeD+RDnerdUpnbsSJ/Ypt1qM0aOdRPfin
hB+HgrclbBijgazn2HC35AYUthWNDCStKp0gtcd7PX49/e94M2nOdKWaikI+
1GjkGiY+0qCn0nJEeA8vFu6OV6wlve97DtfivbXY/c4/fZl0esY1q1H1NZi6
RT3r+DgNHi4ebb6jxI1981Y6j0m/+bAJC69zYnKC4+wf0qW+f/zS5FZhnufB
pp2faPD4zrH8hJIVqJJUVOFOenFk874gSQ58JCJYeoa0+rttVp6FbLgy/PX9
O6T/+/BVoEx3qUFNfnb0HWkPZOs6q/63YZb/5PZ/pL0ztBPjlH43hOlkRW/6
TIOjIe68TLnZhoMTauMqpIUsVCQdJX806BcEMM1Ic45NSaTofm2YFjVb4Uxa
VeOdVqDS5wYpx84aOunAptubKZIfGuyOrIgNIf1Zl72YUHrbcIF7xul/v6f/
MbX6MVOppeH//97eNYld5IiccP3/AY8/O60=
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{3.498745062889016*^9, 3.498746069803316*^9, 
  3.498746344465124*^9, 3.498746675388628*^9, 3.4989781276549354`*^9, 
  3.498979002176935*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"wezmy", " ", "wartosc", " ", "ujemna", " ", "Etot"}], ",", " ", 
    RowBox[{
    "ale", " ", "jeszcze", " ", "bardziej", " ", "bliska", " ", "zera"}]}], 
   " ", "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}, {3.498745127605414*^9, 3.498745155405459*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "kappa"}], "/", "dist"}], " ", "+", 
     RowBox[{"m", "*", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"dist", "*", "omega"}], ")"}], "^", "2"}], "/", "2"}]}]}], 
    ")"}], "/", "3"}]}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, {3.498745189063255*^9, 
   3.498745189783484*^9}, 3.49874626710753*^9}],

Cell[BoxData[
 RowBox[{"-", "1.1052832319206307`"}]], "Output",
 CellChangeTimes->{3.498746344592126*^9, 3.498746675456633*^9, 
  3.498978127737935*^9, 3.4989790021959352`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746273144058*^9}],

Cell[BoxData["0.816496580927726`"], "Output",
 CellChangeTimes->{3.498746344607316*^9, 3.498746675471027*^9, 
  3.4989781277769356`*^9, 3.4989790022089357`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", "0", ",", 
     RowBox[{"2", "*", "Pi"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJw12nc8le//OHAVSpJEZCWSvSOp5IXsdTiHyN7bOfbe8xTJipK2VEYlUuYr
W72bQmUURVGSiMj43p/H4/c7f5zzeD7udV2v6/W67us+j1vYhWrhvpGBgaGT
+Prfr8GW08v15yjA8P8+wV/Dzz8m3Lc2WUMhtWvsmVMPqyWcpu6XY0l6qzES
Oni2irCElFdiUsSIxuVzJpM3CN96FrrJhvRFo9z6bHwRYbk/R+OD3CY1Sosv
kemET900EjkT8VMjhT/IPZzwfdUgWavOWY1s2qZaF8I2aqXPPEjzGifZDbWM
CFe0fBKJ/LCgUa1/lE2RsBGXnk+W25LGLYOXu7j+t//M8rvrP/9pRPEv2s4X
UOBxstDNxog1jezrNz++Jlyp48fol8gAPtdfXqkgXC+5XyC6cwO8sjMbtCEs
prXwo4TECOnjyRQZwsd9lYfvFzCBRsr01pV8CjhP8V3r+sAMfVzs8nmEx0QV
dy+5scDNmuKLJwnX8pxi2HVnKxzrZNYTIlw+OCx9YIYVQrfusriSR4HTN54/
DovcDsVdydV2hD99PUe+2MwO7qxeZtyE61tPmo0l7oB9t/iPJeZSgDR0e3B7
Fwf03H70XSWHAq9MqdFR5lzAuMNdZugsBZRKt53lneeC1CNxt+IJnzvR7Ft/
bhc0NdhZNGdTwKrywt0NQ9xwfyvZi+8MBTSvzX4v8+CFrUYL8XeyKKByd7LB
ioUPRMVrXqsSrueKfLK5gg+K2lHEMJMCV47dukGb5YegUzdCTU5RwOyQv6lj
zB5wryi37qJTIP+9pqmwkBCo5imGHCM8dpT0e/yJECy/6SOLZlCAnfPI3Ybk
vSB6jPHui1Ri+5DOfoceYcihGJWUJFBAVCVff5oiCvOZ51sH4ylgcPGqmGal
KOx52n+Jm/BU3vZdhUz7QV37YlRSLAWESF66Bg/3Q4CoT41SFAUU+kZ/d/GI
Q+b70LfzQUR/928/cHBQEny+8h1aC6TAwawy9bsHpGCCMmXESFjH1HtCKlMK
CtcFk5ipFBiy+qkipS4NhSQu/s8+FJj49qnD4LIMtCndO73iTIGH5UNBXG7y
gAP8ta+cKFAwH3a6oVEefq7qp111pEC7mPW36q0KwOaWJ6hmT4GSTLuNLq4K
INlxhB2sKUBWucf8lksRNOWngytMKPCF+8th1ggl2FYpkbusQpyP617juIYK
VLx+fD5WmQIcbmXtCd4qUOjWfH1ViQLUutlywTwVqLQ61TgvT4FjU0fE7SdU
IGqh6WezJAXW68dV5zIPwt+zlU2tAhTo7p9fifugCtaHWNaTGChgOD+1RT78
MFBWdy7fWCODRuGMGNvVw4AmRz+1rZChNnpz8fTTwxAmddPn718ymAZ759fu
OQITuXlOR2fJMCowouzSeQTaH7HulB8lw52KtFBFbnV4k33Y4HEzGSpnqSdb
6zSg+G28VVIjGfKPq17fM6oB/aRqc/16MnDx0F4eYwIQreHZ/qyWDLe/Xhi5
KgfQvCiSUl5BhhPnKbERiQAW4tmnBs6TQfScDMlbXBP45wQ+vQ0kQ45i6pJ7
iBYIJl9vUKOSYcP35P+06FrAdEjErdiPDIrb1pWEL2mB/r80HktPMqxk4pOx
Li24lM/x7aodGVpahngj+bXBhUlP5oouGar4p0nzrdqg0nX/ldhuMjyrj3up
yKkDL4X8/Ee5yLCkMviPTUoHltpddxdxkCE3dFX/B+gAw4J+7tJWMrCwSJjd
DdABFy3BicxVC5iq1jI1e6oD3o2z6XOjFrBP+KLSUqIuvHmiylJTZgEj17jT
meb0YKmgtSxfzAKSrynU0IcModji5GiWiAWktA1q0n4bwu608aPJeyzg0815
JZstRvDP1sTChdsCnAxu6ykpG8HFBZkns0wWUHq2J2RDphHsKqKtkb+Ygx9T
TcWuo8YQ78G2xnPJHDxeje4VvWQCMipKPWJM5lBN900QVyFBiXNJVfk6CRa4
Dhg5a5JgfuWsh/QyiYhb6H8XTUgw/C5Mk3+GBOkac4u7PUjAJHDQsPEdCYTn
knjFCknwV0Oyl6GcBCs2dU35SyRw/HvIJsaIBA9sm0jPW4jruitVTy6Ygsgo
e++yORma2Hde85g0hUQpK4YgBzKw2tzoHxg0hVrVMunvPmQQ/6j08xyawuau
qcTPyWRg7jNcTKObwp+iLPIPIk8in72yEuc1BSW3BnTko0DDvpBCOXkTOHP5
VtTAFwpEOF2ZdybicrijfNA6wQoyZtT+Kk3qwkOxoc6C8zbwbiTg0GtVTfig
FZAzImQPnWXaN339VeD5SQ4HYyMn6LyjuU+mnAcaA/bY3OJygXyjwlsSWbLo
XpfsQxdyg8FHP+Vc1eRQyVYjxF7cDcaq7rVKT8hhgfWtJSV5N1D/fTGvZq8C
0nLKQ8eOucGZ++efmscrYnf5D1GygxvwC5op6xkqY7dz3WRUiRucZ377civv
YXxup3Pro4A7PIt/YCtnqY1uhnbJL4Q8oC3h3snqM9ooj2pB5RIecES95rhv
tzam9Pp8oCt6AJuKXF3M0eMoxy7FbqTtAYIZhoqi4jrIOirUO+LhAe/knMXq
mPWw9RHXnEGVB/geJPWl/zLE2BR3oz3HPMGoX//+8WZzvMymMC3o7AUMMl9P
fRs1x4ja2528Pl5w7XFXZiWTBWYFNJlxB3tByh3mvGITC+Q2y3fkSPUCbbaj
oawfLTA4S/ff9tteYDpo42qzhYKCjt57RWe9QMA1aLg81Aq9crkE8hK9Yc6c
Y3d9iRUeY/IsKD7lDXd3SrL97LBCvbAl3ht53mDgkdvaxX0CB+xkZWpLvWG/
EGOpbv0JHKwQbO/v8YZWpynugi022L71I8POnT7wzuV2RW6bLb4eP1l98JoP
1JIjHIunbdH22g21ojs+oGKtY9PDY4fVvWfX/1b7wAah9aoCPzuMPCwk/ajN
BzYWDrfP8dgjy8x7VZlxH2B0M98wHuqArbvg2hdxX5C/7C35oswJb1zJvupQ
4QuZ/PL8hi1O6Hhu83xRjS/Mdc+sfuh3Qp5XwvGvG32hcYb8+QizM+a/2sKr
/twXIt91GK64OyOLuYLR6rQvPJbkupAh6YKC2frvdsv7wdG5w69tnriimPDI
hsC7fuDLL1XQO+yKnOdVy1ge+YH4Yf1Jr2VXNH/bl38J/SBBN+7OxAE3nD14
ZXPTaz/4xjDvr3jLDXff6Yl8OucHHqLW63zn3FEj/2HUblV/eJTk8HUx3RNd
tSX01x/7A8N+q/fctzzx4XutA2tP/CGncqqc0u2JbxbMWBZ7/GEh/769LIsX
NmlI7el/7w9elNwRy0wvjOEcCtNb8gf6wsZLF4q80cCPXjB2MAAWH2RVlff4
4oeJitvGFQHwLcPzP+p3X7Tb91JjrDoAWjK/3zZn88Ng6rdu6uMAsPdd+RNg
7oeGzIuUoM4AMOON2rZz0A8pC6n3dn4KgBO+HE9Y5/0xjmJ4/w4HFWqO8ZrU
8gTgzXu3fjzgoYLkTt+B2CMBmLXKk/FAkArcf5QP+CUHoN+sbXCBJBWilj9y
IxcV5ZJ4BUeBCsvMsaQ7e2i4YMfHE+9PBdlOe6doRRpmt7XCaBAVTIJPJdgd
p+GHjICNhyKoMMSbxmnmQ0P/XN63rUlUsFSN9114SMPRckaSYSEVUs51jOs/
pWFc35UPJhepYE/Ou18+TMMCGNukc5UKctnX7a4xBmLTBs37bOVU6Fbf2Ush
ByJfSM0YYzMVhLjqKqQ8A/FArvPLq61UyIOFnRzRgWi95G6h2EUFlvnmK+zX
A9H1WWuu+CsqaEzrR7+dC8TbfpZK/qNUYAzhz00qCsJPHB7jlxhoYES/cW/f
3SCkpPzce3UjDZ58Gql81xGEp2+wS15gpEHyVm+7uLkgLLs7q07bQoOzSS37
H5sF44KOYvdtdhqIiMgOLHkE4y56TZ0jBw1CNbNzTOOC8eHDLxWsnDTI/PVN
UbkyGLfTQzV1uGng8mE5/RdrCD7hCFY7KkCDfR5Nk6/3heBJw7rYy4I0mBvW
onYdCcFaO9KfxT00+J4/Nf3TNwSHh79yZAjTQMIrcWpvSgjGPN2f0CJCA7ax
JxyeF0OwJ0K15vs+GiwZ//os/zwEn5parQmL0YCldT67ejwEtWzlh6TFafDD
4bmXwVoI5v0IZpeSoMGOg48uP5APRYXxgU2bpGgQHn5yJU0/FDmUzU4PE5Y0
+5kU4ByKhh9H58ulaTCd6rw5Ki8UJfi2Lu2TpcEj6t4NlypCUZ4vxfo54a3G
WQL9HaFIp47t85GjQclYrZXIx1DULtinvkzYiseoIulvKG58x5MWK0+DoaTL
vIscYVjp7fxpjvDmLTcvxEuHoXV4iKy9Ag14PX5JCeiEYbJmoeljwtLdlR3/
OYTh/YNjolsUaaBnfts9JyIM7XdAjgHh7FJ9Jp/cMAz67R8SS5hKO3XJqiIM
lbL+VV0nfPkro/SJzjCs1zQSbiTsYPTluu+nMGQ2yC3v/N/xrpMb8pfDsG74
nWwb4ddLnuqvuMLR32s5/T7hHyQd073y4bhSVXo+h7DA+QXpFINwtGK213Il
HO8c//yfazhapoc5ixHuZaZJpseFY6CuQP8Hov2SQl9UxM6Ho8cZ4/gEwiyf
E769exCObwU/HOEhXJ4XfPDyi3D84+k3dYmIR9gIlStyMhyD3sV78xBeUa8K
82CMwJN8HoUJRDw5rKdInkIRaOSia/GBiH+ijfO56MMRuNYgGCFKmMvsldp1
ywjcaNsy4SBDHB8dJDtCi8DZn2+j6cR4MlWtu8hkRuCOfyNcN4jxlsi+8yyr
LAK9FFWyqyRpoCKT8Jc+EoFhW20v5hL5w536qfP6EnG+a3V+gUR+dabzv3/J
FYm7drXLwH4aXNvuk+ZgFImKORV/Koj8zNPqXUOPSHz2lr3BkMjfJlJLqEpS
JDJrKU0OCNHgRNc3TsrjSLzdrLPYQtTD9Lo019rbSJzcUp7Lx0+D+rDl9se/
IrHk2yM2N16inqT4v7hLROETU9PdbbtoYFFkUH6mKAq/+DDcrWSjgejBis62
mih0LlKKiGalQTqLzCDL6yjUzJreospCA3dZ+4xulmhc8mn7HUbU9xH342vb
YqLR14RtVXqeCoIn82Y7i6LxjUnzwYs/qYDMEutna6Mx4I7tk/FvVGh2N5ym
zETjadOWm5uGqWD31COzyDUG1zaNPHRpowKzU8TzlwkxyJg4GVvQSIUwlUNC
3JdikMfnjvHlWioEe2a6vHwXgxdl93w4WEaFFek1jV+msdhbdfmifwYVDgv9
HSv2i8WzkeynSAlU2GI152B5KhblxSzlthDzq8HSfdXZjli0OZoaOO9BBXmG
p0Vd6nH48JtiwyMtKvTZXVcZto1DqmPSITU1Yv7UD9Naj4zDMGnp3ZnyVKiY
kf3m+TAOJ74rvi3jpwLXRtMPHfLxGHCqd/fd2QA44LbnqtbOBAyTSh9uyA8A
6RXOQxf3JqCJvtV9HnoAXLjsVL8kl4Ces1kv9WIDgGMbJz8aJWBFcs3aPvcA
MLFap4anJuD2yzcfdygFQOObHpuIpQS84sfVH93hD6TV0eVNWxJR9eqPVYY6
f/he0d6Zw52IJyu7+6xu+cNMuJLLgwOJ2MewUBd8yh8Uca8Ua0AiBnxNNPAz
9oc15pjBmbFE/Hykf0ayyw8eNll0n5klzrdtYPhKrR+4blZtlWVIQgfOTwrj
1/3gUc99fV/BJAydKtLqi/cDCSf9mE8nkvCsUeN7ExU/6I/u35fzXxJuMnV/
F3HOF+r3GFiMP0jGjI/G7EUKPtBQkePn20aYjRxvx+0DnnYzht/fJOOcV8GD
j/+84THPlO3AbDLmCvFPCHd6Q0f3q2kb+RT8IqQgnXvCG75UbHN8dSsFjzDG
520L8ALn6+btb+tSMEJHbOCemRdIGkoNvehMwWLv5Yidil6QdMCN4/qXFGy7
a/mAbc4TuI0F19uEUnH0zy/3z8GewOfgpXUnLxV3dZ3O2ebtAf7NGhsmr6bi
bU0qT5quB7T7ajLx3iP2d6CpNu3zgIh92x8Z/ZeKPIkVBW7D7lC+/W7j3KY0
DFr/7HXBxB16frJas9HSEAKdGZqF3WAimfxEPTYNE+Kj7+355woTTdNPT55K
Q+3PbmIab11BbV9drs2NNLy+L/EQprrCJe4nTc4Dafj23JhWz7gLtD2enj6p
lo5Hin1uKlxwhpL95aPex9PRJZVVMT7IGQTbf312NUvHd0/3Jp43dIYzb91F
d7un40X3A3lKy04QIT+zIzo7HaOz7J+tWTlB1StficWP6ShxMMw0e7sDKMoN
WhdPpmMt58L+X4/soUZZ67fYXDpSSm5LzLnYw2KMQcIcUwbauv18NVZrB41O
Xuu2Uhm44cQ5SvwJW5DU/ZdyJCADDU688fc/ZQ3+cV/eHw3LQLwe/N5S3hrk
o330ROIy8MdMc+TP3hMgcYVyqygrA0/KehorC5wA+bkf2nLlGRhOOZcnd9sS
yK3xAkbVGShb9+VrjZEl/H7Ekm/4OAN978ffmpumwI/HclW/OjPQsnyx9ZIi
BXQ7u9++/ZSBO82iDZ/WWkDeyr+JDV8zMPQzzwclKwtwWF86vmM6AwNjw7vs
FswBpk31Xv3NQC+LCM8lFXNIrfDWfL+djiW8G1em3czgw30D+jYuOv7QULTk
uWYKg6fc/wjx0jHxXPuF/hETCPiYVDMhQseNb48N7z1hDJoiI6RfynTU8td7
LaVtAE6izOHb1Oj44L7HBdlEfSBygoNFnY7zVXzfh5v1oNI95uj143QU1FBb
Ianpgib/QjHo0zE6zf3+7jAdEL8gy9ZsRMe11+EDhdXHARhGnlpb0NEijjup
U0obTrw2/xFqScfzf59zRntogd5CdE+QNR2zvbMUZq9qgmXmHW9WBzoe+6Sv
LbMbIIK6fUuZEx0/OkvX/Sg9Br+HJ3eLuNLRZqQ8bTe/OiQYNJo/9KSjkOqG
ZiPGwzAYctH2pTcdOabWs+wjDoFhSHdxly8dnzX3/Mf/4yDI9eoxXfanI5vx
m7rTjirwaevf05ZUOnrJfc2ufnMApAyzRL/T6DirsFhUoKMEX+Ty2h2D6Ggm
95738CMFYAryd3kQTMff9vsfGVrKwXCq9Z/xEDr65R498VhNGi6L+YQuhtLR
M9+18qmgBLwcyhn8GkbHRYsttvQN+2E18LxwXTjR34tijAzjwqBxYr+WewQd
Yx/evBx0VhBaVmsO/SZcvS37ZbwVDygKX9/gEElHt7t9M0cFOODPJrELtwjb
qUXWxbUxgW+JGcMbwmd1lNsDL/7WqH3MduAD4e+BGmyeau9bwlmllJ4QLked
Mna19RZ/ecO1NMJTs/uUKzm24fmrM0XShD0yP4uoaXNh+I1wxrvE9W8IlIlb
xfJhLy//cQ7CIWv24pm2ezHC5KetBdF+hmQVYZWIffji/idS0P/6d++bdXG+
GJqeaZQIJPrfTBuXqLkniW2XxaZMifgYE8vWkP9kEO+UXGQl4lcmmrWt96s8
9q9H6pcF0nHkRsecc6Ii8v/+sSBMxH+qQyH00vYDmKPBUB0dQMcPL6xazxcr
4xqbe9R/PnTkYtKurKtRRePta5fbvOjI75XW8xLUsOp9+8YLHnT81POrKue/
w2gwkff4izMdw1flllbH1DH8QoaMjSMdTyYuLJX4a6AqcxJfpR0dux4Jxk1G
Axox2VbNWdGR59z7h8FMWtjvVxE5QaZj4/Jl05EeLRzr/bz1MYmO/ss14ixn
tPF0Us2bVQM6hvlwRhfs0sHO7Af+W44QxzvustIS1ke9elnatYN0lD90mm9k
Th8lqwrceJTo6Jswys/TaYCmzhxVuRJE/ei3eFG8jVBVY6n7DScdN/PakW/e
MUXB093W+kR9T4dMMfVHmmFjavjWC1voyHzol4WOKAkfRq/HvlnJwG1JjRFH
fpMweX9W8sXPGfienhTCdNoCz7L7LS0OZeBb9S/pGb8sMO4Zy5pCfwY2iriV
tVqSMXGmIupoN7F9v3a+lRAFtaRWR5vvZOB0R1PlQJUlHjBkX7ninYGrapsO
CjVb459QJ+ufThlYeLjxQdgeG7zhef6xgHUGiuxx5giMs0H2ktwEQd0MjIOn
UwbqJ7Fuq/ty7d4MvCPXE+z/0BYndqwqRr9Kx6CzbvX5BQ7I+1D41UJHOqqG
Z7YXzzhgbfLmjBMN6bhn/EqBioEjti8s2PSUpqNIWt+ePf8c0THnsnRYZDqW
Oop87qtwwsqlqRtKAulY76z1LZ3BBacH8th+7UhH9alqVwVxFxx8xzFxnikd
7cKbjt4wccGrzyd7H/xMw4YrA7umz7tgte2m/ZsxDZf3MO9flXdF/Xt9bx0d
01C8bkhyytwNG2/skGQgp6Fuas45+VA3DGcOKTurm4ZV6+3W6kVuGGS5uiNZ
Ng1f/thGuT/shreeKPHULKfihp1C9Yvu7uhK47wSlZuKUzNZE199PTBmsZ+s
lpqKWxZH60pPe+CkqfLH7+Gp+Ego1Uuq3AOzFgVdZe1TcbrXUj1j0gM1pK86
S4inItvfnvf5rp5os1OhUP5xCrLsU/OeMfbCTBux9wwDyfjh9haVd/PeaKDy
JjqwKxmr1vySk9l9cPJvXvtgXTKOkXxWfkj6oByLgealomQs2iWytGzvg7Lx
F4tHbZLxAUen4uY2HwysKRgpGEzClZnjakyJvhg5OcPe9CwJTUoObrpZ6Ist
x97kjjYk4Z0QharNVb4YVBrvse9iEjadvxS1/YMvckq08p61S8KSVzaT8vJ+
OEJ+cfHwcCLGqt3on3/uh++bnhmaPU/E3qaXhs2f/NDBQy7YrSkRGzhbSMZz
fvjclvH72ZJE1LPV313H44+TSui0aJ+I0m77x1/Y++NuKF0R+piAA5brbdof
/bEwKb+b9DIB82rntMJ++mOWOWN1SksCpio9qAhc9cf9K1IhC5cTcPBhq0c7
XwA+8ul+xOScgI/WLueNWgSgmfmJ8M9d8fhx0KXh6OMAVP2grmN7Ox7ds6U/
T7YHIG/8IfnhU/H4av+tBI9XAch3C9UYTOPxK+dHzZcTAWjw54Bkz9s4rP3l
89WDg4rL8X+HC4n1tGTnQkchPxUbx9Ttg4riUPvXFs8L+6l4sSQtycguDsfG
I3dzq1FR3eVbRd6XWGSw2E2asaeiUZkbQ2dnLPKkr9o89KCi7u32QqbbscTz
kkaYOZWKYhu7NSr9Y5FVOJ+NKYGKI3E/t1xfjMFP8wNPvS9TMdxQaNfxDzEo
MSkV/f4mFR/l5lbMNcag7Lf/TgtVUfGSWufO+KQYNG8MmZNvpOLAYNsdw+0x
uMhmOeTdT0WPnAJGt9lonOlrzc4fomLw9QMbs95Go1w96XvhGBX76r96bi+O
xsW9uxolf1IxTshM8Y1ENGZ9X3mNG2ho5Bm00W9bNFaN5jceZaTh+JidG/ev
KNRJ0ee7yUxDuZznx0vqotA4Fhs1WWnI8DtNPUwvCpv0tSKucNKwt6zFoVA6
Ct3Ta/tP76Khf1l4Xw97FPpG/EZ3Hhq++pLM5/ouEjVi/7SO89Fwq9i/nU8b
IzGWniyUKUDDF7OHa7WuRmK+nE+O8B4a3ig6rOPtE4nvHgSo8wnT8GPmfWse
s0hkBFJDrAgNA/xcnAcORGLuqx39r/bRkO8ec1HWagSWCGh46YrRMF4i4mvy
WARWBNiLeIvT8HOfCDWzKwLFRRrWoiVouPHsabe3OcTz88lt14KkaDjjHCTw
KSICq6W6XuhJ0/Dg8zOf/jlG4MArrkvsMjRMwc8N4noRaNutK/WUcKju7seu
chHYsZ6TFS5Lw93nIsbu7orArCvbx3jkaPhPeFFz22o4Rr2r0ywn/L1UdjTi
Szhylvp2KsrTkPuAQ8+fZ+H468oUvZzwi8Mhm5IfhCNj70DpbgUaXitkKtpb
HI7N8Y8ORBLeXkDOfpFEHL8yZ/yccM4fo9lMn3C0SerdukuRhuKL9A7iBo9/
tiSlkQjLWJ3gPHI4HFv5tr2PJ4ymlaNSIuFoet+R/zphu4sqqtJbwxF3GHvU
Ex4IN+I7+jsM/UsT33QQfpG3ctb+QxgqfM6J+J91ebZezW4NQ02xGdpjwnfP
TZq9uROGzwp+9F4lfIC189r+vDDcsKX5cRzhI90Jl+jRYfjia/jh/7XHllqt
t+YahlsE2AO4CbuMlV9LNg7Dfb+s/V4T/Wk6EP6AVyUM/da2miQSXs9/F4uC
YVjC3yAuRnhS24wpgjkMGY9/5EQiXhah3CYwE4pT9OcSZoQjjolZ874LxQim
xtheIt4lMcqSm56E4lwHSpoQfsA1IsaSR2zvbx/bS4yXCP+1E/tjQrGu1Eww
mhjPoft1Vhbuofj4+uatz4jxV1iK6hhRDcXVqgdkYyI/xkZcVI8Jh6KknuJw
rCQNfwftjKnaGor8vDOBpUQ+pXTXxD8ZDsHx+Q9cA0T+Fadu13DrCsFjX7NN
R/fTiOeg9dc890MwdSx3alSUhgt5HIH3U0JQZlbmcieRzw8rxWOrZEKwfvG4
0V4i/yc3/2kb4A5BZYdP3z4S9XFSvr+PY0MInuzsO3uOn4bKum8sHvUFo+QE
1XxyNw0LP+pTX8QHo7hYe48YUY+iszefx/oE4xBHyMwVDhp2x9iOqVsG4z2m
+JadO4j2NfMK/5IKxtzvehXD22hoULj5tMBAEP4KE329k6j3ZPkZt7MYhMZF
ZmX7ifnAMPs58twJwmcKX9zkNxL9/5Qy4RT7v/9bHb6ZrVBRijdWdUg0CFP7
bfodv1NxXVOW9n17EPbd8GBdmqBi182NC2xLgSjbLBmXSMxHlZ+PPst8Hoid
povPnd5TMbvQ5QktPJC4vzkN2XdS0XLvDJ3XJRApb4S6Qp9QcVb/gmOfcSC+
XvA0iCPmv5B+MnuaSCD6mNs2OFdTsUqBnvz0BQ1VQwrehJVQUcPkjtcAMW6J
9SIPV2hUFHkvcK2Cm6irnJk1BV8qJnenmuQT/WB2M7Sydqei47cHl1Rqqags
39OWYEPFH+IPLjDzUJFjDyX7qCZxfxh7cI826o9XB92vxm2n4p2ll6ktZ31R
aodeL+fVADTZ4TB9hOqLY6+/L7hfCEDHnRs/95n4YmBbtN29vABcvNAk6sLq
i6Y3V/Xl0gLwnJJuzft0H/zDdn7Pa58AjLea9VNK88b+xP6hDqUALOcILJ0r
9MSdGw+c/tfkjwM1N33+RnqiZY3zn7iH/khX0mTbbeeJgjdPS65W+WP8N5O8
sr2eKD40pzJ12R/ludw0mIh1yX/MlU8yEv2x2Kp+vazdHZd9/x6b0PbHvcaC
f39vdcNTmjWmZzv88OVa7c57v1yx7pzh/ZomP2wcsXFL6XfFYr7D3/tq/TC5
mmoXdM0Vf8Qw8bDf9MMo6wNU5yOu+KnjRKl9mh8av7xRExXogh9vdfwL1SX2
n+A+vuOnE8qu7FIYbffFLRN8S1n9TriDtfbuw0ZfzDU+E7y3xQm3b/7jeqrG
F493a/gnn3XC1K6abJkbvtjhXe+uoeyEDvYL83bJvhihdFRrbc0Bf/nefBiq
SaxXOjrGjxjYoevWZ9PY6IND0tH1Wfx2mBPNKv7ogQ/WX7guyPrTFufJNmFV
d3ywdsfJkfQ8W/QctI4qKvLB1JyWqJRPJ9H5lDeLa4gPGgV93+eWZIO7Ew6k
vpf2wadyqYwFX6zwoPnGvMQL3vju0rmOBV4LND+ysUQ6wgu32rYXk/6ao+Tr
za0NVC/cy6LQ2NVvjkND/vHGnsT6budhH7Z8c2zpdXpNtfLC296OHnPs5hjy
YUixWdkLL1gJ53KwkfA/RsuOmFlPdBwtLn6gaowGRsonzvl6on7OZ3f2bTrI
KCqt0uPigb9PnL9x9PVxXBaO4xm38UBplqKvl88dx5JdSZ4bzT1Q4Wbb+c0i
x/EN831rTQ0PjNj+o9DyiDY2Y9rDV/we6D3BxNEYroniPhZrB/vccUPiHd0E
maMo0/xaXdaAWO8aKQYdiJXH1u2KvzYfdEO3aZE8SeMdoNVZdnJI1wWbrOoy
+hcOgHfXp1j3SCdUCqII5zBqQr5MZYqCuj1ums07Z++iC9d+LP4jz9mgUe8z
7RoRI6hccf7Nd8YKr8can5CZNIWb8bw1ZQ5kTGXIudvzxxSW700MM5LJGGvj
0Ge30Qx81BnaXPXIuOLIdd2N3wx2aF2pF1Mg46nzTH+Xjc1A5YYlPNlAxhbt
rH96VWZg3KxToVdqgYps5sOdyiSQJYnq5k6b44N3l2ibxkigfVdCbjSJhPWt
VuScryR4avWZrTqMhGeWH7DyTRPem3AjxYeEzZvGWQT+kqDg2+wvOQsS7qze
kjfGZg7wve1ToTAJ2+I2psQcMocF5T2DwYlmmCiWTf9LN4cddgFVPgym+E4l
6Fr5HgvwddmbFbfZCNNDEkobRCygvPT9o/BZQ5y3XG7pErOASZ1w0+BBQ1zr
eKb1Vs4ClnuZ7wbeNcSn5bdqmo9ZAHOaCOvpE4b4PiDg2wd7C6Bc43OQu2WA
H/xZBNTPWwDX6/AVCSN9dM/tzeu9aAE+rCwuvir6qLA8Yetx5X/bs2buCumj
smnc54QyC7heQorSmNfD0ATONXqtBaz8+7YrokQPsz9Xdk2/sgADNfkrJr90
cUOgfgIzExmCOgd+RBfqYB/rNwO+LWS4BMfSfJN0MFWgvE6KlQy/mpzn7fx1
cBv/TIomBxnKtTz1dLR10DPBllVPkAw/3innS88cx1W2kclqZTLYpR6XMNE7
jq7vHt/hdiaDXOlVvYG/Wuj+4mvAE1cy5DkMf1r9rIVMl56teXmQQTS0nlPs
pRYamWikVfmSoaCQVSO6VAu3BBmTOcPI0HltPeKQhRb6bu4aOn6KDPKhX8dW
yjXx7cMkbrN7ZPgTc/vLFSfAMK7ahs5qMlTZlU6oGAI63RpvPVxLhlCdnLT/
DgB2yu1d4aknw/xNnsI1ZkCf6dIjZW1kWJflDPt+UgPT0jcv8PaTweWoiojd
H3UcLu6Wtf5LBkVDo6MBMkew40m3ocMyGSQ9l2cX1g/joIvotNMKGcLvfLBI
6j2M9fYjHbYMFIis1GO8GnUYn9efGhPfQgEb77Jbq91qqFxkfZGPhwJca8o7
FDwOIVV13ZymRIHYvWlw8poKPuFbuT19gAJx+srt+yJV8Ofaz51eKhRQ5rEU
mjFTwTovFT7KIQqUh1wxOb2qjLtdg9lZjlFgaZUk9t5aGWcLl3m/6lMg9Mif
W13sB/DCX3qPiiEF3A2CFismlNBMPCwv0Yi4/nfLL3lNShgVkBuxw5QCz97x
nfbwVcIxvt5zrGQKSMYp397TpYhLX0UKw+0ooLb1DNtAnAKeLdKQLrKnwKGN
b4veWirgfM/IaK0DBQzzXtJ7ZRRw82jSo29OFEgUyhIsr5VHzgqF3TLuhIsa
k3v95TA+tmiTlD8Fekf++V8ekUbSBqc/AgEU2HS+4UdbpjR2eBXsYKNSYH+O
SOfUYWn0p3xo/0KjEPXR7adRKIVTbjHyQSFE+3mftG8iSeLhznsOTNEUqP/n
Hv37iRhuX3X7PEoY770xeeQhhiz/hsbrYyhAYbRJSmAVw7zS9+4ucRQwdsIp
Hsv9SJf+Zp+eSIFH1XG3nL/tw4Dl20N6SRRgox2KVsvah7EJSUOMyRQweq3q
wqm0D7lt1pxCUijA/k1C93m0CBZleqdLplOg1W3Y0ppdGIVVDvK+IHxMtKde
o2YvngtZlKNmUIApzmhCwmYvxnLc+HqTToGagsy3DLpC+DzH6tyP0xQIOCUo
m+YvgPzRqaTQTAqcumCsSa3jx5PpO3OWCZ9couvZbuDHaG+bm8tZFPh++d7Q
kQJe/DxYK/E6mwJrK46sgk3cmMy2JnDsLAW8dlTl8G/mxv5DF86XEa5i2EsV
MN+FJ6TLZGg5FLgSt3lMdJwTq5/Ga78k7FPylCQjz4nJXDGfpXIpkHrPuEcl
ciey5BW3DBCuPSJ/ypSNAzVh94pEHgXqXo6W25/YgWS6b2sYYZ2tM5zDWuyo
s09lenM+0Z+tP2ylfLdhXe6CuRFhhgLvC5brWzF6UFTzNOHNh5jqkvNYUDTB
sIWhgAJnuhZnvjYwY0xpnOlBwtbbN7QIkpjQWPW8lxdhbbas15ZfNmHd3hq2
QsK6FxYYzkZsxJdJX9VbCbeKnBf7b9sGrK8xWJokbAxXVbm01luatGc0tp+j
wL/dc8F5Pist1nNTHPKEVR6M1uzKW2qxajKhGhPmfCpeeaFhoWXstZi1B+FN
7xLYRb7MteQWJv8XQ/jppttnyrfNtkSUU//LJizDX1/v7/OjRWp0/sRlwtX+
NXbpDRMtvVr81HLCJymNXte3jbYoc0xy1BD2v1qckdTQ31LV5gL/e1+Yp7yA
nNjQ2fL/3ycW6XtdOlfL3fR/awmhjw==
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{{3.498745176425185*^9, 3.498745197367849*^9}, 
   3.498746069909221*^9, 3.498746344662772*^9, 3.498746675583115*^9, 
   3.4989781278439355`*^9, 3.498979002249935*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "wezmy", " ", "Etot", " ", "dokladnie", " ", "rowna", " ", "zero"}], " ", 
   "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}, {3.498745236716957*^9, 3.498745272352102*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", "0"}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, {3.49874529610408*^9, 
   3.498745299915419*^9}, 3.498746284850759*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{3.49874634481733*^9, 3.498746675649102*^9, 
  3.4989781279249353`*^9, 3.4989790022669353`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746291403051*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.49874634483251*^9, 3.498746675733944*^9, 
  3.4989781279629354`*^9, 3.4989790022809353`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", 
     RowBox[{
      RowBox[{"-", "Pi"}], "/", "1.5"}], ",", 
     RowBox[{"Pi", "/", "1.5"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}, {3.498745551346088*^9, 
  3.498745636973262*^9}, {3.49874571772446*^9, 3.49874576540334*^9}, {
  3.498745802994217*^9, 3.498745814122423*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd2Hc4le0fAHCjzIyMQpKdlT0yz/fYpDgPSko4VtlNkj0OJXklIzsiKkJk
NM4tREPxZiU0KPEqFDLzu/v399dzfa7rua7ne+7zXc8jQQ0ivJgYGBj4GBkY
/l5zJBtaJO5T0XmL3WIdasFoTTxpdSSfisjmfqHFisEobnDVsSqJioSzNe1i
ZILRfxbhbMnHqCiEKutnJRyMVJUNqtQ4qejRmuAvJsZg9DXT/5XYcXeUzJW9
rFochNJjay36VdzQsfSrhUETgeioEHec8vgR9LxPjNEmPgBN6DP2M963R2lL
Lw5/eemLmpsa/K/y70MRzx5vpqUcR2KtPDZ0cTLaPyWWse+aN9rMVqA/4qGK
DF11NMpLPVHS3S+eH98yo/f6Y10vGT1Q7xv3N2Fe0mBz8Qh7v5s7ynAJ1Fev
0gSnLvPNI29ckbD+UePZzYYwXzYTMkZ2QW26OgFTVGNQ5DXTirvtjCY+OwXP
R5jB0m8uNxlPJ/S09WuApIQl5BloMvJZHETUnRot169bg22C1viwkAMK3aOr
4ZJpA/zhk/YBuwkU7OX+IT34ALAIzHkpSdohvnRf6wPrtnCDK2t64JQN4h6n
j9jn2YF0XsAuXQYrlMbcUpGoRIGJxYcFDS5maFAx6lKwMgUYzOTOO9iZodc3
8jmcVCmg1tue+dPYDKnYDE1JalKAdXgmS13ODJGYHPaW61Og9sOV7JafpuiV
A0vgSWsKOBw63SmSZIrehrR6MfpQwEf+cMH9+yaoTkeRUMmjQLCY0cx+bmPU
Vl6s8SWfAtzffnuKMhij2fQnmdcLKSBVYyYwPUdG2z1ejq0XU2CLjOOO1D4y
UthyuqyuggJ1b10qx/LIaMycpP3+Ab5fl+lXuyIZzZDiK3i6KaBqK1fqYw1I
w3ywOLOHApGjD+CIAaDz4/JKom8pcCZp77ytMqDF3FOzkv0UaNyxWGDAB2hD
WG/frmEKfOSU6jKrIyGGBrXStxMUiEpcrDHtNUIuZy5W7/lDAVmz8GMvBQ2R
3o0MuZgNClg47nI5t2KA/tnFbP2WgYBM8tpeqQ8GqOOHx6XTzASQG4aiY8sN
0HnfXo6bbAQ8D6SJHtI3QFnz047tfAT49Wz5eNBdH9VU8LlwyhJwUL0/vaRK
F8EMrVRzNwHbpch5tHRdJOJx/8ZROQLcZGJwheii1mnvB+UKBAycVTuva6yL
Nk3HHlZTIaDqsLH0ZN9e5DRwv2PTXgLmOZeW09d1UMJcqN4DSwLa26w3OA5o
o+YGNddSKwK+vrOS/a6pjShfT/mmWxMQqdZN696hjZ7HUbX8bQg4xHj3R86k
FkrxrSvmtSPAqgnUjBK0UNhBHTaZQzjerqekx4800VDy98MRngQcNhmITFfS
QD96kmusvAhYTQ3lyeHXQFee9S0IeBMQsTZ558aKOioxFDcr9yGAUpPJVNep
jrS5fJmf+hKQRs/cGPdURxU/+qWfBhNwYnBi7Vy+GjJk8xRbvEBAqPS4FoVH
Fb3nibxVFU4AdS0p7fE3FfSOyvDLKwJbxql1qFYFSbav9HRH4uenbRjssFBB
9ssjdVkxBBhsL5KpC1JGk2ws2lM0AnKuGzYYtyghzk1Cq1cSCfB/PbsYnqyE
juW6JqgnEUBS9uRoclRCR2k6U2cvEvDjv6QC7SlFJBiduvw9GZ///X+FLfgV
0WefFUrJPwSIicocaPaSR3m8xyM00whAEC3LqyqPbiUmBLdh92e35/isyCFN
dY+qj1cJoAfLWe1IlUPNQwWXOTIIeEgIHc1o3I3McxvmpK4TkJ++l+jhlEU2
w4PGgUUEND9rs7W7L4lGQmhFPdjZFU2mrjGSyKs7cE7jBgHVXDchyFYSSS+a
ByxgfxvaY3z1Pwlkler7xq+EgMTsuZ3jUhKIV8DLVraMgNk7Bakn5Hchuamw
qVBsvfzvzaXtYqhQjivsBXaO7XrOmLsY4hOYOed7i4C8cEkPz5ydKH5zxmR+
OQEa5nMh5zlFEUcwW8u/t/H5bdMcnP4uhE54fbgjcIcAnq3a9eRLQkgmZe9l
R+ysZJUfWbJCCD29vacXO5nrs4ul63bkuMC30HYX/94l0+uN3YJIyyYm5Q/2
zhiGh9sDBFFuwOgOnUoCNgQCtoSyC6K0xUyRUuz/HPXlDIwFUFkk5caZKgJY
lb6V99znQ1MySVOl2GiXd6WuHR+q9eBU6Md+EeFaXzy9FeXLSl7WuEdAhlhV
U4jMVsQQWTnyCfsmu+PPw8o86I6CzxBHNQEdcWPre324UVnooy51bOv9Io3C
RVyoTswlOwq7n9l9y8etW5Ah3fXJlhoCxg2TNiIX2NBgEPcbJew3/jPC9ups
6FzOt4/7sI317kbIBbGi/kxdtkRsI8383/3fNqPE9Uc7S7AXrLiW78lsRiJs
CZpPsHe+sha5RN2EFEeXveawP2SXXyUPM6Fnl2ui2WsJsH/A0icmzIQiSlvz
xLFjP1VyrTsyIrtEjnf7sKvu3g4OLdygb4q1XHbFvqOooKgn/IfObv1L5DR2
zTXdqPWra3QvSiE1E/vPIPuNJNoyffvW8qQy7HtONmZ2DEv09gKl6npsgklS
VThskW6lofquFTvp2CubsV/z9IM99cw92HPaKbSqgF/0WDwyRrDb08o6wybm
6BayhMs37POVzn8s3GfpTL0Bl39iX1wakFF68J2+ECr4ePWvSXpjzyum6Jbs
+2aY7xMg8tyd80T+BJ3LhFeKEzvhQu8Z9rRxunaVt9NW7KMxQou34z/Ry/Js
U7dhn4g5/jJndpgupfCqQwQ77eZRzU8D/XSqwQijGDZ1TqRXjt5NX42NMxTH
Vj1GvffZpp1+LoQeJoGd7FhfPaldSzdxSWv6a1EHuW1t97JJD8nLy3/vb5hH
PK+uPCFJry3q78J2835/tDf0Jal6z6UoUWy10Vf8kn1vSSHk2jYh7DfDGyzv
v78jXfAP5BT4G39ceKU180dSVeVje27smSjHmnKeMdJYV0E+K/Zpj/hGVtGv
pA+uvJMb+DyYv8RO+8hNkkwRr/Zv7KjVF+eea06TSgUL4r9jX/3aH61EniGd
tXvQ+xn7US5vTqfsHIlf8ZDsAPa0P9/I6t2fpCb5yPMvsM2uKIupasyT0jPs
pauwH5+s/J5L+k2qW1u9UIAdUXy59G37Eqn2glhfCnYoe1TVFpsVUq17a/IJ
7FGb6yxxTuuk1K1Gk47Y0EeKQqN/SCu+Upbkv/kRttoVpc4AYU6iHILY68w2
l1oaGUE8XjZgHeevQRNj19gFJjidlNQzjv1dk8zJSmIGjUWjgmpsRY1PJXbt
m6Aosc/cAPvA6L7jcz2sUCEasCKM7RLhrnW3lA0khrbFLuL6s6q3TDl+nh3W
1sO57mLzn8vi+yrOCQVPk/jZsZnX8r8uBHHDidhy7xFc320JpRFPTHggRyY5
pRp73Kc7vo2TF5aG/Q7aY3eRnolc/cALuxkPvL2E+8VO4QQ5vUQ+GBzWbXXE
TpWhM+6e44P4aBptF3ZbbzqH4BF+qA/hTLyH+9GDkeL2eWUBGGXSCmzB/avX
Mia3tV8QFLgzOu7i/pgsXpJ9U1YY3rNc/OqJPTWPLt76Rxga+j8O7MA2cfot
cHdFGDa5nlKjVRCgbOAs0PxaBG6LsT+xwf02vvPZyalzomB50k/xRikBnK07
ooo7dsHL68d8zLHP8gUKfxISB980/zNTN3F+z79amNsvDqOa4Ywq2JUP/Fq3
NYpDi3/3eEUxjt/0D0/AZQnI+xLd7F1IwJHUMRdbLSkozPjmzIBd8J7DbNxb
Cryc5X9mFxDAlfZRPyxbCrZpvd/UkU/ALv8kj4pVKZB4fa5xWx4+P7xtCzyV
hubKaJP4bHx+jVm7pWxl4VzRIXMubP27vsSLaFnoUogxysgioOhsV+OpWlno
taJyF2fi+l80MekQ3A2l4X7y967h/i7NnRc1vBtG0D7NLDxf+951P1PxlYeq
A/5RPNjWpu4L/+XKgxN15xNaKs4fq/Lkii55WHRpFDx1BfcfUuCKnJoCTBmE
KBhcxs+zDlXWXVKARMXI1/l43qe9ph27QVOC3sFD2zPx/lBr7uB286YK5DfL
Oozi/WL4nxTnp0MqUL9abieD/UVZT/0kmypkF7Ek3QvD+4/tG+drRqpwtv7V
VF0IAb4J35wmbqtC+qPrDnEn8f9xNSC3JU4NbFxWLfbg/ajFuHrEXlsDaqvt
tu/3IMBd0MM+kaIB3lvvFPpR8TzxfMHxyF8D0hljh266ETCxwGIuX6IBs9sV
frK64HhmIyT4eTVhfLb99D+OBDyxq2GVndIEPi4aQ4UZATqcopGOhdowFbZC
3DQloHTG4JVmszZwJUdSC0zwPG45eVGwTxvcnlc0XCHjeftZ2ek9hw7oMe7W
cDMkYHKoJjT0nA480PNXf6aJ61No2mFl/15oZbNxyZPE+bB1O6v8H10Ydp5O
c1uhwIfhffzZHoaw8PAeh/gyBSYpQU6cEYbAciD58OhvClx2tzWOzjQEr5Pe
oYcWKHCnzZUj6LkhNFjn1erPUmB1MDXTTcUIx5M88PILBTKNPyxeWzUCxafM
nWx4n7+VUdpbZg2Q5WHhM1ZEgSz19WOKE2T4lyfwoCZ+XxhiDR/SWSCDfvaX
13H4fSKuAyWZMRvDhPn0YbEcCsCr9hjPXcZAv0R+SU6nwPdP80w1h4yB/2N3
NTmBAi0pjUoJncawmWNKVPI4Bc7e9qZJ3DGBhN7nyZbyFLDl0mr3PmMGd+pN
q28V28His6XHpCFLmJ+TrOZitYNtInJNSeY2cCr3HTyKOQCtptzshtx24BjP
1LX3Xxu4yDq15DxJgSuSU0kstdZw+P2pJUEmBzhuUdP+0MMSWHUiWMa5D4Ja
Qg6fYKcZXFP++JNL0QmKGhX0uqKMgagWH9TNdIZtBvnhZaOGEEiO9rRdOwqd
OkNq//VrwnLXwNPPra5gFEk7s9wiDdEG9aF5zu5Q13BCuGRsji4BNO6OfioI
DVsVvhtSQStzNiT9bE+I6hpgTyIB+vHF68V5L29Q7w3dqvrHCv3+6ScT7nQc
Ss8m6nedskdNNQmjZ+p9ISUrzN/4izNS7z5oIaoRADJW2x5ZcBxBntdOXJk3
CgCfKL+gkypHkLHbQPlb6wDIpfY8Vg87gopkb50qowbAbk879VNbj6LGJm2h
wqsBEMMXZcJj4YIuxDk1ff4VADoZ+yny3a5Igde58svDQLiQ9vC46IorSnNk
5SvrDITRGAmNnXxuaD/TkHpgXyAMtz5cv2Dghoq0SjoFZwLhcmHsOumqG+Ju
ZZNplQwCebVPG6tkd7Reb92RnxwE/4YENDs7uaPPFusPR7KDYOXam9fdge6I
7WXJrExZECz+2kJjznNHG9F2Bm9QENRdDGuh/XZHYX3uFNpiEPywNR0u4aai
3wldeZuYguGnaIbTOxkqOtSyyeEFWzA8tZtbT3Ggol7RznzqtmAYTL/px+VP
RRfOe+tq7QwGuzSVttI4Kjqx2VRhi3QwtCRGKBzKpSKOl8tnJhSC4f++h8D/
AKBF/C4=
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  ImageSize->{155., Automatic},
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{
  3.498745062889016*^9, {3.498745286353226*^9, 3.498745305145572*^9}, 
   3.498745553198754*^9, {3.498745589654595*^9, 3.498745612191046*^9}, 
   3.498745645811014*^9, {3.498745721610733*^9, 3.498745766263813*^9}, 
   3.498745815329627*^9, 3.498746070006539*^9, 3.498746344876716*^9, 
   3.498746675823853*^9, 3.4989781280219355`*^9, 3.4989790023149357`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "wezmy", " ", "mala", " ", "dodatnia", " ", "wartosc", " ", "Etot"}], " ", 
   "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}, {3.498745236716957*^9, 3.498745272352102*^9}, {
  3.498745387753416*^9, 3.498745403942782*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", "1"}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, {3.49874529610408*^9, 
   3.498745299915419*^9}, {3.498745369904933*^9, 3.498745370358263*^9}, 
   3.498746301123456*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.498746344999764*^9, 3.498746675981414*^9, 
  3.498978128105935*^9, 3.4989790023319354`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746307936756*^9}],

Cell[BoxData["1.1408688900104609`"], "Output",
 CellChangeTimes->{3.498746345091053*^9, 3.49874667606546*^9, 
  3.498978128145935*^9, 3.4989790023449354`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", 
     RowBox[{
      RowBox[{"-", "Pi"}], "/", "1.5"}], ",", 
     RowBox[{"Pi", "/", "1.5"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}, {3.498745456274082*^9, 
  3.498745509952921*^9}, {3.498745843204477*^9, 3.498745850028618*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd2Gk4VWsbB3CZiSMyprDNY+ahDPsWZciwV5KQ2aYkm22I5KgIyZTKlHSQ
VFSI0gn7EVFpMBMRKZRMERJ5n/P1/bSu37WeD+tZ617/+16L5E07QGVmYmIS
2MTE9N/xURVpjK+Zir4cSBBPeU9HBzSmpveWUtHmSvWusbd01FhQMVqYTEXu
OYsf9zXT0VajU4GvCCrq+J1xZ+cDOordVpnaO+qLvjyJ+XDuPB01rG9hK+T0
RXtFve3VdOhIsZ/SZkvzRuKSHwP4ckOQaO8Qj4uyB/KmPVK+HROMtPmfC6lb
OKJqrb1eFfOBKHCAfZI9yQZt4tAvM/t4DHnZ9KaNV5miHNXdkssH/ZFeVYpE
cYkGmrAUo8at+CJpy2Rp9JEF/XwfWafG64NyHhb6tWvLwkpmM8d5PS/0/Ol5
B34rHciRyRa21fFARRnZVpZgBEkvXQyKTx1B32HF2Pu5KSxZ5hdwtzgjL4/V
oalSc2CZuXrOOMQJlex88usfQQvIzQ97q2LviLIG7zE3HLQCxfxwyUAZB1To
U5UfLbUfLjzvr93YRaCI7KpjAq9soFj537JT8Xbo45sB4RUVO1itYS9k6bFG
4y+4xz+Z2YPXoMOUW6IFEpFdajhKosAGt/Z5X35zZJIyzSwqQwGxnIJf91nN
kQo3O3eLLAWODdE/rS6boW8OoT4SihQYvz7Fkz1shra4ilDQTgps7Thwa6TM
DNGiZ8NfGlLAJ01h6fE+MxTQuZV78SAF1Ax8omkxe5D8RE9DYxwFjJjdhVtn
AE30DvTsOY/PZ+SHu40CKl/8nf0sgQLzviTDn12ApmyCJRsuUGBE9tuC4hNA
12KFLEvSKeDhx/3rnzhAS5e9/HSuUeB+NNLvFAVkotw3c7iSAvWDeiVeSSbo
olOzRUUVBa6end4+dsIEBf/sdmavpkDn8W5+fwcT5P2Suez+IwrUUhOORUiY
IM6pCZ65pxSIkl9Mqao2Rku+kop6LRQoFy1MLhg1Qsbz9bTO9xQw3+zcpGNs
iCgFb1TXBiggrz13JFrGELFNZf8r+4EC2+Kebm3mMkTdtGor+jAFhDJV2l17
dyOGuCd54xMFKKrmFjdpu5H1tsGOkSm8/zgji+yiXajGKP5V8xoF2joPRy1x
GaBARrDJy3UKPGwVSjo9p496Mkipr/9QAD0yorH26aMt7StP3zER8KPXpUHs
pj4azkrwf8FKQMLTVVkXE32UJGOkG8dDQFZRv4lGqB7yUMq8JCVOQFfK7u0y
wzro1tcHn+axxWrecmW26CCTewJGz7YT8OD23GHmBzqIRT+M01OCgOMcpkIT
sTooc1SbN4NEQI7pxz5E0kF2Pro/GhUIkH7wULCaqo1qTXOigrQJUOyJXE+d
0USbzYc0SToEfI/6OSvXp4lUnhxc6cT2qHnHx2Boola5pAptPQJO397ZsZih
ibr+edD51YAAS4EStSBtTRSXXsasZEKAcMKv8axIDfT7C0l50YKA+P5ntACO
nUjg0ubwq5YE9CnFkBe71dBiXNodXSsCRvaNuJ4pUkNzYsLtdGsCBv7cCLhh
rIZmObaKjtgQcEleKXglVBXl7Zp7kkAQ4KI1WSAwqoz6KwQfCh8g4KTsNRfG
fWVUORqfU4K9b8lCJui0MqpZ/aqIHAio4apr6xBRRrVmUr6TjgQc/hzIXWmn
hNquMPmxuRDAX8L1uq1OARXLK+nFY9e+WIzJTlZA56kbi5tcCRAa61v3PayA
nC7uIa9i14gPOrEtyqOWG/Ncw24E7N1izeGqIo+CAr2v0r0IcL5g+fVQnixi
Veb0bPcnwLAWnQ2PIqE/F1JlFI4SkOYWkNtsTUL1eoprp7EXLhX3CW4noQHD
+e9yxwjI6zJ4WtsghQSK+N4FBhCg0qPWyssqhRZWym+9DSRgx+PogO4dO5C8
xEiz8AkCDgnsPX0FbUcLHm6/3LC3rwudPOSzHeUlRpd8xX7L+PDXyG1xdHxJ
gmcpiIBSg1hRLp1tSDmKmzQWTICm3SvSPWsRFF6ZWLY9BNfjhZb+pGlhVE+2
tnXEzvc3FvLLEEaXXXp7m7CdeF6wyfUKocDlv1my6ATsty0WqfYSRJ5D+VYv
sFM/N9lcYxVEHZZHS39hlzc09cSVbkXf0/TvOofi57+Y9+rwtADKj1t981cY
AYxa+0z+KH609PZF5y7sDa4tT5i386OS3IJZH2yFParMPxu2oKawCyE12HYi
9Z8+sG5B3xxPZNuGEzBR5OtJJ/GiO0bJbiHYuaJxUVR/HsQUump4BZtaEh3m
cm8zio1fMu3DDmrPGrbaxY2uDhqFHojA17+7W8mY4EBrfLm0IGyXaxEbn/PY
0ZB+QuoF7D2NKXvTPrOhjUTVnfXYIRdyFj+fZEVSjfMVPdhHNIl9mY0sKJLP
03EaW7CYZQA2syCrDluWbScJOPb8Y05hwSb0d00fjzr2nO67codJJnQ8Q9LA
DPvNV7Nudi0mJNb+e+EodsubbPeR5HXGW7bc+FPYE8l22TYSawypbjb9i9in
RTTy/q1cZSgVum26i20fkVt4rX+ZYddhL1aLfelgzcm/ApcYK2uSDs+xG3hl
n57bWGQIGHbc78AeK6s4uJK5wCgN8lMdwm6/Kq9Gk//B6N33qW0C+34SoTbx
ZI7x0cQidR77H5nf0ShvhiFWkB28iv1wy0x6deMUgymyJ5w5koDfvEJH705O
MpTKma9zYbu2cfUW8o0zbMelxvmw7/v5leXpjTFMCzQJIWxFj+dPrriNMJhF
9YbFsG/xx5VIcg0y7vqop+7AHtyYShdx6mHEuEm4S2GnhqyUkazfMd7fYyGk
sd9f6dEVHm9i7H35kSqDTZOs4xBlrWScmKvK/8+Lf1bEp/myyRIaMT//W99X
XrLwV0od2crLJJiELbd+z0Dg7ksyu9ISjyS2zHC7DV9iJzl1pbhNHLu4V7JM
eFc/2ZjXolwEO4ed94rk1BA5K3j4vgD2498aWlffjpJ3rxzt4sG2dw2ru1b6
mTy394sYO/bRZd3s4jMT5PxW4twffL/OnBmcK3f+RrZ8XM6zhN33V9vyY61p
soj4Yu137FPZuic+sc2RiedyiZ+wSScvjR2jz5MPHySf7MPuWlCN+DH8gzyc
JVnf8N/6sws2HLU/yffiRgQqsdNbExmXZZfJw60xyUXYZ0k7y0iXVsjcaeSP
57BVHiZa7wn4Td5090hdyH/1SJSq9vSukZ0Ii6ce2HpDlZ+Pmf0hv94dJrEL
u3SqKtcmjQm8Tt47J4v9dxd5Ldd8E+zSKeLkw27K1YqcXN0E3yfuRI3iep92
qP+W5M8CsxM5s9HYSulFG/3ADixee5jcsXtVU47cYeeA5C6JzWRsZg+10OjX
HAA2dnx/8PtYJ7PIKXeYCyQPJbbTsONSzCLSaTzg7Paq0go7lJ4FNF1euOO0
mCKNvcAZaEj85oXc5x0aHTgPDFz5/4gn8oHM1lA1aWxb1RYVf3t+8N9V8WEB
503HqHixZxU/vDWION+M/ePDziBXIQHY5Xi0yRvbxuLADsdBAfCX8ZvNxHl1
S42vwtlfEDLPf1ptxPl43XxO9X6sCAzMzBvGYdfI7nZHn0SAOzE73Aw76XC6
S/deUSgqdO9ENAKeNZh2buIRgzieYfUHOG/X3mqmh+ZsA3pS/TEXnN/Z6Tdf
11buALo4obIZ+7vfsMaUoATkC5M//3sc1/tl1n7JSAm4YtZhKIJ9/accb4aJ
JIwZ28c24f7QRi2jt1ClYO52855vfjgPrlPpGz9IcO2M4HAG9u1C2qQLSRqS
zuVd0sN2s3if9theGnoq5yGWSkBEskhUWLk0eHG7NbD5EmAS4y+2TpUBB/VZ
yQlPAsINaw1t3svCkPCcdiw2k6yHUh2HHMRyV7gLY3dPN0Sq6cmBdzTHiqkH
AZfDeawEL8sBZwXop+F+mFhazr5gIw8L5mfaV50J4AlpvLeMFOCNvYB5PHZi
I00ldlYB5tqqXvNgH8q+bMQpoQhMTGXbxQ/jvJ17IicZrQhbx0yqVQ8RYHSc
97inrhKEOx8tlcH9vfN9daz2XWV4HWla74Hnh9C1kveaV9WgMN5ZioHni/Xl
BO6hVjXwPVMesAN7y3fJlORVNagJsaru3UfAhy8PL3332Al/eAyEjcxxfUgp
iL5RVof12LxvHXh+2Sd/quOspwaIauSHKOL5SODifEXlsiZEjLW7+moRIGn2
yFqDXwtKlt7J39DE+Rx20a5CWQu6rli682tgDzzNr3bTgsKhdzaTqjjPum+k
vm3SAs7J8EUPebzfkONZVhnawHqvpT5MlID+luvsx5V1YbFsJiVDBPcH16TP
RWa6MCTAIJcJE3Bx5aLnhyO6ICx6yH5IkIAL6xu7D6brwmFBXT4dfgK2OZVU
OS3qwg1pl74yLgKW0plPPmjQA/J1+5vPflEgIar0WcVBA7joETV6upcC/P0r
3bQEQ8hv53/s00MB3YYiSWqBIZQ96Y+36qbAK0q5xZFHhlBt38Mh0EmBZNqq
HmXcEIqeVhRnvaHA9FFqs90+Izgdr71wopkCfsV+KrfZjYFT6NPYATzP53L9
ftefZAJqsy8KPidSIPK5HPH6JcC3oG1px/D3Qqaa8FmBPoC2i/7q0/F4Ho+L
4XP5DMBDyCTMnaUAB3ehzfQfgP3vbE9ORFOgbGEiQFbHFK4oWDFl0fD3QN0Z
u94CUxgvSzrww5ECOwSDzhPheyD6puX5ixIU4GSbEZ2RN4eGQ9tYRO3toaox
5s42HwtYcP+34+tuO9isGAkn7lrD7pnrswHvbaDUeaOpz80OdnU21bke2g8U
bsdpCzzHnwBSxFiAFZz8Jl//hccB3p1p5bxAWID/USmYEHaEgju2ea7fzYFf
ZL++qb4TFH3lmD4yYgoNPcdNhCqd4cnUC8twXmPQU913xF7zCHS+NnQ2CdaB
lFWlAC8FDxC19OpMTpMFT77Z/nx1L5hic3gaMTnP+HLktPWHQW+os23++5yK
BqJXR8G5SV+ICB6Y6H4G6GaaSFb4Nn9o/X3mdOXB/YhpOqrXMvcYLMWXtjmv
HUQd5szwvSMQfrTE7qNOu6E8/UWXYZNgoMmntuwRd0c8gY/F1yyCwaWdvdLV
yh1ZyjevSRPB8FfuUqjpLXd0Z+4ftxSfYFj6xlk56e2BBgcuObxLCoZgs+JH
rGGeKFu+20G2Jxi4zh/6FSbgjdhju5bjw0PgkV/jTIqaN1IyD8qdiA2B01fa
qa8svdHD5g9ZjskhUFPBuf4o1hvFLSesONwIAbcu9SX3WW904kb9nQ8vQ+CW
hlarWq8P0v2ZuVmTRIfaHPGSxgUflBmWc91BlQ4tCxUzsfy+SFn0yvw5fTp0
TGhvj7H1Re+lVoHTjg7XY6Qj4lp9UZKNjV5FNB2km75B2LgvulollW6YSIcR
lYfGmWxUJJ+YHdqdSYc7Kk57nMypqF7ZyEr9Lh187mpzbPhQEav6puD5GjoM
PxON7ImjoqHRnNTGRjqMx/AuDRRTUa73ufobb+jwf/9D4H8VFvlD
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{
  3.498745062889016*^9, {3.498745286353226*^9, 3.498745305145572*^9}, 
   3.498745419166413*^9, 3.49874545835095*^9, {3.498745498001857*^9, 
   3.498745512072139*^9}, 3.498745851976731*^9, 3.498746070112791*^9, 
   3.498746345132573*^9, 3.498746676153096*^9, 3.498978128211935*^9, 
   3.4989790023779354`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "wezmy", " ", "nieco", " ", "wieksza", " ", "wartosc", " ", "Etot"}], " ", 
   "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}, {3.498745236716957*^9, 3.498745272352102*^9}, {
  3.498745951953567*^9, 3.498745967114758*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", "2"}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, {3.49874529610408*^9, 
   3.498745299915419*^9}, {3.498745974141631*^9, 3.498745974805883*^9}, 
   3.498746318782763*^9}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{3.498746345255925*^9, 3.498746676217607*^9, 
  3.4989781283109355`*^9, 3.4989790023969355`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746324175039*^9}],

Cell[BoxData["1.2661609883373448`"], "Output",
 CellChangeTimes->{3.498746345271189*^9, 3.498746676232041*^9, 
  3.4989781283499355`*^9, 3.4989790024099355`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", 
     RowBox[{
      RowBox[{"-", "Pi"}], "/", "1.5"}], ",", 
     RowBox[{"Pi", "/", "1.5"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}, {3.49874600926457*^9, 
  3.498746020517887*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd2Hc8V98fB3BbRnZKyohIRlbI+rzN7M8NGSEjO2XkUyIpUVIhRUrITCIr
X/GtPkd8kVWRVXaSTYis/E7//v66j+fj3vvHueec9+t9rqibv6UHHQ0NDQ8t
Dc3fa4JvdZ53mA9aSfDgDOA+jxQj+cU2rHxQT3ac0J8tCspsapXplPFB26Md
W57OUpDgDxqJ44PeKIQiPG/WQkGy6kmDF496I7+x/K2uGApi4QnnaRDzQq8v
s0IZEwVpm1ZyO23zQJsZW/EzvMEozhJxa3i6otgExoCPPkHIf/bJRaMJG9R5
71LW6YP+yIG1gSl92Rw9zf8yUix9GpUsbHhZa+mi+z6So+fYvFFfj8by6LA8
4n30sPm0vgc69OF3lcMkPVKzJ+bsdp5CjCJZ/me2iQMNPTX1DrcrMrX7YUws
KMFQu8O73hcn0Y0FZYu9DzWAxZp7XqjdAa29N2ce0tEB4V1xm2Rme+TIOaHH
t6IHrQ98lmLzbRBfufL5MjVDeHArRV0oyhotPEniaDIzAjupzLNbDpaoJmJp
KHfBGNauWGdZhhJopmukle6CKQwwcJZx7TVHv16IJrq9MoOGwkjf+FhjNKvt
GHmvwhxK1OPb0xwMUevs+znPVAtIciquoWXVQ59OMBz1T7eAUa9/13s2dRHb
4B0HSqYFtFLDm0t/6iJTOWkU+tQC2PRkJL17dVGfxvME33ILEOS4VvgjXxc1
Lc7yDTZZgEd79ld5I120zuJfKvrbArSVuy2+3dBBnzn5htXNyCDsz3HhOhug
HuYy71kLMvCH8SSHbpLQxGBSdMYxMoSRGQynvpFQ5Y+2yDUbMsT6J5V1lpCQ
2z9yOfdcyeD6pda/zYSEbOiajazOk0He/zCLeYQ20vr6XIAjnQyemh0Bpyc1
kdHttFXKEzLUtWdVKX/SRKs6bLFfs8jAVXm/jOaVJqJviKJmPCVDSKFSwZNo
TXTm3PdE9lIytC1MetCIaiJCWDYmpJYMrDX6VE07DTR/4exY+hgZtomplSU3
HkH0TApnX46T4cuf1TnxkiNowy+lvWGSDCaF/wVXPDiCBHRnlcZnyTDAP2Tc
73kEbfet9OVcIYP7g92DRxmPoOJkGcFL9AQo7rV9cl1XDT2W46xP5SGANtif
5tFbFbR3Qn42ipeAXiaOm355KuhEkxXHGT4Clqeit+vEqSCTrmpVdX4Cbkq1
08w7qaA5t2aVRgEC7lZWEx6bhxH94pt95SIEuHWrPHiueRjpph6+Xi9LQA1b
6HzyayUUWxRXe1WOgCwZhk5SnhJK3mW7T/MQAUUhM8qT8UpIjXVd/4U8AR1t
pc56p5QQ3/qoWIwSAUmRk3LMrEpoJUo4iOcIAeVVv8U7bBXRthWnxEg9AkLN
hjgCluSRAUeI3T59ArgPb5/TH5BHWwblBjXY16lCPIKN8sjEXCl5w4CA8ZFN
geZUeXR0dw+vjxEB5uESNOp68qh54luOgDkBgTYxhtnih/B6cpH+cpyAN3N2
Ua/cZNDcWQcTXxsCMmqXBpcVZVAtI13AKvbrS2GlqvQyKI9/6wOfHQF+vNy5
KFsa/Uo9/YN0goDPKww2c6MH0bG+rBSvkwRU0ZsIv/aSQonnBVXGsflHu7qE
1aQQ58DnPi9nAl6JdFKjt0mhuU85hu4uBGxl7N9h9+wAMkn8Fn/MjYAyM30b
/ilJxLlPUeCPBwElJ6Mv/zorgVz01B39PQkwEP1EhJIkkB+5uWAQm959kouW
SwKtveE8/caLgMndlcZ8pfvR/oj6sDM+eD6ttQ4cWxRHLAa71jL8CJCq85wT
CRFDTpetNMlBBPwpYby0HimCUk8/KH+MzTBtXtBjI4KmfRftJrD/PSfHW3lQ
BClT9eSunCNgUO0Bx1yOMGr08WXMCcbf+9Z+le36Qgj9p32+7jwBqzqzjBWR
gijfeamb4QIBwpT0u/cOCiJRpwlvfex/7QIkg9t3o8z6jl012N1abTYa+3aj
62Jny0pCCDgnN8kyXLMLdaqJffEOJSBE0nfuGi0/4pAXVX2M7VB58U/Esx2I
t9Kirg37Q6ZhUfixHSiG/oy/QhgBmV73PoRn8qGeu7tZZ7DL4wLPJOjwIr+r
axMClwg4xDyMUiZ40EbSmTUD7KJCroisuzyocla9KBW7lEFDtHqYG0XW5Dlo
hRMgSEnXW73KhQpfHa85hW2ZMqvEfpALofYZ25vYDZHaLyGHE4V70+xqx7Yz
CXZI0edAttwvR+0vEyCeLhHkFs2GpBPL2i9iX+jRIz3qZEXsz0Z+pmC3qSRx
fN7Piio+rNR+xu6mF4wwr9+GqkqbK3Uj8Hyy1anYMDOhqq/CZQ7YRT7027gc
GNEmy8T3c9jcTmfcml8woJnQhdVM7Jgb4lcNbegRq4H2l0psIjdHhqmADtnt
419qwfYRqjjUsEmLrtBFfF7Cbsn7yWqeS4OEtgRSma8QkDDz9L7Ywhb1Et3x
VAHs86qezsUP/1ATZ18baWDHj5ZQWsfXqdV6C79MsBWii+OdE9ao10N6+uyx
5TRuJy6qrlInVOttKdgGEs1zIjeWqd8zP09dxVZ1ytlVLfeLaleQWXUHuybN
uMG6a5E6rS9fk4JtnL8yMh++QH2jdI0uG1uyeMoybv9PqjlzWnghdk5a/Yb9
ozlqjNAVxQpsX/YjDx9OTFMb8+RE32D3Oj03/qo2Se1KLTCu+/v+1KM+oZgf
1MmGX8+bsLfnRIu4d49S82q5jT5ik7cmvhVIjFDFl7f2dmIXdt2Z/D41QJ2+
XyPbi00JVzUMie2lrmTbXujDTj728QK71Gdq5a93awPYaTElavrxrVRmaYZX
Q9j8lIzMwJZ3VCgXLB7G1n7V3wi8JVTJFpahv6bduS4R+T6J9Ptqm8lf/xuc
CHS/q0kexT5zg9iBhS2L7dqNpITNL6392KYML3U7rnwi8YyKjn3BfsV3x1x7
pIu0mK+l1o2tZx+4WaDXR3rpIlvXjh11r+TL7fYhkv7s3M1W7JJqOT+Da99I
jPRRNxuwbf+Uu28pjZGahMbfIeyKgppf1aPjpO743SpV2Gdv9RIhSVOkoSCB
byXYXOvhT1UNZ0n05SPvn2Krr6YkXQqbJ3WwBE+mYYvffn19oOMnqWZitOcG
9pvcCfu8qCUS3+zbvDBs3uaQM+z9v0iBW8dfnMV2/uzIG3x4hSQq9MLXEnvq
j9JFo7FVUruAt7AeNufzEt8K7XWS8dbH7UrY/i5Jx8QfbJCcdTMecGFnR5S4
bTPaIpnRbij+weuXvVB9UJuggSS/ZeYpbJ17zlWhGzRwfbzk5DvsdRl68WVr
OpBYffbWC/u7aqPsJiMjNDNtMZljS6g8rh1CjMB8NNdDAdtq33v5+jAmiLhT
57GC9+PHZ7Efk38ywzln+4AQbE0Z52GinxWkuQ18bbA13uTP6KewwTfykVAl
bJrwXYSGFTsY0/L0TeJ6cPP2kM2h99uhftoyzwL7Tw+JL/UEFzBnW5H3/73f
Ln9M+BUXtP6m41jH9SeZrvJ3zg5uqPXSaMnGTooyf/ryIzeMeOSaTOJ69mHT
xHXMkBccjrPqmeL6t6peTHqkxA8FP6OB9299LOofUrvLD181Gy2/XMT77zan
X+8sP5w44VThgX2xz99FpGAnbK8aDgjE9bawbOpJvYgAdCzLlDri+p2nkXo1
jGMPTDxJWdiBXdSzHKnptwe+8aQdbaPg/P4UMrv1fg+wKwrKaGJLd+8Qj4ve
C5o7b/ax4zwYDPo52bghBJ76BxquBuI8VNwbOvhLBEYMNJcVsI/s/c5aJioK
nQ2yNsMB+HnXxqFoc1FQqjuer4Xt2s3qo5ArClXL7vemz+L+ZOnfoynW+4Dc
IhkohvOrmktINL9cDFqaerLrTxPguOAZnzAoBo1uNEw+2Ok7cy+EsomDopil
3XNfPN7Co7WWp8ThqpSpogTOQ+o7uZwdvPthzFM1dwnnq3ak//kfQRKgrBxa
G4v9ytSWdyJdAtavJXGKYLPUxThPN0mA1bUxeSN3nC/m8ry/90mCa/rX+3E4
r7e3RBuId0gCdSG3cRHnfbZ3U0O7khScefJim7AtrgdcUrMZy9Lw8/lYdDzu
JxjbEpzvSciARNZz5i3cf0x3HDG9aSMDndIHd321JsDwcn3WpX9kYJmbxfuG
JQHFLNB69bwsDMW9fVSC+5c9ctljLCtyUKFYtayH+6OM4KvsuzblYdyM9Oy+
LgFB69lFWjsU4O6Y+IlRHZzPj9+2u8sqQByv0j9XAI/Xsq+00kkBuLmMjhVp
EXCctY8u4K0C/DDIHO5UJQA5003JXFGE5IziZjdpAoamJ1t/MChDgn7B0+SD
BGxqvr/ycY8yjL0XjHovRcA8sWe+WlkZSPa+WjIHCAg3WWJPclcGFceOgu/i
BNgg+/CTdcrQmVukdkCIgHvLzw4aRR2Gxh+l5pOcBBw2H/W5zqgK0bENJrlD
ZFhyhBtyHOrw7oC9V9EAGdYjbmUtiqsD38Xkp+V9ZFD7boyqNdRh+eBsQnUP
7vcvzbCSfdQhNP7HSulHMjhLJfIm1KnDg8KmVENEhuculQWUMA1YM2rSpcf9
vBDf196+SU3w+9Bpnm9FhuVGGvY7H0nQ70tkbxFkWKP1kv86QYLU9qRH1vg8
8f3V+P5kWoBTO1d9V4zweaPkeNXUbgDu9sOhItpkyEp7N5ZvDtBtxR+zJYnP
A+4miX7lALY1FV8nVy3AyO81rdgVHTjQILqonWwB/aJrqmlCenAtvXyjlGoO
D6zbOGvxuayJPe719H9mIC9jWOHragyXez91rkaZgljt2urlFTPgt4eRaTET
GGQO8e/Wxut8gSnnoKsRlHARZm8MLEF9p163aKAh/HkOFxdNrSGJx6N8Q0Uf
bsxED9KdtQEDr1ouwkEHLPdoXej/aQfFs4Mu7x01IS3ORXEg2gFC0kpbD+1S
Bs79lU4HMk+CDX2a/xNXcdhffHCKk9UV2h21pmunf1K99RQXTg25QVaZx7qN
uzySSvMY2lTwgKfN+rEN2jqoMmvftMV/XnCerTzNOccMLTYYrpZ0+IIjk9Lw
40QbNOu5Idi7dRa0u0ol4zNdUICkapAZRxC8cNlkGH/jgjgzikeiBYOAIXbM
2f+LC7ItnWJqPxAErt/UDdn4XJFVruTObL0gcL6VIBF10xUFciZLy4YGgZfB
UsBwuBvSOvA5jWUqCIqu//OwK94dyfIdSznafw52cMUkzpe5I28rDSGVmXMg
oMLbYNzljvYGTo9pbZ6D8bkjIa/3eiDLACW9R3uDgd4g+1ZFsQciam/43nEO
hoBLnzoShzwR27Dyo0fjwdAzI9i/wOyFBh3eRmitBQP7+s64tENeqO20q/Af
Ngp0Zk8IdUV4obB9A229hygQ464iNSXmjZDkxgb9RQrkzCffrDT3RjxP5qvS
blGgS5QgOi94owiF5A7bdAoY0rYVbmv1RnzPeqLl6igwkL9DQOS3N3KIUOky
7abANGdfym0xH+Tb0rEQO0mB/gFZETLZBwmabahPb1Lg//6HwP8ApxIU4g==

      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{{3.498746000967126*^9, 3.498746023966454*^9}, 
   3.498746070223589*^9, 3.498746345314657*^9, 3.498746676336131*^9, 
   3.4989781284109354`*^9, 3.4989790024439354`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.498745089125189*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Etot", "=."}], ";", 
  RowBox[{"e", "=."}], ";"}]], "Input",
 CellChangeTimes->{{3.498744548017966*^9, 3.498744568889582*^9}, {
  3.498744895344405*^9, 3.498744896370717*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "wezmy", " ", "nieco", " ", "wieksza", " ", "wartosc", " ", "Etot"}], " ", 
   "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.498744673639998*^9, 3.498744688015705*^9}, {
  3.498744734435433*^9, 3.498744736178534*^9}, {3.498744769126688*^9, 
  3.498744804546868*^9}, {3.498745236716957*^9, 3.498745272352102*^9}, {
  3.498745951953567*^9, 3.498745967114758*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Etot", "=", "20"}]], "Input",
 CellChangeTimes->{{3.498652176033926*^9, 3.498652231536919*^9}, 
   3.498744079414316*^9, {3.49874463400905*^9, 3.498744656293548*^9}, {
   3.498744721325155*^9, 3.498744721825232*^9}, {3.49874529610408*^9, 
   3.498745299915419*^9}, {3.498745974141631*^9, 3.498745974805883*^9}, 
   3.498746318782763*^9, 3.498746508874082*^9}],

Cell[BoxData["20"], "Output",
 CellChangeTimes->{3.498746345255925*^9, 3.498746511015741*^9, 
  3.498746676384271*^9, 3.4989781284939356`*^9, 3.4989790024589357`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"e", "=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"1", "+", 
    RowBox[{"2", "*", "Etot", "*", 
     RowBox[{
      RowBox[{"Ltot", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"m", "*", 
        RowBox[{"kappa", "^", "2"}]}], ")"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.49864925995262*^9, 3.498649289813854*^9}, {
   3.498649333396809*^9, 3.49864934323158*^9}, 3.498746324175039*^9}],

Cell[BoxData["2.651723304546313`"], "Output",
 CellChangeTimes->{3.498746345271189*^9, 3.498746528202713*^9, 
  3.498746676398792*^9, 3.4989781285339355`*^9, 3.4989790024719353`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PolarPlot", "[", 
  RowBox[{
   RowBox[{"rho", "[", "fi", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"fi", ",", 
     RowBox[{
      RowBox[{"-", "Pi"}], "/", "1.7"}], ",", 
     RowBox[{"Pi", "/", "1.7"}]}], "}"}], ",", " ", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.498652329713052*^9, 3.498652380576579*^9}, {
  3.498745012104676*^9, 3.498745038553551*^9}, {3.49874600926457*^9, 
  3.498746020517887*^9}, {3.498746547325026*^9, 3.498746573893221*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}}, {}, {{}, {}, 
    {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJxd2Hc4Vu8fB3CEjOyMyH6MrMde8X0+PCehSDZfslVmZO9UViGSksqIJGRU
ZPXckdAQKfI1UlYqo4Uo/e7+/f11rtd1znXd1znXOe/35z6SnsHWPgx0dHS8
9HR0f4+idiSOsgpP5EHSbZP3TkI6zC7zm56e6Ow8vGnWS0ID8l3iNHFPtFSw
x1uWOwm1F35cJa56oOgzTwyWEhKR7+LVO9YV7gh56baXWiag4O1k5c23rugT
g5h4gUoc6qmYninZ54j0shdU0+5HIQdHh+onsVZI+HIdpeq/48j1daAOSd0A
sfC0jT5+6Y9kHsyX0KhSiFvvQrNbgw/ajDKciRRlho2sY5J3uj3QUPY+r9dy
chDclJuvf8EVrR1acFNKUoOVB46JZwOc0M3a5mqng9qw5eeOWzfn7JB13Nvu
7K968I7f9V23hA3iquiM47xuAHsqzv0MWrVCs6ECOq6b/0DdufinFxYtEPnx
o+UoHwDm5Mn0pFpzBG0+XSWSRrDu4LWZNboXSVHl9EWfGUGQ8T2KjvAe9H4y
83TgQWNI94vjlH1kjGwk+Ve7nhjDydanxg7HACXwZeXpqVBhJbD1/YUyA2RP
lyrwK5YKgzMsPBlLuuipjeYFyQ4qsDGWTv+S0EIPngsPcXRSQcKpa0CaTwsd
dbfy/Ik9v+yvb8GkhSIkF1z7uqiQfOlleMW8JjIYaP4Y3EuFYZl/nKMbNJFc
IkNxTj8V+JaKtuYQmgi8U/IOTlBhLIh884SfBpqcsDozuEoFX4Kc/7NNDdmw
d3GXr1HhjPU7u/O3sce1hsN/UsHNzCiAXKyGuL5YyPNvUKGO+en5oJNqKLJ3
67T5JhXmys5nspiroZZvf5jOMxLApTehVDmsig6nkpkucROgyRb+XHiejF5n
M1CseQiItV7VHH1GRoosap3svAQod79tKaojI5JeY1UcHwEvg6u5taLIaOz8
PKO9AAH5HnoB57aS0YEDKnxzIgRoGEwrMMupoNALP46ckiXgkW9NTP8RJaSi
yRYmI0dAvSzPbKaFEmKIHRruwp5X2ttrqa6E+oqibjHuIuBUz4e5/zYUUW07
c3KMIgGkihKRndmKyHUyTWGvKgHxm3NclEYFlHAi+bWPHgE6ueMSC8y7UMAH
T/ev2PcqVwnez/Jof2KOeYI+AdtFDP30B+RRL7nQMG83AQY9M5XnC+WR/aXi
6/cNCdAuNC04qiqP4vfcvDRlRMCb7frCxc5yyIuv36vTjABP5UQG8XoZJOty
V0TTnIAKXqZy70syyAFd2V+GrbF27GR1ogwKW9vccXIfAR0f3n0wtpRBwR7X
ffUtCMjTa2DJ+ERCR8i1U9lWBNhfydVIlCWhAVZRk2F7AlKD9Uhfr0qhPtT3
r44DAa+Cyn9En5ZCKnn/LuZjz/j9qdoSKIU+BUZk2DoSAHL/te00kEIF9E0l
vU4EDAy+J/n8J4myTDm9C10IMDG5/WC/oCSyeyO5MO5BwIdDJ8/0ksTRF/fn
pzU9Ccj9KtOmviyGDGP50jKwB8eeFRa1iqEKx+leTS8CAp8yWJ6yFkMSpbQ9
yd4EhMuFd0cliSLXs9Pa9IcJkEi85sYyLoK+sSZGWGJ3c20er7wpgoxagyUL
sUsLil9YhImgNN/WZvUjBGwhbRgXs4ugeqXHK85HCQgQjvH11xdGJVwH72f7
4+elcZL1wiUhVKzN2jyIbUfxzD3sLYRIjddkBQMIKM4M/mWgKoTGottMr2Cf
IHU5L3cLojpGjaqiQAIer8dnR60KII6W6Y2MYHyerdue1YEfHfTtcvweSoB+
eksWnRAv2iZd90rqOAE0+vXe83M8aO1G7Esr7DHftyPyTTxIJeWjTRX27KmT
Xo72PMiLPj/YOYwA+rf1Oc/yuNG74p38JeEEPKUdXGYlcaKuAfZ73djqwb9M
fTs5UM1y8fACtmbx+K0uTw40pTpfphtBAOOJYM/Mkm3IT0vcowfbwcUigyLO
jppPbD/aH0lASNLkNlkRFnSRcY/uEnbGmTqr2y1b0Y661gyOKALkRH9s6jpv
RfFZ15pNsc9ssTluV8CM8vrkoBn7rsz4jUpBJjS9EanwCltVo2sX0cSIhFus
Li1ii+1Q6Jm0Z0T1fi2bktEEJNfRqUrkb0GbOgarydiuJiqXGrczIP2Uos0C
7J4e6pT7XXrE2PIovQ77u+COjW229OjLsyDFMex4ibtH/c7ToaK8pSvKMThP
7jq5hZI2acN0PJ2A7eYQVfjO6zct4Sop0AZbIJ3ulc31X7TYSR2LSOz9WWOf
QGqDFrp2xbAZW/n2cs5uiTUaddNEtwf7w/LyrjSvVZrn59HWIWyFm7ZswzdW
aJ9PDbl+xe44negcq/yDJjVBjfuDXR3jU95/7DuNiK+S3hZLwHjI50W5u99o
JgmRHCTsRVeS5pj+V9rI6fcHyNji3etseglfaA/E7ET0sV+yt+dffLhMq55y
sbDETiTXMd85t0jLq16944C9lfl+QrTgAu0slF93x/7JSDw3vvaJln7RWCwE
2zSy1Hm86gPt1aJ2ZRR2741XB26rz9Ga9fa2J2KbOZqnnWieobUEhzunYAdM
ab21h2ma6eHOpLPYBSeekFW639Oev1ZXzcXu+MfAbqvlO1rFga4jF/+uN5ho
96jyLU3ZNE7hCnb4NgHB+KYxWjWPU3gx9lx4Rblu1wgtZNzFrAy7Kb2kaeXl
EK21Ne1GBTb789AtTZODtN7C8bO3sA0bMijRi/201/tcflVje31vIGLWntL4
opl+3MaOUyx+xLj9Ma2pdiSiDpuDE/nnkBFtwXnkRD12KvuCBSX/Lm3qEaNA
A/ZKTvHZ1sslNCTgovHXbtkeSoJjqZQax7cjf69Xl2hr6k6roTwPyWL96wiR
UWOTzBaKsfzRjlrs2lGXYn7FToq8ZDDj3/WTat78TufqpVybLHpRhb1x+vcz
uu99lNbedelKbOmIMCvDgJeUC2KJ9Dew81+5MMcdek3Rp1M6UoqtyNAz0Gb1
hnI/iM7xGvZDur2//hiPUvjNNvoLsKf7/F/s0ZqgWFmIDOZh63Z6m3efn6R4
CLkdysbOM1w3X9j2nmIo2xWajn1XyYiTP2WKMmf7RTsO2yKW95hf1CzFq+fO
bBj20SBtx0tf5igcpBzxQGyTnZT8Hr95St3z80ou2HoCBhWqrp8pCz536Wyw
zdsHOo8MLVCGYmZdzbDvfcvqvX5gieJDDbujiS0wPJfDI/mFUqP7pGoXdoP4
n8k7418o08qyymLYafvtFh0uf6UQ5KknTNi3ZEMyy/m+U0aWuUJf4O+BtTri
nSPzKmUpzD6Ahp3pnEriHFylXOAO++82dtoHNdLjojVKRK0Y71lsx4Ccpn/0
1ilbPbMWKdiWA8c7bAN+U+pbPqUoYmfIEAlSepsUAi3lCmBHplne+cb0h5Kl
1s/0CeeBAAc63n+ADjbIcvaZ2Ia3J1Py5ulh+N5SRhi2Ypt7HpsVA/StSGn8
i10Uxdie3MgAOcY3x2WxZ1k5sqNPbgEJNv3LTTivIo80n0gWZQIpE+PKBzgP
Dbb9od9mzQK7zAK4r2IH39llVXmfBTxdfT/GYDfQHCLNxFmhsvbqiiZ2r4xX
Uu5nVmiPGgwqxXn7QmLbU9NUdiiMuzDoifO5Vcey0aqdEzT2CXnqYjvflVWW
JHHBwdYSTw7s8hgzgx8ZXPC9Lrv5Hs7/jQ9tu2yVuGF2htL+529flImPiq5z
w+BkHmc07pfUqYbP3Pm80Or9R94Me/0h+YjcY14IvLX4WAhbQFq4h7LCC6dq
pDMbQ/D9Zi6bRtjzQc1XR6PPxwgYuvH+zbrAdmDg2GIAuL8+M5an+F7ih4mK
qVQW7Ik9IZev9vBD+i4n+/4g3CdK398NrfHDyfGHMYewzRW4JaycBCClp+Vt
GO7D2t2R067CgvCO1Z49A/cpP/3vpY9XhMDDqiDKHPuB242Le58Lwf5Yl8Ns
2H/Uic4bv4WgR1CoN82PgK6snJaAQzvA0sKmMBH3NVoSbxKUEIbbKYNH7XC/
e/epOy1eF4GHMCNWgueHqxDK2lEjBs92D0ybYx86ZxGW9V4MgsN5+L7j+cPO
sTzvkKA4tAsL3iSwR9ituNlOiMNrj9m6CTfcZ3cG7fcqSIA6u5PUCp5fHuU+
VVTYlIBBHp2i+r/zkfG3vUkVUhDCUKlmis21OmYr3yEFm4Sp/IQdAf5Vl26+
HJOCY19ijFiwB1ps1FV4pUEjskDczgb317Y2Mn28NPjE/HN04MDf/Ds5vsOG
BHLGTMI2pvj9MOEsVaGXBUGGNtU3ewnQekIxvikiC8oPPKtdsKWp80kkbVkQ
fjg84GVCgPWm7RzJXxYmze73+xG4nxsb7u95LQtLanZR/wI+76V5nO2WHHgP
9qk/1iHANroA3bfbBR+YNSobSQTsWNMvDGlUgpRu3nQxbI+Ny3EjL5QgNMBk
+LQ0zverrEHEvBLEFCy0W0sR0CmalCm9Uxme5HDUzojjvjSUFNiarAylQ5Er
S8K4H9xTOXwtVeC+sfCDPDzP2874vi2YJ8PU4mgh3zoV3pr+sTe9oAY/KRqH
9+D9wumOaf7U62rgy+MfGYH3E/cUTMR66tXAcpLqM7RChYnfcdsP9qlB59iv
wMxvVHjatnUoYas6HO09sTD8mQo7+6Kq98Wog+0hSc8xvF85ViD6ocVDA3Yf
/vgmG1FBOmQkO0JXC2ziZ35sT6TCwU6eeJNuXZBtkwjxUadCjbKQsWWYATQc
VRzd+cIYtvDK2Y86AaQmRRTNuBsDs9aqpE2sMbiVuWhcfGkEUym12w038Rxm
6Ks1rGEEWdqJcl25e0FkdlIxLAZARdifgdPDHDo8sh6EcFEg/VzlWkm/BfBm
/7h1cdYAVrvQjUyaFRzg1ZkXpNMHKycdIXkBG+hoYPjB6KkNam8M0lVH7WBg
r7zPSzZ1KLZwLRBldYIXxwTtx3LlICmsLMLomCscTF6l67/KAfeiLAsDIj1g
e2hNgvtrMSQfqCJSutsHxte8D6bt342Wjt8trUr3h47pfomLEpaox67s1u7V
UNDyrA5WSrVExTKsEU6CxyH4SMpu2SVLpMca/nBa5zjcqvJO9+48gCYXOQf3
xxwHwfo9ol8fWCG9/blEOWMYhD08m6pXYo2EZtuedqmFQ/m1a8G5Cg6ogvdK
7CxjFHiJDtQ0xzugwMLrtvm6UaCFSpKtXjqg35Zay3KBUTDaX5b2LdERDZ4p
lWQZiYIt/KpPVBedEMN1I8egtmhwKrsm91HCFZW59uxQvx8LsoM3J9zPuaKK
+KVx4ncsCMUuOvAyHEICndXSjdQ4GNBlovh/PoTW9xv/mngVB5dHIulOBrsj
Q/KvzizWBBDRFnxCrndH/Vcbng07J0DGYs2I2Xd3pOL/OvRdDXZdvd1Kggca
9A7XCHVKhBde0V6ljz1Qup9Kt1BDIpz+Xsu6xOWJntXwJroxJ4F1x/6Cj86e
6J/GIO1WzST4v/8h8D88mtO+
      "]]}}, {}, {{{}, {}, {}, {}}, {}}},
  Axes->True,
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", TraditionalForm]},
  AxesOrigin->{0, 0},
  CoordinatesToolOptions:>{"DisplayFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& ), "CopiedValueFunction" -> ({
      Sqrt[Part[#, 1]^2 + Part[#, 2]^2], 
      Mod[
       ArcTan[
        Part[#, 1], 
        Part[#, 2]], 2 Pi]}& )},
  DisplayFunction:>Identity,
  PlotRange->{{Automatic, Automatic}, {Automatic, Automatic}},
  PlotRangeClipping->True,
  PlotRangePadding->Scaled[0.02]]], "Output",
 CellChangeTimes->{{3.498746000967126*^9, 3.498746023966454*^9}, 
   3.498746070223589*^9, 3.498746345314657*^9, {3.49874653208731*^9, 
   3.498746575498054*^9}, 3.498746676503198*^9, 3.498978128584935*^9, 
   3.4989790025049353`*^9}]
}, Open  ]]
},
WindowSize->{640, 750},
WindowMargins->{{342, Automatic}, {Automatic, 0}},
ShowSelection->True,
FrontEndVersion->"7.0 for Microsoft Windows (64-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 355, 6, 72, "Input"],
Cell[903, 28, 111, 1, 52, "Input"],
Cell[1017, 31, 2571, 60, 352, "Input"],
Cell[3591, 93, 92, 1, 31, "Input"],
Cell[3686, 96, 947, 26, 132, "Input"],
Cell[4636, 124, 141, 2, 31, "Input"],
Cell[4780, 128, 1712, 41, 192, "Input"],
Cell[6495, 171, 213, 5, 31, "Input"],
Cell[6711, 178, 389, 13, 31, "Input"],
Cell[7103, 193, 293, 8, 31, "Input"],
Cell[7399, 203, 425, 12, 31, "Input"],
Cell[7827, 217, 280, 8, 31, "Input"],
Cell[CellGroupData[{
Cell[8132, 229, 96, 1, 31, "Input"],
Cell[8231, 232, 229, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8497, 240, 96, 1, 31, "Input"],
Cell[8596, 243, 243, 4, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8876, 252, 93, 1, 31, "Input"],
Cell[8972, 255, 230, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9239, 263, 69, 1, 31, "Input"],
Cell[9311, 266, 256, 3, 30, "Output"]
}, Open  ]],
Cell[9582, 272, 1752, 41, 152, "Input"],
Cell[CellGroupData[{
Cell[11359, 317, 439, 11, 31, "Input"],
Cell[11801, 330, 19109, 324, 375, "Output"]
}, Open  ]],
Cell[30925, 657, 87, 1, 52, "Input"],
Cell[31015, 660, 210, 5, 31, "Input"],
Cell[31228, 667, 427, 9, 52, "Input"],
Cell[CellGroupData[{
Cell[31680, 680, 545, 15, 31, "Input"],
Cell[32228, 697, 177, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32442, 705, 422, 11, 31, "Input"],
Cell[32867, 718, 159, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33063, 725, 439, 11, 31, "Input"],
Cell[33505, 738, 16577, 283, 279, "Output"]
}, Open  ]],
Cell[50097, 1024, 87, 1, 52, "Input"],
Cell[50187, 1027, 210, 5, 31, "Input"],
Cell[50400, 1034, 491, 10, 72, "Input"],
Cell[CellGroupData[{
Cell[50916, 1048, 594, 16, 31, "Input"],
Cell[51513, 1066, 176, 3, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[51726, 1074, 422, 11, 31, "Input"],
Cell[52151, 1087, 160, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52348, 1094, 439, 11, 31, "Input"],
Cell[52790, 1107, 15562, 266, 234, "Output"]
}, Open  ]],
Cell[68367, 1376, 87, 1, 52, "Input"],
Cell[68457, 1379, 210, 5, 31, "Input"],
Cell[68670, 1386, 396, 8, 52, "Input"],
Cell[CellGroupData[{
Cell[69091, 1398, 310, 5, 31, "Input"],
Cell[69404, 1405, 142, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[69583, 1412, 422, 11, 31, "Input"],
Cell[70008, 1425, 142, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70187, 1432, 634, 15, 31, "Input"],
Cell[70824, 1449, 6796, 123, 381, "Output"]
}, Open  ]],
Cell[77635, 1575, 87, 1, 52, "Input"],
Cell[77725, 1578, 210, 5, 31, "Input"],
Cell[77938, 1585, 446, 9, 52, "Input"],
Cell[CellGroupData[{
Cell[78409, 1598, 360, 6, 31, "Input"],
Cell[78772, 1606, 141, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[78950, 1613, 422, 11, 31, "Input"],
Cell[79375, 1626, 158, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[79570, 1633, 587, 14, 31, "Input"],
Cell[80160, 1649, 6698, 121, 447, "Output"]
}, Open  ]],
Cell[86873, 1773, 87, 1, 52, "Input"],
Cell[86963, 1776, 210, 5, 31, "Input"],
Cell[87176, 1783, 446, 9, 52, "Input"],
Cell[CellGroupData[{
Cell[87647, 1796, 360, 6, 31, "Input"],
Cell[88010, 1804, 143, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[88190, 1811, 422, 11, 31, "Input"],
Cell[88615, 1824, 161, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[88813, 1831, 540, 14, 31, "Input"],
Cell[89356, 1847, 6565, 119, 447, "Output"]
}, Open  ]],
Cell[95936, 1969, 87, 1, 52, "Input"],
Cell[96026, 1972, 210, 5, 31, "Input"],
Cell[96239, 1979, 446, 9, 52, "Input"],
Cell[CellGroupData[{
Cell[96710, 1992, 383, 6, 31, "Input"],
Cell[97096, 2000, 166, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[97299, 2007, 422, 11, 31, "Input"],
Cell[97724, 2020, 182, 2, 30, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[97943, 2027, 586, 14, 31, "Input"],
Cell[98532, 2043, 6575, 119, 447, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)