(* 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[    294577,       5885]
NotebookOptionsPosition[    290099,       5732]
NotebookOutlinePosition[    290572,       5751]
CellTagsIndexPosition[    290529,       5748]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Temperature Dependence of Resonance Capture", "Title",
 CellChangeTimes->{{3.414582638799663*^9, 3.4145826507525873`*^9}}],

Cell[CellGroupData[{

Cell["Packages to be Used", "Section",
 FontColor->RGBColor[0, 0, 1]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"General", "::", "spell1"}], "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.414582624075079*^9, 3.4145826262419367`*^9}},
 FontSize->9]
}, Closed]],

Cell[CellGroupData[{

Cell["Data", "Section"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Cons", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      SubscriptBox["k", "B"], "\[Rule]", 
      RowBox[{"1.3806505", " ", 
       RowBox[{"10", "^", 
        RowBox[{"-", "23"}]}]}]}], " ", ",", 
     RowBox[{
      SubscriptBox["T", "room"], "\[Rule]", "293.15"}], ",", 
     RowBox[{"e", "->", 
      RowBox[{"1.60217653", " ", 
       RowBox[{"10", "^", 
        RowBox[{"-", "19"}]}]}]}], ",", 
     RowBox[{
      SubscriptBox["m", "n"], "\[Rule]", 
      RowBox[{"1.674929", " ", 
       RowBox[{"10", "^", 
        RowBox[{"-", "27"}]}]}]}], ",", 
     RowBox[{"clight", "\[Rule]", 
      RowBox[{"2.99792458", " ", 
       RowBox[{"10", "^", "8"}]}]}], ",", 
     RowBox[{"mnMeV", "\[Rule]", "939.56536"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Eth", "[", "T_", "]"}], "=", 
    RowBox[{
     RowBox[{
      SubscriptBox["k", "B"], " ", 
      RowBox[{"T", "/", "e"}]}], "/.", "Cons"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Thermal", " ", "energy", " ", "in", " ", "eV"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"D2Odata", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".001105"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "->", ".03323"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "->", ".4519"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", 
      RowBox[{"4.42", " ", 
       RowBox[{"10", "^", 
        RowBox[{"-", "5"}]}]}]}], " ", ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "\[Rule]", "0"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"C12data", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"\[Rho]", "\[Rule]", ".00160"}], ",", 
      RowBox[{
       SubscriptBox["n", "d"], "\[Rule]", ".08023"}], ",", 
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "s"], "->", ".3811"}], ",", 
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".0002728"}], 
      " ", ",", 
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "f"], "\[Rule]", "0"}], ",", 
      RowBox[{"\[Nu]", "\[Rule]", "0"}]}], "}"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", "  ", 
   RowBox[{"Thermal", " ", "neutron", " ", "values"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Nadata", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".00097"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".02541"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".08131"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".01347"}], " ",
      ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "\[Rule]", "0"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"U235data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".01886"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".04833"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".01588"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", "4.833"}], ",", 
     
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "->", " ", "28.37"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "2.42"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"U238data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".0191"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".04833"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".4301"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".13194"}], ",",
      
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "->", " ", "0"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Pu239data", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"\[Rho]", "\[Rule]", ".0196"}], ",", 
      RowBox[{
       SubscriptBox["n", "d"], "\[Rule]", ".04938"}], ",", 
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".3902"}], ",", 
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", "13.27"}], ",",
       
      RowBox[{
       SubscriptBox["\[CapitalSigma]", "f"], "->", " ", "36.66"}], " ", ",", 
      RowBox[{"\[Nu]", "\[Rule]", "2.98"}]}], "}"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Fast", " ", "neutron", " ", "values"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Nadata", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".00097"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".02541"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".083853"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".000020328"}], 
     " ", ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "\[Rule]", "0"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "0"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"U235data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".01886"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".04833"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".328644"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".0120825"}], 
     ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "->", " ", ".06766"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "2.6"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"U238data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".0191"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".04833"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".33347"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".007732"}], 
     ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "->", " ", ".004591"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "2.6"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Pu239data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Rho]", "\[Rule]", ".0196"}], ",", 
     RowBox[{
      SubscriptBox["n", "d"], "\[Rule]", ".04938"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "s"], "\[Rule]", ".33578"}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "\[Gamma]"], "\[Rule]", ".0128388"}], 
     ",", 
     RowBox[{
      SubscriptBox["\[CapitalSigma]", "f"], "->", " ", ".091353"}], ",", 
     RowBox[{"\[Nu]", "\[Rule]", "2.98"}]}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.4156491584515057`*^9, 3.4156491596444607`*^9}, {
  3.415703768586635*^9, 3.415703772632246*^9}, {3.415703813594019*^9, 
  3.4157038229059677`*^9}, {3.4157038619879103`*^9, 3.415703908075962*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Temperature dependence", "Section"],

Cell[CellGroupData[{

Cell["\<\
Lorentzian resonance in cm frame\
\>", "Subsection",
 CellChangeTimes->{{3.4145882254254*^9, 3.4145882280469723`*^9}, {
  3.4153649569250727`*^9, 3.4153649657706823`*^9}, {3.415539922076974*^9, 
  3.415539930658679*^9}},
 FontSize->18],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Clear", "[", 
  RowBox[{
  "En", ",", "E0", ",", "\[CapitalGamma]", ",", "v", ",", "v0", ",", "vL", 
   ",", "vA", ",", "a", ",", "b", ",", "c", ",", "A", ",", "\[Alpha]"}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Lorentzian", "[", 
   RowBox[{"En_", ",", "E0_", ",", "\[CapitalGamma]_"}], "]"}], "=", 
  RowBox[{
   RowBox[{"\[Alpha]", "^", "2"}], " ", 
   RowBox[{
    RowBox[{"\[CapitalGamma]", "/", 
     RowBox[{"(", 
      RowBox[{"2", " ", "\[Pi]"}], " ", ")"}]}], "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"(", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"\[Alpha]", "^", "2"}], "En"}], "-", "E0"}], ")"}], "^", 
       "2"}], "+", 
      RowBox[{
       RowBox[{"\[CapitalGamma]", "^", "2"}], "/", "4"}]}], " ", ")"}]}], " ",
    
   RowBox[{"2", "/", 
    RowBox[{"(", 
     RowBox[{"1", "+", 
      RowBox[{"2", " ", 
       RowBox[{
        RowBox[{"ArcTan", "[", 
         RowBox[{"2", " ", 
          RowBox[{"E0", "/", "\[CapitalGamma]"}]}], "]"}], "/", "\[Pi]"}]}]}],
      ")"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{"Lorentzian", "[", 
    RowBox[{"En", ",", "E0", ",", "\[CapitalGamma]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"En", ",", "0", ",", 
     RowBox[{"+", "Infinity"}]}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"Im", "[", "E0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"Im", "[", "\[CapitalGamma]", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"Re", "[", "\[CapitalGamma]", "]"}], ">", "0"}], ",", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"\[Alpha]", "^", "2"}], "]"}], ">", "0"}], ",", 
      RowBox[{
       RowBox[{"Im", "[", 
        RowBox[{"\[Alpha]", "^", "2"}], "]"}], "\[Equal]", "0"}]}], "}"}]}]}],
   "]"}]}], "Input",
 CellChangeTimes->{{3.415537893981533*^9, 3.415537913747254*^9}, {
   3.415546312976169*^9, 3.415546347631803*^9}, {3.415546404284828*^9, 
   3.415546424306467*^9}, {3.415789675808961*^9, 3.4157897456999407`*^9}, {
   3.415789793147491*^9, 3.41578993226541*^9}, 3.4157930669430532`*^9, 
   3.415793100577594*^9, {3.415801143532527*^9, 3.415801180064081*^9}, {
   3.41580127790522*^9, 3.415801280296137*^9}, {3.415802234952695*^9, 
   3.415802258655449*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[Alpha]", "2"], " ", "\[CapitalGamma]"}], 
  RowBox[{"\[Pi]", " ", 
   RowBox[{"(", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "E0"}], "+", 
        RowBox[{"En", " ", 
         SuperscriptBox["\[Alpha]", "2"]}]}], ")"}], "2"], "+", 
     FractionBox[
      SuperscriptBox["\[CapitalGamma]", "2"], "4"]}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"1", "+", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"ArcTan", "[", 
        FractionBox[
         RowBox[{"2", " ", "E0"}], "\[CapitalGamma]"], "]"}]}], "\[Pi]"]}], 
    ")"}]}]]], "Output",
 CellChangeTimes->{
  3.41458265773838*^9, 3.4145889720791817`*^9, 3.415364970871722*^9, 
   3.415369447900679*^9, 3.415537915277678*^9, 3.415545478444785*^9, 
   3.415546324083254*^9, {3.4155463909448833`*^9, 3.415546425793028*^9}, 
   3.4156240797505207`*^9, 3.415637323508553*^9, 3.415648021362496*^9, 
   3.415699298882173*^9, 3.4157435968093853`*^9, 3.415789679132766*^9, {
   3.415789728831943*^9, 3.415789760436037*^9}, {3.415789841291545*^9, 
   3.415789887997251*^9}, 3.415789933772833*^9, 3.415801187062581*^9, 
   3.415802259465824*^9, 3.415813528076659*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{
  3.41458265773838*^9, 3.4145889720791817`*^9, 3.415364970871722*^9, 
   3.415369447900679*^9, 3.415537915277678*^9, 3.415545478444785*^9, 
   3.415546324083254*^9, {3.4155463909448833`*^9, 3.415546425793028*^9}, 
   3.4156240797505207`*^9, 3.415637323508553*^9, 3.415648021362496*^9, 
   3.415699298882173*^9, 3.4157435968093853`*^9, 3.415789679132766*^9, {
   3.415789728831943*^9, 3.415789760436037*^9}, {3.415789841291545*^9, 
   3.415789887997251*^9}, 3.415789933772833*^9, 3.415801187062581*^9, 
   3.415802259465824*^9, 3.4158135338008738`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
In the lab frame the target is moving with thermal velocity, vA, and the \
neutron is moving with velocity, vL.
The neutron center of mass velocity is vcm = \[Alpha] (vL - vA), \[Alpha] = \
A/(A+1).
Thus, Ecm = \[Alpha]^2 ( EL  + EA/A - 2 Sqrt[EL EA/A] cos(\[Theta]) ).  \
Integrating over angles gives:\
\>", "Subsection",
 CellChangeTimes->{{3.4155399599149513`*^9, 3.4155401830989237`*^9}, {
  3.4155506883807087`*^9, 3.415550690320909*^9}, {3.415552644252928*^9, 
  3.4155527131694937`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Clear", "[", 
  RowBox[{
  "PhiAvg", ",", "\[Alpha]", ",", "a", ",", "b", ",", "c", ",", "A", ",", 
   "EL", ",", "EA", ",", "E0", ",", "\[CapitalGamma]"}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"abc", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"a", "\[Rule]", 
      RowBox[{
       RowBox[{"\[Alpha]", "^", "2"}], 
       RowBox[{"(", " ", 
        RowBox[{"EL", " ", "+", 
         RowBox[{"EA", "/", "A"}]}], ")"}]}]}], ",", 
     RowBox[{"b", "\[Rule]", 
      RowBox[{"2", " ", 
       RowBox[{"\[Alpha]", "^", "2"}], " ", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"EL", " ", 
         RowBox[{"EA", "/", "A"}]}], "]"}]}]}], ",", 
     RowBox[{"c", "\[Rule]", 
      RowBox[{"\[CapitalGamma]", "/", "2"}]}]}], " ", "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"PhiAvg", "[", 
   RowBox[{"EL_", ",", "EA_", ",", "A_"}], "]"}], "=", 
  RowBox[{"Simplify", "[", "  ", 
   RowBox[{
    RowBox[{"2", 
     RowBox[{
      RowBox[{"\[Alpha]", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"2", " ", 
         RowBox[{
          RowBox[{"ArcTan", "[", 
           RowBox[{"2", " ", 
            RowBox[{"E0", "/", "\[CapitalGamma]"}]}], "]"}], "/", 
          "\[Pi]"}]}]}], ")"}]}], "\[Sigma]0", " ", 
     RowBox[{"\[CapitalGamma]", "/", 
      RowBox[{"(", 
       RowBox[{"2", " ", "\[Pi]"}], " ", ")"}]}], 
     RowBox[{"(", 
      RowBox[{"1", "/", "2"}], ")"}], 
     RowBox[{"Integrate", "[", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{"a", "-", 
             RowBox[{"b", " ", "\[Mu]"}], "-", "E0"}], ")"}], "^", "2"}], "+",
           
          RowBox[{"c", "^", "2"}]}], ")"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Mu]", ",", 
         RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", 
       RowBox[{"Assumptions", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"Im", "[", "a", "]"}], "\[Equal]", "0"}], ",", 
          RowBox[{
           RowBox[{"Im", "[", "E0", "]"}], "\[Equal]", "0"}], ",", 
          RowBox[{
           RowBox[{"Im", "[", "b", "]"}], "\[Equal]", "0"}], ",", 
          RowBox[{
           RowBox[{"Im", "[", "c", "]"}], "\[Equal]", "0"}], ",", 
          RowBox[{
           RowBox[{"Re", "[", "c", "]"}], ">", "0"}], ",", 
          RowBox[{
           RowBox[{"Re", "[", "b", "]"}], ">", "0"}]}], "}"}]}]}], "]"}]}], "/.",
     "abc"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.4155379639772873`*^9, 3.4155379647268953`*^9}, {
   3.4155388656637487`*^9, 3.41553891203726*^9}, {3.4155389424835663`*^9, 
   3.415539031211755*^9}, {3.415539089468835*^9, 3.415539214412753*^9}, {
   3.415539285961134*^9, 3.4155392978339577`*^9}, {3.415539865239097*^9, 
   3.415539867315319*^9}, {3.415540189580212*^9, 3.415540195859633*^9}, {
   3.4155454987732677`*^9, 3.415545558733273*^9}, {3.415545599928729*^9, 
   3.415545735106172*^9}, {3.415550710120329*^9, 3.415550730519535*^9}, {
   3.41562710770224*^9, 3.415627111693184*^9}, {3.4156277208528767`*^9, 
   3.415627723888749*^9}, {3.41562805669567*^9, 3.415628058591188*^9}, 
   3.415789982626305*^9, {3.4157911885773478`*^9, 3.415791192422614*^9}, {
   3.41579123872152*^9, 3.415791269297167*^9}, {3.415791884416909*^9, 
   3.415791986736073*^9}, {3.415792049691592*^9, 3.415792081841432*^9}, {
   3.415793121112297*^9, 3.4157931296125507`*^9}, {3.415793165751495*^9, 
   3.41579316617633*^9}, {3.415801206197096*^9, 3.415801257002728*^9}, {
   3.4158021998476562`*^9, 3.4158022271233187`*^9}, {3.415802344530531*^9, 
   3.415802360856153*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"\[Sigma]0", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"ArcTan", "[", 
      FractionBox[
       RowBox[{"2", " ", 
        RowBox[{"(", 
         RowBox[{"E0", "-", 
          FractionBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"EA", "+", 
              RowBox[{"A", " ", "EL"}], "-", 
              RowBox[{"2", " ", "A", " ", 
               SqrtBox[
                FractionBox[
                 RowBox[{"EA", " ", "EL"}], "A"]]}]}], ")"}], " ", 
            SuperscriptBox["\[Alpha]", "2"]}], "A"]}], ")"}]}], 
       "\[CapitalGamma]"], "]"}], "+", 
     RowBox[{"ArcTan", "[", 
      FractionBox[
       RowBox[{"2", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "E0"}], "+", 
          FractionBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"EA", "+", 
              RowBox[{"A", " ", "EL"}], "+", 
              RowBox[{"2", " ", "A", " ", 
               SqrtBox[
                FractionBox[
                 RowBox[{"EA", " ", "EL"}], "A"]]}]}], ")"}], " ", 
            SuperscriptBox["\[Alpha]", "2"]}], "A"]}], ")"}]}], 
       "\[CapitalGamma]"], "]"}]}], ")"}]}], 
  RowBox[{"2", " ", 
   SqrtBox[
    FractionBox[
     RowBox[{"EA", " ", "EL"}], "A"]], " ", 
   RowBox[{"(", 
    RowBox[{"\[Pi]", "+", 
     RowBox[{"2", " ", 
      RowBox[{"ArcTan", "[", 
       FractionBox[
        RowBox[{"2", " ", "E0"}], "\[CapitalGamma]"], "]"}]}]}], 
    ")"}]}]]], "Output",
 CellChangeTimes->{
  3.414582660958725*^9, 3.414588973372835*^9, 3.415369450470133*^9, 
   3.4155392561954823`*^9, 3.4155392995528727`*^9, 3.415539870402616*^9, 
   3.415540203339308*^9, 3.4155457445619*^9, 3.415550733563284*^9, 
   3.415624093349469*^9, 3.415627116664648*^9, 3.415627758453504*^9, 
   3.415628064232849*^9, 3.41563733431805*^9, 3.4156480321369762`*^9, 
   3.415699310070067*^9, {3.41574360634752*^9, 3.415743610777001*^9}, 
   3.415789985117756*^9, 3.41579213139121*^9, 3.415793187785729*^9, 
   3.415802392037322*^9, 3.415813544593515*^9, 3.415817705740425*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Multiply by the Maxwell-Boltzmann distribution to get the thermal average \
integrand\
\>", "Subsubsection",
 CellChangeTimes->{{3.414588571845256*^9, 3.414588621246562*^9}, {
  3.415540216340914*^9, 3.415540227890996*^9}},
 FontSize->18],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"fMaxwell", "[", 
    RowBox[{"EA_", ",", "kT_"}], "]"}], "=", " ", 
   RowBox[{
    RowBox[{"2", "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Pi]", "^", 
        RowBox[{"(", 
         RowBox[{"1", "/", "2"}], ")"}]}], " ", 
       RowBox[{
        RowBox[{"(", "kT", ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"3", "/", "2"}], ")"}]}]}], ")"}]}], 
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", "EA"}], "/", "kT"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"fMaxwell", "[", 
     RowBox[{"EA", ",", "kT"}], "]"}], " ", 
    RowBox[{"Sqrt", "[", "EA", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"EA", ",", "0", ",", "Infinity"}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Re", "[", "kT", "]"}], ">", "0"}], "}"}]}]}], "]"}]}], "Input",\

 CellChangeTimes->{{3.415539431074314*^9, 3.415539758909307*^9}, {
   3.415539796118031*^9, 3.415539800797398*^9}, {3.415545755224636*^9, 
   3.415545845846096*^9}, {3.415545876648521*^9, 3.415545976989833*^9}, {
   3.415548398541892*^9, 3.415548424259398*^9}, 3.415552278320196*^9, {
   3.415626417602976*^9, 3.415626472849869*^9}, {3.415627820423937*^9, 
   3.415627843445548*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{{3.4155395972316008`*^9, 3.415539617172885*^9}, {
   3.4155396610623827`*^9, 3.41553970543102*^9}, {3.415539747803831*^9, 
   3.415539762052507*^9}, 3.415539801747921*^9, 3.41554023461523*^9, {
   3.415545830593396*^9, 3.415545859422305*^9}, {3.415545896030094*^9, 
   3.4155459191167917`*^9}, {3.415545957611804*^9, 3.4155459803559713`*^9}, 
   3.415548424986824*^9, 3.415624103823975*^9, {3.415626420476886*^9, 
   3.415626477015526*^9}, 3.415627120108316*^9, {3.41562783666926*^9, 
   3.415627844235138*^9}, 3.4156278880123167`*^9, 3.415628070784724*^9, 
   3.415637350642392*^9, 3.415648047152302*^9, 3.415699314922051*^9, 
   3.4157436156314983`*^9, 3.415792172990871*^9, 3.415793195111651*^9, 
   3.415802412021358*^9, 3.41581355872585*^9, 3.415817710522094*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"integrand", "[", 
   RowBox[{"EL_", ",", "EA_", ",", "A_", ",", "kT_"}], "]"}], "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{"Sqrt", "[", "EA", "]"}], " ", 
    RowBox[{"fMaxwell", "[", 
     RowBox[{"EA", ",", "kT"}], "]"}], " ", 
    RowBox[{"PhiAvg", "[", 
     RowBox[{"EL", ",", "EA", ",", "A"}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4155393884489594`*^9, 3.415539417694529*^9}, {
   3.415539785562237*^9, 3.415539819074635*^9}, {3.4155398747344313`*^9, 
   3.41553989078057*^9}, {3.415540454317543*^9, 3.415540455732788*^9}, {
   3.4155459847909393`*^9, 3.415545997852717*^9}, 3.415624117666798*^9, {
   3.4156277441087103`*^9, 3.415627768147955*^9}, {3.415628043436928*^9, 
   3.415628044498403*^9}, {3.415628081708909*^9, 3.415628082426581*^9}, {
   3.4156282082705173`*^9, 3.4156282086759863`*^9}, {3.415792186748392*^9, 
   3.4157921902496634`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["\[ExponentialE]", 
    RowBox[{"-", 
     FractionBox["EA", "kT"]}]], " ", 
   SqrtBox["EA"], " ", "\[Sigma]0", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"ArcTan", "[", 
      FractionBox[
       RowBox[{"2", " ", 
        RowBox[{"(", 
         RowBox[{"E0", "-", 
          FractionBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"EA", "+", 
              RowBox[{"A", " ", "EL"}], "-", 
              RowBox[{"2", " ", "A", " ", 
               SqrtBox[
                FractionBox[
                 RowBox[{"EA", " ", "EL"}], "A"]]}]}], ")"}], " ", 
            SuperscriptBox["\[Alpha]", "2"]}], "A"]}], ")"}]}], 
       "\[CapitalGamma]"], "]"}], "+", 
     RowBox[{"ArcTan", "[", 
      FractionBox[
       RowBox[{"2", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "E0"}], "+", 
          FractionBox[
           RowBox[{
            RowBox[{"(", 
             RowBox[{"EA", "+", 
              RowBox[{"A", " ", "EL"}], "+", 
              RowBox[{"2", " ", "A", " ", 
               SqrtBox[
                FractionBox[
                 RowBox[{"EA", " ", "EL"}], "A"]]}]}], ")"}], " ", 
            SuperscriptBox["\[Alpha]", "2"]}], "A"]}], ")"}]}], 
       "\[CapitalGamma]"], "]"}]}], ")"}]}], 
  RowBox[{
   SqrtBox[
    FractionBox[
     RowBox[{"EA", " ", "EL"}], "A"]], " ", 
   SuperscriptBox["kT", 
    RowBox[{"3", "/", "2"}]], " ", 
   SqrtBox["\[Pi]"], " ", 
   RowBox[{"(", 
    RowBox[{"\[Pi]", "+", 
     RowBox[{"2", " ", 
      RowBox[{"ArcTan", "[", 
       FractionBox[
        RowBox[{"2", " ", "E0"}], "\[CapitalGamma]"], "]"}]}]}], 
    ")"}]}]]], "Output",
 CellChangeTimes->{
  3.414582661083025*^9, 3.4145889734103823`*^9, 3.4153694506982203`*^9, 
   3.415539891943803*^9, 3.415540237779212*^9, 3.415540459202465*^9, 
   3.4155459990312223`*^9, 3.4156241190248413`*^9, 3.415626506723218*^9, 
   3.415627122730712*^9, 3.415627771312151*^9, 3.415627890954702*^9, {
   3.415628073889285*^9, 3.4156280841397667`*^9}, 3.415628210674472*^9, 
   3.415637356002845*^9, 3.415648051095458*^9, 3.415699318870551*^9, 
   3.4157436184433727`*^9, {3.415792176920766*^9, 3.4157921920672483`*^9}, 
   3.415793199013254*^9, 3.415802416227642*^9, 3.4158135588620167`*^9, 
   3.4158177146400557`*^9}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Target and resonance data", "Subsection",
 CellChangeTimes->{{3.415649106810667*^9, 3.415649114515381*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Avalue", "=", "50."}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Target", " ", "mass"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"E0value", "=", "5."}], ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{"location", " ", "of", " ", "resonance", " ", "in", " ", "eV"}], 
   " ", "*)"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalGamma]value", "=", ".1"}], ";", " ", 
  RowBox[{"(*", " ", 
   RowBox[{"width", " ", "of", " ", "resonance", " ", "in", " ", "eV"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"sub", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Sigma]0", "\[Rule]", "1"}], ",", 
     RowBox[{"A", "\[Rule]", "Avalue"}], ",", 
     RowBox[{"E0", "\[Rule]", "E0value"}], ",", 
     RowBox[{"\[CapitalGamma]", "\[Rule]", "\[CapitalGamma]value"}], ",", 
     RowBox[{"\[Alpha]", "\[Rule]", 
      RowBox[{"Avalue", "/", 
       RowBox[{"(", 
        RowBox[{"Avalue", "+", "1"}], ")"}]}]}]}], 
    "}"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Sigma0", "[", "EL_", "]"}], "=", " ", 
   RowBox[{
    RowBox[{"\[Sigma]0", " ", 
     RowBox[{"Lorentzian", "[", 
      RowBox[{"EL", ",", "E0", ",", "\[CapitalGamma]"}], "]"}]}], " ", "/.", 
    "sub"}]}], ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.41564931829716*^9, 3.415649357225047*^9}, {
   3.415651205673894*^9, 3.415651210726491*^9}, {3.415651502841469*^9, 
   3.415651506342553*^9}, {3.415651585417511*^9, 3.415651586910664*^9}, {
   3.4156516700014133`*^9, 3.415651705166987*^9}, 3.415699334889825*^9, 
   3.415699368451103*^9, {3.4157469778064547`*^9, 3.415746987101452*^9}, {
   3.415747190975032*^9, 3.415747191819005*^9}, {3.4157542153649807`*^9, 
   3.415754216194182*^9}, {3.4157555073500834`*^9, 3.415755510875423*^9}, {
   3.415792201102124*^9, 3.415792202347381*^9}, 3.415793388018008*^9, {
   3.415814023922285*^9, 3.415814025239757*^9}, {3.41581421896384*^9, 
   3.415814239440995*^9}, {3.415814532885626*^9, 3.415814532970818*^9}, {
   3.415817398082358*^9, 3.4158174064156723`*^9}, {3.415817451850189*^9, 
   3.415817486023384*^9}, {3.415817719466284*^9, 3.415817730071287*^9}, {
   3.4158183662197742`*^9, 3.415818368321199*^9}, 3.415818405818779*^9, {
   3.4158271643270683`*^9, 3.415827171427669*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"\[Sigma]0", "\[Rule]", "1"}], ",", 
   RowBox[{"A", "\[Rule]", "50.`"}], ",", 
   RowBox[{"E0", "\[Rule]", "5.`"}], ",", 
   RowBox[{"\[CapitalGamma]", "\[Rule]", "0.1`"}], ",", 
   RowBox[{"\[Alpha]", "\[Rule]", "0.9803921568627451`"}]}], "}"}]], "Output",\

 CellChangeTimes->{
  3.415649117853862*^9, {3.415649342903859*^9, 3.415649358557542*^9}, 
   3.415651214425674*^9, 3.4156515088170767`*^9, 3.415651588508519*^9, 
   3.4156516745201883`*^9, 3.415651707120779*^9, 3.4156993414280777`*^9, 
   3.415699375433522*^9, 3.415743623199136*^9, 3.4157470009535027`*^9, 
   3.415747193782446*^9, 3.415754222146678*^9, 3.41575551357796*^9, 
   3.4157922051080837`*^9, 3.415793203160738*^9, 3.4157933889060698`*^9, 
   3.415802422538946*^9, 3.415813565092452*^9, 3.415814032964281*^9, 
   3.4158142669644003`*^9, {3.415814534794771*^9, 3.41581454693646*^9}, 
   3.4158174082589083`*^9, {3.415817458591832*^9, 3.415817487865448*^9}, 
   3.415817614096924*^9, 3.415817730963043*^9, 3.415818371131521*^9, 
   3.415818407025638*^9, 3.4158271810620823`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Lorentzian folded with the Maxwell-Boltzmann distribution", "Subsection",
 FontSize->18],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Tvalue", "=", "300"}], ";", 
  RowBox[{"kTvalue", "=", 
   RowBox[{"Eth", "[", "Tvalue", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ELmin", "=", 
    RowBox[{"Abs", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"E0value", "-", 
         RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
       RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}], "]"}]}], ";", 
   RowBox[{"ELmax", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E0value", "+", 
        RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
      RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}]}], " ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Sigma]Th", "[", 
    RowBox[{"EL_", ",", "T_"}], "]"}], ":=", 
   RowBox[{"Abs", "[", " ", 
    RowBox[{"NIntegrate", "[", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"integrand", "[", 
        RowBox[{"EL", ",", "EA", ",", "Avalue", ",", 
         RowBox[{"Eth", "[", "T", "]"}]}], "]"}], "/.", "sub"}], ",", 
      RowBox[{"{", 
       RowBox[{"EA", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Sigma]Th", "[", 
      RowBox[{"EL", ",", "Tvalue"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"EL", " ", ",", "ELmin", ",", "ELmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"1", ",", "0", ",", "0"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Sigma0", "[", "EL", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"EL", ",", "ELmin", ",", "ELmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"0", ",", "0", ",", "1"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"p1", ",", "p2"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sum", "=", "0."}], ";", 
  RowBox[{"Npts", "=", "400"}], ";", 
  RowBox[{"dEL", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"ELmax", "-", "ELmin"}], ")"}], "/", "Npts"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Do", "[", 
   RowBox[{
    RowBox[{"sum", "=", 
     RowBox[{"sum", "+", 
      RowBox[{"\[Sigma]Th", "[", 
       RowBox[{
        RowBox[{"ELmin", "+", 
         RowBox[{"dEL", " ", "i"}]}], ",", "Tvalue"}], "]"}]}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "Npts"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\< Normalization of thermalized cross section = \>\"", ",", 
   RowBox[{"sum", " ", "dEL"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\< Normalization of bare cross section = \>\"", ",", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"Sigma0", "[", "EL", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"EL", ",", " ", "ELmin", ",", "ELmax"}], "}"}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.41458304075012*^9, 3.414583046995742*^9}, {
   3.415540269451309*^9, 3.415540624050168*^9}, {3.415540679931073*^9, 
   3.415540730786166*^9}, {3.4155462669012136`*^9, 3.415546291121092*^9}, {
   3.415546469183517*^9, 3.415546616453907*^9}, {3.4155468744865217`*^9, 
   3.41554690663778*^9}, {3.415546946880457*^9, 3.4155470219102983`*^9}, {
   3.415548220041874*^9, 3.415548232414732*^9}, {3.415551830444993*^9, 
   3.415551834786635*^9}, {3.415552314298234*^9, 3.41555233042307*^9}, {
   3.415552392765728*^9, 3.415552408651167*^9}, {3.415552483527116*^9, 
   3.415552483620476*^9}, {3.415624242545573*^9, 3.415624249253792*^9}, {
   3.415624324283189*^9, 3.415624325695047*^9}, {3.415624382110117*^9, 
   3.415624486747349*^9}, {3.415624557490179*^9, 3.415624590167341*^9}, {
   3.4156270416295643`*^9, 3.415627062306756*^9}, {3.415627305404591*^9, 
   3.415627394897621*^9}, {3.415627441697896*^9, 3.415627444697424*^9}, {
   3.415627507041526*^9, 3.415627507774624*^9}, {3.4156276556644707`*^9, 
   3.415627683164579*^9}, {3.415627875848173*^9, 3.41562788281986*^9}, {
   3.415627918928083*^9, 3.415627923845487*^9}, {3.415627965624673*^9, 
   3.415627980198*^9}, {3.4156280935263557`*^9, 3.415628095307498*^9}, {
   3.415628133258895*^9, 3.415628140124235*^9}, 3.4156282035532923`*^9, {
   3.415637371947699*^9, 3.415637413656823*^9}, {3.415637454420929*^9, 
   3.415637455674995*^9}, {3.415637485691234*^9, 3.4156375148166723`*^9}, {
   3.415637746090425*^9, 3.415637771655744*^9}, {3.415637816626465*^9, 
   3.41563782461031*^9}, {3.415637903506456*^9, 3.415637959321288*^9}, {
   3.41563803356126*^9, 3.4156380399095783`*^9}, {3.415648150437632*^9, 
   3.415648150738111*^9}, {3.415648537382182*^9, 3.415648602635006*^9}, {
   3.415648748836472*^9, 3.415648778041828*^9}, {3.415649067218238*^9, 
   3.41564909572712*^9}, {3.4156491495220346`*^9, 3.415649171094987*^9}, {
   3.4156492528799677`*^9, 3.415649316152924*^9}, {3.415649362624539*^9, 
   3.415649372933012*^9}, {3.4156512376152487`*^9, 3.4156512709008512`*^9}, {
   3.415651644679956*^9, 3.415651647645138*^9}, {3.415651683623374*^9, 
   3.415651685636747*^9}, {3.4156517192954903`*^9, 3.41565172007684*^9}, {
   3.415699443967512*^9, 3.41569949334918*^9}, {3.415792245060175*^9, 
   3.4157922460490103`*^9}, {3.41579326955651*^9, 3.415793344944862*^9}, {
   3.415814625308722*^9, 3.4158146285053*^9}, {3.415814672483979*^9, 
   3.4158146848174477`*^9}, {3.415817427552219*^9, 3.415817434240233*^9}, {
   3.415817514120131*^9, 3.4158175188378*^9}, {3.415817556688573*^9, 
   3.415817578190528*^9}, {3.4158176220325747`*^9, 3.415817641181723*^9}, {
   3.415818382017745*^9, 3.415818385542156*^9}, {3.4158184236015463`*^9, 
   3.415818451871264*^9}, {3.415827193604924*^9, 3.4158271936900587`*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwV13dczP8fAHDSvFB3n8/dJRkZKTIaiPB6kag0kYaWCqGSlZAoaVGnIS0a
2lNDm9KWkr7tpZ2WCpUy8vv8/rl7PB/3fty95vvuxM0vHz/LtWjRIi7q4f/P
NanKd27HCOP25thMs0oOnBfwrhCPFcbQ2ZdpjvUc4LJspL+nzJ2S4Pq8gwN7
VljEs+OF8dbBoxVfJzkQ++B+U1aiMK4Jqb9nufQJ3Ncv2D6ZJoz6zMWe2Tuf
gDyXzIBFvjBaHp1UPOj1BEJ1xTQ06oTRTO94mraiL3xJd8yN/iOMp8YLduT3
+4G9ZNTkLwk6CnZbZ2xNDAC+rw6PHbXp+PXQ1Tz1hEAgz7eVuNylo9c5cOmV
DgbL3wtl7nF0vB6vWkKMh8KjXJux/fF0NBVr8P45FwoZNzrp05T5F1QMB3jC
YPFUnolZIh1ZFzrYbWvC4MXA9TmFFDq+3zpP7jgZBi21o1tGM+iY4lWQGfo2
DNQimv3U39LxQkRNYHzIc5BRTjWjN9ExsDxGMOdKODypzJQvo2y2utG98344
TKnm8d9spqOWgCJ92ZNwSNUse9XZQseCfztuclLDYbN++6L4dup8gFSp6Hg4
bLjEG3mgh46jvxb7JcpHgIivSe+lMTpqyrSu0c2OgH+dy80rFjPQW05sw/ua
SKi6piRjw8XAwmMiD7WaIuEJzWERuYSBS3OPszu7ImGtQu+LM9wMzHxa+4Vv
MhIO+Wd2/OFlYNYsZyCRHgWuqvq6cksZ+DdAtmy9fhTwZ0eqRLAYuGn5qP3y
0Sigc+S3397CwLzrPlWa66PhlXqd8HZpBhpa9HRc3xYNmrSL3/spB8a7akXv
iQavh+GvNbYxUKjgsssa7WjgcqTtE5dh4IqdvYc97kbDtFWPStUuBoZouvZI
t0dD66HH5qxDDBzff3GWERoD4bP9TzP0GSh2zKWMfiAOYqtSrp42YGDvuKkj
17E4SA25qcVtyMCyjTFbFvTioPCAoMCp0wz8jjVRrKtx0O4m6zhvzEAgzfRr
Y+OAZDubowUDTfg3Khsz4sFj9+rtH20YKH/66clz3+PhioN+1YgLA30Oa3oY
tyXCK6G845wHDAwWfdC7eSQRpmJXdMm7MvCacN4l3vlEsG1q/3bvIQNvu9rr
D4skgbWs8UqWBwPVRI84GRgkwbnxM7aHvBl473xY6qKuJDA0u0SGPmPg9ugF
maHJZAj++eEFBjHwZiVx6vniFGj1kZYaotzkOrn0LJECem++7pcJYaBl3WV3
id0pcFLU7nxlGAOTgnV4r99LAa3G6/k/IhkYttpLV51IhcMqTmbqyQzsYCdI
+KukgQGv8t9/lNUX21srmqSBbZlgSGYKA7Us1WJ/XUuDYAxuEE1j4Mrplttv
w9NgQiHryEg6Az+uST5lOpcGQVKj0m45VP8PWNMZeq9gnHZqvqiEgcq9n4WN
Wemw+L1Y4LVSBtquzeIvXZcOLPd+2U1lDExRrK7dtT0dDi65Yu1TzkDFzDpC
RSUdAn8/6jaqYuC/y0MDxJ10wPF35fO1DGSdFm206E8H/9qtfnJtDAz9cXzK
vSgDGBHrN3dQ1nTIEP9ckwFPrq4ocWlnYNQy/dFD7RnwmM3zvb6Dqoc4PVhh
JgPcznTq2H1mYIBkdXvzlkxwmPGip/QzsO/l0x/fQzPBaNXwE4kJBhIfvTac
8MiC9skuyY+UPdLjuBOfZYF+SUPxjUkG2vjc8mXEZYHu+aKpsikGJqqeEV9R
kQWa6YHaFj8YyFU3uPsf92s4qKwsHDHHQJHMB1Mqbq9BwiaKs2IJgdb8bqIp
z7Lh9rmFPXWU15JGfj8Ss+GjqeGAKzeBkW53nU68zYabxxl7pngILBH0sjEb
zIaq3S59FfwEOibrNmbK58BFLvOd15cTOHJsaeKl1hxIfSbe8XEFgWqrzec4
O/OAy/fuQ1dRAgt+S7SeV8+DU15t2/euJHC8W/CLqUUe/HP0dY0RIzC49pBP
um8eaJsv3ua4hsAlg2K0vsk8+C7de19qI4Hnt6xalZGRDztLIiRcd1Dx7pTZ
wneyEEru+zvtkiFQcVGayZhtIWiCW/Mw5SCd5ts/PQvhXOElNw05Ao1fieZ6
FBfCs+xdQ6xdVHzcnq5B29/AfGJNTIIigazHznr5zLfwxn9uQ90RAtdvnDXS
4SoGteNjjs5HCawYcVAIZRRDs/DnRjkVqj4HxDf9W1cME96lrkGqBGp3bpRk
HS6GNR4+A2fUCbyoYX3E3b0YnB03RE/rEJhXa5WcTLwD5bPa60WNCYw3defq
O1QCexWu2hdTPjD9IdNPtwR2CAa8P2dC4Lnd1q/OWJWAaHrL5UxTAiVW64ef
9imBqT8mhermVD75QgKVHSUQEmCje+88le9yuXOFjqUwWerlPnCFwLYdiR+O
t5XBQGByu9dVAhXkzzwznyqD9gsft8pcI/Cv6p6/IXzlUCbEaLx7ncD3+GSj
+65yCDkdvJZ9k0CbGhDSDyyHwz/i8lQdCWRvmvfcY1IBwevKx1LdCDTZci90
C6MKLMU0dCzcCSw6JGnwQqYKdrCastkeBOrNNnof16mCCoFB5/ueVL1S2w67
+FbBjylukeOPCeSduFRHMN+DRtHhwzO+BLamVI6mSlWDSF5NQoIfgRpyPiFP
NaphIOOkkIk/gQE7t3Gqr1TD7VjL9ooAAr+eoXesK6iGOG9Xu6BnBDrI0DJi
dT7AYqOyMMXnBMomaO5v8KmBGl11rinKexrlfgvl1kCQVuP56BcE2ln8l5/Y
WwPblAZkl0UQ6Kzk/HJhZy2c3sz9/nMUgddcoqKgvxZezynNOMcTuFottW2b
Th1ceFqqUZVJ4G4jz6Sxg/Xw4uXtVLssAodkJ0zKTtfDf+kyQiteE1i3641/
64162Psx/NOFbALddwpvyE6sh6V8jicE8whcFz2Uo8/+D1Id5A0031KvL6iG
tv75D/rcxvJmKf8FLq1C0QZgPY0SDS+i8lnyd3hIoQHupQt3ThYTOPpA6ith
3wDao+MmvqUEMkd8bvT/aIAfp2MsG6sInHXer0BfaATJi0blju8JzFVu/lO2
tgmMHAiJjdUEqppcEqxVaoKyAOcvNz4QeJP3aXyXVxME1hpfZH8kcM3KzVWi
q5pBAVh2hg1UfFwFJe3HW+DuOvc7PV0Evvi6zdiP1g4eetETgZ8J3NhrIP5q
fzv4PX53RqObwIFNs416V9ohdvb30fweAmOuLl5Qa2+Humo7MqCfwJ5Blt9o
RgeIXzNMOTpCYOOA/Cer+10gHX9z7QLl5O822l2FXbCrK8A/a5Ta9xPqKdPz
XaB2tO6W+DiBcS3K7gz7z3Bt5eEjvycI9HhaNqlg3w3lpVu7U6cJtL1xdLTh
Qw98mlM7fnaGQEGrzr8ekz3QvtWqfOUsgWYZyjNhRC9MPItK8vhJYFT4nSWD
Rr3AtmY7mP8i8FAMu2vZdC9cIBfTWYtIHPUq3fFwZz8sO9t46K4giY+Ht1SL
CA6BtqmHxY6lJIrMVB2zUxwCP4N9rv2U7/HEhe+3HgKWZnS56nISh3fxVu6r
G4K1u6+rsOgkrhOznu8N+QKyAqRmGotE+ycWghdUR+D6kipbczaJjWnDtBvO
I5D99w6HKULinymzDVz5I6D4vf/TnRUk5j5T2RglPQqHOzJPqIiRGDEWokgT
GYNTKScMesVJXJXMFPjJ9xWC4vhuB6wjsadl0zoDla/QHlkQcnQ9if8ssPWw
51cwDVzfmbKBxKreSzn5SyfA6t4Pk9ubSKyUIh9tFp2E2zr+lsRWEpfnuD5J
jZuCFzP/XVZWIFFuvemtJYwfoMqd/Wic8h+Rjy6g8AOmGcFx/ntIDMis9c81
+QHHdph29+4lsaUgF88k/4C5i2OaTvupemUIn/BQm4YTPUu25RwisVAzgGuR
3wzwf5Afl9Qg8aCrZvnDY3NwLSLQaqk5ibq7PIV2ayzAdFN34DvKgVuJi20P
FsBeUKrc3oLE2y1SoSH5C+BgXyDea0lieq7vueuS/8BJvaf99XkS73Au277i
WoQec1IapjYkHv803netdxGGab+RyXAg0bPHvmdZJBeucuM1O3+LxOiGLy6D
hVz4okDLR+w2ifI1fouyWrkwUqJv1O0OiX0X5v/JCC/B2D+80aedSDTad6w/
yGkJvorXZvG4klgbKP/lsD43li0a+KXHIZHzJ4R5d44HU9WMNix+QqIdyTz9
VogXgwIaNRIpd+wwdOiS4EUbyfKIP74kktOxKzNO8iJLO/ZIZACJxd/tlMJT
efFCuJXfaDCJm5YU5HUa8uHJkZ4C/xASZTc/6H9xmQ8PyBkM7gslsaiyNUHW
lQ8ZlaoKnDAS39kfWlGTzIeFE5s/y4WTGPZHf1XGbz5cfuCr1N1oEt+uisxO
4vDjnPvZExIxVP6rdhXJRPJjf32XYx3lrY5S250y+DH37Mc68TgSH842pZxv
5Mcz3mk3KhJIfHZN4+IlpgCqt2wKv5xI4oGu5kHGRgHcJR5RJZJEIp3U/3BH
XgBpr5+svJRMYldv5h3OcQHM6rzyTiiNxCU8LfdSHwngC4nR0VzKg4/y1HSC
BdDDzpw0f0XNx1FF/vRYATTmPnE+K53Eur5nx2KKBZBXWn6ZQRaJzLEc+olv
Ajh1I3kX12sS13745VT3RwA7ijaYJVG+vGj3Yh5+Gr46wcz8m02i0NAh9+BV
NDS4M6sflUfN5xNDHmNlGjp68exQySdxWsS92kqThi+CSb4Jys+lXHW36tGw
L0f29Z5C6vPU0m26z9OQu/Lg427K6ov4DGsu03BTs7bFwzdU/uceVNvdpKH1
tC29/i01X0k3ye6HNPRZ4jRsX0T1t3WLRfxjKh6Gd5FYMYnOsR5rJf1pOLMj
ycbqHbV/6aVTquE0FMH8w8tLSLx1xIszFk3DvVrvV2ZRlvhen3UokYZONl/e
/yslsSnxerVwFg0jHGcjYspI3Hm2x9g9l4alj3gcjpWTKJY3/SajkIZDIaTW
N8qE6JNSv2Ia8ieul3hWQeJfwYdqm8pouDlP9u++ShIvOIXJ2FXSUL3qYGMf
5YZvzwxtqmlo26Kd5FFF4tPzcomra2noO2Tqsu09tf8rN/G71tEwa8bWoJHy
Q5qU6fN6GjZzO+24XU3iguJE2MUGGs4R3nxrP5B4+POWjLFGGoquD/tcTllY
IshvbTMN98kmvb5UQ6L0r1F5gRYamh7Mf0yvpfbP6rNfDGUX7fcWOZSl+rcl
zFCONm3da/yR2u9joTfnKVfYfqEvqaP6XzywKJ3yyN3Z4XjKvns/qolRXurN
U6z5icT27ytNFKnP2xZGPpumPOjjd2R5Ew21k9bbhtST6CM6y3hKxXs1X1YZ
/6Pu62r+pjoqn4D3B8WGKBu4e0a+pfLNadX+8aiBxK9HZZ+coerR/sW0WqaR
un+3ZxfkUvX6M2sb2UL5962ig2VUPVfzOjncbaL6uevrIVeq3sj01lrfTKJ5
+ae+Oaof5hvCJN5TLpIcO7ruDQ0fyiX9tW0h8bNYYuiiPBpW67xPym8lccUd
g9LONBqOm7W6mLWRaErEz7VT87Dc7osBbzuJ5y20u31iaHjCh4f/eAeJh55O
RYmE0ND+Odn9k7Lt1jVzE9S8BSWvz37eSeKPDUsn73jTsLP6oOVIF4l82w3q
8+/RcKFNW5HzmcTxA/uO2TjQUHzElLGzm0Ql87YHLXY0PMfnVHy/h/o+idhp
3WhGw0mlfDGRfmpfTQ4cu7WfhsO3tvPWUH5q2lC4Vo6GPWnRk/cGSKz2Tra0
l6Thp5Wcki+DJDZ3KOpuY9Aw44eFVc4wiYoyFy529wpgolSbzqURygVLJfia
BDDKVFNxzSiJyvuFLvVWCqD/B4Xl7mPU/WHMNZeWLIDXXy7LOjVB4m6evsCK
qwJo3ebynDZJYnen06IACwG0FJpze0tZq1xz+YqTAqh7p89A4huJwady29Sp
+2nXiZx/Mz+oeZVcKnL5G3Ufcp059nSe2l8NUWtlU36cUmiWV/1FIgutmYPH
+HHY9tjqv5SFEz65qijwY1v7zm+Wf0h0Uh7nNRHmx/wM2jP5fyRuu3Roa8cb
PnQ0z+pr4GYiz66eG93L+XBxCd9thjATm57VLD7J4UHzpckm+ZTPDP30rbvJ
g2WntJXM6UwM0s19vN2MB93GgpZmMJjY45HzLnwHDy5lbg7XYTJx96nplac+
cSPLSr3siSgT3+z5yftOkBull/suF9rIxMhyg0cxN7lQz1D0peBeJq7L6JM1
8FsAnw8bQt5RdkrfZnfQaAEq9m33vanIxPGEyyWqEgsgv+bwvf59TLxW+tu4
KP8vMAZsThcAE4vz4y+aDP6BGpt3xCVlJqbV5vmfPvQb0NnKtVqbiUqnNbbU
is/Bze9X79zXYWLDdfnn0d9+QqrF3au7jjNx3u6td+67n7DqiJ9Z1AkmVru2
yQWZ/4RftMJ9t08xsffqI1/luFnIChCakTJi4sAro+R3+2ZAKj77rOd5Jv6N
0y+MDfkOQZ/1aPusmGi7YXUb59p34GfOp05QXvREiOuN+ncYclacO3GRiZ8f
pmf3//sG0aeLvVbbMJG2ipZpceEbrBaqTsu8ysS1Uukx0qpT4K186eS5a0yc
itbO1JOYgj93ls6LXGficucx+4olU9AxrHnI6QYTKz7MnO5ymoSgkoZGFQcq
ngdJ12uuTQDDvmu+6y4Tc7eWbA+1HwfnZKcXT5yYyD2xXuzwnnGY6lujpHSP
ikfz+qTtnzGo0zJ/HH+fidIvb1uouI7BY6kvq689YOKv5SoOf4JHga9rSonf
k4nItV13snMY/ihx+8gGMNFlQNYjtXAQGn0at7ZRvvytoNzbdxBS2qJr7z1l
YuWdsyk15wbB5LLystpAJvrs2ixynjEIRSFuj88HU/UP2uH812YAXL7xPwp7
QZ1/1scruasfDPe1bVYKZ+IGmbBVecL9IOeeUD1CWUHDiBUx1gcDYmq03ZFM
NHSs32ET1QfKKo88618yUd7eMvg70Qf84cs8eBOo88fWVK/m7YWeka5NKZQ1
mvadKh3sgTz51MoTiUw8tflUY2Z5D1ys1uSLTKL6VyoiY+/WAx9mOG6KqVS+
n24duyzYAz7qjId2mUzkRARXnlL+DMQvlkv7GyYe8cye0+XugJJLGvKct0z8
p2e1TOlTO9h1PRhSKmJi8oJnTXZYO9QUT6mlFDOx4LJIvfTudnBzrybul1L9
FTOQ2Hq1DeZIp5iN75n4TdPysv6/Foh1y9Jrp8wdazgo3NACJ+dGBTjVTNz7
02+3YlwLpHfo2c59YOLoxKx/j3YLXIyS2f3hIxM7hx6nqyc0Q+f2wSq7Rqr/
b4TpcTZN4BW58s7GJiY2qqRkTx1tAgXi+NZ2ymEvrkwEr2uCp7NvfZVamPiY
TO7gbWsEzTdBhqx2at5CaGlb1RrhnZr6eEE3Exmran3r9jVA7NlMIb4xJgp2
zOUNG9XD2zGRN2mUl8k1RIkp1UPTFaeL+uNM3Kq52dNLqh64nVTLE75S9RQX
b6X//ARnArtvq09R8Z6aDn8T8AnEKgWHfGeY2Nq7UUk4qw7kNa74751l4qbs
d8Wb3evgWEMz9lOeIl1+xRrWwe3uyFD5OSZ6xSxyKuaqg5afu3VafjHR7sUy
q9N6H8FP0rJQbBELpR9Hr9aj10Ji6vsL5ZQzon/+ZI3WwDv57WzbxSy0SHdC
qdIamMJfV4q4WHhS9EDCQfsa0DB4ssmch4VRTQNbNvd8AH6vQr84GgsrGrY5
bCiphrtj5AVZJgsPjxVbe+ZVQc0uelkLZY9TBfvSQqtAzGXZmrssFmqo09qW
OVVBvghvUyWbhVZfFu2+plQFc8o/0ViUhSG2xyXqPlVS/y/aRDzWsFBH98GE
7kwFXNB9UdUlxcJlpy1X7zcvh9yIkPUPNrPwLT43c1ErB/7xQCfJLSyUsp44
+E22HGJdOPLXpFnovdHmhBp3OQyk3Y/g387CuMRYwe74MjClWTjIybMwd35E
8uB8KegWbZL0PMDCAK8vepysErgt4/V0H7BwuPfxtvzwEoh4Oc41RVm8OS9I
4FEJjLqnd506yMLrOcvLSfMSuK+t6L/hMAsfmtXyNNBLILFXc6FIlYWXhK+8
ir35Dv4tsW+cPcnCr7OZSs8Mi2GDfevBRF0WPrqv/yVftRjUhvemGZ9ioaNw
au28QjEE1PzzLNNjYZudmG0Vuxg2P/UEP0MWmvkmjx1OKYKTEs8Tt5qx8KjI
0SXbM95CgkrZfctLLJSb0JOQ7C4EhTrdEUFrFn4JeHBWubwQKnS/6GRSfrRx
Y7NbUiH0mdPWc9mykCf3udwNh0IQu6tdHm7Hwm0qckZzRCFw0jsFOm6wcIcW
57KJTgHcEJ3x07nPwqI0t/H5/jzgjnT7PUe5WlKSI/4xD/w2iVhGOLPQMNG4
5mxuHqTKK+6cdGFhsUf3USXvPBjWdG55/JCaF+3pbHWFPDj9YJlYlRcL/zNK
8rH1z4WD4xtjDjxl4b6jOSIpRjkQG7xq7TTl405Z5xxVcoB2lBmaEEi9v2hp
4jn5HPgvktuXDGIhiB3RerQ0Byz1+h1HQ1jYP5eeL/Q2G9xKIk4GRrCQ2+g1
vXNjNnwIWsk9kUTNH7tnXIvvNew4Qtx7mczCZIGL0hKzWRDwg/ZLP4WFNl2b
/9AHs8BI89dkaSoLxWoEc6VKs2BiSVtHcDoL8/bsjdG9nwXClwMzlXNY2LlY
4+b9hUzQVaZbPC9h4QtDm+17BTJhHZ9C25FSaj7u+Z+99TsDpqpMtKYoh6Wm
d1R+zQAv9WRFpXIWrpXz3+L3XwYUnzhKjlSysPVgQ1dbWAZsMb9bLl9L1feQ
3p518hmw+N7IppoWFtZ4qgVw26TDRxR+caOV+rWzWl+W60w6hHHtJte0sXCz
6XnHpbrpsMvNddHVdhbeP1r54cj+dLD2XtPG7mKh14DDGe3l6dASqutl0cfC
5gjhdfEZryAl9934r3EWLv56RXB2ySvIgHvCkV9ZiNw/Zff+SIPcin3yRydY
+Mld703Y5zQoacy54z/JwqYOh7QXr9OgdSqFJv2dhQ30PNLfMg24pUIkjX9S
+/We1ShRkQqGwVctixaz0TvmrqfasxTgvb2+8wGLjY0y5nfdrJOgd9NRhgKb
jdYfdpvUGyZBYdNFlXHKoXbz8VtVk+CaTEbWyRVsHB0NWiMqkQQ9w9Rqi7Ex
ARXoxT2JUKBncqBSnE39nxJ3kTZIBLudwRG0rWzUnW0Z+q6bAB1Ty8/6KrHx
m0js+O/bcbA280zgnsNspD8/sWLSJg4sb2RV9lJ+fGtc7atZHIzPG2yWPcLG
kMB4L9rROFjgipn8T4WN0ktz11YScbCOqXib1GRjh98z7gNpsWC15zznmT4b
FU0zheq/xsC0c1Hecxs2To1WPpR1j4bllxz+nrBlY77y1+bEO9EgeVLmIO0y
G3Nth+el7aLBWOJllb0dG5+f0/XQNIiGimq3Fq1rbHTQ1+DwSEdDMKE5s/gW
G529eQr3NLyEAzGdO865sjFm9S+LK5IvwbNqPm5bKBtdT3hEcI1HQqe9anAv
5Zp3vws/9EXCto3BXgFhbBRR4/CEtUXCf84KNr+eszFvQvivTmUkiCnelK2I
YCN/wr3XjJeRkJoyXWgcy0YrgfaxntOR0Og3Wf84narnJJy6918ErDEe/D1a
yUYhx7230qLDoX2N5bBvFRuNwqdnrgaGQ2BfX6PCezZa3tyyDj3CYalVT4pb
NRtN4xLO/rwUDr+udpiuq2WjcKrBy1a5cGj2qC81+I+NEVZtNiblL+Bx1ptH
VZ1s7Nfv7dw3+RzmBQNXxn5jo68xfwHHKAy07D41Dq4WQS+t5GObDINhZKLO
YIOqCE4VDc+KzgSCa8H9wz5XRfDQjjvBm74GQOj+girHcBGM8qpaVbHEH2bb
8kUuvhfBzWfzDly94Qu1Ww5zdc6KIH1Joui1GQ7k7+e60kv5fOoj69WU47WK
u4coL5wpKno/zQHXa/vefKO8tUzs4lrKewvkHfh+iqDkq9SPtd+p82obJ2Qp
T/+SbJGc4oDLBd52T8oflhwZ7B7lgN2dMhUO5QOLT595RNnE2yUngPLEhnef
d1Hek74QEE6Zq/R+l/cIB6bmZrVeUyYCa0YVhzlg7PGlvIfyC9snh4MGOaAW
EiM/RFn1YVaTEmWFZIuXY5RndKSsJgc4QH7qvj9L+fkUt+8RytXsVsWlcyL4
xunB9+k+DuRKBSbSKd+8J+0TQTlWkdo2yisWNm5Rp+xs+mlWnPKYHu3Cy14O
2F7xObuJsuPGGQEtykYP1BulKVue3JH8q4eK5ylNSfb/57sKtGIpK8RVpe+m
LFTwbFqH8v8A9L512w==
      "]]}}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwl2nk4Vd/XAHBJkkrce869lzRKIjKkSGotSYZMJfNUhgrRoEEZMiWURBKS
ITORMWMZMlNKZIoMDaa4fQ0pDb/zPu8/9z6f5+5z9j57r7XO3s9zN9mcPWrP
ycHBwaI+/u+7NUfN/VoKP/7LjtTrMA2FUytC6jel8mOP1baxdNtQ4LTrEGii
3C26abOvcyjsEbRNZ6bz4+ilIc8DPqGQ6ufdWZjJj8NyejsW00PB26Rcevop
P77ezcGj9zMU5DllP9mW8eMS95ahZt278NBQWEenjR955n7ppv+6C1/zPEqS
f/NjmGSOUtn9cLi87fH0r60CGFknleGtHQHLv7nd9tAXwIJKPUF1xUggTvXU
+HoK4P6+vOSJhCiwW/xbezNNABWD5vadgodwq8R5Yl+6AK41P2Kecfgh5F/6
IDBLGTZufD5j/BCWsEutjmcK4IBsl1PiuYcQ9+nigmK2APa11KmfePwQul6N
bx/PF8D2RNfUq8tjQSvhfbj2CwFsKNKQP9gdC7JqOccFOgXwu5H/NsHIOLjb
UCBfS/nghjN+p5PjgK1ZynPlvQDKfDH5Wp0fBzm6tbkfugRQoPHVdHhbHEiY
9HKk9wqg03nbVU9XxMMWJ+7E/YMCmPuePdx1PR5YYVZDThMCSB6u+cZSSoB/
H/hs6pfQ8LmRvMJTSIRGV1VZZ04a/vo5uzZbPRHu8rpxEEtp+OywckiuXiJs
VByKO8FFw9Jdk+KN1olw4F5B329uGrZ6LX5Q8k4Ef00Tw52raFh0lb9BriYR
eJ4laiQwaBjtKz0QrvUYBELlpa9tp2E+efeFuUsS5Gq38UtL0nDjGrHxYrck
0OV1/G+E8lebhxbCfkkQfCO+SGcHDaMu2ez5+yAJOD14lTfJ0lBQfrxla00S
zJ4e1GjcTcOpv5oadFYydB+4bcM4QPWvb9P0sDkZ4udH7ueb0PDFprq0dO1U
SG3MvmBuSkPlU0flbpmkQk7MFT0uMxpKn5jZc8UuFSr2r1xhZE7DkZvvR508
UqE3QM7jpyUNRU8ODDc+SQWC6WODttTzr5i2MVuTBoEK66VfO9NQZe2gidiH
NDjvZtI45kvDewoXeX9FZkDumtKjoX40zI1LX3I4JQPYqYL98v40XLh49Eda
QQa4dPZ+v36Dhh0QTQ94kwFn5CzXMgJpyFjPf6uDNxNOTp5wORBCrY+cbki9
TyaYHXciHj6goe94Qr7t1SyI/tESh1E0XPnx7we9oCzoviMp/oXyoKZ1t2Z0
Fhg//7ZPNoaGvRP8vqdKs+CY0LlTDbE0NNfNsD/8Kwv0Oi6WzSTS0I2/e+Cl
1xM4qOF1XPsJDa0H3suyw7PBlFvtzz/KQYryHUqp2eBSuzKmIJuGQxKJOlGl
2RCN0e+EntLQ1ala7PZgNkwpFh4ay6PhX0lvGTPpHIgSH5cMKKbhxd7IzpQ3
OTDJa/SzsoaGncH/RM0FcmFJk3Ck60sahnhxK6gL5wLj5oicWC0Nx4gvdvvF
ckFl6fkzd+pomDLnaqC7LxciF299tGikoaK7mti8Qy7gZHXdz1c0TLerMHet
y4V7r6TCd/bQsD/OU+bBjTygJYhI9FHW136eKhOWB3cvCNb49tLwJrd+TGds
HtxmLvvvbR9lLdruA4V5EHDiw5FzAzQUufD3W9pIHrjNBQtkj1DxpGIsUnAw
HyzWjd7dOkVDofg2TWX+Auid7t/2mvKJ3T7ldusKwKTmXdWlaep5NH7WxEgU
gOGpSnYtm5p/99l+abUC0M2L1LedoaGds6gRh3sBqKip8Scs0DDTV1z78FgB
bHV+HCq4lI5Cx42Ost4WwrWTf/e0Ub61/oGX5VAhvLY2++TPRcf52AtPi74X
wpWjtD3sZXRspNvfj6cVQaOC73A9Dx3d5BvKLQyLwJHTZtdFPjqyHy76PBgo
gpwHm/peC9Jx34pLxrMcxcAZ5nnDX4iOvWf6F7roxWAU3COttJaOx/ZvaG7f
Wgz/PML8U4TpaChUfY9Ppxj0bZbs8NhAx717Lvx8HFMM/0kOeYuL0vHwkcDa
E0olsKsmYau/DB2NSi5csAophRrve167ZelYku0wsCapFHQh4P0o5cOmLa8/
lJTCyQqnAJ2ddGx2i55v+lQKD57t/sLYTccdrba35/eVwc/M1pSMvXTc7rTZ
+eVcGTy/t7Cl7RAd198c/xfpUQFaRyc8fNTpeMGSL18osgLe8w907NSgY/l/
Gv5VTytgKuSlf5QmHcUzOt1dRypgQ+CdTye06RhZkuzmqf0cfDy2JM8eoeOK
TkUvU9EXoGavLyJkSUda65V5vtlKUFK8cLmKskZYcowoTxXIrIxoOmlFx/R9
e42NGFUglNd1tsCajkonTkmw5aqA/duqQtuGjsmCIT08zlUQE+FseP0UHbvf
vjaY/FQF0y+Db346T8fvW1wS3T9Vw6fIJ73BF+i4p+KlktxcNfQ6vJaSdaVj
3L5PQhzcNVC7htbheZGOmhPF5sNiNRBjHr2ReYWOR882Z9adqYGDM2mlmh50
pHNeEJderIHozXUTOQF0tAoyiRzaXgt2wjpHbG/S0SYr+NQdlVqQYXQ+YwbS
UcEhc8LSuBbqV3z28Q6iY/y/2J0n/Wphhs3FOnqbjj8UAv+I9NeCTuXBg3Nh
dPy623A66kEdLLGojd37iI6xBJ9QqUQDtBpqc7IpH78S0Wam1gBReh2nkuPo
2J85SpM43gA7VD/JrU6g40ED2cGjkQ1gLsHVNPCYjv+ZlW3q4WqEogXVOZ90
OlqqDN9mTTSCw/2XOo0FVLxvGz2S8bEZ4pKu5ZwrpO5fQuon/22G9jzZNYJF
dEyIiDDvWNcCSq/j3zg8o2P9kVPunRYtsGq5h8HKUur6uJe6Dv0tkOMmb6r7
go5OIqCvPtYKM+Ypdh2NdDwjuH4kS7oNtjla1Hk0UeNVXtWobdQGFm70raLN
VLyPdj2S92yD2gifr5da6Lit0Ly1r6UNIl9ZOjJf09H+gmuFzu43oAiMc2bv
6Hj1/fEwHHgDnptvug/209GvVGW4QKcdAo2TpyIHqHi7F7963+l2CL9dfULn
Ix1LOb5s/ufbDqnzi+plg3QcdG+9I1PaDm3N54iIESoeLw57m217B5tczbLV
x+g4OZlJt6V1gGT6lY1/Ke/MfsjbItMBu/sj7hWO05Fv3oj3il4HaKm3Xd00
ScezM+5qcSEd4Lr24KHFKWo8nn47cXUn1L2U+pgzS0eePZXMu0Lv4c2C1lH7
OTr+NtmUFbvvPfRKna5bO0+tf8ZnqcXj72HqweOswB9Ufp18M52X/h6YZ5hu
Nr+o+nEktXSzchc4EEsEGBwEtnj7/eC42A2r7TsOeK4k8PHc0xt0mT7Qtw60
lVlF4IG1RbzPrPog3FTZf4RybcWyiVchfcDQTa7T5CMwsONmre23PtiocFGD
IUBg/FTMcYP8DyC3gtB9yiAw+kxXaMGxAbi4tNHFhklgoSG9wfnWADz74x5K
sgjsuTXjWFIzAHv/G3njLkiNR7W7hHvnRzjYV2CgIUzger9t5darB8Eo28B0
aBOBuW9frAKuIYhKW34tYjOBUhK7kmK3DEFvYnmMugiBq45KTBqoDYF1pMiH
7C0EFiv8ftRxcwhOX5+xuiZGYFpK2vQBgWG4duSeHV2KQMYDGF2hMAIVh9Vv
1FMWx6rQFMsR+Ku2mHJ1B4GzXgplef4j4Kdk+/WjNIFiR4ZlN3eMwG2RnY5P
5Agsczvnd8/tE8TNtZ9VUyTwh9TnG/o9n0GT69mtScrCRqzBNK4vMEuLTru3
h8DNXAM73GW+wGEZ649DSgTqc6+ujA36AguOE7pe+whMEj45+VflKxgMLt1R
fIC631ys02zrKPyd+qJppUqg4a6UvWyOMcj802S/7CCBnsJ0RuCuMeBYe/eR
gRqBHXarZO4njEGOofBqtjqB+77Oip72GgeeFvnJbToEft2Te3+H+SQU9jB5
3lD+HteWy4yZBOvRXyJXdAkMMjf7Wd8zCc+WVZvX6RHYJMnU9Df/Brao03Li
KNWfStOp9tNT8KLIPjPGmMB1xwsvdRuxwaFWo17FhMAj9+YeHvRkA/Fu+/Ao
ZSV9niTpZDY4sb8LKpgRGLPxkm3wDBtY272C3lkQ+KXs8dqQqO/gmhB5epUN
gRE9HC5lf/6D2c6PkdWUlxqId6+XnYHLK8XrLtsS6PD3YJWq/Qy4XS7fNGRH
YKlu1IMPbTPgpT3YW3SKQLtHavVPcmbhn4/4CqfTBI4KhjKvjM6Cd/EFhY0O
VLwMqmStF5kDv83LIoIdCSSc8/aTD+cgcEFcx9qZQPK3tZdA5Dys2OHqTrgQ
6Kcb3r67ax6CbSsymihrpf65aSb4A0Je63LLn6PiQUpu2dOEHxCe7FrJ40r1
73n55+qKBYjVfy6b70bgZLisElt8EdYFcB8/dZXA38dztu/zXIS4cr07wtcI
vDlmMdD9dhEStw6PB7gTGO7ycDv9+m9I/c2dbO5F4MPdr8J3jf8BMTn9dv7r
BCaItu2J0fgL6aeiOeopq0hdDfVI/wtZ7dstZXwIdPywq639zD/ITddnLPMn
0O38Ys8efw6UHYg+WEb5jmRBZWkpB+bTRy6cvUEgV1KzpwCbA4u8LrX1BBC4
R/eZgqH1Eiw7FhOYE0TgxXJhqwANTqzl+PTLOJRAluz9tWoaXJijZbFlyV0C
+8aMirYEcmFURIdOJmXbH88lxhu50HlbXcLvMAJFZYUrd+gsQ4Z+6qHECAJF
vo/ytx/nRof40+Hj0QQq/Gk/W1LEg8fGBsvvxRDYlrtC4AnXCty/0/Sz8kMC
uxyMlz84tgJpDZqKobFU+5vJJfo/VmDFlMTAzngCv1XkN+3UWIl8+7+JeyYT
iPm32x+TfLhw095gawp1/VFHqXWufDjytt+jjfLxFZecLd/yYYn967ZNaQRu
8PD3o99dgydCnl6qzyBQW27H4ZJgfiz8cL56zVMC984dbj66loZxW8fHSygb
WE1P/adDw8BzNoRNLoEBV4O5o7xpaMllcKowj0CXTXm9rZ9pyC0pv9q0kEAf
rlm51zl0NHWfN3lcSqDA6uo8e1kSPYKXyWiUUfXncTHvWSsS46KJ5VOU6UOh
083BJA4XyxXtqSDQlB54T22IRK4GldsfKXP9OJuouZKBYu/1bW88p/JlaeKo
ozwDz8y6CLx9QeVH++CGNj8G3lnqNXq5ksAMDgvmjwwG5tJCKoWrCBy+LufI
3cbAOZks59PV1PyZ2W1PJZnIwrKDfDUEGht+uyWswEQlvaa1hZRJ1Ycqh4yZ
6OX8tenfSwIrLUzNfe8xMcFjPiGllsC1d/i/ezxl4stby9wO1xGYx7nZcLGJ
iV9iCL3vlD91p2b9N8xEnkyRrQ/qCcx53xtp8ouJEqVyf5QbqPZ+p1VE+Fmo
3ajSMUzZapPFac0tLHTp0s8KbKTqkVeWXc1uFoZ9sfbd0URg3CLKB6izsHDO
xbSD8t9oN98QIxa+5/KSudZMxW/B36dvbVm4QA9ZvrGFmu8Bjl7tsywUEokd
qKMsnvp48KcbC5XlsoqcWgmUX55xpsObhdYqZbcFXhE4t+bIbHcAC331m2yL
KSc+Ua/kusXCZOtuJcvXBEoeKxszDGFhvctXgaVtBD6RuHayhfKY5/xoOmUm
S4rD/jYLV4Usq9J9Q+BWvWBTkSAW7oglHsxS/vgycXK5Pwv1s0RcYt5Sv7Nk
VpCeLLxQJqeG7QR23uBs17nIwogmFeEvlAsOnNvzzIGFxd36M7feEXj+Ofv3
UUsW9n61bpbtIFC9vyteXI+Fv+ddErsoO7JU7RFYuJ7by82zk3pfvFDVjJNi
IZIheiLvCVSVP/ZTS4iFNltitzZRLqaJteouY+GNnVl/XLqo+rNusCR/ionN
R5qyyroJPG1/K7W6gomTx7t9j/dQ9aZV0MYnkYl8576acvdS4+WRi2r2Z6LB
nWU8R/uoenznxZSYOhM/NKvYjfUTuMRE8PvJCgb+7dHfGzpA5Z9C84eGSAZu
GrOm7fpIxbts5SjjHANPLveq8h4kMFYwWqpoAwOnVcuEWSMEhrUELTlzkcTR
q9LcrZTZy1wej6iSOPg0efr6J+p+peGyUTQS36wNrfn6mcD5u0EpP7MJzJ+x
PV08SqBQyjdxY2ofdzFpdaHRFIGpgqwHv2RpeKbH9xHvNHX9Fl5Fid8CaLdm
IeAF5TrmzhGXegE0dB823fqdQF6XDXyW5gK426D439wMlX+jvsOlN/lxgfPE
4fs/CXTV7N/WNLga2Yrv5TV/EVhhu9hX4bUaR10Or/9D+YV1J+fjdauxp3fX
d7vf1PtoiSK/idUqLMvnfSD/j8A/js9K5yd40cOmcPgdF4n/uHWiPZV48GKU
eOvNZVTeq94r3j22HM+8jivay00i69tF4+XRy9FSKSgwaTlJ1UtpqclFbtxP
s9rhykvi8mjZ/Umvl+GSmuXXaPwkyr6KeimZuhRtVj2xKqPsU6/+d9xhKdYa
6avaCJD4yyn4ls+OpRgwEbUqn5rHO9fXqe0o48RVpET8EZJEOwWx3gN9S5Bx
Wrv2rhCJgol6GpIHONAtn52uuJak6lnoCS0aB/b8jggZpMyr98ZCrv4fxIb1
G8msI7HwoXW75+p/sKncZaxtA7WOSv5ZLVl/QJIvjG+NKIlRA0bN/CKLcMdk
18wzykN+e7+v7/oF0497uqy2kgiiP34fuP0L8hU2J+aIkeglV+Jf/vMnKJwo
2KkrQWJgnEPawuACHCjsNAmRpsb/TmlWpH8ekv5e3bdLhsSvBuYC6QnzwKW5
fnM/5UBuu0tL7eehof/khKQcicH5voZy7DnQXr7g2SpP4jNr4UN0/jkwNhNK
WqlEonbD4SIR1xm407Ilppqy/VnaeNj+GahXlg67spfEBanYbS95Z0B+w8Hr
I8okmssGGZ9I+Q9on5zNy4FE0kvtu+zId2h1rqY7qZF49a0a3duLDUs/tvBu
OkRipJu0XNYxNijpv1/SRfnSwFlNx+1syJCbmD6gQWKl1Mt++aJpCPhBtAod
JjGoewP9n9cUoM9p/2Z9Eo/2hf339sAkXPnvgrv3ERIdenPW+PybgBxbzwu7
j5J4LeOVO+35BKw7FH78sQGJWSuuhbxSmoBfvBXK14xItAjuuaiC41AYsWZO
3IJExa6PPsrnRmFimdDkR8pSCoeSru4Zhc1XtozctyQxd3VJsP3SUQgz3fN2
iTWJFbJxNfujv4LLetvsnhMkvhQw3pXd+gXE05/ZB50iMe+8YY2B/meIGjDm
VT5Nou5lpc8btnwGHvJnzhRlQ7GB4LMLn+CLz94FA0cSe8PaksMef4Jk86rg
9c4kvhOt6X3yewTWr2l+WnCBxJICdlZx/TCEqDkdO+lKYkfwysVPicPw233V
TxZVJ5KaCt2ueQ5D36juAa9LJG6pcBLp3jUMUTXvOjTcSHSjOd/enTUEtMv9
P/s9SeRZeSfEMG0QfJ54xd31IlFZpRZtAgeBPbxBVfU6iTMZH/RrHAahTc/m
dro3icKi53ndpAbhtvjX9a5+VPy2l0/tPv8RlvezVXmCSHzdRhcqK+2H36pc
d+QiSKwBq5892b3QcadDqodyl2Cm/ir/XsjuSX51/T6J+6R2HIwz6wWrs2qr
X0VS6ytWE/mBpxcqYwJun4qm4nOc+TnZoQd8v/Pcio0j8VY6a7v33m4wU+6R
UI2n6uC2lAppZjfsvJnRPEb5t4kve8dMF3wS1uJVSCRx56b2Z6wnXaCmcSvo
bRKJMUNWT3w2dQFP/OpA7gwSH0TwrftP8D0MjvWLZVM+lH1/QOFXJ5TK5zQY
ZJKYf8HAqaq3ExybdZcnZpFY5Wq/6kFsJ7TMhQbszSGxv/7or0MinXBHm3bj
XAGJNifOaq9T7gD6L4Zv73MS1bbT7ZwftEONk4586AsSXR0Wyqs92uFcv98X
1UoqHiWjrXVt2qG1iq2VXUWt99P+9Xw72iHgZjPd+yWJNw8KcTAb38IC4ZUi
2kT1X5/09M2yt5AaUGjcS1nLZKWS/dQbOLYwviK0mURLRw0U73oDeX3GLgst
1Pie5zXty3gDjo9lFVpek7gxPt/xlt4b+CD9ufFcB4lCXmYGKRptEJy41l20
k8SmsvZTuhvaQJF+VKqX8uTnt4Ia86/h/vyLMNUuEn/cP79RIuU16D6PMmP0
UvlmtLE5nPs1VGtpT5Z/JBFlhaZ0u1oh1b5gzfIJEp1Uvs2nZTXDiwnW86eU
b195Vb89pBk6z3s5mkySaM16v3nEpRm4vDTrMr5R9dlO+9m4XDOciPx4TZtN
7QPPLVlz5HkTCDes/BI2R2LRgsdV4b5GCN9mVyHMwcB9T2U05hUaIDOnyaGO
commlpnIhgaolpdmuixhYHnhlaor3A3Axl/nKzkZiM1PVO531oOO6V0xm2UM
zPzhcyTmUj3wBFeEp/Ey8PYuJ96gijrwnCAc5EgG3uufNuK0q4XW3QK1XZTF
DVXC2Hq1IOy7eoMng9p3qmSv/b23FspY3J0NTAaa9F0HR6IWFtR+oKUQA7uy
Op/F1r+kznc9rEDqvS5ekpm2XeYlOBjGNfaLM3Bprqx2BKMGShJiRPwkGKgY
z32ncVkN8ExGem3bzsArFjqL9LlqSPUNlXeVZCAHVj2dflcNn556J/BIM7Ax
7OXPovBqsOa1ddtJ7YMtJTRXmBPVYFgpti1oPwM/P288xiFWBddkg+8rA7Xv
7ZOKimFWQULSJCeb8uPVljlaPFUwfjOv30iFgQMqbOuud5Xgrb/33paD1P1J
Nx1Ti0rIHNL9W6nJQE3j8CUCl17Av6WXO+aPMVAlKDPg2MsK2HK5WyXTkIFK
7VKZZzIqQGtU6amlEQO/dp8PSAutgIjWf0G1xgzM2HzL1d+iAiTuB0G4GQPz
g67c6PhRDse2PsqUOs7A1Vn7LWPkyyFDo9bbzomBe7d0eyvXlYJim+HYyjMM
3J64UeNZbinUG349UkDZ3WLpaY3YUhi24RXhdGFgU8Esq8i1FIQ99eviqX2V
v98ekScipRCa92FF3yUGMvz3DbNvlsAlobnwI94M9OCW3DdyvBi4EgMWFyg7
mOToiugVQ7gYyy7Bh4Gzh8MvXdlXDDnye3dN+zKQsNG7ZS5UDKO6Pl23bzCw
/87l6DWdz8Dcb7VwYzADtWJ9ePh0n4HKpGjK/vsM/N2lNpusUwSp0es2zlKu
+vYn4bNyEfCqkw8zqH1g55PMqT2SRdCeyBVGRDHwkrbEjPDKIrAzHvEYj2Eg
Z5O7aVNTIQTUJByLTGBg6yPBtY+1CqElai3XVBYD06ojz4abFoDMIfr1pCfU
erSNGvEfLoCIGd5fJtkM7AY+tzjlArDQ/TX9MoeBBvSa5RMbCmBqaU9fdB4D
9Xln6g58zQf+s5EFasUMXHdl6b1ut3wwVBOwfVTDQFVO27NcaXmwebliz6GX
VH48p4vJP8wDdqOVHpvyy/XHss6G5kGw9pO9qnUMfJVA9q1wy4MqA3VirIGB
to+ltzVp5sF2G886+VcMZKft1js5nQtLro+JtXYxcHhDof2zA7mQXVI9+WuS
yk9ryZ1Jq3MgH67zJ35jYLzdlhXOi9lQUq8srz7FwAdJOmbqY9lQ01Hsfm+a
2jfHrt8gVZcN3exsXsn/GGihHLGm3DMbuMRjtln+YODbDd/ftH1/AmbRF+wq
lzDx6+uWZIvPWcB9TeSDH4OJ3QznHYbzGTAkpk5TZDLxaHqD6OUvGVDR6agx
SfmMs0hsyvsMcJXNLzwmyMSUG0mV0sUZMDiKIaLCTKQpfluccsuAcmOr/Q2b
mKhQf/Us7790OLcrOoFXionXmb4+qmQ69LH57MNUmbiYrlV57WQqbCw4Ebnn
IBOj7Dh3RZulgt2lwoYhyhpFhxnVuqkw+dNUQu4QExfg24ykQir85UyZbtdg
YjWzhOciTypsJvdeI3SZSN64JFn8JAVO7zkV+sCEiTtVq73eLCbDrE9l6SNn
6hySL5jsnJ8EfE5ufwxcmPiATGp9nZYE247JqvCeZWJ7mvINxUdJYLk1qfHy
OSbuPr3gIxqYBPXNAV16rkw8FpmcZWeVBNF03bklV5m4T+H2j7OrkmB/ygeZ
k9S5p10rI9jL5TEENf5M2/GQiQPDR6wrNRLhw2XN6CHKpY3qDElMhB2i0cER
sUzcc8FBIl4hEdp9FJ1/PaLaJwoeidqaCMJ7r8jVJzAx8sr+8U6uRMjJnq2w
TGViy4GPs241CdARPv32dh4TX9jV7dx1IAE2WH5eHG9gInvGui4lOQ5+roxc
m/qdiUs/KyvUhTwEvXNvOj6vZ+H0m5jnHIVRMDbVZrpFk4UtPndapM9Egn+5
98E7F6hzsILmdN+1CHi4r7zRI56FX2UOBX5PDof5njKWYxMLjRRFd4qKhsGr
7Qc5P8xT58xi7WJ1l1Ao28d5fojy/CXDmjnnUEjXq/r4hfJTj7jXSZT9XZWf
f6d8SOv7+L8zoaBULu+2/AcLWTLRciVOVHst0Sk5ymdcJ/9scwgFXwfu3iDK
vEdc/VfahcI591qNUMoxt9yKS21DwSrEtziCcu/+5xOnKO/J+xsRT9lwjsuk
1iYU2AvzekWUD0/eBs8ToWAZ+LVukDKf4VWtaatQ0IpJkf9CWXqNa9gjyopP
bJMmKG/pzeo5TJl489F7nvL9wF7ndMtQaGZ27121wEKH8O0pJyxCoUQ8MlOA
cs66dfNrKKfupbKB8pIAHY0X5qHgY/1mfhPlkneObCHKLufv2ItRlnhifKjJ
LBQs/LQ7JCkXqPs9ukJZ6z6vqhzlhvhvc6KUFdMa8xQo8+X8//8t/geKTveg

      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{4.5, 0},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.415552332181415*^9, 3.415552436215261*^9, {3.415552490787064*^9, 
   3.415552515134419*^9}, 3.4156242601178303`*^9, 3.415624334390402*^9, 
   3.415624492524732*^9, 3.415624594894004*^9, 3.41562651118985*^9, 
   3.415627128048079*^9, {3.4156273060211763`*^9, 3.4156273301363497`*^9}, 
   3.415627445811562*^9, 3.4156275086091413`*^9, 3.415627791398924*^9, 
   3.4156278931687593`*^9, 3.4156279245676413`*^9, 3.4156279820035887`*^9, 
   3.4156281438551197`*^9, 3.415628213236376*^9, 3.41563003895278*^9, {
   3.415637363498624*^9, 3.415637414561982*^9}, 3.415637457026574*^9, {
   3.415637488464286*^9, 3.415637515793867*^9}, {3.4156377469134007`*^9, 
   3.415637772844616*^9}, {3.415637818504286*^9, 3.4156378255236673`*^9}, 
   3.4156379051470127`*^9, 3.4156379611517153`*^9, 3.4156380414072123`*^9, 
   3.415639403781006*^9, 3.415648058975222*^9, 3.415648151491064*^9, 
   3.4156486078197727`*^9, {3.415648768603374*^9, 3.4156487789673767`*^9}, 
   3.4156512218784437`*^9, 3.415651274992085*^9, 3.415651513540532*^9, 
   3.415651596635828*^9, {3.4156516509046383`*^9, 3.415651723369739*^9}, 
   3.4156993491200113`*^9, 3.4156993807727823`*^9, 3.415699526819964*^9, 
   3.415743628924259*^9, 3.415747006437543*^9, 3.415747198796173*^9, 
   3.415754229357196*^9, 3.415755518902605*^9, 3.415792220954727*^9, 
   3.415793217291051*^9, 3.415793280986595*^9, 3.415793313068069*^9, 
   3.415793350922172*^9, 3.415793397158599*^9, 3.415802428210343*^9, 
   3.4158135858626957`*^9, 3.4158140382856627`*^9, 3.4158142714681273`*^9, 
   3.415814538240279*^9, 3.4158145705640507`*^9, 3.415814631308566*^9, {
   3.415814679836301*^9, 3.415814688540975*^9}, {3.4158174128431587`*^9, 
   3.415817493127925*^9}, 3.415817530603777*^9, 3.4158175817632523`*^9, {
   3.41581763166354*^9, 3.415817644988777*^9}, 3.415817735389923*^9, {
   3.415818377120822*^9, 3.4158184552655287`*^9}, {3.415827186623049*^9, 
   3.41582719648604*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\" Normalization of thermalized cross section = \"\>", 
   "\[InvisibleSpace]", "0.9711267544722355`"}],
  SequenceForm[
  " Normalization of thermalized cross section = ", 0.9711267544722355],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.415648611352868*^9, {3.4156487721175327`*^9, 3.41564878246176*^9}, 
   3.41565122402354*^9, 3.4156512764376173`*^9, 3.415651514985532*^9, 
   3.415651598079349*^9, {3.415651652348106*^9, 3.4156517248205833`*^9}, 
   3.415699350599358*^9, 3.415699382219385*^9, 3.415699528260882*^9, 
   3.4157436303814993`*^9, 3.415747007887122*^9, 3.415747200252823*^9, 
   3.415754230799534*^9, 3.415755520355027*^9, 3.415792224434722*^9, 
   3.415793220798267*^9, 3.415793284506914*^9, 3.415793316590784*^9, 
   3.4157933544845963`*^9, 3.4157934008114758`*^9, 3.4158024296840754`*^9, 
   3.415813587331286*^9, 3.415814039702729*^9, 3.415814272929015*^9, 
   3.415814539694717*^9, 3.4158145720166407`*^9, 3.41581463276154*^9, {
   3.4158146812905416`*^9, 3.4158146899866667`*^9}, {3.415817414356678*^9, 
   3.415817494571671*^9}, 3.4158175327714157`*^9, 3.4158175834109287`*^9, {
   3.4158176330347157`*^9, 3.415817647092968*^9}, 3.41581773684445*^9, {
   3.41581837865797*^9, 3.415818456713772*^9}, {3.4158271880910397`*^9, 
   3.415827197936419*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\" Normalization of bare cross section = \"\>", 
   "\[InvisibleSpace]", "0.9712870972924003`"}],
  SequenceForm[" Normalization of bare cross section = ", 0.9712870972924003],
  
  Editable->False]], "Print",
 CellChangeTimes->{
  3.415648611352868*^9, {3.4156487721175327`*^9, 3.41564878246176*^9}, 
   3.41565122402354*^9, 3.4156512764376173`*^9, 3.415651514985532*^9, 
   3.415651598079349*^9, {3.415651652348106*^9, 3.4156517248205833`*^9}, 
   3.415699350599358*^9, 3.415699382219385*^9, 3.415699528260882*^9, 
   3.4157436303814993`*^9, 3.415747007887122*^9, 3.415747200252823*^9, 
   3.415754230799534*^9, 3.415755520355027*^9, 3.415792224434722*^9, 
   3.415793220798267*^9, 3.415793284506914*^9, 3.415793316590784*^9, 
   3.4157933544845963`*^9, 3.4157934008114758`*^9, 3.4158024296840754`*^9, 
   3.415813587331286*^9, 3.415814039702729*^9, 3.415814272929015*^9, 
   3.415814539694717*^9, 3.4158145720166407`*^9, 3.41581463276154*^9, {
   3.4158146812905416`*^9, 3.4158146899866667`*^9}, {3.415817414356678*^9, 
   3.415817494571671*^9}, 3.4158175327714157`*^9, 3.4158175834109287`*^9, {
   3.4158176330347157`*^9, 3.415817647092968*^9}, 3.41581773684445*^9, {
   3.41581837865797*^9, 3.415818456713772*^9}, {3.4158271880910397`*^9, 
   3.4158271979786463`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Lorentzian fit", "Subsubsection"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Enmin", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "-", " ", "\[CapitalGamma]value"}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}]}], ";", 
  RowBox[{"Enmax", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "+", " ", "\[CapitalGamma]value"}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], " ", "/.", "sub"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Npts", "=", "100"}], ";", 
  RowBox[{"dEn", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"Enmax", "-", "Enmin"}], ")"}], "/", "Npts"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Sigma]Thtable", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Enmin", "+", 
        RowBox[{"i", " ", "dEn"}]}], ",", 
       RowBox[{"1", "/", 
        RowBox[{"\[Sigma]Th", "[", 
         RowBox[{
          RowBox[{"Enmin", "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"i", "-", "1"}], ")"}], " ", "dEn"}]}], ",", "Tvalue"}], 
         "]"}]}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "Npts"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fitfun", "=", 
   RowBox[{"Fit", "[", 
    RowBox[{"\[Sigma]Thtable", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "x", ",", 
       RowBox[{"x", "^", "2"}]}], "}"}], ",", "x"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Normfit", "=", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"\[Alpha]", "^", "4"}], "/", 
      RowBox[{"fitfun", "[", 
       RowBox[{"[", "3", "]"}], "]"}]}], "/.", "sub"}], "/.", 
    RowBox[{"x", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"E0fit", "=", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"-", " ", "Normfit"}], " ", 
      RowBox[{
       RowBox[{"fitfun", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "/", 
       RowBox[{"(", 
        RowBox[{"2", " ", 
         RowBox[{"\[Alpha]", "^", "2"}]}], ")"}]}]}], "/.", "sub"}], "/.", 
    RowBox[{"x", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Gfit", "=", 
   RowBox[{"2", " ", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"Normfit", "  ", 
       RowBox[{"fitfun", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", 
      RowBox[{"E0fit", "^", "2"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{
  "\"\< Overall Norm: \>\"", ",", "Normfit", ",", "\"\< (fit)  \>\"", ",", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Alpha]", "^", "2"}], " ", 
     RowBox[{"\[CapitalGamma]", "/", 
      RowBox[{"(", 
       RowBox[{"2", " ", "\[Pi]"}], ")"}]}]}], "/.", "sub"}], ",", 
   "\"\< (T=0)\>\""}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{
  "\"\< Resonance energy: \>\"", ",", "E0fit", ",", "\"\< (fit) \>\"", ",", 
   "E0value", ",", "\"\< (T=0)\>\""}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{
  "\"\< Width: \>\"", ",", "Gfit", ",", "\"\< (fit) \>\"", ",", 
   "\[CapitalGamma]value", ",", "\"\< (T=0)\>\""}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Sigma]Thfit", "[", "x_", "]"}], "=", 
   RowBox[{
    RowBox[{"Normfit", "/", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"\[Alpha]", "^", "2"}], " ", "x"}], "-", "E0fit"}], ")"}], 
        "^", "2"}], "+", 
       RowBox[{
        RowBox[{"Gfit", "^", "2"}], "/", "4"}]}], ")"}]}], "/.", "sub"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Enmin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E0value", "-", 
        RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
      RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}], "]"}]}], ";", 
  RowBox[{"Enmax", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "+", 
       RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Sigma]Thfit", "[", "En", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"1", ",", "0", ",", "0"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Sigma]Th", "[", 
      RowBox[{"En", ",", "Tvalue"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"0", ",", "0", ",", "1"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p3", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Sigma0", "[", "En", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", " ", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"0", ",", "1", ",", "0"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"{", 
   RowBox[{"p1", ",", "p2", ",", "p3"}], "}"}], "]"}], "\[IndentingNewLine]", 

 RowBox[{
  RowBox[{"(*", 
   RowBox[{"ListPlot", "[", "Ltable", "]"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4145831080620327`*^9, 3.414583111420779*^9}, {
   3.4156284847994833`*^9, 3.415628689161697*^9}, {3.415628814841309*^9, 
   3.415628851975069*^9}, {3.415629996135886*^9, 3.415629999478713*^9}, {
   3.4156380787406607`*^9, 3.415638084962056*^9}, 3.4156513364709063`*^9, {
   3.415651744733471*^9, 3.4156517474591217`*^9}, {3.415699538085273*^9, 
   3.415699590210579*^9}, {3.415699713490834*^9, 3.415699726126894*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\" Overall Norm: \"\>", "\[InvisibleSpace]", 
   "0.24300861788030445`", "\[InvisibleSpace]", "\<\" (fit)  \"\>", 
   "\[InvisibleSpace]", "0.06576650540987412`", 
   "\[InvisibleSpace]", "\<\" (T=0)\"\>"}],
  SequenceForm[
  " Overall Norm: ", 0.24300861788030445`, " (fit)  ", 0.06576650540987412, 
   " (T=0)"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.415630015613813*^9, 3.4156300534375*^9, 3.4156380196513977`*^9, {
   3.41563806821159*^9, 3.415638086158373*^9}, 3.415648103492948*^9, 
   3.41564816859966*^9, 3.415648627407845*^9, 3.4156487937323513`*^9, 
   3.415651295540291*^9, 3.415651340202121*^9, 3.415651519445709*^9, {
   3.4156517297723703`*^9, 3.415651748501418*^9}, 3.41569939305716*^9, 
   3.415699594309807*^9, 3.4156997366655416`*^9, 3.415743635125722*^9, 
   3.41574701208885*^9, 3.415754236464196*^9, 3.4157555249494534`*^9, 
   3.415802434518076*^9, 3.415813590464177*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\" Resonance energy: \"\>", "\[InvisibleSpace]", 
   "10.000938088387876`", "\[InvisibleSpace]", "\<\" (fit) \"\>", 
   "\[InvisibleSpace]", "10.`", "\[InvisibleSpace]", "\<\" (T=0)\"\>"}],
  SequenceForm[
  " Resonance energy: ", 10.000938088387876`, " (fit) ", 10., " (T=0)"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.415630015613813*^9, 3.4156300534375*^9, 3.4156380196513977`*^9, {
   3.41563806821159*^9, 3.415638086158373*^9}, 3.415648103492948*^9, 
   3.41564816859966*^9, 3.415648627407845*^9, 3.4156487937323513`*^9, 
   3.415651295540291*^9, 3.415651340202121*^9, 3.415651519445709*^9, {
   3.4156517297723703`*^9, 3.415651748501418*^9}, 3.41569939305716*^9, 
   3.415699594309807*^9, 3.4156997366655416`*^9, 3.415743635125722*^9, 
   3.41574701208885*^9, 3.415754236464196*^9, 3.4157555249494534`*^9, 
   3.415802434518076*^9, 3.415813590466893*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\" Width: \"\>", "\[InvisibleSpace]", "1.353120435839107`", 
   "\[InvisibleSpace]", "\<\" (fit) \"\>", "\[InvisibleSpace]", "0.5`", 
   "\[InvisibleSpace]", "\<\" (T=0)\"\>"}],
  SequenceForm[" Width: ", 1.353120435839107, " (fit) ", 0.5, " (T=0)"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.415630015613813*^9, 3.4156300534375*^9, 3.4156380196513977`*^9, {
   3.41563806821159*^9, 3.415638086158373*^9}, 3.415648103492948*^9, 
   3.41564816859966*^9, 3.415648627407845*^9, 3.4156487937323513`*^9, 
   3.415651295540291*^9, 3.415651340202121*^9, 3.415651519445709*^9, {
   3.4156517297723703`*^9, 3.415651748501418*^9}, 3.41569939305716*^9, 
   3.415699594309807*^9, 3.4156997366655416`*^9, 3.415743635125722*^9, 
   3.41574701208885*^9, 3.415754236464196*^9, 3.4157555249494534`*^9, 
   3.415802434518076*^9, 3.4158135904698563`*^9}]
}, Open  ]],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwV13k0VW0XAHBTigopL1EZ7r3udQdjiZS9DaUoKTQRJUWaZEoqkWQMKRkj
QzLPyixDZJYpynBFKaRIVIa+8/111m+dZ5+z936evdY50lZXDp/l4uDgKOPk
4Pj/1e7hax0VFTHUUyOtbfD3ArmIq6Ic28Vwx0XK+rchXvAlZvNE8w4x9Hfm
Zn947AVnk51Dz2mLYfaURPefBC84VUwZjjwkhuf5vM2PlnrBkSFvdw57MXS5
UWAVNekFWow9xS0ZYji+7JK8/9BdEKupk7eV3YgDQy8uN5C8ge85B1VYXhxT
j4ctior7wOd9xyg6ZAncul/pwU1FP5BO6emJE9mEOU1DTru9AyCu5H7YVYnN
2Bz577jV5SBYFWh1olFwC65lnxXY3RgChevsmRFikug4L9ol0hAKp+Q6evcK
SKF/hVy5tkQYnGr50scSkkJHPSmXQjLhq//eCwtL4ax7tIO8POFi5kC/iBRe
sObYr6JFeJ/PsP1mKTzvMC3uZxMGlud3jkcwpfBag6+QXWEYWKQmL3zRl0KW
HV+v3JHHYCZ3Y4ufjxQmNIaf/psRDrWaGXzhflIoQFmtIVkcDkyTgdmkACnU
vmagf+B1OCy6Q9OrYCms/vld/vVAOMR2cbv+DpfCLbkXNlAEI+Dj7cC3NilS
+LuB/2GmYwTY9Ty5u6dBCi31YsYoupHg6vVqnHu1NPbldalzrYgGi66lkak1
0nije1y2QTgadCkaA30C0viG9sU+Sioa1tUXtmULS2N70E1d853RkMafXmAm
Lo00HY7VFk7R0P8gzD1fThpnYj33nPkcDVpP7YTP7JXG1NdqWa86YoC/fMOO
am9pdGToLx9piwUbsyL/RB9pVKrY01o8FAu1f8w+3PWTxn6vGnu5H7Hgrpp4
U+++NPKYFp9WWBcHM9lKlc2PpPHf7Qfn3hnHwfsEQ913idLIfb29MfFDHKT5
+hpOVkljzSdNrjduT8HAZNHqv2VphGrXIc70eBCyLnkixSGDDZqThsOF8dDj
6NpL55LB2toPQY2v4sHq4ewB4JXB/bEns0p64sG1Y1LNVkAG20hdecLcCZBk
NCBQskUG634cx4yTCbC4v6LkJMig7hXanX7xRMja7SH8zEMGb07HSckVJcFM
iL5Z4B0ZvK9e4DdfmwTb+jckOd6VwZbSUs+3b5OgzCF1m7avDD5OLzdJmUiC
hriOY0PBMrj43n7FiOQzGP1DjtsYJ4OHLP6+ifB/BuJZDYygchk0IUtUXzif
DBZ/Hjo5V8pg0PzAy/pryRCva1FuXiWDrekqG5TuJQPtw4wh47UMDuppaTET
k2Hbqk0Ob5plcLgnW7VuIBmMrC4XcfXLYNIbKWa46XPw+W/97mt/ZXA+e81B
qlEKKHfOX/NZlEHvTLTMPJ0Cg8H9aeHLRDxsfrPbMQVU+ZKFirhIaOHZwpv1
OAU+L6j1/+YnYZhN1S+LwRTYM2zh6CZBwinBQIHjV1OBJyMt/tZOEu4QOhlv
l5wGObbBXfc1SXhqurTrckkamFOcVsYiCX9/DG30bU2Dgthdlyp0SXh4bOW1
5bk0OPOgXe3ffhKavCpSmN2bDtUuc20eJ0n4fNIjmzydDp5a2ste7iQUUM70
9DTOhCSD+77+HiQ8oPfhi5BdJtSZ9go/uEPCenkj42KPTFhjd5kae4+EsZlh
laZZmfA4NMqoKIiEj1Rw+BV/FmSMzCRMxpKwbGnDZY83WdDjnahnWknCidkq
++PMHPgdPPXWrIqEB5/0T5iq54BElLq5VQ0R733A7PSeHDiV1WZ/pZ7IP9b7
8rNTOTDRsxDp10bC8YcNb0fDcoCDZjJZPkT0Y3nfWRPuXKA38jyQ5SDj6qcW
b56M50KjoLkcPxcZCwasxRN+58IF0/yqb9xk/PM1Q+AFbx5ksE/PFKwkY+32
4m1rSXkgP19hrCNIxmNCDybczPNAmXz9v1NbyBj/ybLfpiMPdrh/i47cSUbV
yPvt8k358L5Gd+stTTIKpSw+N/mQDzf4YppPIRm5HXqo9yfyoeyR/hJVl4w7
FDUDVdcUgGb6c4sXBmTslyocvmRYANq9p6U7T5DRbtcGGZ+eAtin3JO85joZ
h6Tq8xXnCuHwVz25jTfImF4itOI6/wswf1qcTrlFxulWntQPW17AZYEnOZqe
ZMwP3f3gh94LCB0/U2LvR7hTOOhk9AtivqdbuiLJiHcyA+30XoLd+rW/YkrI
mEX6mPOwsAgcG92vpZaRMUN5Xe6tliK46fnjT2EFsb7BbcjnUxEEfe9cbq0m
o/W5yYglkWLIa4paxdlERo3m/4ycXIvhrxdt09kPZMw2nX1YoVMC/rM6OqxF
MpanSGQe+1UK68aqOv9bJmOKxWf9dsEyiOgDaw4OCirlR+rZ0csguULDu5Ob
ghI/jZR4T5VBla9K/fXVFKx9ah8/0lwGvzeRDOrEKRghY/3ZOrMczu7mMrZU
p+ChvOLzwgGVMLndY2SvBgWPe2KxdlolONCXHZV3UVBDaSo1pqESPAT/Plqh
RcGxmgPbbXheQUzfdE/6XgrusTy6P0L3FXRcGj7x+ygFf4tdiNvQ8Arw8Sur
By4UfDid23VupAoU/x4bSHSlYBp9h825+SqQtJg++sKNggr6C75eq6thmSJz
4IM7Ee/DqyG6tRpKC722U3yIetL1M1PvVYNq9541xeEUXA9X1JOVaoC5oaWQ
XURBZbf7XWef1sIm17OKP0soqB+7oHu9qBZW9y+lrSgn6r96rSKnvRbGE+Wf
0qsoyJ+ynBrI+Rqeq4T4OzUQ+Tu+cFu2fg0yxsaWq95T0LPimtaCch2Ihfat
Ulyg4EHjnIyvk/XgrtP7Xm2JgpdjN7pY872B0dmeDK1/FKw+mzSwSvYN5Bzt
MjLmlsV2Nw6rNss3sG9zW4TLallUmxVRO979BtxSamnlErJ4nKlxJ/hNA/RX
ZO/V3ymLm00kp+fbm0DbPkvcWFMWv5i3PF/7swlSpDMnzVAWrW9edzEWaQan
u2khl3Rl0W/LD2H7E82wVv9Zb8h+WfypICA387kZNLujbN+Zy6KjyJdoTb5W
iJvw9rW+JUt8D4mG/VFvhzu9Ibb6t2Wx7IbDxbb97XDudfReRU8ifmxiodSy
HeTjclct3pVFVq/M8T7vdqg4POD7MFAWn0f+8KJ0tsNg8Va/6mhZTLgazFrh
8Ba2+I74SZXI4qdtXqF1rzoghqwV0D8ni/fXK65Vj+2GL3fOmU78lsX4vh2q
NkXdsG04QPLvX1mU1zvcktLRDa1PevJF/8nihBb7g/XKHuAQvdh/eCUVOZ1N
OF9c7QFrvnBWgygVq8sbHlvueweMqcm2QjUqRmnG5Bly9oHrfuGo2h1UtL1x
xXvF5j6oTdtu3bmTijxrFD5VqPWBuc2d3z+QisK2N0+o2vdBEFtUmrmPiuM7
fjXIDPXBTIe2Q8JxKq5xcFHYWfYeSl5GbAh2o+Lhq9QtGy/1w91DjZNuN6m4
8Pl23UvvfjCcWKg9507FUaPMHs3YfhjZYumseYeK0hpXA7ja+mHNPUrPNz8q
psRJ2e9RGABLk7zw/ZFUvMlbcGby6wDQp0auqEVTUXvIghLPMQizPiJ7yU+I
enN3iCiIDoJvqevvhadUvCsR4dehMwh5MnA8PYWKdY/GvkpEDwLPdJM4fxEV
96ospvfiELT6L838Kqbi2O2pFykmQxBBVmgaLqVicdHggpbtELCOhd4orqTi
Yu4BbA4agiOVx/pt66m4x75cm9w3BFLH/QtMGqgoJn86RWZiCCZmSgOxiYoa
Agz3j4tD4CEruUusjYrHhkSPvRVkQ+r9T0/qe6h4fXDMd4USG5yooi75vVR0
2t9lu6TGBs2qvYZx76l4xVz6UAeyoWM2fdllkIoeKsrX1xmxYdHMwVL2MxGv
cYjed5ENdb8Stwt/oeIn9wq9TEc2hAR3Cy5/pWLSTy8hMzc2UGrUXnV/o+KZ
FqcN2j5sMKT/k/L+RcWnJlanVOLYELwu3Mhmnoqe8ZsfZiWxof23vMe+P1Ts
+7VJiiuNDYfrLYbWLlExcunEGqUCNjzMmhP4sUzFonduHauL2dAVFqTZwUFD
jukXr8rL2XDEuuLJYx4aXp0a+Pusjg3hBkdaXHlpKG2zsvJDIxt6lacWT6yi
IZfEYPRYKxtOcG4xl1xDww/Pz6R69LAh+kthAKcADbn3TVivfc+G/rYDpSOC
NJzYNGJ4ZYANlrG3xFPW0zA8OXswd4QNT71F9P1FaNhycS0p9DMbPl7MvH5R
lIYqqosde76ygWSyO9VwIw27hki1bRNssNYY6FWUoGHI8bpW5hQbnsk4r1q/
mYalea9bTv5gw2e+tWq/ttBwtiY11mqGDdTpJJt3UjR84fRMWmOWDba9O8OL
ZWi4OjpafewXsZ+VXXXRZBqKOJGaz82zYTz54twtWRquKgisKv7NBkYQj+wp
GvE+5ggP+w8bLjrHmGrTaWi/asar7y8bMs23epOZxP3oPQrJC2yY0mku4JWn
oX+d8t/di2xQYFiPflGgIeVb6/siwvbCi+ublGi47vZCwzLh3D8PdTJVaOgt
bFUjusSGGTbDMXgbDa2FXF5zE1Z5U5NwdTsNa6ad3lYT652yzTqM1WmYKHBu
woRw4eOfnKoaRP7LG8VLiPfP3QpQEttFw04fOZsZIr/tZ0mn/2rSMJM6O7hE
5O+6vzSkH2moSYoN6SfqK1YxflWhTUNzZYGgIKL+v+IT35/q0nB96fjE+jk2
aHB5SXrtoeFw/OOqi0T/bn4VP3h2Lw03ZrPkIon+lrfnuevp0zDA3wMiif4v
v9TPkttPw3uGQuRLxP543HNbO3WQ6H83jIQQ+1l1SXhX+yEaOjyW5R8m9pvL
NO1injENmVHN4StG2XCX9L7J5SgNX0rYxzQT5+U1v8PCseM03HDS1OkKcZ54
Z/gYGmY01G0pcRojzpvvK3X/fxY03OlJbz7axoaG52+Lh08R52miM9yoiQ38
wee/1ljRsP8W3Vmyng2BJyP3+p6jobjJsg1WsOHB39+86y7TkIdlWf2YmJcd
xulSm+xpuE192KWDmKfR9JM7qA40LHPOuz0Uywa1k9WXdrnQ8KfnJuUroWwY
rAzoOu9OQxnTh7Sr14h8xDSnnDxouPY0xU3Hng1KV3+s9LhDQy0L8bRxW6I+
GVONx/doGM/BE158nA107y0J1UE0NFAzLLBRJ+ZxoL20JYSGm6tv/Z5SZMMt
Va/u3lAaHu2KUzpII+Z7bGzV98c0fBKwPTvmPzZc08+9IhFHw3m1Y+VzH4eg
VlB3l1M2DScH/XcEaw3BZdu5I7dzCW/L+azDHALRqhR7/3waqqsLfSr7bwjs
HNYmPX1Jw+qS0CSOiUEQ6u7hb6mkoXM114cfIYNwMup8L6Wdht8vrvD27ByA
OXKw47sfNMznPc35XKMfthQfiK6coaGFoUTnnEw/7D6wpub5LA0V+eI/KfP3
wyNnP2HX3zR0s/6m5Nj7AZTqvHLFOOSwqMzoJcPpA1y0cf1+XFAOXT3xs0Lq
exhJPX2hnyWH7ioBns8F+qBdfuuZYTs5dGlX3eT6pRtObV3rPHNRDkWXDR7w
tXbDD/XP97ivyOEsS8o0K78bhHZHpFEc5TCKz0xR5XY3GJktTp+/IYdB3HBN
RbQb2nxqPWYC5PDc2KGtRgZd0Mo2ieXOlMP69YWUh687wOIzK2dDthzG7/3a
KpjVAVMTvNWUXDlcJtdkJj/uAIH5ok96hcTzhzjsSLYdYLh2MyuwXA4XyOvF
zNZ0QIv6p9INrXJYdm9HcPvxt9Ac6txH+S6H39Ho09Kqduj8VKDNmpbDvUc5
023G2+C92mz61p9yqLJJj7Gypg3GBh3cdeblcPjbkwK2UxtwMexJp//JYYoy
p87C+1bYXmN3KUaQjmvUB7ULc1og4acl53pFOvIk6idqeRLfS3vi7MSV6Tjs
Hm44YtkE2ZGDndJb6WhjpsJVoNkEZXAyWVGNjub2XFXDi43wLuCEgSHSMbX+
5s7zNxphDdk0zM+IjveaJhYWPBrgmom+HMdVOn7+l5BelVEPix/Dno460PH+
4uv2ytB68Lg6LNrgRMeKfVJKw6714B/kyhvqSsfZNY0zgbvrIeZN8gjJg46h
G1yfebHr4NVO7li9YDp+1Weq12+qg1WUMuHgTDq+KZIOdM6shcD8lf6O2XTM
m/PRuRdeC+u0jTmP5dJxHYg3FnjWgrjl+HepQjqajGc9sTtSC8xIsZa8MjpO
OAd1beWsBaO1zj49TXSsT1d/LmleAxG/WIubx+mYP7STX022Grx/PTzJmKRj
A++EGUW4Ghx+/alQm6Ijn/LAuMJyFez/9drTeIaOagonVsT0VAHHr5Mr/f7S
8dRon0STbxXYzt4XnuVnoOehqA0nZ16B+s9vtEYGA59fn7utmVcJsj+N/d+x
GMj3ze5wWUglbPhZPDGqwMAT5/zen79SCVMz3ln/VBiord3QeZhVCQkzm7du
28nAmJ/KS8bpFcA/cwCe7mdgqOj3g2+zy+H9j0xTl0sM/CCpZlPVWwq7mVKv
T11hYET17m93K0ohxyZ0q8FVBk5P/Dxol1QK9wZd1ks5MzBq9Hx/pH0pKDXD
24abDNx9VvLmPf5S8H3+dv/mQAYaCbgp2emWgKrFL53aNAbWlE/rCLQVQXyk
TX5WBgP9JL7xZxcVwZruPpnILAYm8P5ovpBQBB8NKjgv5zEwUZwUqeZSBPfV
fSpFS4j8eb9Wi0sWweiGjRoXGhh4fqm3YsD5JYQ27VQW/sJA+Q+kN4aqL8DJ
ROTPzFcGUr2cnO+QXoDpwLfKzgkG8mdOXmgWegFi32MPhH1nEP8L/XkJE4UQ
J8x5Xmyege4tE3KmCYWQcfx13BZeJoadlbv8Yn0h1I0dWMsgMzH74/mZyBUF
kGIv27ValomCfHPVZr/ywe/PctQklYm7XP7IK37Kh/2rc2hZDCbu/Fp+Reh1
PnTKC+sqqTDxWuj34nTvfGA797ht12aiQYkodTt/Pvzhtvyic4qJdC9rvXnJ
PNBdPhR5wIqJ6q+MG1KF8iDot67+UWsmTmg6plzmzAOZb/QMO1smVoXwCymM
5IJ+z9yVB/ZM3LvTJcDneS7EpATND3owcd1QCJeHSi5oHqjgdYtnYmjYNd/B
ozngq5f70iuRiSGfppPT9HOgQyvJ9v4zJtZGV664tysHzqn6Nz5NZeJv2zsM
O1IOBEseCarPZWL1Ng3rpqFsYE9/ExGpZqJqF/UHSzQb7jzeTMn5yETzVyIc
KVGZoN+8afTxKBMbT1yN0vPOBGGuTYm3PjOxz+Og9dKVTEi4JC5tMM7E2ebC
dem7M6FaV3Tz52kmrhj4l7Z6JgO4fgqJbOJk4eoFB2qDYQbcOcizwleShTVW
ijddRdNB35u79rI0CxWbSLEpPOkgXMrlZUpioUlya+TEdBokyHJykqgsLBD4
Xp7anAbVS0uL5fIsPLChYwuPVxpwZczP/tzFwpQvEyObfqaC56qJUUtzFo6d
OOSjzU6BvGyeoCsWLGzXdCzY2ZoCH49s2X77FAstF9Mb9pWlgM6zQ36x1iyc
edK4JyYiBVZoF7EGLrKQ8t2tRfJwCvjd8HY5cYuFRfSBmIA3z+HBlOQq01gW
vi2mhRyuTob4LhOGPpuFtQWUwJnRJDB66fZR9SORz3BM0Md3SbAc+TSSNMrC
8dW5LoONSXDCcnLl0hgLtVtfV/7JSYIN416j2d9ZSP7lm1d2Kwnu/ct7IvKP
hRIkRS3bjUlwkb5OaHizPJYmefAtHEuE7R4ts9fM5FEvc3Dhylw8mOcc5zY7
KY+XngXKMCbj4Q770zpNS3nsjd3KnhqOh2Zcluc5I4/vwlkLwS3xcIpT4XyI
nTzuvHA/c/uzePC9EzKQcl0eHXKzN4ociYf3d43r+sLlcWP8rERT+VO46dcX
rtElj/fPRv+o2BMH/AyPEL0eeYwrjtJ/sDUOwptl/Yx75TGkkEl1lImDPEEn
twv98nhLj0Q6uxwLX8METkaPyuOZTvo+jpexcCRBR2ZhVh677L5m8TBiQaE0
M6NERAHr/zIvdG9+Ah8nPavUjyrgXKhvnSkjGqYatr6ROa6AOg/krvtvjoY/
yWOtq80UkPI79EaDYDQInT7QP2ChgH9aruudm40i/ufF5j3PKeDhqcOjGyui
ILI8m9norICL20tiTxlHgVHQQPiJRwr49NWndjnvSKhQUL98460CjvXpLZVy
R4CyO6/l/U4FNAy+Y7dnPhySmzsPxnUroIqTw+n34+EQdP6yUk2fArpxqtZt
6QgHi8TEWf6PCrh9VN1KLD4clv4TuBE1o4ABEeeq6jActtyerunnU8QLjzwP
vL33GCzO1bUv0BVRsveSagMjDEy2smf/MRTxbUzpm3jpMNjH9VeMh6WIivO6
QZ6iYbA1lnV6jYIi8m+AiCPcYcDX82h6k4oiPoz49s30wyMo2H1GWFNDEWvl
a74mBj4CfgqHiYeBImoadH8zmn0IhSPq73guKuL6xJzhxJ5QqG/fGpedooi3
lj902i+FwPj4rbRr44o41r58KjAwCOqdDK8NU5XQoF949ERjANzM/KL4n7US
xpQf+KZ/2Q+0179l2T5TwtvdlXvOK/lAzRecNPuohO+shIwWzbzBZGL520Wy
Mu4bN+haFe0F7Lm2TifCo7tN3XuivOAiV3zxzf/f15eUSyJ8d6OOdyBhDc1f
HkC4QM9nUwbhdzYp6i6RXiCSJGgwQdj9Aq1tJNwLesylUs5TlHHvPSf/ykde
YGU7ff8q4aMRAfvvE55yrHa8Tlirz0rQjDBvgDX4Eb7ceyh87qEXqBan9KQQ
jjmtlMMi/FhEacUXwoEGP1fEPPACGWmuie+EP95Z3XWecBazs32esHrGSMJ2
wnU6zjErZZWx5nWndmeIF8xdLVahEv6VmhfAT/jOLf+NCoTXr2y07A32AgE/
s3+qhDu3x29NJiz7dKlxD2HhFcFsLcJ56a05hoRnI5+9FCSs+TLu8RHCBj9v
BA8EeUFjtf1NC8KGa4Vt0wkfbdWyOkeYNGqjdZ3w/wD8a+b+
      "]]}}, {{}, {}, 
    {RGBColor[0, 0, 1], LineBox[CompressedData["
1:eJwVlHk0lQsXxg3djJlypRJneN9jzFXJVOyNoqSuSulGSKgk0ZGMJ1OmopOU
ech8JUMhTpJZFAodY6JBJYrqSojP98dee/3Ws9Zeez97rYfqcO6AEx8PD0/r
cv2/u9xoMt6yRRazVagPjvM5g3KCxxoebVlU0S+lnRZyho8pGz4/05PFODfj
b0wJZ3DKvRDrbCSLWqM6zTfkncG+ihxN3C+LUfFWlxb0nOHw68ssHndZbHjS
r+rp6QyGqiZV7YWyqPflkp/ouDPINjSrn2KsxchxWsfJkZMglMejKKW+DlXn
JIP3/+UCY7uPkMbEemwUCtYsyXcFaj6Xm/6nHCbXdKVxPc5BOif6psf6DRi9
oTU81NcDBK86HG0Tl8d9s+3FbGcmKNQXXJ9cLY9xcvwl8eeYoDX77YmErDw2
23lJ5ngz4YRjiLYVRR6HN/lKdEcxoWZ7tvQ7DXnU3mx2x7+YCczJ9x2/98vj
rVuOnSJzTBjed9p4U5w8am3qMDGN84RySXe1BFkFnKVN6XKGLwDnTahUvpwC
pn2gLVh/ugA19xNnH1AUEO5ZJP/x4wI0H2po7FVSwPNMCttdyAu4iTLH1ugq
4O89/2zO1fSCWeqj6Fv/KGC/9hqP+KteoLdZ+GtckgKGFE00PTO6CPbKXX27
xCgY+UJh8VKbN9i3f+zfKEHBcXcSA7nL7LE0ICVFQT1z/WORb5a5Su3V0J8U
XD2jYF4xt8y7w0fdN1CQohj700/VB+xObx9PUKOgOSgF18T4gO2/ufMfzZa5
XbTJwcYXrJX95CPDKbhi1clkJxl/aDQoFIqPpKAdPxwjGP6gZvnqR/YVCmr7
f0qe2uoPCyx4WnuNgsN/ZHaUHvKHtB5+79l4Cmoa6bX03fSHN5euvjiZT8Ga
6AfYIhsALtzUUJNWCuqHfTmaqswC75DacX4RKkbO/HUmQz8QbHt+v/0iSkXO
vLxAh0kg7CC3veoXo+KGqIJb/BaBINlS3lksRcXZBoOKMIdAKBC+U2a9joor
1FbzfQkPhKHrN1n3lalYm77/wO6eQDDMcJE6sYuK2d83qvoyg0BpOk9knxkV
UzoF4qoCgkDM6P0KXXMqfqrg3cYXEQQDb+1+illQcUyTx6M4JQg8FQ8Ncayo
OH/rWa1vcxDkFkGu1EkqcvfO7H69PhiEH0nr1V9e3me/W59WRzCctK6Mygqn
otjCeIxffzA0/rIeDI2kIkuVv6DtXTCwtLL8TaOpmPuhwChiPhi+FW96/CyO
iqKaksxg5RAYyNy3ozeLigWPE+tZ4SFQEBGxb6KOit21ez3n94SCgKJaensD
Fd/57f2cdyQUHJs6vxY1UbG/3DrB3ikU5PnWXPdopaKeNPPlDCsUYv1yumee
U3HE+FBe5b1Q8HWrP8I7QsWr3UbOGXKXYY/lgoPMIhXbaAPO/nOXQcKRk0rh
oeHwfg2OkHAYcJnefSp8NMy25v6ZvzYMHG782AsrabjKRnevrG4YeHdN6JwS
oyGf7s3EVu8wyLZ4JcaRp6FTJAvC58NgwbyGcwxoSJe7NqYlHgH1Nv7/nTSk
4e1fbP311AgId9XTOG9Mwwn5RJTeEgFSVytywkxpeOije+cOqwhQelrELv6b
hoc3+9ZmZUTAod3pzrz2NDy+5YvtRe1IKNoZKJUTuDxvbxvX2ysKvrHNrK8G
01BUYti4JzIKtg5JZzNDaShYv71tZ2oUVJ//d6tRBA1faHQWn2mMgtb0riOv
r9HwDXtTa/PqK/DuF5G+Np2GtrIfwgcqrsC6olbVmEc01Ll9Lf6KeDTY/rrh
eeExDeum7uz/QUbD7R22j2zqlv0QbVLy3h4NSoPf9qk20XCF0/GxF6ejYaug
3Pknz2i44La1rqspGiwc3Cr5hmjYbZhyRTMkBsJlVu+8OEfDqfdLxz2E2bC5
++fF8AUa9uVV+2fKsmH42lBB/CINT/Flun5gsEFLKFeiko+O8+c/n7przIax
eZ2hWWE6+pU77lBnscFk1Jbpu56OrDAoU/zJhhWFBbcDttOxq5AqvHXuOpSc
utYTbUDHbwXzqltFYsGG9BRIQzquLJj6tE8uFsrS9M/W7KCjlpu6zHODWDhx
/bnOkjkd7RqsD4tfjoV6r5nOwGN0tFSurCuTuQFBhkaLISw6ntR1UTlqEQfZ
e6IjogLpWETVmstyjFvOvz6p68F0NNK6ZyvmEweiLm6KaWF0HHg0ImGSGQe3
YpMsKmPoqNu5wWf3TBwUvv2WOZFGx/HKOYux7JvAvZxleugxHY1tZ34NysfD
7LUvL6zr6Cg4JnJ4UDse1ifp2jg00NEltSB05f54sC/qdD/XQkedvV9XvQmJ
h8/c+cTITjq+sSoW5nyOBx4ly4lHr+kYrOe4yrIxAVTaVlxn8BDI2Pm2+l5K
ErSJ2ygL8xF4JEo7fr4qCc4cul83yU9g0Qz/O+/eJCgcOf6tTIBAnyW/11Kr
k0H9Z81BY3ECFWts5rquJsNmwkfGXp7AlZz5tt9RKdB1+nmRMYXAw8Vueh13
UoBZrGiqSCOwT318YfBZCpTpcS9+IQkMExmZLBVPBa39W/r9NxJIkTIW+ich
FfRYk8mJ2wnU9rwXx6hIg4GGHZoBBgRCVm/w2sE08BNKeWaPBE5XC7pY8qRD
dZzZb8UdBD5p40jl7EkHgzt5thV7CIwnJzua36eDUd9xavdRAvXlaMnH4jJg
92ZurqgPgZ7fp5SLft6GA59Mldf6Eah6KzrpomQm2GRU3SEDCJztDenzU80E
N7HUEoMgAreJNzhq22dC7PgJjnskgW82iTdHP8tczq/p9p5EAh217FiHS7Pg
7T8nLEaTCUwzOpJG6cyCCYmXXZOpBIaKPZTZNJkFi6zK3pWZBIp++L5dSSUb
aNaXRnUKCOzNeCwelpcNLqtX/ZfCIXC93M/C6pIcYLaxLv5bvezvWJ6CRXcO
+AdN/SqvIVCqeE+G9kwOxHztXuyoJ7A+qSpypX4u3HuaJMj7lEC76d6Kmo5c
mAtRknMaJNDcfFDDiz8fon4YG29cWP6/NHPgXFUBSH6o65ZZJPC+jHWAzfsC
SOgHRx4eEi1dLc6USN6B3Jptl7v5SVQWYm2KPXMH6iK2tPiIkFjmY7rqKFEI
s3L0Pc3rSLxh2UMfK7gLl8RvDxbLkWg9xNrd+OourORTcE2UJ9HT1rZ4nWQR
rP6wLtqVRuIZ9lK5pk8RqJdKdUqpkHjVddpB/+9icNrJd9BOl8R1VBPhrsgS
mNAOfLtrG4lWYzHTpnklcF5lkblZn8QKtZ5d/E0lECg+F/eHIYkZVn1l5ryl
kNI/zb2zi0Sdj2dztgWUQtfZ0aOzViQ+CT86LR1yD/BWrcN1LxKRoTnJeFoG
GnNHXmV5k+g1Obvy+X9loGA7bVXhS2JEtVbSG0o5LJK0vYMsEl+3rL9k6l0O
D8tDtMlwEgsF1PwV1SpA66WJaFU8ibynhyx/5T8Ahu7rsKeJJCbb5+unDjyA
P1Mv8g4nkyjFq9P2QLQSfjjmz/BlkOjYJihgeb4SSn8IjZrnk5jWNHoKjKtA
Tbq9fKSSRI6GnOz1BQ7IeTtpfOeQGMerwfm99SGIDP0u+OMRiV2P7Y5PnnsI
41nqGSp1JIKhZfaR9w8hbws7yrOVxJG/mEPnudVAO3jQTnBg+T6bLs19XTUg
G9svqDFPYtD2V/Vnk+qAZdw3oPObxEVKirvR4zp494NbaLhE4jaDu2vt3tVB
iVWPxUF+BgaYXRPy+Ksedm/oTPASYeCzC18z55/Ug29+o9Kj9QxUMpEbaRZq
hKGa4l1m2xl4+INq1Y36ZjByL1p30GBZn5rqY081Qz717oQ1MjBfoHeEK98C
nqEF7LM7GMiKpErq+7fAKrOcPrY5A3cav9pwVu8JGLxMOtVrw0C3pd93VRpb
If3z5QjHAAbm7nHVGBFuh+A+9imzSwzM1iLIFv12cG5K3qURxEAZfcHht+7t
oJ5eKrgQykAR14WK6t52qDnwKuLGVQbORSxuGM7rgOEqzcj6ZAYmnK7xvq/0
HOQj3kZSOAyMUWtYin/zAlIIwytDMww0/LVROiboJXwMdj70eZaByTdrGWdv
v4Sto1cU5uYYyHXFl+Z1L6EjlXt/zRIDhR2/S4vwcoFnjevQAQFF3CwzqesW
yIUz389O1a5RxFlZpwfZAb3wP4feolI=
      "]], LineBox[CompressedData["
1:eJwV0Hk0lI0XB3BLFNlSSmVJhWaeeWaGvJaSe0MqUkqbXeWdt2RfE0oY2bey
RbaEomzhDSUVJVQipSgmFVLCSMuL3/P74557Pufc8z33XpWj7vv+FhIQEHhJ
1f97ydJ94dXGNPwp/3dtQfAreM7UOjboTEM1n49zRTtegqOWpO+kCw0/+awp
+UZ7Cd/1PkUIu9PQ/eTQmk2LX4LMtvTrqt40dLE2+zTY0Q0WNv9NnAikYToR
Psa17IZn5x+GTMbQ8DnvIHHFsQueDuzPFr5B5eX/tpnidoL9J7J8WRkNLWd/
+g8f74RvX0Tvq1bQsFIxQ3rKrBOkZv79uL2ahj/2Zu81le2E3ZKKZOwdGsq+
b1wlkfscOvQ+1i97SsPiq//sswt5Bu3Jvr2q4zTMFT9obERvh66PtwzJCRqK
f+kPaJtpgze6/BKtKWreYD43obkNPr/zOmM0Q8Nh9WDNjiNtIER4rDsyT0Pd
HIUmtctPQOeBs2uWNB1533H2onIr5E85CC5l01HX3/vG6/0tUGyS47xKk462
1Wl6NvQWKMt416WiRcdK0XfDygIt0AB2hWxdOpad1n3uXtoMr2KszXYjHd0i
hFTpi5pBYv2BlCgLOvan8TeKtD8A//2mNAFPOmJ0oWh3eBP8x0vJHfKiY4eD
f3OnUxOEeA6uaPWh4/0Sp2bRbU0QHX9KNPkUHYcO87/KiTZB1uPCD+tC6Pi7
Q93rYfQ9uKcvnL09gY5aOf41/vqNYNy2W46RRMfk+azDG4UbodXqUqzMBTqe
Lk18e+jJXXjhpxHYm0rHmrR7Q71Wd2Go3OHwyWw6Lo8+6FgTfAcWqTbIJtyg
Y2R6Ls/9dT3EVi2M9i6j44t5JZuIq/WwxNBS8HAFHcXkk27MedXDKofR8TXV
1L5if37LStcDI0O+o7KBjoml55sbzOrAQtL3fE8bHW1WRGkU9fwL6dPkf4qj
dFznZd1ixqoB7vQFO2KMjmOVu7bsXVADXtO/7up+o6PhDruimt5q2DXdfM5y
ko7L4mILs7nVIDBttzDqNx0fae6MfT54C47z42T54gTGBS+XUiupggP8KW9B
SQIXSXhw3nCrwJBv9VJKmsDV45uEhxyrQIGvmk5bSmDABxcOf0UVPJ+6o+iw
mkCXh+1uM1GVoDf1dcMTgsDvbWN9/KAKUJuyjH5FEti81kqzzLYClk3d/jLE
IrDDMcPrln4FfJvk3pzfSGCERgr78mw55E8qav2lT+BrOa5mT2g5iE+aQ+4u
Aovu791Up1UGb77fOODnSuAbz8zvm6ZKYBtjTbOjO4G2i5dwtreUQPk/yVpm
ngR6SXyfq84ogYh3fkvX+BLYUhDBqNlaAhrt0NkaROCdca979NTrEFnUuUsx
lsD+31FVhMU1mOIZNSyMJ/BpVsKhHxuugb1SDTGZQOCmweuLlgtdA62US+KP
LhC4I0T/zN7qYngfeuyxRyaBpSl2gVeUi0Hbftro4XUCl2k2xkkJF0Fexj9V
N0sJNIpqyH0+UAgSL3vXZtwkcN8exeV9dwuBZ3ZX0K2SwMqvPLf2wEKI0zvf
uKKOQHPdSDmBX1dhaNnKzSdbCbS2qC9SEbgKyW36mrLDBK5ITDRYtOUK+OyX
+zU5QmCtcW2B4borcKD/a2PXFwIf/pFZXSF2BeTHs81TxgnUfaRksuFVPuTI
Cp6QnyFQk3tcy8snH0qtmnOURBlYW6cxL/FvHsTzLnPmFzLweOeYmUVeHrif
9CMHxBhYoG1pfy86DzSD1evzJBnox3ZUU7XPg9u5UT3r5RgYWB5ePy2aBy2f
zSWJ9QxEHV2+lWMuFHuodS9WY6B8xA9TvV25EPVr7tKYOgONSngqoJsLuxaX
b7hJMNDKpPN9g0wudDFljTU2MjBMvMj3a0QODPj2nNYxZKC+2JNAlkM2/BJ2
GDZyZGA3f+RaBWaB8dzeDPOjDFz/rnPNCqUsiP9pbHrIiYHDdvueNv7JhLVf
6aXOxxkYe/3V2te1mWDa88M9yYOB9G5fL9mNmZBVHD/zLoSBVY1pmSp/XQID
87uip/MY+DpIxzH/WDpEbq+oDbtCuXcPuW1HOrzYWnA87ioDzz7yXL2STAeO
dvST3GsMtF+RHLxvJg0SlA/GP6pg4GCSeepcXBoMTHyVk7vPQJ8788w3D1Ih
NFVRtZzHwM6dzQMdu1PAtF1hKHWIgcKKvrFbNqWArJDCleBPDHTi8+bbVVMg
33WVitkoA0U+mDBNZy/CfeMVip8mGBit+uNIzY2LIDQlI6cgSKKOesIqpWUX
IXTPApFIZRJFV0qrxkwngylX+KGbCol/HnsfDBpKBtl6obAD60hceWdWKKIr
GfLVBAXXqZM4wFwpxqtIhvuzs//dYZLobp/MrXVPBqHSGf7UFhL9UvcXlY0n
wblFX4YcbEn09mlLLZxPhMqyBfHu9iRuvd04rPA9EXgHlXTOOpL4TPt147WB
RDC6ujcq24nEA5I+JqL3E0HE8F+y34VEK8+tT5u5iRAVyPWzDibxiIW1qN6S
REj6przoQDaJp4T5WWWu8XA/Ra/CKZdE+muFc1mW8TCpb2ntk0/iSdGjqaV6
8WAZE1F6oZDEOkVL079E40Fuw1eLFzdJvL2wtXowNw4yjtRl7GkkUTOxM8a9
LxbyuvcTpgMkbot6k/nQOQYsak/ztHkknr1UcET+YAzMZeRmrBsiMei8bHz4
1hiwdhhbOPuZuperffqqfAwsGw0bKhsnkTYxai7yOBoi5isvy81T/7Ks+3s1
EQ0u9CUyg4pMrJTeGXhVMApWSeo86lBm4q+bE9p3xyOhddz2TJ0KE7Pewsho
fySoVRePXVBlolWxds3Zukh4D9hqQjLRNE0s+4JPJOw74B5aqs9E7sfCp6Zj
50EnpIPvb8NEhsrjqI0jEWBbbiVsY8fEM2PfOgZ7IyB04OMSAwcmjpe18nOe
REA7zjEXHGNiemNCgllpBDgKsk4kOjPRXtuUa+cRAZGhif3FAUxUsf3nkNUf
LrwJt2zpTWOix/mXdcoKXJi/9b67IYOJOwRtyTwpLqz/ePJDTiYTW4NO81mC
XHDbFi7AyWHi7nFBXd9P4SAkWr1pspCJIx2NOVkV4UCcX14mXsNEkwUZVXtM
wyEoqjdtczcTo8075l6Eh4E4EZK4vYe6N+Ce01v/MEhrV4uyfM3Eu8PrZSac
w6BS2uf0yT4mqic9HDPaGwYjKVJ2mUNMnNbeVc9VCoOD+UZr//CZeCFJZJFz
fSiw6m+U1smxsEBuPtZg7hzwxs416R1ioau+F9FQEALfWrUer7Vi4bXAdz8N
M0PgV+Hnp4ttWKjlflT8RVIIyBwx7+u3Z2Fuk+2LRWdDwOCl/Mw5DgsLB6U+
vbIKgYw7ZYwnvlS+zGjblEwIWMT3p1lfZGFHqUh3s9hZuMvScwvsZOEKl+Rt
qcLBoHlG1CGui4UvpD43TU0GQWF7156clyw0YdU9c+AFQfwJN40HvSxcE55s
fawpCOyvXOGL81j4eOiB+fKQIJhdLhV4aZKF7kZ9wooCQaB0duJBnxgbFVxF
zm0RDwR7TsvzP3Q2PrunMZJsFAD7tQb48wQbjUnl31baAbBT6Lf8ApKN4h9C
UghaAGhlk0ckWGzs+bjW4JtUAIj1XJxQ2MhGjnr36u9vTsGtbcdkDTazUUI3
4GifzykQVxXYH2LGxrcR2iMS5f5Q/UHv1QIXNhqInMwu3uoHzs+sQmco0/rl
3wbp+IFSfQA56srGgIoTGw6TfhCRfDvsmTsbq2I936xb6QeHt+qxL3mz8Sbf
YEhqwhd+Z+tGsQLZ2NC6a4TM9wUDa53N1tFs3HcqYZmTmC88eq6VU1bMRrfo
tALnUW8YHQ2+7j/KxnaPnNU5LZ7wyGe3/6C6BlrItL4fWOoBQTeG2cudNJDJ
vwQj213BcGknefyqBgrl94Vcu+gMD4ZxzIangfSgnYO/jh6H/V/mvrqs18Tw
V009t805MPDjWZcP5Z8WMsZnKLsI5d0OotwjuafSkHL4SiNuLOUXWy4nte3i
wK3t5xVKKdsPB9r0m3FArkDa7Atl7+ktCoKmHOixXVN8QlUTl7tU7N1pwoGj
xyfiPCk3nJ3rk6T8zfu+dwDlfK3NJ15s44BojBNEUY53DQu3oax9u7inmPK5
+N0trsYcSJXTEBmmvEm5zS/ZkANrVYS+jFMW/yS95BDlm4yu5zOUpWeNSldT
bjHyzVqoponOC09/LNjKgR+etzeqU5YQUDxWixwIDY5eyaIsVv92QRBlqSib
eW3KIZejipCyWu7sExPKr5srx1uBA5UlT8t3U3YTZKbEUzaozUk9SDnQIn2z
JeUn9z2C7CmvujnBW0H50NOtRzmU18ttiukz4MD/ALlm9V0=
      "]]}}, {{}, {}, 
    {RGBColor[0, 1, 0], LineBox[CompressedData["
1:eJwt2nk0Vd/bAHBTJSpFJPNwR9c1pG80eZ5QytCESsoUilSGSJTMmcfMmSNT
xpQhZKxQypRkJpXCNSQU/e671vvPOeuzztnn2fvZz95r7bWOpPmN05ZsLCws
gszL/91toprUlZQEcTZUwvtwniVQ4+x3sCgL4qPxj1NbSyzh20PRH237BPEH
xSNwqNwSLLOcIq3UBDGQ8djMr9kSTCuII/GnBPFLR2vtyoglnBnydWexE0Tt
Dec5PHZawSHakYq3+YL4bey9NneQFQg2NMtdIe3Epz3Z51JcLsPGxyxkXjkh
nOXhXBgosoaJY+eI6gRhrCwWmLrCYwuS2T09KfwiWBx7rkT5+XVIqQyJthcW
xc8qYW4f1tkDZ7D5+RYeMWxK/551sNsBxOtzI6b4xJA1Wz/z2YAD7Fmae71V
UAx9pWPP/zfhAJcsvJXPSojhxt+71tR/O0DNgUfbxxXE0HRptTJrpyM4Tn15
t3pKDKv2LT32M3GEwePW6ooPxLDHf5NLz4wjlG2zk40TFMfDHArRf6WdoHLU
hzdbRBx/fTN7OEx3gprS+KXnEuI4Ub43qlXZCZoNGho/UsQxnfHO+pm2E/TE
C1zcsVccj90etx92dIIlyeqQGENx3NbXXBfU5AT7dnHNPEgQR+X8A5k3rjuD
KbWj9+gWCTyt/+UdV+8tMH377RN9qwT+l3bzl8AY0/b/+nh5JVCopPCVzDTT
FbID/fwS+IktMsiawwVMj90fsROVQInpEkMNRRcwsT4wGScrgV0v1O2cgl3A
OCfrzzctCYyR6tRfOXobjKhuYgH3JXDL+f3f9fpdoVE1f2NsgAR6by0hZU26
gqz+wMKjIAmMDZHrZl92hb/u0PoyTALHOf6WjfC7QXIXu8tSrARC2Lshykk3
GL0X/OFytgQOrP35MNPsBjY9ST5H3kjgaJFudlX1HXDxfjnJzi2Jvjaj/1l0
uoNx1+rY9CZJlCNNbu8YcwcN4v6BT1sk8ZJzyzbdBXfY9qqsvZBXElXUph7a
8N+DXK68p0ZCkqjBV+F9+Nw96I+Idi+lSmLF+sfU+8P34FCqDe+lo5J4P6NW
68NdD+Cq3r6v3lcSPdNmOMYKPeGyUXlgxn1J1DwrRH1d4QmNy0affQIkUV3f
I7aswRPc92Tc0QyRxM9e61xyPnrCXKFibdsDSSyzWtVd/OcJfenHNT5mSGKd
hnaO7GkvyPX3P/6zThIfctW9PbXqBRvIsilvGyRx7yElih6nN1g0tc8UNDHH
O3mOx5jPG8TYdkTYv2H2X+VHYTjVGyLdMjsX30sy52/cweaMN7herz/HOiyJ
a5eC/iwXe4O2/l9zgTVJrLlMd+i384GtFpVJEixS+LVo8w9Fdx/ocXTplWGT
wldEX+sHQT5gHrWgC+ul8In4+iqPLB9w6fipcmWLFLYeWwhf6veBRycHtlSK
SSHpXJTNTR1f+KtTU3kRmN8z1loV3eMH9Rfu/Lp8SAqvjlCIyof94L7tPgUH
dSk8H56UY6LvB7zBzzL9NKWQ4BsvPOTgB5TWgvDCE0yXv5mLL/IDg2MpVqym
UsixFHaBV/4+FBz24M30kELXUEOLJyr+MBeuZRTsJYXVjoMHnx31h//6tz9y
9JHC7rtnl9+d84cXDjn/qflLYfruZWuF2/7wJqXj3FCYFNaX6/uUVfrD+DIh
ZWeKFD7NefLoMwaAUMEbWmi1FBKfRNkrnw8E4+Wom061UrhNnqvrsm0gpGkY
V1+ok8KC6HVrue6BQPk8d5zWJIUpKl08NhmB8B+niMPrNim09Qjc8G4qEE6a
Xy9n65dCq+nOeg2/ILgvwHf41ooUZh3lFy5tDoZdnb9v3f8rhVL/SH0KA8Ew
GNafG7vGbC99Yez1fDDs2Zi1tZxNGn2Fx24dlQyBiT8q/Utc0mh4omJ34N0Q
ODJi7OgqLI0C9xosD+4PBY783LS7B6RxtL6uOqUtDIquhHWFqErjvo5DOi8m
wuAC8eaGZJTG1zeFbt5hCYenyQev1WhIY6XIavB3oXC4FPFe5Z+ONA5FhSX6
nwyHeufFdo+L0nhxbxl1riYcPA+prXm7S6PW4C8WyImAR9oh/oEe0rh/t9Ll
g+URzP2qlzfCSxq3BZt1HH8VAZtsrpOT/aSxyfWPct6XCIiJTDhZHiqNN7+O
ftCVjoT8sbn0n8nSiGa0p2fTI6HHN0PToFYaZ/24BG4XRcFS2PQHozppDDjg
milcHwXCCXsvmDdIo1jzztbuzigwLWi3u/FKGieUajwifkfBj54/8QHt0mie
OI40fAAsFP2f1UPSKNX1Im9LzwOQaeGIILEQkMXt7/oh/hho4blA5WIjoJDv
0PpL1Bi4alBaN8VOwF9VU3UsB2Mgf9hs7ukGAtqNjogGWsSA3O8aPXUeAvr/
FbP/9DQGdhFuC5iKEfDF1UvzUedjYZ/7VGL8AQKensyk9lXFQV+Dxu67qgSc
DFJ/utARB24bH7aZIgE/OHDbyk/GwYsHWqtkDQI2/pA6tX5nPKjmPTZ+pk1A
t8cNEegSD2q9ZpKd5wlIZDOvFz2YAMd29WRtuk3AzpMVr06OJsLp75rUnW4E
nF2o2KPJ8hAupFbkEe8S8D1FIOma2EO4viWpSNWTgGlZoV+1jR5C5OSlSrsA
Ah58fr3L8OND5v4z+7YrnoCXWkaLdn5KAhu+zb8eVhIwxXDIWY49FRxb3G/l
vCAg38exshj+VLjjyVguqyEgt85RezFKKoTOdK69qyegS1ura4huKpS0JnCy
thLw0zK9azY+FVa8KSKWnwm42c9WXmBvGgQuqKvT/xJQvWd0bmtUOmz7Wtcp
sEbAwI1f81yz0yHuE1iwsBDxWgbH1o016ZBVs9+3k52If71/EKO/p0Odv9Kr
29xE5Ht8usNRLQOWRKS1m4WIGNGnkBL3JwMsD7PpmewlYnGAL6dUaCb8VPYY
O7qfiFRNI/bu7ExwkFlz3HWQiP6ckesrGjLBg2flwbpDRNykOs4qsJIJDz/N
9uQdJeLme4MiitZZ0HFt5PzSWSJaceSqGJ16DBjz0jzCmYgrcw+SdXVzQGHl
3ECGCxHVBB0U/trngLjx7NlnrkSMMtlxYzkmB9aIUrqf3YnIuc7BdmwkB6rK
vJWJ94l4zlLGKtQtF/Z0H9lUEcscb/tQJ29tHpD2Dvm1xhORISm5Qfh7HvAn
3WIdTCSipk3EEU++fFiwyF5kSyWi11zhQVvrfChe2Diik01EC0q50uedT0B2
+9uy4XIispnkiBJCC0DExVJhvpKIsu/nLWNrCoC7fzV3XTURG8ZzZt1mCmAy
Qy5Vpo6IsUan1JJPF8JjpfDAm2+I+K+Z47YyTxFI6emZcPYRce5p3OemiSIQ
jPzEqfCHiPUeV3MzW0rAXb23T2WViDFpH3VGJktgfKEn/9A/IhICvDZd5y6F
orNdJ/XYSfif9NYt93RL4Zhoe5wzNwkD6Z1ySd2l4JrdSKkWJmHN0yTam8Wn
MGzYsNIkSsIZoYFnD0XK4Ah3fds7cRLuyW7Y3qVWBrzXa+2GpUmohEsZjuFl
kKdUWc4uS0LBfI7WhwrPoL+m8KjWARJKPLXIG/Z9Dmp2BUJ6qiRk83bUSi55
DtmST34aIQk3Znd8/Tj0HG765IZf0yDhsnQYb/P+ctisldkbrkPCbT7WPi+X
y0G1O+HKxwsk1B732/c5rBIe+cXvGzYmoZnC4OzepkrgUonb9N2UhIlrYoq7
/1RCT0J00YoFCZ9G6V/msK6C66bhS6LXSLiz7VGXgc4LSPnh629xl4Sfve99
nlWqAa/e8Cta95g2nI34fKMGrJoSjyp4ktCfhSfdLr8G5FKKOf/6kHCy5nJ8
FrkWak4P+EcFk9AjSslUg+8lDFbsDqhPJOFsb8gSEuqgPguss5NI+GYxYYP8
0TrIitI6FppCQrW8Uoa/bR1cv2a68XwGCT0Pd1f4PauDNYnggNlcEvZvMYnp
PVEPYv5jARKVJNRXqPpkFt8AbE7T1utfMOfnfIvBQGMDTJgtH/tZTcKTvzNa
RxgNULCfh6u8jjlfZon8oVqNoDqzL/DEGxIaWgidKuZogotnIgPde0n4R7dW
gh7TDKieZGPRR8IF7mlZrdZmIChka2n1k1Bmi0SSAOsr+LmxhktgmIQUi7vv
Xl9/BXeqvwc++UpC0e95+8ZOv4aHhENB/YskLGno50rHFvjmZWXwY4mEj3K7
h/Lvt8B/I0HiKyskDB+zbqtrb4F3ST2lO/6RMOu+dLbIpVZg2WHbf3oDGffo
ctmKxbaBjlNEltlGMg7M5dmwfm2D+M5ndnbcZLSi/ndYQOUt7ApnXRfKQ0aJ
f5vlfgy8BYuNsfQ3O8jI6X7kwaBKOxRdfvH7404yXsw/FP4zsh1Wm0bqJoTJ
iBsGFc5Pt0OMl+wZDgkyClcvP3im/x5e/6lzBwoZk88Q0iIFPwBt+md7mQoZ
P17RTrnc3AEuOrwJjfvI2F/ylLHyrwMac5UtOg+Q8fAAX7rTvk64cNlriYFk
hMmlFoWSTggd3iEpe4yMxslCunJPuqBP9eCPfdpk1K9OuWA73QXkJPOyY7pk
nEoP2F+i0A0vDZ8cu3yKjBds7S9YlXfDXIeaQ7ohGSu/zbl96+oBVcUrB4qN
yCjZFco7IPYRgsJC1r+8SMY4f76Nf60/AkGnN2HAjIxbloSzWth74UzTtQZB
GzLeq1E4RjzyCSqfx20PcyXjmtth3rCN/eBzquWn6x0yptGSXveZ9cPxH38a
rdzJaOZWe9Ciqh/GxEycVL3IaP+yMnK/wwBs8iP2TAWQMeOzpYj6zCCY6JfE
6sQzv3do7bPi9hGQmR67oZJIxgZL6YIGrRFYuM9/lJBERm2DFEcWrxHwr3JZ
+pNKRlaRg8dvz41AiRQY5mWTcelfQNjpz6PAMdsqxFVOxpeCdVq7WsfhXeDq
3K8KMm4rM7J+x/0F4gjyrSNVZNw1/+HEvM4XoJ+LdKuoJWOVkg5n4/svcKb2
XP+VV2QUW3kcMTowATkhX5Je9TDzdWBULWfDd7hJ3uFc2ktGN4+K5ke630G1
7ujxlD4yqr1blHwc9R06FvLWnAfJ+DNuvV6o+CT8NXIwIU0w6zG4ONXmvx9w
XOafhO8vMjbu6mvv0J6CsG2xJy//JuOP6Qqn7YFT8H5JzuPYMhmHfIwMrF5N
welXxkObV5n1+zX1zazqNEQVLG5hrJHxzYqBy/rb09AVHarawUJBjvru3JHi
aThjUZMUw0HBxZ9Ws6niMxCrfeaty3oKzlt9eWGiPwO9u6b/nuekoFJ9YnTC
/Rk4zyp2QXwTBUUWfynv/z4Did/Kgli3UFD7cJN7DTcD+tt1q8Z4KLju+E+H
ZUkGmCTfFcrmo+DHHJnbidoMSPXl1wrkp+CByD8Cc8YMGLV9ctt2BwVbWF/+
PGTPAGn9wznHd1KwYt/W295eDLDYP9CrIEzBha2lEYWRDMiUcuLkE6Vg2lSV
TW0aAyY2blb5JUbBG8pXk/ILGUCefXT5owQFHZ40ytm/YMCV3gOxFVIU9L7T
Hr3hNQNyaruaEwkUvKRys86mgwGTWbaLd0kU7DabEo35zABaKAfJlELB/W9q
433HGGDr9NBATYaClcMKuf9NMuDJhd2+BFkKOv9WPJExw4Bp9ban6+Uo6H+2
XKBtngHyNIvxb/IUvH/RTbxgkQF2vH/5WhUpaNDBK3x0iQHFy1HqT5QoWMtv
fyKe6blhmmPYfxR0Xy2LTfvNAKXXDen2yhRszz1nYfaLATcLjTr09jLztfa4
8tMsA8pi5ln37KfguefpSpumGLB4N0hR8CAFs4qoKysTDFC2lDZbUaVgPlv4
6+whBnP/qArvRwqeJYV+EfzIgAolvZc1ahSsvposceItA1aEfsykalBw6VDL
2xP1DNjP5i3ufYSC3J8ZL8WfMeDOd6ETlkcpmJAXfrQxmwHV70vcNbWY+d0p
RlFPYMDac60Cqg4FP5jPRaQHMsDDz3Xz9AkKFq3afhC9woC6a7wH359i1suZ
0RVNAwawGeTaluhR8DYb+9EbhxjgI93X6nyWWY+DNte6BRjg/3Jv4D9jCnre
6VftvzYDbx5/qBgxpeBhsvOTvr0zwBVm/b3BnIIDPs0+xhwzEHwx/qi/FQVn
s/xaMXoaIlaW1m+7TkEdy4uFt7OnYJ9enoSIHQU3zI+0mltOwXjexX1kB2a9
zz9HVckpULlYf+2gMwX/6je+JcX8hMHaoC5rd+b4dURuStz6ATK+Yun1ocx6
E7F4tEL6Dl0D76vehlOQoXatPab/G9zd493dG0lBsrj41pOR3+D916+cMzHM
/hgXzuuufYVbWsU3hFMo2K/+Qqz20wQ08mgcvFnIjHdc+Utu1jhcv7J45l4x
BW9GdUtNmI/Djrpsu8BS5vrpPq3+WnwcbBw2P0p9TkEWZ45zEYljsLW7h+tt
LQV900RrIxJH4WKCdS/xPbNe6uWt1tUOw/p54TmFDgoOipwi9/oNQ5H2O+4D
XRQcpRcI2p4YBo7VXXC6l4JhmhpceU1DkGvyN9N9mFnP2bnu8sRBWCSEOX5k
MNfHfPrlO9TPIFahm1g7x8wnPcxHoacPDutuani8QEHR4lO5PD598MApgNdl
iYLpth0EtbFPoNjsXSzIQkWT7tpKocJesL3sMmPIQ8XdKnTn+Ns98ODPnh1q
26iYH5NUfEW5B6rCfqnK8FHReW+pguxiN3CVO4SuCFBxqDhH3su5G7I5r9ET
xai4lnT9BZdXF4zlmF3tp1Ox9uXmwMd1HcClKhHVKE/F4YuUZmpIB+zqGKzM
V6QihUeJ39awAzxWjLju/kdF4+vbWYkLH0BM+0yO6EEqzrE/qhDY9QHO/9T6
aqxDxVwR1leiie3g5bGRR/M4FU/3DYR4arZDzvbXe+RPUpEhEfJjy693sHTg
sN8/PSqW+gpdj9J7B9EhQEw1ouJIYklNvehbeC+3+9KIDRVHM1b6Z0ZbwHT3
Zqc5Wyraiwld1UhoAcbeCT/2G1Q8xjYYnXW6BbYejsslOlJRR8hgt2fzGzhp
9HfW2o2KW6uKbHgrX0P7/UaPuSAqHiqdU1trawbjkKRI9lAqam4vapGJaobp
SOdH28OpWEhwH5Q+3wxbkimv9zygojdr+eS1ySY4/jSYx+0hFYMEGBMc25vg
3bB+MvsTKq7st9gX5N8AxhP0ou2FVLQj6wDH+QaY/rG+nlhMxcuREYU7ZRtg
y+/yL5plVOzHiWyejno4vlmUHlxNxWumKq5jhHp4u/dL1fZ3zHyuNqnSxl/C
Bah5S3xPxW+JJ8tP1byEnxqxQ3s6qKjsGLTGiHsJm04eYzfsoaKKlZ9V5ImX
oGP1RCtpkIrmysseRj610Bbp9Ik4Q8W8gl7zB87V0PnlqRp9loqS/Xs/3d5b
DX0qC3m756n4/MCqEP/qC/g66OCu/puKk86Lp839XgAbzU7a7B8zXy2hFjxp
VaDcYHPtIY8MKn0ez49ZqABVgdyejG0yuEJnzyysqwAN6++QxyeDHz7ppJ0L
q4DTPFd4K3fIoPkel+DTshVga2T5/KO4DPbubjhXaFsO6fMmrHwKMqiXvJJv
t455fjiSYiO0SwaLLp0deNz3DArjBzsld8ugy5KJ1umiZ/ACLmYpqMhgamlj
5sqFZ/Ax6Lz2cZRB7zhj/RtVZbCJYBAdcJL53PZLU17QU7ilr0VlsZfB6KZe
zwy7Evg7Gp067iCDvps0G2r0SsDDfmTHm5syaL8mKiSiXAKBoS7rI11kMGJn
4p/4tWJ4+DprTNpDBj87ib0LCS+GlwfYkzXDZLChv0e0vLYINFqP88tGyODH
R7vTk7OK4I1hQvDWKBnMyJZUfBFSBB3Oim6fYmSw+mTMavyFIhgvMjl3NVkG
aVNRF1KWC4GT+II37IkMni3wqToWUgDBpRsCHQuZ8T6ufrtmVgDb1PRYzxXL
oFDFo6vT/xWAkMnkjESZDG6+oXWPZZh5Ho0XfFvyQgYL4MOqrsoTOLnZ6X5P
qwwu//JwqFvKg7hf9L+ikzJYbHHEQ7k0B3x/RV2k/ZRBrsLFQ4ahOeDwa7lG
ZZoZ78yr7tfWOaDzq8lTb04G4/37hEMlc4Dl18UNASsyeHXK/EzCg2y4shDC
u8BFQ56WqCFu/8dgsDDvyLqZhtqO0evyLz8GtQXD7i08NEx7pd4eovkYRBaI
cVQ+5vsWP2P4OB/D+/lqURNhGvaqhzvwB2fB3vkpSguNhiuy7pyWKZlAmtcL
/EinYdRqvVGHVyZsn6/4MS5Pw/oOlcu3rDJhes634J8SDfvNE65byWVC+pzo
7v8O0FCpe1jgzctHwDWnC6k6NLyU+oUhO5MBfYwnBs7XaNiwcweP+I10OCwr
0WR6g4YpmiJ/Np1PZ54vIndr29PwGJXDX+xwOvgNOvNJONHQsTfsaJ5wOii2
wYc3d2g4G7o4Y96SBv6PP+iIBtOQa7Kle6NcGuwx/qXemMuMX0Jf1d+UCmnx
l0sL8ml4u1pyj+FiCmzq/iQVX0BDy9iwdtWuFBjVrmG9XkLDg42dZ/LCUyBk
7/3aHZU0nGE9ELPKnQLj23fuv/qGhk7+VHPqpmSIbD2wi/cbDUu5Ce5f5R/C
TX3+5bnvNLwWNaZB3f4QDAamajt/0PDE4aXwqKVEEJxJ1o2eoSFf6CPP4bpE
SOFltRb8TcPDpt1X9M4kQr5hU4rYellcizhSesEvAZq/6m6mEWQxb0uPlRBr
PGTbkbq4SbI4F9GQNvo9DgKW1xJ+kmUx7Mzm8ZaOONDhLqIU0GRx+u5g9Mqj
OOiU49VQVJJF46VN195rxcGwU4+rsposJvWecCM+jIVldpNv6qay+HkhYXCj
fgxorJ2K1zWXxb0a1asvVWMgdElD66wFsz9H2zRjqDEgNSWTb3NFFjfKybsl
/IsGrZ7FGxF2skimD55kfxIND7NDfw96yKJqYd0IeXM0qOrWrHdNk8Xt+lv/
XByJAn/N4ufeGbLo/FRZwbk9CjoOPboSkimL/ff4LLOqo8BqT2BLao4schrv
IBxPiIIw8TOhr4plcWHHJ6tn+lEwPDvFz18vixbfQtKa30aCV4wosWhUFv13
D6n1t0eAVpvIeMy4LEKEf05+XQTwsolk3J2QxWjjDeMhpRGQfk1IUntSFisy
1l31io2Aeo0dohOzspjOUvxdwTQC2Oa38ouw0lHzxYeIkIVw8DrBsc5fnI5V
2o9tkynhoOXL3nhdko4GE+wNhiLhwFvF5m0gTce56j8rUlvDIZ3EyipNpuMv
yzO3q2fCoH519W+1HB1T0kaj5sLCgC3/98L8QTpynm7af7o7FDw5f4ybXKDj
enaLlUGHECgp5Ai9YUzH5B15JuVGITB6Rkz5nikdpR8qf83XCAH1zFMByRZ0
5P47+HRCIATWqZXTB2zp2KO+rz7uRTAEuPk6n79LxxK2pxadm4IhYlqc0yCZ
ji0vM5ZjawOhPnpvsUUqHXccbG4OyAuEuQN652+m0zHzwBJrZGwg6AX55Udl
Mce7r+bglxuBwE+ZOtlRQEe5Zt89UZKBEG9WGX+ilo6nhx/rzPoHQFqXPk1r
mI4OI26Bo1b+cPK56+ieUWZ+hCa/hZzxh7X41HjpcTr6rnTFHDviD+dNfm5Y
/UrH3amdgbNEf9g+6T1eOEPHgWd891gn7oPfv5Ik/n90NIystOKxvg+2Mtu2
jojKYXFUuvk5Vz9Q9ni7cMtIDl+sBq7La/KBC0WG7EYX5fDQ/dyM7RU+4DX8
ZZuqiRyaD/3LDMj3gTZck+O4JIfKi+mp0VE+YMoqbx1uI4enLybNK5n7gL9X
+ED2bTl8LUPzcGf1gT4fveZPsXJ4kBQUP6fpDXcCPsXu75LD0Q02Dq2znsBF
8wjX7JHDozbetM4JT4htIwXo9cohXS7q+/hnTyjhuel6tV8ObQncGqRmT/ge
veVi4rgcLl8NlRJP9IQz6epSfxbkmPvpt0NDRzxBvupJfiW/PLJLUW20Mj1g
9Kdn3d6z8qi2WF/K/9sdpt/sfi1lKI9tSX/35Ey4w3LW13fcRvKY9ruv9liP
O2w10+0fMJZHPoHnecVl7qDaLfjb00oebewU7ISd3CG+ulC2xUkek3V6N2gv
3oWToQOx5x/IY2vVUPNL1rtQI7/3utsHefQ1u1Vgu88NdrmvNwnplMf0NKJm
Dc0Nsto6T6R0y6O0ycZnwqJuEGp9XbHhkzzeNG0aW/7nCsYZGQtco/JYcItL
Qq3ZFVYFtrglzMljrGKT4W99VxC7N9vQv1EBbYTz2NNdb4OxVfP7PzIKaPDd
WSBy/Bbo7x5e+EdTQJHZKtvnvbfgGNuKIAddAe3DBcO+tN2C3cl0s03yCrjO
y4b/Qtkt2NjzYFZESQGThVhViv1uwdPDl3hV9yug1o48xwWZW8BFZNH30FbA
846VLZm3naFsbO9HDlsFrIq3mOeSc4JX73enFGYz46k857zC7QiTk3dzb00q
4G2xgHlWPnt4dfP4rRGyIi4+v3dpee063HnyTUHAQhFVDx4ru3vCFtT4PtCv
ZCqiZF/Xo4kma2j4hj+NRhVxxfHONqPSy6D/Y23KlrALE2/wCF+psIThxfbO
m0wLj9x1kmfali2t4g7TYiVv3i+WW4LPTnXfYKYbTZcDfJl+qnlfJJ9pK4fy
9RnPLYH/EY/2D6bxVjF1sMwSei5IZFsTd+Gr0Bv8BiWWYH5lNsSe6ZR9nkEi
TE871jveZpry7AbbeLElrA+ygACmTbnG5hyY3lOR3ZPNdHGO1GB4kSXE8Cuu
+8a0xkfWwbYCS5CSZPsxw/T8pUGzB0wXyHa+/810mEvshBHTzepODzeQmO3J
KXOTTyxh0b5Cicy0kgQH/0amve4G7pRnevV6Y8r7fEvYEmD0bw/TY5JmtDim
SamrLUeYlkqXP0xmuiTvXdFxpr2szDqn8yxB9XlKzBmmBe7eMH/GdEu93R1j
pq+OGM7dZfrsu0PmVkzP////Fv8DfDVtYA==
      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{6., 0},
  PlotRange->{All, All},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.414582666104485*^9, 3.4145831138761263`*^9, 3.4145889788149567`*^9, 
   3.415369462033064*^9, 3.4156286396810303`*^9, {3.41562867657683*^9, 
   3.415628696859037*^9}, 3.415628861273712*^9, {3.41563000212436*^9, 
   3.415630015237604*^9}, 3.415630055242258*^9, 3.415638024209326*^9, {
   3.415638069366228*^9, 3.415638088314986*^9}, 3.415648105690448*^9, 
   3.415648170690754*^9, 3.415648629512608*^9, 3.4156487957942467`*^9, 
   3.4156512998857393`*^9, 3.415651342104991*^9, 3.415651521430991*^9, {
   3.4156517317110357`*^9, 3.415651750576062*^9}, 3.415699395155546*^9, 
   3.4156995963713207`*^9, 3.415699738745739*^9, 3.4157436372348003`*^9, 
   3.415747014285604*^9, 3.415754238571681*^9, 3.4157555270376577`*^9, 
   3.415802436706757*^9, 3.415813592644659*^9}]
}, Open  ]],

Cell["\<\
Note: The actual resonance shape (blue) appears to drop off faster than the \
Lorentzian fit (red).  Let's try a Gaussian form.\
\>", "Text",
 CellChangeTimes->{{3.415628885531952*^9, 3.415628902207514*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Gaussian fit", "Subsubsection"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Enmin", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "-", 
       RowBox[{"1.", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}]}], ";", 
  RowBox[{"Enmax", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "+", 
       RowBox[{"1.", "\[CapitalGamma]value"}]}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], " ", "/.", "sub"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Npts", "=", "100"}], ";", 
  RowBox[{"dEn", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"Enmax", "-", "Enmin"}], ")"}], "/", "Npts"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Log\[Sigma]Thtable", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Enmin", "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"i", "-", "1"}], ")"}], " ", "dEn"}]}], ",", 
       RowBox[{"Log", "[", 
        RowBox[{"\[Sigma]Th", "[", 
         RowBox[{
          RowBox[{"Enmin", "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"i", "-", "1"}], ")"}], " ", "dEn"}]}], ",", "Tvalue"}], 
         "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "Npts"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fitfun", "=", 
   RowBox[{"Fit", "[", 
    RowBox[{"Log\[Sigma]Thtable", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "x", ",", 
       RowBox[{"x", "^", "2"}]}], "}"}], ",", "x"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dE2fit", "=", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"-", "1"}], "/", 
      RowBox[{"fitfun", "[", 
       RowBox[{"[", "3", "]"}], "]"}]}], "/.", "sub"}], "/.", 
    RowBox[{"x", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"E0fit", "=", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"fitfun", "[", 
       RowBox[{"[", "2", "]"}], "]"}], " ", 
      RowBox[{"dE2fit", "/", 
       RowBox[{"(", "2", " ", ")"}]}]}], "/.", "sub"}], "/.", 
    RowBox[{"x", "\[Rule]", "1"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"\[Sigma]ThfitNorm", "=", 
  RowBox[{"Exp", "[", 
   RowBox[{
    RowBox[{"fitfun", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "+", 
    RowBox[{
     RowBox[{"E0fit", "^", "2"}], "/", "dE2fit"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Sigma]ThfitNorm", "=", 
   RowBox[{
    RowBox[{"2", 
     RowBox[{
      RowBox[{"\[Sigma]0", " ", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"\[Pi]", " ", "dE2fit"}], "]"}]}], "/", 
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{"Erf", "[", 
         RowBox[{"E0", "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"\[Alpha]", "^", "2"}], 
            RowBox[{"Sqrt", "[", "dE2fit", "]"}]}], ")"}]}], "]"}]}], 
       ")"}]}]}], "/.", "sub"}]}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"analytic", " ", "normalization"}], " ", 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Sigma]Thfit", "[", "En_", "]"}], "=", 
   RowBox[{
    RowBox[{"\[Sigma]ThfitNorm", " ", 
     RowBox[{"Exp", "[", 
      RowBox[{
       RowBox[{"-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"En", "-", 
           RowBox[{"E0", "/", 
            RowBox[{"\[Alpha]", "^", "2"}]}]}], ")"}], "^", "2"}]}], "/", 
       "dE2fit"}], "]"}]}], "/.", "sub"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Enmin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E0value", "-", 
        RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
      RowBox[{"\[Alpha]", "^", "2"}]}], "/.", "sub"}], "]"}]}], ";", 
  RowBox[{"Enmax", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E0value", "+", 
       RowBox[{"10", " ", "\[CapitalGamma]value"}]}], ")"}], "/", 
     RowBox[{"\[Alpha]", "^", "2"}]}], " ", "/.", "sub"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p1", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Sigma]Thfit", "[", "En", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", " ", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"1", ",", "0", ",", "0"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p2", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"\[Sigma]Th", "[", 
      RowBox[{"En", ",", "Tvalue"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"0", ",", "0", ",", "1"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"p3", "=", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Sigma0", "[", "En", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"En", ",", "Enmin", ",", "Enmax"}], "}"}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"RGBColor", "[", 
       RowBox[{"0", ",", "1", ",", "0"}], "]"}]}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"p1", ",", "p2", ",", "p3"}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.414583139164957*^9, 3.414583141747271*^9}, {
   3.415628915502582*^9, 3.415629156628426*^9}, {3.415629245870842*^9, 
   3.415629282940509*^9}, {3.4156297644212313`*^9, 3.415629770644369*^9}, {
   3.415638136219746*^9, 3.415638141778121*^9}, {3.4156381955569077`*^9, 
   3.4156382068421917`*^9}, {3.415639332701733*^9, 3.415639339770822*^9}, {
   3.415639375621611*^9, 3.41563938126714*^9}, {3.415639442749649*^9, 
   3.415639446043022*^9}, {3.4156395161518507`*^9, 3.41563959718076*^9}, {
   3.415648126326679*^9, 3.41564812684352*^9}, {3.4156486479014053`*^9, 
   3.4156486595790577`*^9}, {3.415648824666739*^9, 3.415648883736212*^9}, 
   3.415651350732438*^9, {3.415651550779072*^9, 3.415651560080173*^9}, {
   3.4156517625629463`*^9, 3.415651765199005*^9}, {3.415699813614644*^9, 
   3.415699834452005*^9}, {3.4157016473739223`*^9, 3.4157017200988703`*^9}, {
   3.415701857228829*^9, 3.415701861442897*^9}, {3.4157019064641743`*^9, 
   3.415701923991572*^9}, {3.4157019572242203`*^9, 3.4157019609258423`*^9}, {
   3.415702080234495*^9, 3.415702085713682*^9}, {3.415702193970887*^9, 
   3.415702204144465*^9}, {3.41570472569508*^9, 3.415704749626606*^9}, {
   3.415754265007907*^9, 3.4157542875577717`*^9}, {3.4158184842782393`*^9, 
   3.415818489283526*^9}}],

Cell[BoxData["3.345955228421154`"], "Output",
 CellChangeTimes->{
  3.414582668103998*^9, 3.414583144378921*^9, 3.4145889796702967`*^9, 
   3.415369463788773*^9, 3.415629168178162*^9, {3.4156292504671307`*^9, 
   3.4156292845590076`*^9}, 3.415629773056148*^9, 3.4156300790364523`*^9, 
   3.415638147417362*^9, 3.415638209138191*^9, 3.415639346974599*^9, 
   3.4156393843489037`*^9, {3.4156394230956306`*^9, 3.415639448362474*^9}, 
   3.415639532861706*^9, 3.415639599195346*^9, {3.4156481185773706`*^9, 
   3.415648128158286*^9}, 3.415648176480714*^9, 3.415648662827572*^9, {
   3.4156488048547153`*^9, 3.415648852115075*^9}, 3.415648887626216*^9, 
   3.415651354115673*^9, 3.415651527599185*^9, 3.415651563355934*^9, 
   3.415651768450889*^9, 3.4156994228912983`*^9, 3.415699858165496*^9, 
   3.4157017275315437`*^9, 3.415701862559807*^9, {3.415701911554407*^9, 
   3.415701926972662*^9}, 3.415701963324387*^9, 3.4157020972979727`*^9, 
   3.415704754321129*^9, 3.415743641887897*^9, 3.415747018647704*^9, 
   3.4157472059374638`*^9, {3.4157542442809668`*^9, 3.4157542889412622`*^9}, 
   3.4157555324701557`*^9, 3.415802440572233*^9, 3.415813596164612*^9, 
   3.4158140461829157`*^9, 3.415814280836804*^9, 3.415814637532493*^9, 
   3.415814694835346*^9, 3.4158177430309677`*^9, {3.415818468565727*^9, 
   3.415818490534966*^9}, 3.41582720286579*^9}],

Cell[BoxData["4.051489729158674`"], "Output",
 CellChangeTimes->{
  3.414582668103998*^9, 3.414583144378921*^9, 3.4145889796702967`*^9, 
   3.415369463788773*^9, 3.415629168178162*^9, {3.4156292504671307`*^9, 
   3.4156292845590076`*^9}, 3.415629773056148*^9, 3.4156300790364523`*^9, 
   3.415638147417362*^9, 3.415638209138191*^9, 3.415639346974599*^9, 
   3.4156393843489037`*^9, {3.4156394230956306`*^9, 3.415639448362474*^9}, 
   3.415639532861706*^9, 3.415639599195346*^9, {3.4156481185773706`*^9, 
   3.415648128158286*^9}, 3.415648176480714*^9, 3.415648662827572*^9, {
   3.4156488048547153`*^9, 3.415648852115075*^9}, 3.415648887626216*^9, 
   3.415651354115673*^9, 3.415651527599185*^9, 3.415651563355934*^9, 
   3.415651768450889*^9, 3.4156994228912983`*^9, 3.415699858165496*^9, 
   3.4157017275315437`*^9, 3.415701862559807*^9, {3.415701911554407*^9, 
   3.415701926972662*^9}, 3.415701963324387*^9, 3.4157020972979727`*^9, 
   3.415704754321129*^9, 3.415743641887897*^9, 3.415747018647704*^9, 
   3.4157472059374638`*^9, {3.4157542442809668`*^9, 3.4157542889412622`*^9}, 
   3.4157555324701557`*^9, 3.415802440572233*^9, 3.415813596164612*^9, 
   3.4158140461829157`*^9, 3.415814280836804*^9, 3.415814637532493*^9, 
   3.415814694835346*^9, 3.4158177430309677`*^9, {3.415818468565727*^9, 
   3.415818490534966*^9}, 3.4158272028991337`*^9}],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    {RGBColor[1, 0, 0], LineBox[CompressedData["
1:eJwU13c8ld8fAHBZyb6TSl8ripKZrPucT4UUycgoFRmJComyk5CMooQ0UMjO
ClkhW1HZowjZq0iK8ju/v7zer3vd55zP+YzziFo7G9kxMzExabExMf3/77sc
LW+vFH7QSbI0rnVjV7ffFFEvmsoPp3yZ8wW4OdSZbTtITdhvdl6NYKrapK66
2SZNII0fFngyrr4I4lFPvenfWZjBD/Kv9+cF2VLU/c3LZOdf8oOIWgJ/dYGI
uhKz/KhNKT9YyZR9MIxWUX9kInT0aBs/bJq8cOtYnYP6eJ5PSfIaP9grive+
9fyofnXns/k/kiQoaEGjZmdlNDbOeoT7GJBgT2x13zF0U2PByShFzpAE/27B
7jeegRo9C7sqR7Hfxgg2/W4I0khb+jKvZ0yCjcm/pVuTQzR01jSNt5mSICDe
+Kr5rbsaoVykLZUWJOh4XLeuWPlQg1c6/cX6ORKcHSskRTHlaSxnBlQV2JOA
0fN+bRtfvsYXmVO99udJMGOkOErsKdB4Kc/H9cGBBMG3/fP1nr3SMFB1d0q8
iNezbS4r8kaZxgL0f3Rxxutl+nb0h0e1RpTO/r37XUlw91Z5idyeeg0Fgxdx
JHcStPD1NN8Xe6fxyYxn7es1EviJbVgJ6/io8cyiSuaIJwnut4bvPUPq0Lhi
ecUy34sEPtY/Jn9Ed2lo2khGbfEhwd+YHXIj9r0aVPvemgBf/PtzJRvb7g9o
FN2fmvrkh9er1Bh7dGRQ40TVKlncnwSXzkTyFjd91Vib4Va/coMEL0b0PG92
j2g83fyfzdsAEkTd5xDhFx3TAG3ZMEogCRpPrt7mzJ7QGHGFApsgEkTvMWYa
CJjWCE4w7C8IxvENjn9fkTGnsfOdNQtrCAnYCoo/jNO/azSvXNl1/DYJZo/n
7hTkWtS4JBFknBxKAlJYXXz4vp8alWUR9nphJHii8fh5i+2yBr9RjPcS9iTP
ooLfs18a+X4vnmtGkEDFj0JuN/qjwULLLZ7BtrDc7H2mYVXjeGZJS/QdElxo
HRmWM/ir8au7aXH0LgmKvG4/DTFkYug4fdoYEYn3W7ptvGzfBkY8a//WvVH4
/Je9qBelmRka8jMHg+7h3y+q9V3ex8qIaFgyk7lPgi1Po9H+42yML6f/XujE
br75d2TqOjvjRihvtMQDEmww3tiZwbGJ8UlEIO09tgezHGeQPSdje7FwuXsM
CQj/ws2HerkY9SNyo3WxJEjkj5K895uXIeClunIpjgRXDUXcxEz4Gef5D3DT
H5KAd7pWdYsqicHFMFY6F08ClrYW97OSFMapdgsd3kckyHsv+ZtFnsrIcbA9
VYT9a/O+zunDNIZBjHsg+xMSUIWrkmLSBRhJu/3icrAtzB9dWfkuyPhRE5xl
+pQElucrHb7obmE8mI/tSE0gwZL2eNT8vm2MsaDECf1EElR+7yDT3/3HUBFK
X1vGPjiepjugJsLo0yndfugZCWQ4L9LX2cUYuwZrVOax7V5oRwavizF83Fv0
Yp+T4FFwZJ4K23aG8LMBt/FkHL/ARxtH1SQZLirfQu6mkCD3wt5cs4s7GDWt
s4/3pZKAI0C3Yzp7J8N29V/trRckuBd536nXdhcjrOTSNCONBM6nz64dWdvN
yHcfIC1hS+ifi2TN3MPYsPD6jFUGCZ4qhFxulZJn7MzeGUTPJMGBqYupO3cq
MAwcYzPfYQdvKbn7V0mR8XTUbUUlmwTvzewUrXz3MuqSRv6bx5be4n3yZaky
Y+aMkVZKDglOzl43/LtJhaHWK3uPlEuC1n2sF1OH1RhnY56WNGAzCwnVyzhp
MEKMeQZ980jwQfXDBA8/weh+P7VrKp8EDy+/nZh/D4x/oSeMEgtwvlny5EPt
foakTqOHaSEJdPTrnzo1HmC41aTU1bzC+QAzpr5/NRmPrlNnPIpIcMJygWOj
nDajRuMmWbaYBNecejqqrhxi8BdbWT4qIYHQ+arf8QpHGPvc2oIMX+N4p5A+
N7zUZZyRJ7I2lpIg3mL2JDpwlJGdufX3lTJ8XgmJrN9MDBgd50OFpctJoNy0
dVv+AUPGmsRvrSHs+5ns1lkaRowjiV339Cpx/RYXX6BYHGdcPq31mvkN7l8b
nn3lvWHCiNtSOFiCnV5xZdmk2JQxHh21W6KaBFCxymNpfoLBa8Rk3I+tX0Zi
4a47yVDmc/aMqsH1PKDS2K91ihF4W6/+71sS3NpqTjYIsWRkapfNFNSSwI1p
0GBnqhXjE4s0xbEOnw9nd2dM9FmGqN9Gq656EsRN5769+tCGoaN+NTi8gQQ2
EUFuFnm2DJeV0awDjSQw76qdfNBjx6h0rfmd00SCn3767jbm5xnfZOVF7JpJ
IPBOMKYy34HBPZugvbWFBHqL1yqYRC4wLOx97996R4IvJwu2Ht/nxAjYPvOa
8R7nc2gw04MRZ0bG15NDi9iHxOT677S7MFYsVGSs2khQ47a071a5K0N48wtj
+gfcz2TULcl1VxjaXTSvd9hLTMVxC/1ujGiDxXqVTyTQ7io/V7bvGqOMx3p2
Drvmdtk5PT8PxnDzB0pKOwnKjoed8eryZMhr5ViROkkgtpPX/8obH0ZkQ4FS
Lfa9+WDmkiN+jIXDrzmudeF8sDbq/jR1nZGjX5s70I3Ps5ZVZqnoBoPnY1Pg
3R4SbNbrvBWfHsC4aNxmfqAXn79RZeO/tJsMafM+prQ+HB9aqWT6hyDG7b7B
jpP9JNDoPjNr/ieYMXHqWxrPAAkeRPA3xiuGMF6cXTC48hnXp83ANf6eUAb7
6M/tkl/wfl1fG61phjPOnVtd6cGeDBYzma6NYGy/wJ5EDJHgyA63zCH/SMbN
WS7379jxf0KTLK2iGMPOpMPJX3E9vz9WfE33HiPJTej7phESUBSUS7q0oxlM
v0TryrHtW92DrM0fMKw8dzx0HiWBQv++ZZJHDOM/P4X9nd9wfzwnQq0ZjWP4
MqnQQsZI4FiSwX9UMZ4xEMCYVBvH8/kUpzLp/iPGo1s6UYkTJKhYH7/UdPsp
488mfTvjSfy8yKzuo6RExslwY1X2KZx/HCe0FSYSGYJRZ75emCbBy4zP5XeK
nzGuUWxf/TeD+2U6j/fOvOeMrgcOtz9if7c4GX+jOJnxIN5Nft8cCZ6PCney
z6YyloS82KawA9R3FUZtS2McT7je+3ge1wNZ62D86XQGOTn0BvN3EqTMZG3R
5MtiXJaMNHmFvaCeWhIRkM34kPZA6vwPEljvfpyUyvGScTc78eP7Rbyf6ywx
wTvzGHOyqSn+S7h+OHcMZHPmM/TzMz0Vf+L84BI/QKzkM7iLi0QfLuP5t3jM
ZfhXIeOCavlP3V84/y5nLu/iKmI0l1U3/cPuDWOvO7K7mHG76t1l298kkHs5
K9MU/5oxceCTlsAfPL/cDtx/OlnK0Knr3tyMnWKa8Xzv4XIGe/Nwtewa7r9n
5NhWGW8YdnoTD4axHSgT7Fcdqxi1rbMOD/6SYNCR9z4ZqhkB7Suk1X94/pz5
QSvhesv4avLvW/Y6nse9qpEfmGsZ0MNSasVEhlNim7J2s9cx1gd4res3kEF8
775cQaUGRuOVg/KXmMmQ7ViYbqzbyIjk9GCispBhs5urtoNzE0NE5evTs6xk
OH3XJ9h7qIUx0Upz4mDD/19vqFO45z0jz+4I4yV2tb6dD4S3Mg7cL+hfYyeD
kswTQ57yDwxO6YmM5xvJQO55OT9/8yPjU5WQ1xEOMpRflBxbPP6JYTMXJBi3
Cf++bs0cG7WDsTuodJzgJAMtsufRD+ZOxtLW+aJv2OKkfRqVa52MwMPmJorc
ZDAICm8dI/cw9IbCt/dh+5Txp+jK9TKo16oX/XnIMPMjJyfhZB8jJVn6Xisv
GULFrC3O9QwwnNQtz7rzkeFuzadHJJkvDOVP9+WE+MkwZurfuS1qkFG/vtbq
QCJDrLRqFpfWV8adGPmn/GQyVHR7nVDnH2aYypy7VIy9evj8WMToMGP8ZBsX
G5UM18W3MKtnjDJyv7P0Z2IXfBhwdXn8jeERopJhRCNDj8Hecb+HYwyOoiSd
RDoZvh19s0YqmmB80OsSOCRAhsg/NIuHXZOMuBHO8VnsoOKCF+ss0wwpkluQ
+mYy6La3GmXfmmX8eJF2fBh70tjJWHJwjlFKfBa/vYUMIX82s30XW2Acuahd
07WVDLdOm1mMJ31nkFm8o3yFyPDVL0pFyekHo+/hS6vt28hQRosPCDu4yLjQ
ILh++T8y5I0KB85y/mSEi02r84iSQfCEx/ANwd+M46+FuQqx1QvNFjap/mEI
GRzvOylGBlKdsKHGuVVGtm+FR5o4GRbSa859nPjLuEr9cejYdjLoR8cHWKB1
BpEpKbCMPcW68YmzERPR2nP31UFJMhwNz5EpbNtAxDjXBk5h+wYdLBuNZyYs
2X8bR+0gg+y1YL45JxZiQdH6x5edZPj+PKuRT5GNIN1VkvXaRQZ7qapgwX2b
iFy9Nn7Z3WQoYk7vFbnESehzOv4YwZa/T3AW5nERoUEJr47uwecZ9CNxlwMv
sfOgWiyzLBnSY+zjvT/zEfVMnR7F2KMlGaKxefwEsw+nhqg8GeiPd3v9MiIT
iaop27qw3bmcS00VKATxC62HKpDhV+sTWUEhKuHt6l6zpEgGk4Vtr9P46cQW
Of7kdCUc72d+LR5CAkTJbEbQmb1k0PnPLShbSZBYOj+k06hMBtaPu+1uhW8h
7kl6S/vuI8PnJeXi0++3EnKjNG4FFTK03nsdvSq0jbhkqdv2SJUMwiURVzSm
hQnubWO5BmpkGJkNr01+JUJk9PnfY1Mnw+7qLd6efqLEmEmRibMGfp4q4/L2
3eJEIMVw33YG/r1qf/2PpO2E2MdpwV7s/mcTP35vkCAs9UQHDiAc3676lkSW
HcS/TeUVv7AHFCyUttF3Eo8bTBOygAxVrt0a5L1SRM+BcGv6AVwP2zTZk+J3
EVeZdmi2YN+QTHpiP7yboFZWS/gfJINkwm22Dap7CEPVXxOTmmQAas/rF6Fy
xMJyVPNTLTLY3k7YrKQvT9wp3J1lrE2GFg7BiYgtCkSzrLVzxSFcz2mb3qq0
KRLnZ9cMXHXIUDgd9t2iRIlgz4xV2HEY15/00u7vmXuJg5Lvf0YeIYPAV8Ee
ncJ9RLmQio/tUby/xJ/d2WYaxMm+T6c36+N8NIztOPCMQazEXkKt2L2r6ss3
VwliL+U5s4oB/n6JeZ/eHBAdHxijM9i/FydCRvL3E653euqSDMnA3jxqKOV/
gMjZxHuby5gMZ0XY5w8qaxJ6DWmOVdilt6BkRFiLmAo8qOd+HNer9NXeRoo2
IcnkyTdoQob4qmxOKSEdImF55EG+ORmyFHPdXpfrEamN2a4WJ8hgVfAq6e78
USIn/tox1pNkEJW/Gta35xhRTnBtMrUgg8eyyH92UwZELX/nt3/YsXcYARpP
DYmW4ac1L06RIUKD2c38hBHRF6zg8/s0ru+WgobAMWPiq/ma+bMzZLigoBMj
+/o4MSFdv1fXkgzLoh+YjkebEL9aT8w/tiLD+GmvSks7M4IqcMMabMhgfDTt
fKyPBbF18giaxN4rY/ub98kpQqyMKnTPlgzdgxL27E2nCXnLtM4ROzKw8Bmj
RDErQkXeNT/8HI7nraUR1RYrAlg07u61x/nE/4QI9DpL6L9o07l1ngyml1zl
9s5ZEyaeDyXkHMhQUgs8OwpsiFO6Nsy92Ed4pU4MXbclLswvl0lfIENKdpvw
Z/lzRMi+/2RbL5Gh48zeC7wbHIm7mya4rjmRYX+CdIAKzwUipj9vQtgZx99C
g/5V/CKRcl3r2WUXMmzIWZ3RcXUisgz5rm+5TIaLXsT+wQxnolC81+It9t8R
McbHay5ETcNFGu0KGcRehO6oa7tMND1U/lGBPdC6Htke5Ep8uMDUds6NDIZ7
SC0LB64Qg3zRISXuZJgIqvF2aHcj1szKVi08yKCScmTsjIkHwSId1MPqSQbv
2m3/vuh6Epxr+q+ysUNSne27dL2IzYnDTuteOH/2rmfGOPgQIq5ZumneuN7m
J9y2hfgSOzWv7jT0wfNDJiQ5Ns+PUJ7YNPzMlwydn8/qBN73J4jS9kpdPzIc
3tw7PaN6g9AKf/JoCVtq9+Rs+tgN4ricvIm2P863M8JPSMdvEpc9zBsnA/B6
8ndwpgYFE7l8r43u3sT9cK+k9d6zt4iF1M2flQJxP2jjLR3UDCGcOvu+Xw8i
g9quVtFciVAi56K6j0QwGa4ZfFi/LB5GzLI8Zm/BLufISQmWDicuKpzeSg8h
w/E1XZF/JneIzKaKlDJs812RyhTvu8SU1X9yZ2+TYdPDoFFnk0jC4e6gZlYo
GZRpn7ZTYqKIczNnnQ5EkOHyyOWWC1bRRMrNmpVx7HgXvVV/8QfE6BbxmxF3
yHDs/O+q1NkHhM3h0dieu2S46clGEY+LJZ4NaYr5ReJ+6Rlz859HHDF8LSVL
PIoMij6j7ktnHxJWqeeqne7hfi16oN/o6CMigdGgS8U3qT3snONGBo+JLx07
ul5jP9rK2/r41BPiNMvkFOsDMtiwZy2euJdAnLS6QH0US4aNLqc30ohnxMNf
LU8hjgw8ux7Uh00+I3ru7JYaw26S/Xt486PnhFnFLEM+Hse/L7XlGCmFiDmu
39iF/eA47eP5rhSiazrHyOcRrq97CXUez1KJ41tc7BsekyF17/4jrgZpRHTe
h+8Xn5ChuN5au18hnWjXkfchP8X9YapNSm1bBmF07Ufk6QRcH0fPLc1zZhHH
OtxKF5PwPFM7HLtbKpe4c6FT8+EzMvxXIS4FlblEK7NyG/Ec34dOZIkHmucR
evK/Rm4nk+Env1StwZN8IrzRzEk2hQz/Dm+rUtYpIFosS1Y6sP8cuDi+8XcB
ceSOJ4/IC3z/Mtyoue/SK+K2RF9sHXZ3dYLRulwR0ViuJnYhjQyO6lHmwatF
xKHpVeWidGwrSeaSlBJCU8fPSi+LDB81uu6cNSknTrBr/V3HFtX5dXP9QgXh
VMsVX5BNBhdZCf2F8EriITxs3/KSDCemfX4+u1ZF5PyzdGnF7taf0y8RrCZq
yyV5AnLJwPucjbfiTTUxp1KoPZmH539mcvwlkbcEyy+vkcf5+D6atOIS2veW
EHy139+gANd3ddeAx6Na4qB82+viQnzfag4j1yjUE3FSU7uDi8n4/Ynz0oHK
JiJ7PLdJtQTnW9XdH1uKmomalGvnZrH9dTjyZotaiFlRtsTjpTj/Lke2ve14
TzAPtWhwlJGBfy1mwOlHKyH49F5vGfan0cZ2eZ4PxIEtIhTxCrz/+6V39LQ+
EmY9Yy+7sMMf55i+m/xIXIzJ1gutxPV+quGPbPQnIpasHvz9DRn2UcYuXllt
J2Y4TX+/qSFDxhTnC6Yt3cSGJqGYK29xPhZJL6QvdhP0WyMKO2rJcIhv+76N
nT3EfpbLF+/UkWHLqCDv9Ms+wrR636YD9Xi+/hdW/De9n7h4/V/KT+yw4KKC
gawBImY1bPBUIxlq1OKRfOsXIvO1kQ9/E75vEsfeaM0OElXXNm+uxV71bNdV
CBsiphdTjXa1kEEv4ek4T+1XAmaq636/x/3VasfTt56jhElGiHV2K643tr5K
T+VvhOP5Y+tWbTifjwucffX7G/FgdECl8QMZzO5oaD6KHicynj3v9P6I4zVe
MldycYKosnJ0lf2E699cleOg3iQxNfArM6adDM+Ng7NixaaJ9fhKnSMdZOC6
mpUzunmGoJ4I+vYX++2h9l3PN88SRCdZ2K4Lz0ObWoNIxXni/nuZe4q9+H6S
paP+auY7QU4Ul+7Hvs2xrgg3fxCRrptrAvrwfbxua1Hwf4tEuADbj4/9ZFD9
1BnYY7NEcE79ue05QIbgJrMbsrw/idvlC6Kin3E9x0q691b8JILPDhi6fMHn
X5iT3iTzi2BV+jQpMIjruyzjuu3cLyKAvfHGG2wzKvOt1cIV4npmQR7vVzIM
jupkdRz/Q3j8DCVlj+D7jI39Q8/Rv8Ryg3/68VE8H8aNRxT7/hFu8Vf3r2Gr
zFNkvbrXicvI+rLuGJ5f3tm/65eY0ALJfNMi9iOVa+7nr29ATqNHk+LHyfCG
HuitxcuMHG+rfpycwPX889YrOVUWNGEh6xA1SYbXpSsW/p0s6NweiQ2qU2Sg
pk373LnKimza+eVDpvG8fX9K8EoDGzq1bSJSco4MjWsJD5zSOFDf/Oedrdi7
eTakljttQuY17VXu8/h9c+Ntvr9qnMjE/s1C7QK+L+ZG72UZ40Idqq9CLn7H
+fp+gEm8jhsZcWeKUH/g/KE+N5zM4EH6eTEGNos4f9mbN+8K40PvboZPcC6R
wamn9TDbQX6kaxrgn4/tPsx0fHctP9JZvZS7YRnPHz3vSq5PJLRfS4s/cQW/
H5T3y+0UpqIqAfW0Q7/xPAtOXLSroCJiSg7msU8aHdqsZUVDaneFXIhVMlQu
37h/oZCOlHoX2wb+4vz6pfOfZPxmlJ85aR/4D89LN00Z5bNbkJzf4PqudXy/
bu4cZ9+zFcmIt8h6M1Hg8O+rW8J7hZDkpWd3N7NQoMs0Y9s5QgR5nfun2ob9
8MSx/G91IqjV8uRoICsFpkljxrGGouiaEVl1gY0CEqC3zHtVDL3TdRpJYafA
kYFfqnQ+cSSi1RxhsZEC1flaWz9liaPGfQHD9RwUWGo4sGtpeTsSkv8S7rOJ
Ao9sHRafPZdALtJq+xQ4KUCrEoiLNZFEm7f9CHvCRYHIPNfHc007kCOz9V43
Xgo0SIT9aV+RRpWrFYNSfBRw+5snu6tjFyL/3Bw6iK3b7zghVLQblY1//HKE
RAH2ylCjyDt7EP9XmdtMZArUmS7buAfJItu+24pF2P6eCqOf98sh7vf7Q0Sp
FDiYKTKbqCuPrOqfKHRj7yj5136tSx4Vvvk9EE6jQLYnR8QlOwV0Oj9PfoVO
gX0vvO7vilREObGi/a2bKUB3iV4zE1JGzFG+QYFbKCAlYeUZ3KKMTEN7ZdW2
UuCDumdSx/V9aN0nKjBFiAJqDHE7qxUVZGC9YY+PMAX8FE6mnJLRQMkWp3vk
RSiwPs9hmb6JgVaOvw4Yxy4x+mR+bIaBkg5d7jYSo8DkTVUmszqEfuz+6i8l
QQGLmdWjed/2I21Jxq5B7FmDotY3zgdQvPDDzmhJCjh4RDx+/fcAOkg2lGba
SQHHw8FfeiU0UfSvqvYuaQq0kNfnBtK10cSCkF/4LgqM9v7wNzQ9hDSmPHYe
2E0BdUONE3EbddDogJxvtgwF7COODm7zPIz21iRKBspRQMa0s6b0kR6q8b/v
pyxPAUH6ioOb01Gkj4K7JrAXRU8EDWrro3PlF4KPKlLA6vc39xxWAxRbpDxG
V6bAoajKrvszhmi7uxTRhK22+ZveeycjlKsoFOO9D+9fEj3wXTJCjbkbtIdU
KJC8RcksnvM4+p3xLiVdnQLy482XtW1NUbDDm78WGhRQ6t3K8YrVDJF35pvw
MvB6cwJ5zdLNkHRKLNsVggIXc2N6n6ybI4sEazuN/RQw9d2VKj9igSZOm1TM
YQeMR8oopp1CbkI6tKQDeP0Fy7o6rqdR+EOZOjZNCpwseC0sRbdEFfdXtrdp
UyBwTDdO0uAsOmI07XPjEP6+wKXQ4sazqIv/S4eiDgWuHDb8YappjeYi3gbG
Hcb59zgaOLVskHDIndGzehR4HWsqau5kh7K0b2hQj1Jgj5eHv/OGc0iFzS26
Htvwn0OlfNw5ZBBwQnPXMQq8OhTYxf3RHt3w2Z68ZIjr+UvI5khXR8SrLrCW
akQBvtgP250kLqD435uOnzCmwHbTyJKcgQso/+o8S+VxCjh9FlUpN7+ERlxK
bW6ZUeCUzo2lH5ouyEU2u0zVnAJc12ibT+S5oLXZBMoMtnH16TELkcuIeiHo
rcFJChi8vIZc2V2Rlp2B+JbTFOCd1mRELF1BaiquV6uwhz5mc8Q4uSE5ruim
c2dwvf7rkzeYdkNb8rqdCywpcOC7+RDfrDtaWDtTrmdNgQSfKNpRigf61ubP
t4idU3g9xjnZA/U9e2b90Abnl3ChdYmKJ6rTGeMYs6XAt0Oje9YcvVB89CWT
6/Y4Pk2nftJXfNBd+7tpEudxP/LKPKOU7IsC1fJWW7C3SKv0/jX2Q85DS0mC
jhTw7JLQ/Vx5HWnu9pnLu4j7nW2K28zZG2j+beit0csU2CDtmcDXGohGY7L6
Ql1xvztXwM6jHYT6HFpl5K/gfuhfpnmzKgjV8pE7fN0o0DGq1c1REYziLR6K
CFyjwDbRV76NnSHo7p6yKxXYUoM1If3Wt1Hghs/1Nh4U+MUsb9r04zZyfiFy
KdeTAsFDQg0OW8OQ5uKL14d9KHDGZe8G1agIpFrfxL2ArVuX3fBI+Q7a83Da
MsaXAjWsM8mHBu+gzUiWfcSPAsL69gMnt0WiudAiI58bOB9zuFybI6LQQ7G6
6ZxgCpAnurQKmqKRrdBRQ5tbFMgKu3Rm/eADJEfvLBIIwfnorWzUWfUA1W/6
dsP/NgUsp82/H62JQYsLrIJG4RQ4fm61qu5bHKqcCvdlj8D9j2KvWuD+EIWO
UkdKseOPi56N4IhHIj3bc8TvUuCqjV24xr5H6OgbTc2fURTINRo9wlz6BAm+
fpeefo8CJu8vLrZZP0Wj+cf5ztzH8d98+c0KbwLySrXtq4+mwEa/ZL0WlUT0
IiLQJS6WAs+DuDZ35iahDadqH6s/ocCN8GnukL/J6J2JHvMC9oiZtZ2zSwqK
O9Zhn/wU538+S2rQtxS05+CoAk8ijo98RJd1VyqykGZt+vKMAj+qvhZ/bU9D
ktvD99x/jvN5qNN0/4l09GMbNfpQMp4XWX0Pt3xNRyGk7Za5KRTw+Bbnc3Yl
A71aOfjzRhquN473f05qZyP/Hy0nldMpMNgZrp40lI30ZoyrprBrApVH831z
0PCgTZhxJgUSd+Y62VS9RHwNN0UlcigAV8AqRikPOTx4e7SxgAI9YU//7JYq
RE+fe+W4FFLgY2eSt2tIIfqUJ8+3+RUFhPK7h9FkIVJrTfjgUITX/7WhKfjl
K8S90ceY6zUFpN95iZkeK0ZAUywswOZ21Mm7UFaM3MSnqKdKKbBJU4JLUaoE
fUbmXVllFPjpuunLzKbXKMdD6YR+JZ4Pr/b6vxgqRcPB06+XsatETzz+drIM
0R8825LwhgKb/+TxMXeXoet5/APzVXierU+zzXaUI4OpmTNRbykQ4TKeumei
Ei1apNh2NGLL7VBqLalGOx1P1fk04f3c2lF8kL0GnfKgSEo0U4B17llGz/Ea
VBt9Y9y9hQJ/J37dv/SzBsW8P+0o0IrvC+Wssk80a5EKorucbKdAXNE207lt
Deji0fcfmDvwftPftw/6NaBEi0D5TGxdlx1aEl8b0CaPHz9WOymw6sb+6HV6
I+rLbXV/0kMBzZP13fa6zchX7Jb30GcKPJBXTR9qf49CzJLnYr5QQMRU0PGt
SSu6F1599ugg7v9bu6qu9rai1OXVQ6VDuF8tXP3v7Fgbamt2oUaPUODJ0pTM
Wa6PSPTKyexDkzg/PHV52x3b0e60ayL/sP84XHluWtaOlD9H3y+cosCCfOzF
Me4OdORQm6foDK5Hc/8AxYIOdGWrpvbqHAWKP1Nqq7m7UN1bmcGcJQpYf3Go
j/vZgz6sHDGy+0mB78/+bR8360V9Mufrti5ToCAsxM22rBfNxT7LDPlFAZvD
q3uXgvqQwEUBD+s/FBBI3hafJDWAHKgbSHQmKvwY43NZyBpEbof/C3yHbX5+
L+MuaQj5+an/CthAhXJpTfOjR4fQvQn3z/PMVDD4fPtETf0QKiufSm9io8I3
PpH7L6q/Ih67jgO+XFSID/x6mdIzggwsQ2zkuKlgpP9kV7HQKLp3QiNwBPv0
g4S1jdajiK6fXHeYlwpaO30MCuZGkcg+Nx06iQrPw+mlwD+GFDZR9V/SqeAZ
f70/1W0CubE0OlkLUKHr7BybRe0EKvrrfZcmSAX2274tb6mTSP3HyAfvzVR4
ppJhUFg8iTT7C4x1hKhwaa1IQJB9GplmG5/4KorXfzR3C0/dLIp7sdErWgzv
jyNPWlRkDvUllcUfEqdCRs5AxrT3HLKMER/I3k4FTtJAyPO98+j89cUzXjuo
8PphKfN78wXkZXjfliJDBYdHwat5376jpz8/OWupUCFLwDTEef9PdJi1KGwG
WzJZNnPP1Z9oifzwxX1VKtQ49crszfyJdOUsB7+qYT/u8KXTl9GK47S+H4MK
tn85ioXml5HxEMue4gNUWFDNf/EmewVxtCjN7DxKhUrj1W0T5mvoSmLMeW5r
KvB5H8mr2coMS52DMdXYjwxH/UYIZrjKJVV31YYKLm/s6TzWzOBxtUz0qy0V
NpnZmlxLYwY/vaG+V/ZUMFt9Gxq5lwVCVqSOWl6iwo2KtAaKASs8NqiQz/eg
woBGvc/qdXbYFsxuZe9JBbmItByfZ+zwtOzYHSEvKgSpbtOerGWHJMnhqWBv
Ktw8OJPpv2kjpK6xJ1v4UYFSkH4yPGoj5KYZ0NkCqVCraj6x9zEH1DKN/jG7
S4Udb8wcK6I5IefIqe0bIvF6tqbklOdyQlx0x9EM7L6d6xvuvOOESzvrEtei
qJC84PD8OQsX0A1StZOiqdAziuIDnbnAIeH8vamHVJhX+3fdicENxyeHyu7H
4/180g6XM+UGQvHEN41HeP/tNO0SJ24gNxxWufuYChO/zn6gJXBD+Zz0F8UE
Kux7pzq3tMYNvMSslG8yFVg+KD5IyeKBlVt2xpIpVLh+/bPQpWoeGPn42acN
29TBft98Jw+U2LW2ib7A8TBQeCjyjwfORrx0r0/H8d+p8kJSlxf0unckOGdQ
4SdTvUHeaV5QFk1sFMykwnvpDWnfnXmB81Xk1gtZVEDMNylR93mhcOByNd9L
Kvx2f3NQvZsXnkpOTZVgd07vz7MZ54UQF2uqdS4VopuvX2L84oXTrMb2hXlU
uM3b/Z6Zzgfsu5V4ThRSITyE/qDhGB8suGcpM7+iwmL25baLp/mg/812q0zs
HcOP7rxx5INcY1rB3yIqzFXe/nYhkA9OeC+bP3tNhcfpFyM6CvjAJ5RNTqeU
CkxXagaeV/LB04fUjXPYaYWqgyJNfDBcrPBKtZwKX3RkcsU+8wFrw/7wQeym
U4e+pI7xwY4uA5ugCip4N9f7dc/zwcUlJ9LHSiqwPhXzfMrED3dY/CauvqFC
fXnu4xJWfsglR7wRqqKC5YV7l9I4+OGnXOal89U43oLNDzn4+UEQSjV5a6hw
bT3ttQeFH9SONW0txF4+t8BXSecHv0vjTetvqTDjd+lXgxA/JPosJ6bUUuHz
sR6zcGF+eBvG5qFbRwX1K4IREmL8MBZPPfYdO29gZ8D97fzAkSEuGVuP89nQ
SrxLkh+kXyv81WigQvVl1zOzO/lBr3F/xzD25pcuYr3S/ODUbZAZ0kgFJ/kW
y4e7+SFqzDJgTxMVZtMMOfbs4YfCn04nOrCzeZT4YmX5oYvVT86rmQqNeiwO
HXL8sEKJ2CjSQoXB4TyuMXl+2CL++Esd9rJN9FyLAj9oKGS+uvCOCkkLZzmD
FfnBcn9pOOk9FYLjlE5QlfghwKDJphh7W8btrqvYyZY9aqdbqZAzv/96Nna9
0ziJpY0Kd3uTDEuwJ32XJ9Kwwx4YHo3F5o5gq9L/gE05feUI9p7H1NglbKdf
ivUf8fMMMsWd4j9SgeqRflgG27VUQQs+UeHFriXm03h90U37hcawBauzV6zw
+ot7DBbD2vF+Uj8pqOL99Y1bNst3UIE/obr0G47H2rJTUjf2C1B/dE6GH/5j
9/Pw7cT5f+7sUNkufgBaxDHxLipsDy2K/SbFD9bbH0s2YffzbfswsoMfghQz
/zp1U0HYu+thkQQ/NBs2ZZb2UOGc939SfSL8MGPVE2DVSwWRVBb6rv/4gddl
/AR7HxWOtI0NG2zlB+M7bBxG/VRodlb32krjh6tPqIO/sJt4O65WkfghLku8
6MkAFSzMi4rVePlhoHm/7eRnKvxnYuL5kp0f/vUaqN/9QgX/e4PKGcz8IDpp
Sd47SIWt5TG2P9f44NxGvyr/ISrIb3nUqTbNB/MHS4UER6hw5v2m+v63fDDh
Kcv+DvvC8ruxwyV8MPQyef76KBWUBi3tzbL44MPWuzXj33C8TZmPGt3ng/xF
m/PFE/j8PB1aOk7xQYZUr+GFSSqs/QxcPKLPB88s9dWFp/B67rw7ZYv44H6L
Cu+taSpIOzeI3hThA7fnPIWmc1S4p5V8yf8LL1zsDXjCOY/jv5DkL/KeF2z5
VoIrsRu1VW6eL+MFE+/hE5LfqSDB4jX+Iwb3K+Pi9Z+LVMi176lgOcILK8xn
dR/8xvFqfbK5IJEHFlS6lA7/oYLNXPGb7aE8MOGk+99fbOeitNbDV3igt2/v
d9s1KoxpB1jkafJAaT5nrNI6/vxGuyDPKDf4WBcOt7PSoG3w23whjRs21Gz0
IvPT4PbybRceXU6w5s46U4p98qvI5dQ9nFBranDQmkSDxXNB3TQyJwRPx3Hn
k2kgScyxePdsAm6adIIhjQbVWczK7jabgH5erzZyCw2ehvF/HHDlgN28Ubx8
EjSgIcmErb7sYHZyy3MuNRqsTg+XG19ghjst2+OrsZnDTrFs0GGGeg3ZqGvq
NGhZPLhQtJ0ZlIQ1r49o0OD6jjBu+LIByKOXLMoQDbaeydr0xXADvLtUTbmg
RYNdW/iEV9SZAG6cD2w2oMEor1Sa0f41dO2Hq7e/IQ0IGOgoXV9FOTa+rspG
NPjXUjJtVbmKtmnfs3pmTAON4WbZCxqr6A9nuYaXKQ3k/WetH6j/QYXRfD+l
TtHghmZ+3pZDK0gqrcjutj0NtO+cly2KW0JxX8w4Nc7jeFozrzw5u4Q4aL9z
5rCDbth96pVeQmM31FeMHWmg6PE05FjlIkq2qAr97xJ+vp1WAjH5A/3H1/yy
wJUGQt9HWryOfUcRWheOn7tCg7ePwwc+bPuO1ry5fwu60aBZtOWw/8wC6p/Q
P+DnTgNOetH6dOgCiqtp79DxoEGg8Inu5dh5RL76+fdnXxyfFzEuEpUz6EaW
39NIPxrstMvjFzo3gxaGhQ8evE4D1b08mo28M6jtmHV4mj8NvoWf3LnfehqF
S43/d+UmDeoluAbySFNo4+eFgxy3aSCTJeFtfW8crR1kvaMQTYP1b6TVq4Ij
qONOh0wv9oKLVNvfvmGU3Zv8/voDfP6CL/MEng6jM85aPO9jaKBnaHqaWXIY
vYkPDrd/SIP+mzztzzS+ooDvHGGPn9Ig3GS1cuOmQXRSo1f6YAINJJjn1UYT
viDFW+nNk9j3/NosP+37gkaFjnDuS6JBXQpt+rPjZ6SlE3b743MaNKRG77z9
uR9xJPCEsKfTQJa33eTLYg8amvy8Ixs7vmIqRuxRD3qtlNNgnEED0cUCEauD
PcixWX9jUiZej29ubFVMN2r5eTdYPYcGnbGKtpM6XeiOHjnIpYAG5XvPnbfv
bUfnYoe3CxTSINRq6KnX3XZEDOfXVmCLc+oEXtRuR/PXjFm5imjg7NoQ/q3o
EzJ4Hn0zpYQGX3Nym9eSPiLKH3pAXwUNPBpWVG8QbajmwlGlu5U0sIoI+Puy
tBW5fL45dvANjscFx/Z0lVb0rmrhSHYVjrdut/6s6nsUfKuZ4v8W59eDTm41
wxa0QvVLkWiiQQA57bBkWwNKDS4068O28vK6yHSuAR1fmdp0t5kGb754U87+
q0d5/WZOKy00aDeLVE9RrEeOz+T3tbTSIMEnnanuZS0akP3W6NJBg4eO3zg/
TVWj0KSt3hKdNLiaprG0M64aqVCMZPqwI9SIsjrtavRguTLqYDcNUq/luFa/
qEL6FXEn6X04HmqtUdkLlWhtTxt3C3ao+1bXkTOVKDOR7c31fhpIDY+Hnm+r
QByBruKTAzRAbaXkg8XlqPqI3kzZIA02x3uTnJNKkUt5wFOXIRr0Cpb9MBYt
Rf/teW0g8ZUGQ9vTK/Sfv0ZeJMlXd4ZpMOXyysExuwQp9az7nv1GA8uxL/e9
e4tQql0B38ZpGlTyvWegsAJUOS1Y8RJ7t6z23m7ZAtR52c/RfIYGjkp5v751
5iNWv8N16bM0yPz3FR2XzEdnYwa99Bbw804e+TPTnYs8hLR2/sSuVE8TLLyV
iyKfZXQ++U6DKnR9rEUlF1XmuMvN/8DxSv2eXsR4iYQauMaiftIgT40IW1nM
QkpHL99XW6aBbkDYqxj7LKTb3gUj2ELn5SLufs5EXoNJj5RWaECpfhTu/CED
df/aZ9j9hwafD64Gb/iYhu7ttC0XYqIDW83YX0/lFJSR0+RQh81cuL5zrSUZ
VSvJCjhtoIO83+YhC5tktAB/Lr9hpoOZvXZ4bOxzdPRE5A5rNjqo3xONYNCe
IdvBnx2c7HRYDTfaN1mShHzsLAIKsG/8qs1XP5OEMi5LfmbloEMV6TSjLDsR
cYSW33vBSYcDPxOjxyufImF+MTDgosP8ZOK+w8JPkXLMrdkV7P0NWQvRAU+Q
3TMjnSM8dOALCXpVZvgY1ZRO/J3ho4NoS7ooOyke+U5THRRodDgkNB/AzBmD
3imTaruxJwxPuArHPEBCATzC+MUNmt+ZD9tuf4BKBdk7GwTooP+88pTWoWi0
ovULTm/B7nTbeD7lHtKJWnzEvJUOtFc2Tu3776HYgfnlF9hWrFn3xIaikPKV
iawfQnQo7hv8rCsWhd/vegVDhOkQzdNv9+XaXVQz3XlFRoQOUQPj/Srv7yDy
vk+tn7CZL59LHJS4g3LfNwduE6ND2umK4J6hcDT3p2whfzsdtAr+HV3wDkUO
Jk8bP0vRwei7iLuoejAqSYwXvylNB3K7itt0ZRDimInx27mLDncMipCCVhBK
DbirdGU3HVCpfgfHiUA0+tI/kUOWDpynqEUpjwOQ0qrPajb2yN9NT0QUA9BN
bU9TYzk6HL0pNFbdcgOJfb7M/VQex/MVT/AxlhvIktPGQ1GJDhJhkzRe2nWU
Y2LZ3oP9aSj9yyFrP/Qv0WKP3146zOWXsCsX+KLH+46PNirTYTxNeWLC0gf1
2WoZnFGlQ3iM42jzsCeSzt2fyaJGh/bOjDsNup7Ic5XBno6de0zw9Y9iDyR4
T7l8UZ0OsueelMQ8vIZM3uzYeZugg9COY29d/d2Rl3zoAw1Eh4wJDsnwZTeU
+HyGeQH7qskrRq+zG5q6lffZdD8dNh18X9p04QryN1C/v12TDjf17n1SvHUZ
pdQ8YerBvvzr8FVB2mXUrMTkFKZFh8rRDOWXKS6ItrlO54c2/v//6hqdnzij
jK/6/94cpsN3qStN+r4XUZtx3oUrR+iguKckQ7D6Alqqo/RK6tJBTFzT6Pym
Cwhl9BRE6NHBMG4gXfGFA+pwtXawOIbjF1f87gPVHv0Zre3iNaBD4J32R1u9
ziFhsx2aNdieLInkPSN2yFF95j8pIzpkqR14Ylxli9ZZrnYsH6cDPzny4alk
a7T9as/+DBM61LCJ/nwsYY2OTKi9PG1Kh2THLTtYMs6i6Hfrt2vN6GD5q7Xi
V5kVkn5wG907SYfCwcCn0mxn0LGNM1laFnTg7lB3l1Y+jdw89bf8xqZpRzy3
uXgKvTlN+Wl5mg71tS852SdO4p71JEPGig7KfE3bWgTNkWfcusBX7KXymzyG
LmboKad1UPRZOnAZvind+c4UTc5JWq1a436y/HLYNtIE+RXn0prs6HDR+pvf
cR1jlCxFuelzjg6MN/wn0t4YoaZH7guy9nR468CrEKZmhKj+ai0x5+nQ31l2
OFXNEKXr1PrbXsD5EXmer+2tPlJpM5nkukiH37nvLjv8PIrqTcYNC7Dv8qhu
Ftp9FA1bc4ozO9EhT+omm3GaLhLyNahLcMGfHzrmP96jgzJYhmUOXaaDn41a
urCcDlINvRIzh92tamAXFXYImcbG2DOu0MH3Pz6JCn1tdDdvYFO/O67X7CfH
olk00X8qTq4BV+mwrc+tYNHzIMqqXO+TukYHjaEN8HHxAGpsEcvy8KCDg4Pm
tfzv+xHz2Hl9ujcd5NxVGT0UQJEX/xRVYEPzQvvQIYSEF8OE7Xzo8KFLc2v7
JQJpML1cKPClw8KMuJpLswZy3/LznqE/Hez28zj316gi1qTg1RXs7rgdCrY8
qujeDkHbxBt00DnjR0k7pYJylNT3zgfg/jzEZjXNsQ9N6N/oDg+igw0o7NJ7
oISudpJBKZgO3k0yXA+YlRDbqeS0fuybzufJlVcUkZhDg6d0CB1SSaqJ+mcV
kMVNHqHGUDoQG55R+8/JoalNCYHOYXQwfrHh0/avssgjUm6WHk6H+2E3hCf+
7EHRT4wr7SLoUEpnHG4GGdRa/NCKJRI/j9qaUkWSRqeIXU0Z2Ohg05EoYyk0
XVsubxRFB+tcv+jb8TsRx6dB5qR7dPDKDv/puG8H2j8jkUI8wPPFsZ7PNX07
Sn24TWQJu5B34roF63bEeYj2KD2GDrH3w3bJ2IqjT0msUdQ4nF/xZWx/5cWQ
8rE1rmZs428b7IsTRNGjtcXg6w/pUERfOv2JJIpszUZ8puJxf0r2581iFkGN
rP2/Eh7h9R6yf23BLIx2539yNXlMhyTLnsrdG/9DP7lrzlc9wf03Za36+A4h
FFyTeDwmkQ4NAgeGB8sE0ZRzXJtuEp4HMVP9574LIP1tkUc2PKODmtE7JSEZ
ASTgcR0uPsf5tDYkdKqIhnwkrpWJJtOhbJ733R0WGhr65KTcjZ2/PZETTKko
TebM7gOpdCgp/mXqzUVB3P0mL35hj4k6lhddIqPLIUfFsl/QoVxRzqi9g4TU
RjQEBdPpQFqyOf1fHj9qidvKOpeJ999HLXfS4EFy2pTrz7Po8PoNd5DVIW4U
vcj5xzwb70curXL9BBc6pf9n/m0OHTZUj2w482ATqlr97uj5kg5dDTN9WRUc
SCJ98tueXNzPGl5nN85sRHMsvf0P83A98Cb0xp5hR0Z5H0yP5eN+MpQtOZ7I
horONH5kLaCDvUrFfvVJVnT9dXGDcyEdJHd2cqmGsSB+55gCrWK8f5EW1zq+
DchN6M6eVez7vfILF3yZUE9TUHpuCR1G9zvohg+uE0nb3RO2luJ4G2xOY6/9
S7B/urjlI/alrlwf9HaNcLxu+yC4DH8euNWxrnmVUOgzDvteTof33b9PFH7/
TcTc0mV/UYHvJ79c1hL4fxOrSgdvnKrE80rd5PNn5RWi9q7CtYY3dLB4tnzM
IW6ZMNEi2TypocMtvaXJ5vRFQmyjSq/2W9yvf54XvcWxSCw0njm2gP0lW+Sb
6sUfRKhelvrBOjqseJxJm9v/nTDnbc+bwS6nRG1ML1ggJD783hFTT4fOrXeS
m6UWiCrjQ9TJBjr8ipZX9dOaI+7QnELvNdLBdO9bevKRWcKiO5pJownfV7je
pcSZzRDLJ4dn7jTj+HAt39gQMkXssvatU3qP+13S8tNijXFiRTxZ/Qv2Hjtz
ym6vMaL+W3PerVY6nB9EX99UfCOsHASf9rXRwT1vbem26SixZxeiBn7A869u
oskuc4RYnbELlflIB94JFtYsthEixqXgqv8nOhjskKn40PSVsFXom5Fqp8O+
1d/1ZxS/EvJLTDbt2Mv95WpFSUPEu2v6xyQ7cb849HH2hOAgseH65I533Tjf
fG8IvtvZT7QC/1P3Hjq4qFDPNkAf8Zh5H1W4F8877TiaqlUvoRwcyOTah/O/
rT3tQWE3waqTeXVrPx1OR+6v45vsIj5u+jRTi73G+Tj3v+1dxMUI4V6Bz3je
vX2ktyu7g1A7pn2sGvs8Se1f7Z92YiPpUp3jF9xPzModd+m3E8+iS/MqBvH6
JU/ap7J9IrofmYTaDNOBImufJWvXRiSf9mHiGaHDUO3HHYnGrYSr8POrRdh9
S67hPLrvCe7nC9abvuH1t++9V2ncQvTZCvTmY6dZnrO6ZttMvJAkjp0aw/e3
98JGAr5NxIGMMPWX43QQZkoN4qlrIPgv5ueZT+D6zHzDM7lYT3yW6d3BPEmH
ikk/qSKpesIjT5JqMkWHttYXrw6m1xLZJdUzf2bwfg9kNln9qCby0XX+pFl8