(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     34806,        784]
NotebookOptionsPosition[     33248,        723]
NotebookOutlinePosition[     33775,        746]
CellTagsIndexPosition[     33690,        741]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Solve Wave Speed ", "Section",
 CellChangeTimes->{{3.405144865341703*^9, 3.405144874435622*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"TravelWaveEquation", "=", " ", 
  RowBox[{"ForAll", "[", 
   RowBox[{"U", ",", 
    RowBox[{
     RowBox[{"a", "*", "s"}], "\[Equal]", 
     RowBox[{
      RowBox[{
       RowBox[{"a", "^", "2"}], "*", "D", "*", 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         RowBox[{"2", " ", "U"}]}], ")"}]}], "+", 
      RowBox[{"A", " ", 
       RowBox[{"(", 
        RowBox[{"U", "-", "\[Alpha]"}], ")"}]}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.405140355729174*^9, 3.405140410841956*^9}, {
  3.405140471803191*^9, 3.405140477867083*^9}, {3.405140560710653*^9, 
  3.40514061869506*^9}, {3.405140855599902*^9, 3.405140894359185*^9}}],

Cell[BoxData[
 RowBox[{
  SubscriptBox["\[ForAll]", "U"], 
  RowBox[{
   RowBox[{"a", " ", "s"}], "\[Equal]", 
   RowBox[{
    RowBox[{
     SuperscriptBox["a", "2"], " ", "D", " ", 
     RowBox[{"(", 
      RowBox[{"1", "-", 
       RowBox[{"2", " ", "U"}]}], ")"}]}], "+", 
    RowBox[{"A", " ", 
     RowBox[{"(", 
      RowBox[{"U", "-", "\[Alpha]"}], ")"}]}]}]}]}]], "Output",
 CellChangeTimes->{{3.405140866594586*^9, 3.40514089501035*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"?", "Resolve"}]], "Input",
 CellChangeTimes->{{3.40514087342153*^9, 3.405140874911013*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox["\<\"\\!\\(\\*RowBox[{\\\"Resolve\\\", \\\"[\\\", \
StyleBox[\\\"expr\\\", \\\"TI\\\"], \\\"]\\\"}]\\) attempts to resolve \\!\\(\
\\*StyleBox[\\\"expr\\\", \\\"TI\\\"]\\) into a form that eliminates ForAll \
and Exists quantifiers. \\n\\!\\(\\*RowBox[{\\\"Resolve\\\", \\\"[\\\", \
RowBox[{StyleBox[\\\"expr\\\", \\\"TI\\\"], \\\",\\\", StyleBox[\\\"dom\\\", \
\\\"TI\\\"]}], \\\"]\\\"}]\\) works over the domain \\!\\(\\*StyleBox[\\\"dom\
\\\", \\\"TI\\\"]\\). Common choices of \\!\\(\\*StyleBox[\\\"dom\\\", \\\"TI\
\\\"]\\) are Complexes, Reals and Booleans. \"\>", "MSG"], " ", 
  ButtonBox[
   StyleBox["\[RightSkeleton]", "SR"],
   Active->True,
   BaseStyle->"Link",
   ButtonData->"paclet:ref/Resolve"]}]], "Print", "PrintUsage",
 CellChangeTimes->{3.405140875761231*^9},
 CellTags->"Info3405144475-5436130"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ResolvedEqn", "=", 
  RowBox[{"Resolve", "[", " ", 
   RowBox[{"ForAll", "[", 
    RowBox[{"U", ",", "TravelWaveEquation"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.405140660008911*^9, 3.405140667329401*^9}, {
  3.405140704589896*^9, 3.405140724034268*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"A", "-", 
    RowBox[{"2", " ", 
     SuperscriptBox["a", "2"], " ", "D"}]}], "\[Equal]", "0"}], "&&", 
  RowBox[{
   RowBox[{
    RowBox[{
     SuperscriptBox["a", "2"], " ", "D"}], "-", 
    RowBox[{"a", " ", "s"}], "-", 
    RowBox[{"A", " ", "\[Alpha]"}]}], "\[Equal]", "0"}]}]], "Output",
 CellChangeTimes->{
  3.405140668333696*^9, {3.405140709896691*^9, 3.405140725029342*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Soln\[UnderBracket]as", "=", 
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{"ResolvedEqn", ",", " ", 
     RowBox[{"{", 
      RowBox[{"a", ",", "s"}], "}"}]}], "]"}], "//", "Simplify"}]}]], "Input",\

 CellChangeTimes->{{3.405140730090506*^9, 3.40514074613774*^9}, {
  3.405140787691478*^9, 3.405140826577784*^9}, {3.405140971616575*^9, 
  3.405140983875491*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"s", "\[Rule]", 
      FractionBox[
       RowBox[{
        SqrtBox["A"], " ", 
        SqrtBox["D"], " ", 
        RowBox[{"(", 
         RowBox[{"1", "-", 
          RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], 
       SqrtBox["2"]]}], ",", 
     RowBox[{"a", "\[Rule]", 
      FractionBox[
       SqrtBox["A"], 
       RowBox[{
        SqrtBox["2"], " ", 
        SqrtBox["D"]}]]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"s", "\[Rule]", 
      FractionBox[
       RowBox[{
        SqrtBox["A"], " ", 
        SqrtBox["D"], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          RowBox[{"2", " ", "\[Alpha]"}]}], ")"}]}], 
       SqrtBox["2"]]}], ",", 
     RowBox[{"a", "\[Rule]", 
      RowBox[{"-", 
       FractionBox[
        SqrtBox["A"], 
        RowBox[{
         SqrtBox["2"], " ", 
         SqrtBox["D"]}]]}]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.405140735591216*^9, 3.405140747050545*^9}, {
   3.405140790362928*^9, 3.405140827041141*^9}, 3.405140984939232*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Quadrature for U(\[Xi])", "Section",
 CellChangeTimes->{{3.405144828833828*^9, 3.405144839673553*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"UIntegral", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"Integrate", "[", 
     RowBox[{
      RowBox[{"1", "/", 
       RowBox[{"(", 
        RowBox[{"a", "*", "U", "*", 
         RowBox[{"(", 
          RowBox[{"1", "-", "U"}], ")"}]}], ")"}]}], ",", "U"}], "]"}], "/.", 
    
    RowBox[{"First", "[", "Soln\[UnderBracket]as", "]"}]}], "//", 
   "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.405141820449045*^9, 3.405141883032774*^9}, {
  3.405145182096433*^9, 3.405145183553097*^9}, {3.405145233248428*^9, 
  3.405145235069809*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SqrtBox["2"], " ", 
   SqrtBox["D"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"Log", "[", 
       RowBox[{
        RowBox[{"-", "1"}], "+", "U"}], "]"}]}], "+", 
     RowBox[{"Log", "[", "U", "]"}]}], ")"}]}], 
  SqrtBox["A"]]], "Output",
 CellChangeTimes->{{3.405141851241543*^9, 3.405141885040107*^9}, 
   3.405144359660694*^9, 3.405145183803035*^9, 3.405145235727779*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"USoln", "=", 
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{"UIntegral", "\[Equal]", " ", 
      RowBox[{"\[Xi]", "+", "const"}]}], ",", " ", "U"}], "]"}], "//", 
   "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.405141905091769*^9, 3.405141911162948*^9}, {
  3.40514194234222*^9, 3.405141996777429*^9}, {3.405144343183703*^9, 
  3.405144351942557*^9}, {3.405144384937002*^9, 3.405144386810107*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{"U", "\[Rule]", 
    FractionBox[
     SuperscriptBox["\[ExponentialE]", 
      FractionBox[
       RowBox[{
        SqrtBox["A"], " ", 
        RowBox[{"(", 
         RowBox[{"const", "+", "\[Xi]"}], ")"}]}], 
       RowBox[{
        SqrtBox["2"], " ", 
        SqrtBox["D"]}]]], 
     RowBox[{
      RowBox[{"-", "1"}], "+", 
      SuperscriptBox["\[ExponentialE]", 
       FractionBox[
        RowBox[{
         SqrtBox["A"], " ", 
         RowBox[{"(", 
          RowBox[{"const", "+", "\[Xi]"}], ")"}]}], 
        RowBox[{
         SqrtBox["2"], " ", 
         SqrtBox["D"]}]]]}]]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{{3.405141977495708*^9, 3.405141997551266*^9}, {
   3.405144344861907*^9, 3.405144387988831*^9}, 3.405145186630204*^9, 
   3.405145238607252*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"USolnWithConstantChoice", "=", 
  RowBox[{
   RowBox[{"Simplify", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"U", "/.", "USoln"}], "/.", " ", 
      RowBox[{"const", "\[Rule]", " ", 
       RowBox[{"I", "*", "Pi", "*", 
        RowBox[{"Sqrt", "[", 
         RowBox[{"2", "*", 
          RowBox[{"D", "/", "A"}]}], "]"}]}]}]}], ",", 
     RowBox[{"Assumptions", "\[Rule]", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"A", ">", "0"}], ",", " ", 
        RowBox[{"D", ">", "0"}]}], "}"}]}]}], "]"}], "[", 
   RowBox[{"[", "1", "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.405146174857265*^9, 3.405146247257248*^9}, {
  3.405147670127761*^9, 3.405147676065358*^9}}],

Cell[BoxData[
 FractionBox[
  SuperscriptBox["\[ExponentialE]", 
   FractionBox[
    RowBox[{
     SqrtBox[
      FractionBox["A", "D"]], " ", "\[Xi]"}], 
    SqrtBox["2"]]], 
  RowBox[{"1", "+", 
   SuperscriptBox["\[ExponentialE]", 
    FractionBox[
     RowBox[{
      SqrtBox[
       FractionBox["A", "D"]], " ", "\[Xi]"}], 
     SqrtBox["2"]]]}]]], "Output",
 CellChangeTimes->{{3.40514617786807*^9, 3.405146184242688*^9}, 
   3.405146216151054*^9, 3.405146248289501*^9, 3.405146289423126*^9, {
   3.405147671730965*^9, 3.405147677190397*^9}}]
}, Open  ]],

Cell["\<\
This solution goes from 0 to 1:\
\>", "Text",
 CellChangeTimes->{{3.405144903423289*^9, 3.405144920481616*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"Exp", "[", "x", "]"}], "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"Exp", "[", "x", "]"}]}], ")"}]}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"x", ",", 
     RowBox[{"-", "10"}], ",", "10"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{3.405145540443712*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwt2mk0lk8bAHA72T32/VnIUpYkKnJNG4kkCcmSJRWSvz2EKKmQXaQUhVAi
FZJnypKtEqFFqKxZHlsI5b2d836ZOb8z17lnrrnmzLk/DMX5rMUJFiYmps1E
s9afHz4ws7pKwx/1JC9vfssJQyFiv7+u0PA13xM3Hs9wgplg38LzPzT8ILGt
/JQEF5C3ef89O0vD0zUXo665csHrq0lc/UM0zN21vvH+KhdwqXfL4bc03Gll
c/TCTm5I8TtuEnGThr8c+TqfPs8Lf7lUzOxu0PAhm20bdBT4wO3WtPnWVBoe
yhHZKmPBB7oNkVZT8TT8MOninVeP+OCL+H2n45E07KNgvDXMgx/IL0YD0Wka
LtFYr6c4KwBFzL73mHRp2Pbm1ezd2UKw9XAhR8ZmGqa/Pnbd7qUQNNz7fmqT
JmGXlsGCr0LQb2Su5qxCw7l5m9Z/kyCBWLzas1ppGnalND95mUqCSKmRN9FM
NMzlYjj9KEMYbDbbj/E2U3HL0i+HmVZRGLqYcuB+AxVn3X8rqT4jCr5dLSU7
aqlY/WKW+wNxMYg/t93Pq5qKX99iunDcRQxq6eJ/3z+iYkbdy/iev2KgbtrO
n5xMxew9BvxnQQLYTxhpSdpTcaPsMbsdi1LQ7rDJeuwoFT9U/C2ygyYNt22k
Q19aUXGh2bujF8ykQceUUe9kTsWuRVM7Pt+XBrfN6UcLd1Oxl7NyfLu1DDQy
D0foqVLxsR69q2ZNshB7O/q9/QIFJ57kme30IYPNDe85jTkKDs89dlAgjgwK
SbaSLNMULHw44JxPPhleXlJ3yftFwU1tMzuf95Bh0rNrfvIbBR+NquLcZUwB
c731chF1FMxRGjV6Q5UKwp/qPXMSKXiQ5m1lJKgAeMPdmIvxFBz4SbUlR0MB
PMND77ldo+CNOSxI46ACNKzf3KN6ifhe32O7nngFCPa/Y/IkiIJ5uUZ7rwsq
wg9SiGqdIwUvn5eMuCu7HsoOaI4MqlFwQegx2afHlcHxLg9boyoFV6TSbN9F
KwPv7yH5QiUKjuxpk6M9VAa3rFvWXhQKPm1uVx6xpAxSY9xv5kUo+IlCqaPY
DRWIjBm8z7lCxnf2P5Jm6lMFi9pMF5VmMtaItGBJuK0GCh8Het80kPH7VwbX
m9rUYH5A3datlow/DDjttGFVhwyOWvPcajJWlHK0DjytDv3GYztkS8g4TjTM
rHWbBpxp05MQSiHjpjfNToX+mhDzrad10Z6Ml4uOUCVMtMB2cv2+dFsyvtSx
66mMuxZsXPWu3WJNxk+NjXqOXNGC92T2Kh9zMvb+e87kZKMWiLmq5Y/vImNB
v5lEQePNkPvr/IV+JTJeqeTtfmupDSGMHZ4JP+Wxs3+17K0bOpDlpDuCvsnj
3HDB7XbVOlDdoek63SWPpVaVew/368C/5zQ7i2Z5LMQfS5tU1oXIiHWmIqXy
+HqD87VP1bpwVahzw40weWzMZzjeydgKmZs9x7Il5bFaf+NDqUt68OL+iZPm
wvKYa989ltzHetAj7viTiU8emyw+2+jdowdyK4e+OjHJ45IxpfIv2vqQU6fb
QhuWw/fDwrNXRvThwRHWovxyOTx9pSy3y8EAKgIz3UsOElfGmQOOAq8RUF3/
u/9pnxzmhBmmll4Eceb7+pl3yeE9Lzr7ni8jcFaZP2KpLYf/eYvkr9fZCTxf
D+38IyGHKTn7njs93AkOBuskdv+QxWxWKrmZ93cBK3tQfZefLObe8Nf3Fn0P
eE2bMTF5yeIVPWqExs898Ombop7KSVncpmxz/h/HXih+1vE45KgsvpPxatzy
4F6wOqWRRTaQxYb1RpsZP/ZCQcuQjzunLO4y8zmyVcwIzJKOkP9lyODnrbur
9R4bQ7CWgLFRsgx+4xHa5f/VGPLaG/9LiJXBlxRNsxc59sM/kn4tNVwGC85c
UfzisB8eJlPdjFxlcNSk03UpkgnwpjKKEtRl8Ltp3v7VaFNoTr+iS6uVxuNu
zOHMjw7CvO7u457V0jglW4XT7OdBoH5aiXn6VBqbaJ0kBfOaQ7CE9xejAml8
BZu/zN9lDioZVuc946Qx/SY5W6XUHGIyabVPraVx97LRNE/6ITC69dJs37gU
Ll8yOHgl5jBMxfWssA1JYalpTy+1/MOQGbZc+KpPClPdA68zNxyGccdtnHod
UlgiuoJVj80SkqhP6WpVUnizTQXHcqQl9OYXawrHSOFwr9a42StHIKDsJqmX
JoXH11E/MmFrIOdW4UxZKfxpdG8zDFpDU/JnL2txKRzImitcyW0DMv7iLW3c
UvhHSjevnZUNvNZNjqqdlsSDXGOlZ6ZsgP/l1bkCuiQ+XGZ2iqRpC3lvznX5
HpPErSOynqL9dpChFbD7gaUk7oy7vkdDyB5ib/k87j0giWdC5fyzdtmDj6/H
VWMkiRv4rYV98+xhh7y9gbyiJKYN3B0l+zjAxwCU1zwpgTfIh+eu8h4HlvWc
AZRICZyT9y70X7cT/E5g/WkVIoHjW/J7LOadYGR59WCsnwROtHtU/EXEGd59
WFRZcJPAdbbb/1yycIaM87++tZpIYOGidSOWbc6g2fl2b5CoBK4+t0WB1uUC
DhdTRN8XiONPhsyaknxu8A7mOYNyxPGBdF1NQU03gGXrJXKWOBbTSQg2OOwG
ZB/pPp/r4vibyR/jXZlu8MMxp0AsQBwP+9wNnFI9CW56j/UcdotjHrG3c1us
T4HXTIvTxDcxzEt/prz0zR16H6pZpnWL4eOoNCmOywMOnr5uCB/EsEEu6Llu
9gDNfosNiXViOL20bLLnigfMvP08p10khi/kiiVybveEgAfDl0ODxHDM5V/n
7uafgbDjrI94hMVwhFFuwkV7b3j83ij+I68YDpfoGBsL8IYfBrFetzjE8NhX
SlxwgjcYyopqaCyJYlO6VTJ/nTfwf1n/2OK7KD6057GthNp/cPvw/rKMElHs
xFznPsrhA3TDxGdKpqJ4E8nnr/FXX5h52pk2tVcUF85ZzYgu+YKColRgJYji
Uw+zvIQk/SCGNVfXZLMo3iW3iXzD2g/McXmFl5Qobt/3xHe8yw/6t32qejoi
go96fCyc6fUHJjV5+u5LIjiab+rMTbEgCBXb6XApTARf8g8/LgNBsLDq/Lch
UARzGmeyfz8ZBJPteXr73UXwte36SrpVQfA1SL3C/KAI7ss5oE5zPgfP6naU
2kuK4Dd2FZ8d6cFwxs4uN/ChMNYPmTpMLzsPI3vDdlXmCeN9nLUOtwbOg6vG
ne9L2cL4puX36HdiYWDLMiAfliSMTwncOrQxNAz2FXhkXQwUxgsJArMD+8NB
YS44NXGnMN7QmKLkEhoBX65lXC7qJGFDNP1nw4FI6KgnFx9/T8IK/uV+zacj
oXU1v020iYTN22uMUqMjocb3mWRENQl3H65jK8WRkHOso8gyl4Qb9ZctuHSj
wH0DX9tfbxI+dQUf8N14EZabIyTMeUn4rej20K360TDHxrWDnYOEf1yudRl3
iIZJg+tOVatC+FS3z5fOC9HQX3arUGFWCF+P7qjY3xgNdZlV+n8+C2FRNr6N
9daXIdZ97nhOvhCOc9jPZnghBmS4Tz2Y2yWETzL/2PWZ6Ro41LTGUvWFcO+W
YKuz6tfgrs8mb/MtQjhKztfykN01WP91SadYSQj/COhaXqq4BhrFcXUuvEL4
omtx+ox/LOwye9LX3imIXY/nef/4Fwenk/6Jlp4SxPx+RfraBQlQbOjyp9dJ
EHPZfz6kX5cAk0tveniPCWJmZ7Emn/4E8HFJzD11QBA73zvP4SmZCKHaCpvI
WoL4FxjOpMclgtzD5BtyCwJYLTnnvkR4EnhsUWX/KiOA6xWevpyISQEf9sr8
AAEB/ND00sWKwhQI7jTaT2IRwMcsGPTK1hS46n/iuvEwP15v3T5wXCgVCp/e
kawo5ccKbvu0TG6lwugWcfVUQ36M0iNffqWnwUkdNuuD3nz4hoX1DVetDPDi
SP7zy5kP95jZPHlwLAMCuihZ0Uf4sGlCa7HmpQyIDoDvL7fz4Sputyb2zxmQ
9yzYYyM7HzbpzSrsjsqEQZ3p8HWZvFh0+nar79BNcNHtLait5cEztsgj8Mdt
cCvQtTJ8zoNNT3qy2Atlg7tkImtTIQ9OcepKvIyywXt5t8O7RB58sgeJPMzO
hlD6A+HPjjxYi72EL0fzDqQY+YdNLnNjklHYreXFO9BgzWspqc2NWw9PnJbI
y4GmxhPMWUrcONqq+vr21znQuo3+SE6aG9+1Fbp8qTcH2qV91imwcONPvI21
aeK50NfXXaPetg6Ph1r6FV/LhcWTuSp7PNfh+/zH3DpD7oFq0PZVr3tcOH3n
Zb6S83mAv5mp/kzjwi+9Iw9q5+SB1W6XI9ZXuPAEj7nMTEMeRPLHFoEXF86z
m0iQEMqHz/e+WQtu5cKeDWcTJQvyIeZDxOPSFk7cZOpkpvW9AIZU3zjNzXBg
tT9LH5gjiyA04WvsqUEOPFnueFvnURGQ5hnPe7o5cMe2146PvhSBwSsJ/vpq
DtyU+JF6SrsYUq3dq1KjOXBQWbwL/1gx7LnIJ6wrxYEts7su+514BHe/Ha47
t5Mdd+XyVqY2PgatAa2lj9rs2Pb3Z1Xv74+h7peQpqYyO17/4Ppbz6XHMLTw
/uYQPzvObGfseruxFFRJpr6He9jwUPeIgUZyKZQZ7qGqBbLhKpkXJpluZVBX
oh3RX8yKpTM1VJVUyuHIM+FnendYMUlgr12qYTkMVc+MpSWzYjN++ltV13Lg
an5sfSCYFX9r/CjSml0OBwY2alQZseItFuX0cYmn0CWh2Jv8gwVvsHrL3UJ6
BsMXRPWNJFhw96scbsr6CuA+X8JI5mbBU669KSUmFaAeZJzbv8KMz1sdd/T5
rwICvM5zB39nxopn4FrqywrgODb4qbiQGYvXbBW3sakEBe2n/qQdzFiY70b+
SnoVOA9ZlnxzYsJRZi2VW3a+hOjvk86qlky4Rql+3ZMzL6GwJ0Ys0JAJw/of
Gz0yX8JMe/V5wQ1MeMl85Hns3EtgCTe+ET2ySjdZWDyZWVQDsGn0wyG+VfqJ
A/q/LpEwVKUq7x2x+kvnc3K5YZr7CriODbsw6/+lX7MUn5iofgVW5LxIKcpf
uvl7pZBHXa9gppCGTcdW6AP4aGoc92tQwXJ6peEr9LiirQv1vq/hxi+RTcEF
y3S0rJuoaVoL/sAsy7v0h67rzZzlKVoPr9mwnmLfH3r7749xkpvrQbA5zNag
7g+94+7tqWHzeii2XEn3jv9DT3zwdudobD0MnF4gdVL/0D1myhPDORrAImWC
67bJIj2+qYiDg+kNaI5+ntO4NU/XMHmSWcfbBH5TairiF+fpqgn7ps5sbIKK
hQv2/9zn6Wf6pf23mDbBTg7Vhtat8/Qojzlp+dgmsKAFp5/++Jv+5vKwnj1f
M/jbS23P5flN3+s40hot2gJVH2zDxYJn6SXHZGMS9N/Cnqqv6/5aTdH1/gl8
2DXTBo8K7297ozdFF7wQERfC+QEkbp49nUieoge8491cL/MBuhY4kkw+MuhK
d0L4rxl9AMUQK8FbS5P0M4WHWvpufYDXF37zIqMJeqCkUX7OgXZYid/MEd0/
Su/RlTqlWtkB3g9KFoWEf9LLa66MG1Z3gZUb7W57aBNd/Pq+jyIdPfCm9fPV
zYcboeDiiDfXxe+wXe1VmXn6TxAslIkVwQOw5WNXWk3TT1AR9O5MaBkAzZDx
4I0rP+GjvO0Tzu4BWN8kvpvLaQA+LqawDk4MgPAJrw666iDcyNH/bCc9CONZ
0nPqNUMQtlGhbqv/IGTzBmzhGxyFrxdfjDopDAFLgGiPq+QUxIQ9SXvrNQy2
xwvbK7SnYI8c8ow4Nwxl+6GJ13wKrm388kL14jA4y59+Vh49Ba9/ipafyBiG
143VCWxzU+ARxBvqVzsMUdKue++9m4aBon7vadERYHtV9uhn1Cz0nh6tPFg+
AnZFRvd1s2fB/siOpdGXI1Ce2nPzWtUs9GVX7Qt5MwKu7hxXNk/Pwte/JdxR
X0agXtjWJcphDqY+bWFRYhqFaDdmCdq233CgIE4CG48CJ6/5BefJeXhWGiKv
2zUKNkq2WGfdAvSnfRm/2jsKhbtcV7kVFqDUrJLSNTQKZueCzj85ugBvjNsa
LRZGIW3ozjm2+gXQ+voODYr/gvWvpv/Lu7kIxpHi1d7Wv2BvQLLzr31LsOFV
c7N22y9IS7x1t8ZlCcpuq+wK7v4FI8X5/UlhS/CCrh7yvPcXxP584aBXvgQx
gi93Sk78go/mA7ax8suw3/R3QvC6MXDdsOWw+vwyzLzfPfAcjUFUf9dun9y/
0FKz2XDHgzHorTQ5E1T9F/R6Am+Jl4zB9mScFtb5FyKOGhmPlY/B1N7C0auc
/8DqWk1JGB4Dh6Lz13M9/gHf3Mgmq24iPkDh68fNq5BrHTTRwjoOf9i3VUo7
MyH5rAS3GetxyPOZfXvbiwltnc3jZ7IfB8u+hz8oIUzI8R97NqfzODyuoPIq
pzChszwbYtk8x8HNg9dRu4EJGVbGU56Fj8PHtj62A6rM6N1e07iV+0T8zWjz
sBlm9Pb22+Btk+PQIhw9LMDEghidN1XHZsZh6NqlsLt8LMhD1/JM2sI4yIRe
LK5TZkEyTyPrPzNNQIxdJBePIwsiz4wGiAhPwHHZMJzewoKWb3yPd9SZAMFs
f83Se6yI4zxiDQ6agA3i/m92lbGi8pH799xDJ8Dwup/DRzorSsylFhyJmIDQ
cN/YhS+s6Lb15KRUzASMOP43skOIDaVeXa/ukj4BmHzmTvN5NpT86fVfRvkE
nM1xFRywYkehJlty1o9PQPy5jKx4F3akVB47qsuYgEfm75S3ebOjQ1mri3tm
JmDyn+7O+CvsqHO6JspicQI8j/H4bK1mR64Vj9r02CbBXaSsI5bCgTT6CsJM
pCfBLXo1bcs4B4ouZ9WJMJyEaHttWv8iB5q4pPtC2HgS8rRPl1xl50RM+NHt
XJNJGPzR/qZPjhNd4DROqjCfBFfIW7xyiBNd5ble8Nx2EpwXTW17n3Gi1pdX
T/F7TYKDe6bs5QtcyIH0TU0/aRKWP697Nx3Hha7+GFx3NGUSbhifC7PL5EJZ
e6XHfNIm4YOyTd+mJ1zI7RGt9VbmJOwZFr3bM8CFHrbF/teWMwkqrokKWvvW
oe7ALZeKyyZhzj5mwzd+bmSgVj9t0zYJCW8XvhpJcyNGkFqEbPskqO04GVum
xI1Gc+su9XUQ+coYTlxG3OjkzQ5Wu+5J+PyFtUTLhxuJmawo6fZOQo11uFZM
JzfafthkNGhsEq4cCti2OYsHWckKOiE2BiS1qnoWP+BBhug26Qs7A7KM+m4r
PudBUbMHDf7jZMDjHfvYJNp5kJjsLOkGNwM+qUi9X+HkRXHd3BKNggxQZqlx
bfDlRUxu5MpyGQY0lLEnHDXhQzpqpQWftBjQplb1ut2GD/WKZgYpazPgS4HX
bxM3PnRhZ51mwBYGTN7utjW4wIdOspx349nKALFrDxRpz/mQwqUlA9kdDDjh
cuDFOJUfDRmesvhuyAA2kdSh8CV+dCXx+51VGwaQJsP4SZwCyIQteXr1KAMo
jad07gkLIE7aPP5nywCDUP3oxo0CSMohaGDBjgHnBn4qCjkKIDfj48GfjhPz
l2udyKkVQN1fDSosTzFgXHXkP2dDQdR3iD/HKIABc+LfLhWZCqJgD1KDdCAD
ltnaM+YsBJGccGfBJGHuvhevoh0F0ffy0dKEcwxQSr4uVBQkiGpfCrA3hDLA
eVmndLZQELVaHNyxEMmA7pZLUxcFhFCc06HR8HgGpBjc2ntETAhJ7kv+veM6
AyxKyzMVZYXQ8DaOxiXC79J+7mlQFULlzLZvvBOJ/XXamcFpKIROHxdwM0th
wLOFlZ1XQ4XQ4p/e1g8ZDEij+SVfHxVCrVsXfrPfY4Bl2rURxykhpLh46vgd
wqR1uTs0F4RQ1bmLzdvuMyB+8sNwGxsJXQzecu90HgOiKzX0SWQS+vLGLvNl
AQMCDo4NpFiRkM5q9pTuQwYcCXHWyawloeKscHrEU2K+WyIcC00kdJsW+/43
4S56Q+fhNhKaTjBXcX/GACu2Df5830jowfzRGvPnDLCJm30SMU9ChS+iYkUr
GWCbfXHTSRVhdNfx+b2wagYcr83bqH1dGG0Q2fWh8zUD7gzarCSkCqPM/t1i
OrUM+M7J0zpxUxhVa/TXpRJ2Nj3rkV8gjHwfXnQ+XMcAl07dQqnXwkgwOKji
TT0D3IYblZh/C6MGn+Sk1EYGePKMUd/ZiqB5p+VPLW8ZMHauio3HSQRV+z3x
lHzHAPeRK0NGJ0WQ0THDK26ET9YrF77yFUFVV985rhJ2Dnfb9DRWBI3dGUvd
0MaAo7PfIatGBJkpSdK92hlg+PWTnTtVFE09CE+/18WAeuMCg3xlUXQ/5UDi
AOE9FYHkAXVR5OiXq6zQzYCdqWID9nqiqNTcQCuHsP5BS49DlqLoUrDzq4xP
DNCqfX9ua7Qo+jq33izwCwPkixrSOH6JIvr08XHGN6K+2RWVzVOi6OPUcVmV
XgZcTS7siV8QRXyseNKJ8FxIPFWcXQw5FIr4txNuNrV6pEQWQyHtmpSSPqK+
k4P1+6zEkIwlDrX+zoC3m9h/X30thvK39eqdGmAAy/oF8YNNYsi1aU9GGuGt
UqPbhdvE0C5/9gd1hHNY3kbc/CaGvP/c+C4/yICgjhTe4kUxtC/md107YQV/
BYW36uKo9sfeV6rDDAip2H1Y4KY4mkqLf5s3SuRn1zbleFccZTvEfW8i/GvV
Lv5xvjg6tz23ZpywrVFA06FycdR4+3Tnpl8M2N5VAClvxZHc1cywSsKLs3wb
JZkk0DuuJ5z0MSIfjU9stBMSqEbmqsz9CQaod7jm+HpIIH7XCMYrwkMB01D3
nwR6VWGa2kv4SA13yIkwCaQzHmErPskA7QM7pvPTJZDfuen2aMKz7jnfNjZL
INvLbs/tGQzwzvN8pqMhiXiKtotOTDGg9710Kr+OJLqsJvSSY5oBB/40+w7p
S6K8O4e2kQmrmqpsStsviayXdB5bEB6YHiz6fUIS/TcXnf2UsPUOx7vlWZJo
csTkte8McV99PBSrxS2FfLlp0+2zDCj+y+TBLSSFHjzVLh4mLK302PiHuBQK
HODcs0J48Rw/Z5KiFBI//p1bcY4BZfItF6aRFBL43bbsR1jRY0/g40ApNH3k
Rxb/bwbwsug6qw9KoaBn4baq8wyoDGQ9PDQmhTi+WLduJ+w28X737RkpJPir
gmpCGH86pcjPJI3UElrPexD2LckcmZCSRm5VyjOFhL/Y/fN6aE748NYlhQUG
FDyvD9lYLY3EU3Y/5Fgk9lctyXPgtTQ6I3ySVYQwS66DfVaTNBrp27mbQtgu
fsGAt1sa2TeVRugRFjihwjI+LY3OWoTOexEOIMXFFCnJoOIylokPhPeesUxV
TZZBDmmaczF/GOCX9ILpbIYM2qB9gieZ8L3nVM8n2TLIT/Qb6RZhFpbpXfrF
MujvhcWZx4RxWtzUgXoZ1PQ5cqmb8HbcYPLfggxKj0nSoy0R50V4G0vFMVn0
ZrmUVkLYfuudMytOski5m9P7OeFYe87P6JQs+lbiUUInPJbfWdLkJ4sUhx9y
thEu0Pex/xonixrOG5OnCNPciir+0WXRw6vfWdSXGSBRKXt2r4IcEop5KHCX
8D+XfGmhDXIo183Gv4DwAP+mxp5NcggpR7wrIfzYdQ/FD+TQ4deZx2oIGwl6
dOTYyiHFpN0/vxAOOFmhy5Qoh8a7N1QJrTCgU/Qw04u/cqh3dx1fIOEq3FMc
zSaP3vz0Wnee8B0Pt6MWPPLIV2tpKZKw56tzZaMS8uiK5gI9njDbmTuuEtry
qKTCueM+Ye26iUZ/d3n0S9RVvZ1wkk9Moma3PLrbU0BV+MuApdbiiphv8qjO
LyhfmbCz0oe+/p/yyFJDQVGN8OYeSfXEKXnUQtrBo0P4V86RwUU2MpodKooy
IhzLO3zhqBwZcVxdVT9N+F3fuhdSh8iIT4dSkU+YhS/rnaElGbFeYj1TRFhn
u/oPH2syqjrELVNC+HaKxboWezISsi489YzwWeOb1iHuZNR2efPLOsKCTzbO
fY0io8Qrusf6CVtEH1S79YyMbnpk5wn/Y0De0SnjykoyUrV8dFeM8J+NiW6d
1WRkPSR3Q5LwnY72bL5aMpLO2BkoT3iCbEUKe09GM80Wq6qEY17YLTiMkFFz
7EP9nYRrpk6/IktRkAPZM9OdMKmOu1dfloJOV9bwnyF8Ir1oyYZMQf7Fu8PP
EuY1mNicuJ6CdJRXLP0IH73mnce8mYJ2aTZ3nCc8pxh47YcJBT1X7KAlEFY9
dvHIvfMUhEhZm8vW4v97RmK7QIzXNrqWr63n8sh7l4sUNLjBMOkZ4cEnpvtp
1yjI5YvJj6q1fHhFIecGBS1FqlnXEhatua985wkF+biHZ3QQZqK8Wb45SkH/
3nAXTBPW0P1TsTxOQcMacV6zhB0ObAg4NkVBNvHCWr8Jvzh3fUpqgYJCdNDj
RcIB7daDGWxUxJ6ae2WV8FjUyLt0eSrKLmku510l/g+G1uUkH6GiBoMaXkXC
DZZXsiptqIh1ydx1PeFnr7nS+45RkVVOY6US4bTbnLEbnKkod8zXQZXwESv2
gFovKrLJFL2uQfhjPdP+2ctUlJZrlrCNcPu9hWmLKiqq4v6hZEr4NSlwPOgl
FY02S+w6QLgsYn7oNqYi1aoUWzPCSXa/v/5qoCKXFbdoc8KHRGbrIzuoKJ4J
N1kSbouazCgbp6KTM3HM9oTfOQ/uJMnTkAdV3MiTcIbwfE8IlYbUMo13nSHs
WscRNKhIQ6dn2vS8CC8rKj+q2EhD0jdGVbwJK4+6Szlup6Hv8wkzvoQvnJ2a
KTpCQ57JUltDCGuHLufujaWhotideTGEVzfyQsl1GlKZ3PLfFcLN32S+SCTT
kN6WD9uvEnYCA8HxDBrSvsnefI1wPMuF0KR8Gjq2eu9TPOHhGI4jfa9piPlC
UUUK4cxUIY5zf2jolkFUUjbhH9GvKcmEFbrOpK9ZJch3x0PCF9RsstZcYdvp
951wo/TGvDV3yWf8NF6iocH5tso1CxeSX0kt09DOQv7va46t0QitXiEsY6tx
Z60+Jf3pXYQby/W11yx9N/HJFOFYPvltay6Mmv2l8JeGelZ/7Fpz477nR+MI
M4ZPWq2ZrcNA1+EfDVU1nji/ZtM6hkUQ4Ry3PZFrTn56xyuJ8Eod9fKaqTdY
8hoIC9X2JqwZ2TeIaKzSELn+8L01x5gFahoTnnLTfLDm96Bs6kLYtIn30Zod
qFej0gk/eV/3fM15wnrZpYTrfe9Ur3mSbbyqhbD3p5BXa9aZz+oaJPzsr1XD
mv//Pgbd/LmpZc3/A+X7FjI=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  PlotRange->{{-10, 10}, {0., 0.9999546021127687}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{3.405145540921424*^9, 3.40514631960702*^9},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztXQeYVTX2f8y8KSAgiIhdUbEXdC1Y6SAgAiIiKNJEBnQVEV1QyqqALmsB
FVyUMoxIRxFwFKUOnaFXwaENSGcowlDkJZvfuSdv8pg3c5Mn++233/+f73t3
5t6bX3KSm5ycnHOS1G/dJaXdX1t36dC29dW1O7fulNKh7WtX13qls3oUXyQQ
iFsZCBRZc3UA/8tAQF/ywmpczCf/yQel6E+ibNmypaxUqZJkeuR9992HZ0F6
nSxOnDghu3fvLitUqCDvv/9++tujRw9x8uRJocIwjjZhwgRZpUoVee2118qH
HnpI3njjjbJdu3Zix44diEYZ1+G4CIXETz0nhCVGlJLyVG90YuLIkSN4FE+v
i4pQKCSffvppWbRoUbFq1SqKvnTpUpmcnCybNWuG28BQjjt8+HBZpEgR2ahR
I3H8+HGKu3fvXnnTTTfJSy65RGRnZweM8hazAeRV5TkhlQt/Af7EUSuLF7m5
ufpLBOLoGi9mz55NlKma52dBitOmTRs8FxkZGUM56rBhw2TZsmVBuBlVjBs3
jpJo1aqVUeigGR/RNMSIjsep54BCy8LK9u3bE0SREFGC0aNHUxYdOnRIMx7j
C+FxmTJlqGAN+J3x3KYENtlyCQL0icX69evxU00c13hx7Ngx3XUDQb/c7rjj
Dspt7dq1EbmtXLmScrvzzjv743GxKIkneFeKf+rUqXADVI3Ra9CPc4J4icSQ
kXrpWwUWREVUAbVl9Qt4V68VWFdBqVKlKLf9+/dH5LZ7927KrXTp0l4VREnc
rIKdO3dS/BIlSqDTRVbBWS99q8CCKNt2HBcXR0mpTxeR1NGjRymp+Ph4zZVX
rFhB/IIDRRsyZAhFa9CAWnS4wyJwfLPDGtGtmrsFdX+2oL///jslpd5zQWVK
Sop8/vnnRU5ODkVZvHixLFeunDz//PPFmjVrjIIWo2x84v+pkhrk2ZZUZUtJ
7du3LyKpPXv26NZhlFSMGjUK1Mtq1arJm2++GcxQbNy4USdrDrSFxU87B2TZ
llCzgNWrV0ckpe4pqYoVKw7jqKB42rRpdBfR53Rx1Ev8P+oc5Ok6dowdOzYi
KXWvmbhJ/uDBgwsmX720It8izyhCDsQizbg5VU9yyBtL1eePGEvRNHgsNQoh
u3XrVmAh8BL/jzxnOUcWRTZt2lReeeWV0RP0hCAl7ICDi8zMTM28aOBo3rw5
bgNmWZSwJG+44QY1eChBUxEEtsvBK1/ts8rHACN++PsMPXcknvX9wHohMKen
p4t+/fqJ06dPI5qWeCEU9uzZU15//fUk8Sr6ZK9evSJF8XjZuXNnqmxqWXxF
JY8cOTKvvI/yl0GhAEC0BG6MKuLZhf2TpHE5KeeOXs7m1Mf7p93/R3KJ5ElK
4WnkclyKeI9Ke1Cxa9cRPBWZS3YELqBrfGjypHVixvRf6e4875mOQ3fBKPH4
mVi4YLuYOSMLPzWO4hof+vmnTWLUyBV0F/SeUXpTJq8X301ap/LANT70Vdoy
0e/92WLkV8tV58A1Xn75xRI5+F+L8VM9E9ckOfCzhfLDf2ZQGt3f/FG83vl7
3Mu3e01Xcf6urgmy+1s/ydfUc/176cXvVGovtPsmXjZ7erSK0KTx17J+vVRZ
veqX9KvyyGAVQT2hCHGIqhJTUKo+lSh+8cgnjshIAIGhzz5dIFJTl1GRUQFx
4pcNe1Uhf/11v5qUnaCfCudH+VBXRXl2tWW8/2vYEnT1xuBzmNVFdA2qYel0
+Gvt23csTonuh9XjzVkH8DkVa94BWBy+cjy11ZJi/LjVaABy6JBM0bv3DPF+
31nUSLqotta501Q0H9mm9QTZtMkoamZVK3+Bv3HU+krIds9/I3uoJgoImq5q
UUhOjBu7KoRuMTdjK2W9e9dRxSD/QPaB8v+DH+6/hf0PNZho8byR+VI1tp8W
hw7lou14rUZxS7CFyd+tExMnrBGpw5fKTz9ZIN7s9qN4qeMk2TFlknz6qVGy
ZvUhlN5TTb6Wz7eZKHt2/wlNgpoDWhnSQXp79hxVs+oz5f4Hv8Z/C1tIKxhN
1wtJxMI4go62aeN+sWDBNjF1ygahxiL5yYD51Elf7DCJhosa1b5E95WvvjI1
/IXQYTPmbAktX74TaQ2OktWYKM/GFkztAvpTioTBM2dC4uCB42Lduj1KQN4s
xo5ZBarEG6+nywaPp8laNYYSRaDw/fdo5BTz520FI/s0Sg50ZdlgOf0pGc5m
797f1axnNwovPx+0SCKL5irpxk+MpOT7vT9HTBi/mipJhQ99kl9Pf4qrme8p
SlYxNarPv744WdapPUy2aTVBYvD8ZuIasXTpDuo9Krznk+wm+nNemOqdOw6L
n6ZtQlLy5Zcmy7qPDqfhftjQpaHRo1eoifbJQG+fNLfRn2Lijz/OUEdDZ23X
dqJ8rG6qTHnhW/rU+ACnTv2BbHv5JLeT/iQL6Eh+3bQflUZ9GqPAKy9PkUO+
XAKyevgks4eTAVUrVvxGglDLFuNk2zYTiCKwGRW6+SRzgJPZoSpq1NcriOu0
V2PTxx/NE+vX76VKVOF1n2RyOBkEJTLS0AdKuK7F9u2H8Oo1n2SO0p9EsWXL
QWpkaFxvdZsmpispUg3Cr/rAjzEcRM+fv012/duP8tnmY0Tf3rNCe/cefcUH
nstwhN9+OyI6vTJF1qsznOri+PFTf/WBnzTgaK6KORNb6NtnFpr4iz7w0wZc
8XHIEbJZ09Hyb6//IDZu3NfBB37GgKNNfD91g3ym2RjZ+92ZkFna+8DFWXDV
7+Szz4yV/T+eh9bYzgdO/xcx0jh69KQcNHCRatUjSI5XJWqL13lJmMh4A6la
IsosWz43jsY2FVrb5J5gpDFv7lZiTmh66svhUUubNJLo/yA+n/zX54vQgMBX
rbBFGbtu7R7ZquV4+c9+GRjvW9hj9bf/g/Ju3XI8cT4VnrFJ4zwjDVVv1HcU
UwidOEHlb2aTRgkuQ05OrlASqlCSqpqmn3HAUv4hsMLPBy2kvguWqUJTmzTO
5/wRFi/Olo0afCXU7NAKW5qxiuVT21HcDIygiT3Wqzs1Nkr0/J49fibBX4XG
NmmU4fxVS0enw6CIduSA9fJX44IcPmwpZgd6etjIJo2yRt1NUrxDDdBi+7Yc
K2w5xqan/0IMc9s20mk3cMEiZGRskQ0fTxNZWQessJcwFrgnGn6lhYf6LlgE
1UfR3pGvFfYyxqpxBfK1rud6LlgECEP/eH+22Lf3dyvsFYxVIgQ0DtTPVKjj
goVaDe1DSUrgzISNzlWv8hAYv2gY4LG8tk1ujAVGdn3jB8XLFzpjERSFssUz
Y8WUKeutsOUZm519iFoSZEMVarpgERYt3C4aNUwLZWXtt8Jey9hePX4S3d+a
RmOgCtVdsPgyvd+ZqSS4zFiwNOqokUPMmrXZCluBsaq3Ul1t2XwQyVR1wSL8
oHq9Eh3B7aywN3iRZJ/eM0mwxl0Ve2CQGoWaloi5c7ci/8ouWIQ1a3YTuzhw
4JgV9iYPS93m9dfShWeeFA87YKkjgKWrhuyMZYGSWKuSxaywtzAWYiDEEO68
D7pgMZQoCRyTUGcswu7dRyFDYhixwt7GWLA3sCeu5/tdsKhnsEYojjU2OnvT
CITNWQeoCyiGapXbHV4kElHVxJjuKtkDg5AqIdyL1at2If/7XLAIqniYIIHP
WGHvZIK7eJ+T7u61BwYxA6P2d/gwDXX3uGDB2DAthaoqBiwxtqeafI1ZshX2
L14k8eMPG9Ff6e5ue6DXiCC/ciP6iwsWYfmyncSglBxrheVSybfenEYyFO7u
sgd6LBw9TXFG5H+nCxZfp98/5tAXigELxSFN0/bvP2aF5fZKkxSIFrjzBFhf
jzo23MratWvLRx55BHdi8uTJiXkmRpmYmAgzKNs1rUgJolHLBvVHkLpHBU+c
t/Nz80hKlAkJCWTihFXyyJEj5LhQtmxZtt6SVTIYDIp58+a5kYYwInUZBIM8
0qyc0rzHybJTp06yY8eOOivRu3dvWbNmzUgvGzLDklvDRx995E4h9BX16gwL
LVmczRTaO50FUWeyevXq4sMPP0SMIFNkuA0FgjYUMesW/3hvtoBRD3elmBxr
BzDv8eHDh+WDDz4oBgwYcJ591t7Mn2cSSIYNcza+V8n8GL47devWJd8g9T9N
7RPCmbeNkicChkzomYw8LdygvDy13uCUbNiwIdqrYrTEqJNcSu6lcQYWKczZ
/agwXJSK8WM4AqKfXHDBBXANtMqd5Qrx7TdrZbeuP9JdSc7awneomPE4LS2N
HtepU8ch6yApXDCUb9y4D8mUOKvBFeL64+WeENESkpKS1PBzxq4LPmiki1Ed
2jaDBGv3nUTZunVranUcKJrqlto5CD1UU2HVFU3K0CcaNkiDCMGU2XvjxMu7
7rqLmDp9RWK58Wq420Zw8Farb8Uisxg+LFN+8M8MuitOVxcvmWQxdepUsHHZ
okULJXqRvl527dqVaP7iiy8caAlCuafHTaRTnLOw8mPx2PX8+fNlvXr1VA0q
7qWq6rbbbhNDhw4lV1EnUlCMTz9ZQGZiFW5yxJKgj/Hz6NGTVtjK/G2h2oaq
H3c32gODYK805eRZzQ0uWBT29S7pmLbGghWLFm2X7dpOjAmLGRWUi0pAtMLy
fJqMPHA5wd319kBvnoEvo74Q8q/ggkVFwYqifrFgw61C0WCFZaUF/GyowLi7
zh4YFGvVvB5SMLcK0o1En/tpBIrYMWUSFIhhhG1uCLANdnplSkxYDJaNG30F
tYsVlvVQShpcGuZh19gD8xgOvAhUiObvUSAWFaUyhYU4FixZ5NDpVAewwrKy
jyp3zuzNdHe1A5BAembuAAxiRk3Koe2eDvsqFyxqCTz065ErYsGS1hI2C9U0
rLCs+IVHEHxB6O5KByCMI2RjdQR6VgqM6fAaUeEKFyxqaeBnCzEYx4KlyS5k
bPWlrLCsWieV0ty5W+nucgcgTcm1cOkA9GYDjz82gpwfVLjMBYtagkSfOnxp
GBudmRkICKHw+gELtMqNzR2w2ZPFDHeXOgDFpG/XkpeZI5CVuE1Ha4XSJS5Y
1A0sEnCliAFLQ3gHxfZjwUI3BLat+IMVlg1ZNJIO+TKT7i62B5KvJ02qeFSL
5qBVIBYVBTXhiuW/xYIVo0etJM+cWLBKroVTCswCVli2NJLhcfrPv9LdRQ5A
Goxh+XcEet0UlMIVSoWyLljUMFyjeGB0xYplS3dqadIZi6YIvejmrANWWLYl
k68PtG+4u9AeGBRbt+bIJxXn5aZYxgWLikp54VtwmViwJBFivqCx0XmhicjJ
yYWTFeRPq9zY1A9WLzLmbKG7aL59BQHFtGkb5Zs8TjgAofk4BU8hzJ9AeGkX
LOq1x1s/YZCKBUtdresbP8SERccB3YqHW2HZEUP26vkz2b5xV8oBKGbNzApr
9B2AwTClLPtGc5gvEIsafuftGdp84Yol74LOr06NCXvy5B/Uhg8ePG6FZVcb
Gh2hE8NdSRugJ0R4HniFaOE9T6sErYGn6PC+ePWVqWFVeK1atcxVSfKee+5B
NCuNhCbDR/XODl8RWnfo4cDFDxw45umJEtS/BzQV0G/SUwcygqaaHWC9Ls7Q
sodLzWp2wUG+p3iVkk+0yio3N1eTEijqQIKpRzdIMNXovGRY64S2wVM18Kfz
NRTmZr6GvpzzTYTOSV5zzTWyTJkb5YMPPEy7KqjGQcrDONy51LmpKzczNlTl
4Yyhpb/00ktFSsonobFjVjJPyqWMnQpraKnNPA1VeThPhNTUVFm+fAWR0n58
/jyjK+t9cjLU4RE5nTlzhjanuL5C/dDu3UdCx48fdy+doQg38zT04NwrgiF8
3oSEpFCf3iNjy8tQe5t5GVrvcF4wQLVokaI++OUh9bHd8zL022ZehnrbzCt0
6NBhWaxYWfHMMy2d8zI11iYnMBTWAe3B+UK7D/TqR8dskv2U0WyvCGVmZqsm
UTWmLLxhpUAls84CHkyP1f1cXnRROZTGKosWXAPgf7BFYMWdCq5YckhnodgZ
G56fbT9khW3JXxJrZ6Bao0pyAMKFg5YbOAI9LwMYjthb0opJaqzwdChwSYoF
S6pSKKJjwUKJjQmpmnZYYdn9mpxlhg31atjKvMVALCSAE6Ar0Ksl6BiwTkKF
BFcsjMjsFuGKFb/8so/mOxobfb7SlusG3iKKSg9vk5UGYh4HvzRXoG4BrkDP
KwzNlhWk8S5YVCoWGI1IXRYLFnoMrG2MCYvuBuXqb78dscK24xrGahjoQXAX
5wAkfRqWNjkCE+EyT2rDfZ7a0IvBcdpHiU/iZt/ZWh1mE58UoUpot46vJyKY
VBcU3/s/KN/++3RakOUTj1TWWAfmFw9GLJ4eFxpP8FT4xx82+sUj37C33pzm
Gw9qSpV31Hi8pgZl0O+9Fy/kixSE42fACEYFR8ZTFYz3qHD9caLGE7yykJ1g
/eLRfMQmHhoeGuDZ8czNAODIumjR9vyFMCKJMaNXwmRUUCQSk7HYi+wIrFHK
LzebkUE+3PKWL9tpFVmkjaBVlnaRFT8j5Zbib4VGljrA11pxUbvI4NGfD1pk
FZlGAizkLiByUkRk1Aqsnux/Yo5EhYMwumF9uQsIPphkeXcCYch46smvYQUs
CBTJGrV4jSWcWijz6twOSBIkVjs7Aj2hDOY/LNjV3aWlJVawyYZ1wdGxBDW9
ceAsymOaa25hYReGQxssD/Q8Bhv1+pwdUKxauYv2HXAE5hN26WULS+xZwq4r
1hN1vDbrjDWEXSusdluCuUcLu/TmWTug2KCE3WdZ2HUAsrCrxEheBEUvn3HA
fjF4CSwDsWDJ7Nhc5W2L5XkSCkoFDhe2uSVQ8a2wCSMf0EmXGW/qMrlzJZML
LzR61113HfxACQanry1btiBetHmPduCEbxJMsLx6gV42i0KdlYqzhKnixBO4
ncnKlSvTzpa//kr+6VCVyVtuuUXN1S+Cb1xUHacmDyShM6xbu6cQ8hxUn0mm
6pNIxGZviNOnT58IbjVmzBhKsm3btn4UYriEj0AhFNpoRnVvbNy4MUVdtmxZ
BEXZ2dkU97LLLvOhiGxfYGBgCAVRZKEz1RSVL1+eoubk5ERQhIbHKsdoFBXT
n6HVBFI8SN3+n45CjoUmVZODRoeohw4diiBHtTCKq1pbIeTQWAIHwkLJsVCy
EjkUK5AkJk2apDMiWj799FOKqPqh18JJaNLfCAtUMBSyJEmwpvZUmApY7RSb
lJREzxT/0J8UPTROXnjhhQqieEG4IqKpBwznYtJmsRG+QMos1LQeZUXlFVdc
QVGxH1r//v1l/fr1af/i1q1bi82bNxMrK/h7eZ8WjpfsMFwgSRbaXN2CYPVB
1MWLF/NWIwVnvXPnYZqKQygoKGsL5a4eMvv27UtRBw0aFNF4oWDmLfGiVUZx
phy+L5mZO+iO3jyVnxwb/a8mR30U6k9Vq1aNIGfGjBmwi/mQI0aPWkF7NhVC
jq+eOJmjDRgwQNcM1Qac/GvUqCHvvffewskIQgCmJXG8kpVeNolCCUIh6uRk
Ixp2Xr777rvl7bffjqaD7Qv17tyFUSLYSZM9B/JREk1OMLC0qw07sThj9+79
PWJG7IfVaw/ga7Zgwba8j/ikHZDcHN7rO9sVGAy72/A2MvSysSUW8eFpOc/Q
JThgydIMi3Ms2IMHj5NiCQomG6xe0AOL+mz2DKU3T9gBYcaHOd8V6PnPglLe
S4heNrLEIj6WDmPHrhiw5C1hMut8WIKaCCiu6tYeHkL92uTGDiuiy2vfi/Tv
N+RVT0NLICzPPXv87Ar0fEkNL0l62cASK7zFi7SnVgxY7JCiLU3OWHK/8rQa
Vlh2tSK9CyZb4Yp63ArozSwMT0B6Wd8SK9gTkDdEcMUSE4WfWyxY3cXhhGOD
Zec9+FLR5gLhinrMDogdEGiXREeg1xRpuYSh/alniRVs9FSycCxYsnJhhbWe
YPhh2YmUzHgRU4G6dkB4qZOawhEYpO+JhQS82phe1rHEIj4WeKgho2AsQQ0E
bUiE1htDblSf0E2s9/QMvlj2Yib5GOuiw9XzqBUwCDfSfNVT2xKL+PCOm/Tt
2liwVFEmP3PBgp81MyrKD3u5UVEfmAJrLSugrqhU2sRPE1zTEov42GLJrCgH
rFiyOJtsCbFgdUWx5soXyysraNPHiIqqYQUMwqWNKor3IqOX1S2xgo1ipvTs
gCX7ksnSXLB6hMSeWjbYq7mizCVO9KaaFTBviRNvyEQvq1piBTvt8W4nrljS
gcIbW09q82EJyiu29KqzvCJW8cmKgbTzh14K6AD0xjbYt9gZgl5WtsQiYF5q
GhFdsNjOj9e6WmF52R9t6RkxRD1iBQyaaw/DBD9sidVDFLYjjQFL5gYogLGS
zAbLCzqxeyd93XBhH7ID0rQNe0Q4Aj1XDpgOTZPBg5ZYBExW4TcVC/bEidNk
aNuy+aAVltcH02iDZhEu7ANWQG/BEZoEtubUBN9vYM+B1WDBggWkisORYrAa
4C9OrVm8mNZVR1NyaOoE7ZGTqTdxKJC6GK0Gs2bNomfY04JP89J7XNDq+jlz
5kRVGBrkhZkru4BEIe/PWA3IzRVxFKkRarOZM2dSkpUrV45G4c2cxaCBC6Xe
9IPeVMpPnovJgJW/+TTiOTk5FFdVeMHk5AmBhg8P75gSSZGDyaB48eKF7mdR
smTJwigSbPI31QjRKHKwGmDHCESdOHFiBEXjx4/XmnofiqgxYfIPhl0QRbaG
g0BRsWHDBhxoB7ud2LRpEw/MG8mOV7FiRTo7Ls92oNvOkC+XRKpc77WnwjQc
eIVNJntckyZNSKMP1gEomrpiKwXxgVu5TrAVEerD5FL3RKHG3lggsrKy5K23
3oqtIujwPzC3Ro0aEXzgwIHEwQojSXj2ONHn3ZkhvSdPNJIcjAVYIYCoinOF
Gw4OyUK8yy+/3K+SVCC9VvWqQ0J79hwpkCJ7GwKYIEVVjSpMkXnInPo/GkW3
c5JY9R3hUHF3fnIcbAhh86Bqw2FytLHy4osvxjBQMDneWAcDC2/pSC//ko8i
XzMCZ4GtTSjnlJSUQBxtLxKUXbp0IQp506bCqNGiAjYQKIwaBAtTAj7fAw88
QO1N8UdaX5Samqo3tPKjBKp8CD1s8omgJJr8UJGxUCrBT8fUvt1piYWgBYlf
peGMRWAf87AGzQ/Lm8eRJyb2vJG6WVa0Anry81kuqLz/oD9WizLmZMEBS+6e
WN+v9RB+2LsZC601tn40HS9us8Ri0pbS/tsIHbQtFgG0Yn+Lw4dzrbC89yIZ
zbEcN/yBbrUCej0csznTjeMWSywCVlJCc4hZYVQsQXlPS9LE6m0jKcrNPllV
4qywOSCE7q1b89y+brLEImAnanib6km1CxadDlYtOKLaYB9gLLQVsHzyboD0
8gZLLAIqCnM5rch1wertS3kvBl/sQ0Y9wyK22dt3jl5WsMQi8FaZYSucH9ab
hxODibRlXGcFDJJGAEszeIt+enmtJRYBk2z25rfCVmEsNsp7vs3ECCNeeUss
4mPTYpOxOWApz7ZtxocmTlhthfW0XsHQjuxDOIoltHVL3mh+lSUWgyIOYwDN
2hToh63BNGOqjeMjTJPyFZZY4WnTaOG07ggOWGqI2KKDl4n7YmsxFvMJjO6m
feGys7EErWXkhkMjoIzg7WHyI87K7VHG4rAGuK9nGd3uEkssAsyPYFF62PDD
1mUstieA24o5SyhniUXAxvUYYnewYtYP+xhjsRMOFLqmWFDWEoswZcp6cqLL
Ym9aP6w+8xNqMfR21evD+ZaxxCIsWridWPoSZqt+2IaMXb3a21HMHGVLW2G9
VTPwGoalAOZpUyw53yeNJzh/KP5xpA4kZyjM7LF5q3ZwmMTLL03WhwVQpBI+
aTxp1B3WPsGxDm3VBvsUYxFf4UJfpS0P6eH7PCusRztE0Ibqm7+rOIhpISrq
k4bnehhUE8RjAuIVZGDs+G6Dbc60z5+3lZxkYfjAvgIISVZYj3Zof3AoEziY
2WYTfNJ41kgDq6NQ76DBnHPE+6TxHJcBPAXHCbzY8TtyjoyKJWgrXer52+jM
Hyi/tT6miE9urQyKcewQbFzQe2FzL2xn4gNvc1al9f94nj61BcPlGR/48wYc
bhjwP4ENfdnSnaiz0z7wFww4JHgc+IRzc2ACUgL2SR94igHHkh3IBpg0Deg/
H+dF5frAOxpw7KPWSxGPQ8EUHJ39mA/8Jbp6s2aYmiF2qrLTeKLEZr2x1mGf
ZF42kkEzgwUZLAOUYOmjqlY8z/FJpjMnA0qwJxgOAYN3P46Mg1DK3G+fTzJd
DGqwIAt+4DgNCacxffjBXLFyJfnQ7vFJpisnAxEediuc3aaaFZ2EgCPh0FJV
2OmTzJt0LUrUoDqhAoBQg06CEzgHfrYwNP3nTaifwHaftHrStRjV8S+/7ANZ
AqtDceoflo/DQxH7zfNOnVk+yb1D1+KUHE6DxTkP7/WdLZ97dizmhbSVPLKA
fI3T05iDrPdJti8nq0sMloC1FrA44jQxNZej+Rg+6/dTN9DpTKoyVxvJet49
ibJly5ayUqVKEsnjd9999+FZYfYTn+3fo2l2PqBrSWp3aBzQM8FrA0fhoQkr
6YwOVMTWUWkjluGw4ND69Xu0/L8kP+E8cJ7UhGt7SUGmFZtd4qNR/gldS9GW
RBgjcOIa1IgdU76VTzYeKWvXHErVjeYPxgTLEKpdH9gJG7camzKMEpAHWsEH
uudfqWGxjXw00v9F14u4leRSB8e5kJhcwcSADQHRttFialQbAtGL1prj9Fp9
FK1qsPpYYxQntGnTPkpLhXSbImlNqcU+89GKkEbX0iQVQUuambmDDpqDHRWm
JPQl7OWF0bORakN1Hx3Gh+yOoK4PqQwHG6l2FfBObA6Gi4XBS3Uc77DmeDpK
kR0SAqNcimaxZ320on1D17JUm9nZh6hxKXJow2MsTcS0DnN2WDzgCIvGhS6C
YuEISDBKHHRZ5ZHBdKpktapfyqZPjaIVqXRwdRydW51A9YNE9A+1QAZaOsU4
jo/nXrhgmxoR1vEB3qpKElAj8aiSd1xqw2IP/TzLjrchZTGKcPr0mdCuXUdC
+A74QReFDoSzV9D6MjK2eMcuJ9DGrZBg+EhwqOdUXjgMPFG++84MNGzxaqep
4qUXv6N9UxWX8Y7yjqcNKnBWqldL3iHgYD6oQX1CMw4Bj5OVHx4cxCNVznp1
h8t6ikshehDXJKSWRPUbLyCKoPfjyPHiuCYSKfECk/mUF76hu2Rci3l1L19S
rZUNoYEL+Hv87fV08UG/jIhnfd6dIVAe4p0uH8HiUIFoTXIXs1aMQvokbPUL
5eQcDyn+Ry1C5QjXFHwZ75Bzr/dg1czZJ8NP+3Gj4jfbwHPUWKWaWSIf8T5+
3Cox6LMFdJfsPYP6j3omDnxP8o59F13fSKc9KHB3AT979+3pqDvzWaj/x3PF
++/NOnf1ZFgLrfZO8c3I4pCDc5ORxXkG5yYji1MLrDIyh3MIFno4j1wRZHMU
gUt2WNRi7n4XmZ3daQNOxYOCD4to0tPTRb9+/WgrPBW0scznOAGrrKIJivki
UVht+yBQ5N+Iv8WP\
\>"]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1392, 969},
WindowMargins->{{-21, Automatic}, {Automatic, 0}},
ShowSelection->True,
Magnification->2,
FrontEndVersion->"6.0 for Linux x86 (32-bit) (April 20, 2007)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{
 "Info3405144475-5436130"->{
  Cell[1985, 70, 853, 16, 166, "Print",
   CellTags->"Info3405144475-5436130"]}
 }
*)
(*CellTagsIndex
CellTagsIndex->{
 {"Info3405144475-5436130", 33586, 735}
 }
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[590, 23, 102, 1, 143, "Section"],
Cell[CellGroupData[{
Cell[717, 28, 664, 17, 89, "Input"],
Cell[1384, 47, 446, 14, 89, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1867, 66, 115, 2, 89, "Input"],
Cell[1985, 70, 853, 16, 166, "Print",
 CellTags->"Info3405144475-5436130"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2875, 91, 289, 6, 89, "Input"],
Cell[3167, 99, 431, 13, 89, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3635, 117, 391, 10, 89, "Input"],
Cell[4029, 129, 1128, 40, 125, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[5206, 175, 108, 1, 143, "Section"],
Cell[CellGroupData[{
Cell[5339, 180, 559, 16, 89, "Input"],
Cell[5901, 198, 444, 14, 125, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6382, 217, 441, 10, 89, "Input"],
Cell[6826, 229, 833, 27, 174, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7696, 261, 707, 19, 168, "Input"],
Cell[8406, 282, 548, 17, 193, "Output"]
}, Open  ]],
Cell[8969, 302, 121, 3, 59, "Text"],
Cell[CellGroupData[{
Cell[9115, 309, 350, 11, 89, "Input"],
Cell[9468, 322, 23752, 397, 525, 11902, 201, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)