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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    103973,       2073]
NotebookOptionsPosition[    101350,       2026]
NotebookOutlinePosition[    101692,       2041]
CellTagsIndexPosition[    101649,       2038]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"mbstar", "=", "5415.4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mbstar0", "=", "5711.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"B0M", "=", "5279.5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"BsM", "=", "5366.7"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"BplusM", "=", "5279.3"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"LbM", "=", "5619.4"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"D0M", "=", "1864.8"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"DstM", "=", "2010"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"piM", "=", "139.6"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"JpsiM", "=", "3096.9"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Psi2sM", "=", "3685.6"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"kaonM", "=", "493.7"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"KsM", "=", "497.6"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"phiM", "=", "1019.5"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"rhoM", "=", "775.26"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"rhowidth", "=", "149.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"omegaM", "=", "782.65"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"omegawidth", "=", " ", "8.49"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"muonM", " ", "=", " ", "105.7"}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"squarkM", " ", "=", " ", "95.0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"bquarkM", " ", "=", " ", "4180.0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"cquarkM", " ", "=", " ", "1275.0"}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Bplustau", " ", "=", " ", 
    RowBox[{"1.638", "*", 
     RowBox[{"10", "^", 
      RowBox[{"(", 
       RowBox[{"-", "12"}], ")"}]}]}]}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"C7eff", " ", "=", " ", 
    RowBox[{"-", "0.306"}]}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"C9eff", " ", "=", " ", "4.211"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"C10eff", " ", "=", " ", 
    RowBox[{"-", "4.103"}]}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"GF", " ", "=", " ", 
   RowBox[{"1.1663787", "*", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "5"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"alphaew", " ", "=", " ", 
   RowBox[{"1.0", "/", "137.0"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Vts", " ", "=", " ", "0.0394"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Vtb", " ", "=", " ", "1.019"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"numberofdecays", "=", " ", "5404696"}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"b0", "=", " ", 
   RowBox[{"{", 
    RowBox[{"0.292", " ", ",", "0.281", ",", " ", "0.150"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"bplus", "=", " ", 
   RowBox[{"{", 
    RowBox[{"0.466", ",", 
     RowBox[{"-", "0.885"}], ",", 
     RowBox[{"-", "0.213"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"bT", "=", " ", 
    RowBox[{"{", 
     RowBox[{"0.460", " ", ",", 
      RowBox[{"-", "1.089"}], " ", ",", 
      RowBox[{"-", "1.114"}]}], "}"}]}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"NRBR", "=", " ", 
   RowBox[{"4.37", "*", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "7"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"NRauc", "=", " ", "0.00133"}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"jpsim", "=", " ", "3096.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"jpsiw", "=", "0.09"}], " ", ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"jpsip", "=", 
   RowBox[{"-", "1.5"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"jpsia", "=", "184.39"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"jpsiBR", "=", " ", 
   RowBox[{"6.02", "*", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "5"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"jpsiauc", "=", " ", "0.2126825758464027"}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"psi2sm", "=", " ", "3686.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"psi2sw", "=", 
   RowBox[{"-", "0.3"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"psi2sp", " ", "=", " ", 
   RowBox[{"-", "1.5"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"psi2sa", " ", "=", " ", "23.07"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"psi2sBR", "=", " ", 
   RowBox[{"4.97", "*", 
    RowBox[{"10", "^", 
     RowBox[{"(", 
      RowBox[{"-", "6"}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"psi2sauc", "=", " ", 
    RowBox[{"2.802257483178487", "*", 
     RowBox[{"10", "^", 
      RowBox[{"(", 
       RowBox[{"-", "10"}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"xmin", " ", "=", " ", 
   RowBox[{"2", "*", "muonM"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"xmax", " ", "=", " ", 
   RowBox[{"mB", "-", "mK", "-", "0.1"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.768042511440446*^9, 3.7680425886548386`*^9}, {
  3.768042655722806*^9, 3.7680427829496317`*^9}, {3.768043469101739*^9, 
  3.7680434767802153`*^9}, {3.768043520665967*^9, 3.7680435461338577`*^9}, {
  3.7680438318670406`*^9, 3.7680438326758633`*^9}, {3.76823577028961*^9, 
  3.768235831700453*^9}, {3.768236772772351*^9, 3.7682368423044157`*^9}, {
  3.76823944287317*^9, 3.768239442937003*^9}, {3.768244201338833*^9, 
  3.768244216523238*^9}, {3.7685482113105383`*^9, 3.768548213647256*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"037a651f-7204-4dd6-865b-a90153e54791"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"mB", " ", "=", " ", "BplusM"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mK", " ", "=", " ", "kaonM"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tpos", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"mB", "-", "mK"}], ")"}], "^", "2"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tzero", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"mB", "+", "mK"}], ")"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Sqrt", "[", "mB", "]"}], "-", 
       RowBox[{"Sqrt", "[", "mK", "]"}]}], ")"}], "^", "2"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"z", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"Sqrt", "[", 
      RowBox[{"tpos", "-", 
       RowBox[{"x", "^", "2"}]}], "]"}], "-", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"tpos", "-", "tzero"}], "]"}]}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"Sqrt", "[", 
      RowBox[{"tpos", "-", 
       RowBox[{"x", "^", "2"}]}], "]"}], "+", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"tpos", "-", "tzero"}], "]"}]}], 
    ")"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Pplus", "[", 
   RowBox[{"x_", ",", " ", "y_"}], "]"}], ":=", 
  RowBox[{"1", "-", 
   RowBox[{
    RowBox[{"x", "^", "2"}], "/", 
    RowBox[{"y", "^", "2"}]}]}]}], "\n", 
 RowBox[{
  RowBox[{"fT", "[", "x_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"Pplus", "[", 
     RowBox[{"x", ",", "mbstar"}], "]"}]}], "*", 
   RowBox[{
    UnderoverscriptBox["\[Sum]", 
     RowBox[{"j", "=", "0"}], "2"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"bT", "[", 
       RowBox[{"[", 
        RowBox[{"j", "+", "1"}], "]"}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"z", "[", "x", "]"}], "^", "j"}], "-", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"-", "1"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"j", "-", "3"}], ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{"j", "/", "3"}], ")"}], "*", 
         RowBox[{
          RowBox[{"z", "[", "x", "]"}], "^", "3"}]}]}], ")"}]}], 
     ")"}]}]}]}], "\n", 
 RowBox[{
  RowBox[{"fplus", "[", "x_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"Pplus", "[", 
     RowBox[{"x", ",", "mbstar"}], "]"}]}], "*", 
   RowBox[{
    UnderoverscriptBox["\[Sum]", 
     RowBox[{"j", "=", "0"}], "2"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"bplus", "[", 
       RowBox[{"[", 
        RowBox[{"j", "+", "1"}], "]"}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"z", "[", "x", "]"}], "^", "j"}], "-", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"-", "1"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"j", "-", "3"}], ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{"j", "/", "3"}], ")"}], "*", 
         RowBox[{
          RowBox[{"z", "[", "x", "]"}], "^", "3"}]}]}], ")"}]}], 
     ")"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"f0", "[", "x_", "]"}], ":=", " ", 
  RowBox[{
   RowBox[{"1", "/", 
    RowBox[{"Pplus", "[", 
     RowBox[{"x", ",", "mbstar0"}], "]"}]}], "*", 
   RowBox[{
    UnderoverscriptBox["\[Sum]", 
     RowBox[{"j", "=", "0"}], "2"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"b0", "[", 
       RowBox[{"[", 
        RowBox[{"j", "+", "1"}], "]"}], "]"}], 
      RowBox[{
       RowBox[{"z", "[", "x", "]"}], "^", "j"}]}], ")"}]}]}]}]}], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQAWIQnfZN/eoPoTeO3QEaYDpgeqTBTyDdcOyrNYjmm7bJ
EUTzTHPzA9ECH6xfhwm/ceQ4Pg9MB81RFwkH0m0S15VAtNadED0QrcTpaAii
l57ytAXRecxFLiB63XZWfxDd5yUNpg++vB8Eog98nBoJoi9Y2aeBaI5v73JB
dOxH0SIQ/S9yRzuITj98aQ2ITjvrtR5EF+yq2QGifyn/2gmir1+YujwKSLNN
1lsDondq9u4B0UFmVftAtPJ7hTdvgLS/7ee3IHre43fhb4G0VH99LIh2+Xyl
AETf2BdUBKIrt2iWg2gha8VqEF1zvWM7iPYQeAymAUy5mnM=
  "],
 CellLabel->"In[51]:=",ExpressionUUID->"8d2adb61-5946-4313-bd6e-7ea714537656"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"f0", "[", "x", "]"}], ",", 
     RowBox[{"fplus", "[", "x", "]"}], ",", 
     RowBox[{"fT", "[", "x", "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "xmin", ",", "xmax"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "5"}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7682376421465816`*^9, 3.7682376439527864`*^9}, {
   3.768544739336843*^9, 3.768544758258254*^9}, {3.7685447936366663`*^9, 
   3.7685447959544687`*^9}, {3.768546782832716*^9, 3.768546786192734*^9}, 
   3.768546820592787*^9, {3.768546870700826*^9, 3.7685468872954226`*^9}},
 CellLabel->"In[60]:=",ExpressionUUID->"1bb2a22e-ee7e-457c-9fa5-3ad0259e7bb1"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwV1nk01N8bB3DKUgop2rRYskUlikTdx65vKAolZEn27JWKRAgtFBKSnaxz
P6mEEqFSWmwz1sxnECFFUil+9/fHnDmvcz53ztznvp/nfiQdvc2dF3BxcQ2S
z/+/e2Xj45o93WFRyOoEh1YmClAbasd33OHi8puzAcRLDTTFkxrdwXqe7/IV
4t3Og7l2Eh5w2n9wpow4OVujaqLVA3r/RHpwtTHRYQl6UETTC7wbvoflEw9e
vbO+iMcH7vjvOyPUwUQMUT1YL+4DZlyX5mWIL6R9dbyx3Qe4wpsK9hCLFusW
+Nv5wLqoJDsvYp034ypaFT4QGC9v9Y44Y7GOUbOHLxRfHL6QxGQim8gvfhMt
fqDfXZaBOpmod1RgVG/UD2SNNmbZEtuZKTqlLPQH0ZSrTeeJj6/zPKS/wx9Q
snNcBbEjNa6amugPXtwq31S7mMil79sPgyMB0JGtza3SzUR+ajOn03sDIUZ1
2nd/LxNFDS08bzJ8FpxfXWBNsZnI3UlWim/+LDRFwfxqmolM+41ePxMLgoAw
U6O9xKJd11Yp6wbBBjVT7WjirOZV5cvTg2Djdq73Ehwmqi5XHGceOgfq8ges
LQeY6FvYoeOONeehqL1Ho3uIiY5uzNY9kxQC+2uq3I3HmOhkhu7Hp8Uh8EAl
ySiI2E9ywI7nRQh4u1Ub5hHHSMuci58IgegFk/e5xpmoUi6fUWR0EQaMFece
Ea9VLlrf/+ciBO6W+awwwUReGZ3nc8pCoYR7pEF+kom2YLaFoGgY3CmbFo79
xUTG0SUV+9ZEQHXdEZ2di1ioxqR6dvRLFEiwSypWr2GhD68ix3/9ioHdgir7
v2xhod+vRMrqfl0DNWm79VN6LFTV5S0f5BgHIpomH5ptWGhfcfnBxqF4mD0s
FXIwiIW4fqx1MzC7BV313o1BSSykG+IwklOeAH0Lb/mFP2Shdg5v7S61JBhy
6TBb2MVC8TVvrszfvQ0jff/UmudY6O7pifxZ6TtwaIHYfK1CJ9rhM7W1IysF
TisPhclZdaLpT6vzOPxpMCLITDt3uRPxKqYFyMbehVtyHYZ8TzoRGgPVtB/p
UK4+tdx+vBM512fmSO/NAOskuzfdCl3om+otJ9e5DEicfz+7060LxdklpjTi
TNiRMGPBW9CFmKcNglFAFrwOqEq2/taFVP9bxPqokA0KXeG+gru6UfveypOt
3dlwd5mRiHdYN9ooKZ45kZYDL/QjVn5ldqOQZ1cPyZjlgjb/r2V+Sj3o8r4S
HCiaBzr9YYpnrvWgD9NPdXe25YF3qZ/DucketN32k4lQVD4o/9AxX2PZi3ZT
ZQsX7CsAiUpv64OvetGRurjtwnMFUBhy3Wt+Vx+qmnK4FVt2HxLKWrY6lvSh
xuV0hpRPIaQapR2NVfmEhsb8ZFuli8BfS/Y/94efUGk5ZZU+VASm6l9TzWT6
Uaz6Qu7wtGKw26l/2PtUPyrp1u3rMC2B1B/1jfCkH53xUnXRX1YKyr0mFXmC
bGRZqeraXF8Kn7/utvprx0a9P2hRnytl4CpWGjVTwUYW8rV/YBMDDH/7KUaK
0KhRdU1+RzADjG7uSu93oNH775Kqxy4y4Py6JjNVJxqxyjbXfAplgHZAWlLE
CRp9UdRiDoczYMe/yWZ5FxoJbTrOPxvNgDYPtVdunjSyEs112ZjEgGhXqci2
0zQa/qEs71rKAA3ryT3asTRa8sjo/u8+BlzJcUsVoGgUk6jteqKfATarc29s
fEAjgcDdcu/ZDKgSaKJUy2m0aIdSbs4AWc+1MfjoIxrxUMKZB74wIISzdCCl
kkZzJR3JudMM2CVpMztZR6PJ3BNXzJZiYK1w1/rcQiPfCFvDakHiQ+Gne1tp
9O2EJZ+cMIas9303WtpoNCFtdPmvCIbcF92Hn3TQaCxTMTR/FQbn/sn2s100
Gro7eeafFAbfbWuEnrJp1JUY6nJfA8O77ePpuhM0SrNg8whqYvjuVlSw5BuN
7MR0sny0MMAtPvUWYnbCwl51hEHLZ0mr9SSpx61I80Y9DKtMSjKOTtPoZ/zV
PZwDGB6WpKsLzNKowmy8y8AMw92+/PoK4nMipmcLzTFclxdMdv5Lo39xwg98
LTDMxizmrvxH6hF3U57rGIbS+gNOJlwctOJ68or1JzEodTy5tJaXg9pMfjNC
XTCIJD9zKiJOErQ2HXDFMCOQrbubj4PWXhOPLvLA4CFoMWLGz0FSV9PnNHwx
9N/OQq6LOWhg/3zqXT8MPuy9l4eJ85bYa3AHYNAZWJB+UoCDFGIlA16dxrDY
pUjTZgkHKcfkjFheIOs9V0YpCXLQ5D7eqMpgDCeFhnpSiMsXn9y04SKGW7O1
3XxCHKQeLXd88BIGadcLT5nEe68UtvlHYRB6y/CwWcZBXEZL/JhXMEii+ROP
iev4PYU1YzD8agtlCotwkH7Ulv8WXMPgNzwn9ITYOJJRE3eT/J6jtejYcg4a
yXa9uTwBw+2W7Bn1FRwUWSvhnJCIwUGdyg4lfv43TiA5GcPaEUXVRaIcZLt2
X++aFAwc/V4eQ+I/6tyM1FQMltN958OId/r7WtxLx7DEVmPnBHFLnIK8ZAaG
lUZczhJiHORTyv6TlYnBU3xK1JRY8O2d5k3ZGIJ1wzTPEheOmGXk5WCIlOF+
n05syC/gL5+HIYWR21BLPLCpTr8wH0PMWJQYm/iSzrnVSvcxNERmlv8l3mCv
MlpSiOGzTXOy6EoOqgr+8nRbMYb1vC8fyxMfSc2KwyUYlA9v5tcgnq6wdlIt
w4DYfhf0iW92LFd7yCDnFTW71JR424+mReoUhnr6X4k58VuR8O6KBxjcPsQf
PkTstk2zdPdD4oHi6QPEfCZTodWPMDyjJSONiLPdiw7trcBgqJrxR4sYrjjJ
Pn+CwfzQA8MtxL254r+1q4j16hzWEJ970frmRTWGwhmhPdzEq9ix6frPMJx/
/LB5gOyvfE7X92UNef5qMX89sfm6v7r7ajE8v1fX9/96TWiUr3xTh6Ei9Z5B
IPFVK88R43pSL74XKobEmwM3Vb9rwJCYGp4kSvzyZs/1gy8xOGUIOfeS8zrB
SHBoeYXBRM0lPpOY+53xjsNNGBwLnIUciNNHefg73mCIfvjonTgxSzawuPMd
qW+YDCeM5CVQb8vFYx8wPHGpFt9GvNxx0Kz3I6l3SZJ9B8mb8V2Lmf42DBsC
eTrE/p+/SsEmxw6Sp6NJIyUkn5GshrQBJgabipPvgfj5CjWd4S4M+w2T2y1J
vm23fxV178HA7+zI7BQm+TPN+zzai0GcV/OUFfGOmJXXJvoxTDWmiOuQfvmY
/+64L42BocUSKSP9daohUmWKg2Hzp+s2K4nvc80wfw5h8Na2Lm4l/Wmwoazw
7DAGjR7zGwrEA5ouwX9GMFgpxoQFkX7ecIYpNTeGwVSz1Zyf9PvN8cdevFMY
yv+5p7HI/MhVsvf3+oEhFuo/LyGu8FgU1D6N4Zt/QZQGDwd9GjlyOfcXhoFF
WZfDF3CQ0tDvFP05DAvDPmR/mKMRksnKKJnHYOzs29BG5pf5if/yxLgpCDvV
Ut5C5lsQnYIHF1Ig4TBQXvWH3E99mq8iFlPkfRd2wE8ada4faB4XoKBq1Y2a
ZWR+jtlcbbVYSgFrvP9S1xSNVnT39MkIU1BoNGJ99DuN7JnB0w2iFATf3rlm
wRiN/ryvkeKToOBBzuj1sU80EhR2lT8lSYG/R4W3QR+NJEyXbe2QoqBpFG1M
7qGR4Vt7jTwZCkRfPAuX6qRRwivugwaKFKxjjJ7/8JFG+fyFFqVKFESve/x1
8D2NKg3Mj63cSoGMdY7VZDOZ//VZJ4eUKXg4Np808ZpGW2t1gyPVKJj42hB8
tpZG2lxjYV/VKdBsHbDUqaHRYZRwxVKDAisPuW0LntLo/NOBW7JaFNh6W1k5
VdDo1ZOIwkZtCvbYf6/OL6XJ+/NLJp8JBW9pkaDVd8h+xA9WJ5tS0FjhKdyV
SO6TcFbG5oMU/BvI2Xb1Jo2ELUbcTA9RcH/4nmgtuY+rZwT+Jh6lIDyBx/Jk
MI1c7G/1yx2joPy6lkJuEI2WvxZveGJDgUi7vE9LII1cU5Su9x6nIHBL5/YZ
LxqJ7TGVkDlJziN21UoHWxo9z+3geexCweLATdJLjtLIU+j4iJEb+b/S9QqZ
h2lU98mb8vSk4ES0RMON/TQ6FRavV+5Hgd288U6vXWR/L9tcdS9S4NXxcaJF
kEb+yrYm7aEU8P+Tp4QX0WjDncHtLmGk3t0mcioLaBTo/vNPTAQFdYJx3rI/
2UhKcPW1llgKDFuuFzzvYaML5sewQzIFO+6J8PrmsJFcFSdx6g4FTmV+Zf1p
bNQi7XkuIpUCXkeX16qJbKQwfUH3fjoFU3vhfnQEG7XZ8sprZpB6nZp7GBfM
RqGN15Y2Z1JwY9VqvcuBbNRxO739Ww4F9ybj8wyd2egSl1xlWB4F9vmnKAVb
NlJyK0sXLaBAI2xweKMFGzE/qofn3adg+G5Et64JG4Xtfu6yq4iCFXqPJmr0
2WhLtpFxUzHJq+ayB/172Yi15KOyTSlZrzU6rKfBRpcDjop9LaNg3w2esTtq
bLStl/37Iqag+U1IZ7sOG/0PQpAgrA==
       "]]},
     Annotation[#, "Charting`Private`Tag$1132#1"]& ], 
    TagBox[
     {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVlXk0lXsXxwkps1LduDJLcb3mCP12ToXCJQ2akExR1zEPRZIpQ5kvQmS+
OOd5nlflGjJUSpxjPEhLpsxjOOEk3vP+sdden7XXXvv73WuvtWUcPM477eDg
4Khkx//zoGJiAu2OG6SNlLpqz9OQj84EA89wg08Z2c+fsFngjL5EWrMbWBnO
/TXF5uNO44W20u6gXVLZm71AQ+n5ejWL3e5grFnRKrxEQxekR8dF9e/COVd9
K5EVGhqPy5As4ybDw3LRX482aAgTOwWSEmTAJDnMl9h8P2vB4ak6GSjyGjo3
WDQkVk4q8bYlgwz3wx7dnzRk1DqvYVBFhqmN8pNrmzSUu9vIhObuCVNjQeMJ
HHR0PXLGa7HLC8I1NJ8e20VHg7N8s6dmvcBmv/7nRjbbWinfyuTyhl8j5f1m
u+nI7vc71qe1vOGc+9aaIx8dORDzms9SvcH+FuXfbAE6cvm6tHrGxgdqs6Uu
qIvSkZfOml/OoC9IJ4ejEXE6iprgumc+FQCqgsVrzmp05HZLUXbndgAcc1Oi
7VGnI4thk5Y3+wLB6fsoqZ7NYgPxB9RIgTAnpLtPXJOOXtAOVO7JCQQvEsOk
T5uOaiuV5/usg+DChWfbHvp0tBRmbedQfw9YNav8EsZ0dEUqn+SfFgJSxwzb
r9rTkXMuqbOuPASKCyO0lW+y9cl8s+V+GwKlk2bHN9kcI6cQlLgYAgfxVKnc
W3RUfbgYKzN5ABcy315ddKYjcbUyyWHWA4h2ss4ru0tHd3M/3yughsL+sW9D
z+7R0R/4yEVBsTCo2j7UaZtBR2aPK6pMD0aAdXRSWvUoHdWb1/6cnYkCuXpm
sY1GO+r4GDm/vh4DyYKj3pDajjY+ilKb1uOBr5zrG4OzA9UMeCgFOiRA1Epv
31P9DmRaXmnZPJEIl9w4V8wediCOVfHbZ6ySodFddzSopQORQm5OF1SmwPH9
K0FGBzoRY4ynUVcnDRiRNGtOn06UWN8avZ39N9xD2XxTXZ0o22+x+KdcBhwM
FpF/QupCWuQV1d4XmdDRwvF7W20XYg79VjTGmwVpTa/fyOt2Ix7lLB/F2GyY
fDUnL/apG6E50MxazQHFNzcbqbY9yOldXoHciVxYNPDvvifEQEuaybdct3LB
Po3zZFArAyXYpmY243nwUm1WZDa2F/X5nQlGPi9A+XZrMwe5D2me3dXfeSQf
mOIeuMqVfsQ4Ue3c/SUf3LZCyD3nPyMpGYm8xawC0Gs0nv95bwCFvImzVrAq
BGHV9mFm+RcUblqB+4oVAdMySiv54CDqYNaRtHuK4IDN09/Lir4i9RtD5kJR
xRDZYjjhJT2MjhNUrh2mJVB+48M4/cAIsmlKUBfeKoFk/9OCLw6OopqVm8mx
1FJosrsckqA5hpr3jObKkv8Bw6yJv+Njv6GJOS/Fbrky2FYiDTv8GEeUSuJy
zkQZNBdf8fn1chLFHuPifJRVDn5+ZT1DKdOo4gvpa69FBUjxvGSSqLPI/66m
y2kRCuiGqDIuXVxAl6o1XWnvKGDQ4CyYnb2EwH9i+csHChyLtw9yL1hCR7Uy
gmc+UeCM1LXaE2VLaKtiK4m3kwKaxP2ulaolVJzXUndykAIpcUx6es8SWo+2
2/uKSYGNhLorXQLfUebl2Dc5ilRIDpXKPPzwOxpcHRUjR1OB+eGUR2LQMjLz
crpOjqXCxEjLQnrYMqpdnMwnP6FCzIV876KYZZQ5M6fhmUKFgNTKou7MZXRp
+IelVy4VWtQnnfxrllFHK1+8TxW7X5m+Y+bXMnr3QoMncIoK+v9s+pZFraCL
So0skMdAtFx35AS2isKivayKFTCQ80gduPzvKsKm5EoED2PA0SX62r9pFQmU
Rl4cOIKB+uCcTlvPKnp3xAzz/g8GAfZCci0bq0hDuc+p6DgGDp90Lh88zURC
qvMd/FYYcObatp0dYaJmzYPFvcEY9DQ1cPapr6H27zKa1x5g8Lr0bOZnwzXU
Tz1aPxSKgaeYy+SE6RqaUTbom3qEgXxKO0PGYQ0Jydvx/nyMwVPfqgjx5DV0
WazQRSoNA/3suiS5tTU0taqm5ErBIDCKP37zwzr6TuhVzlExYKn5Tnky1hGL
bASeOAYl0b++rIyuI4F5a5ugSna9RbBdcmsdqY/7RcdVY7Cedv49v/YGCmLU
TmHNGLSpVl61K9xA/K9MSje+YqBBVTB0SWChmNSTro7DGPimUdwTcliIz/f4
4fYRDIrvcr1oLmehXVoqhQXf2P6mLdyvtrAQNyGc9+cMBpsOC9dkuH6y76U3
vZCJQXPywyCDwJ9oudAx2koAh81J9TYV8ibyjLhhXCuIg5Ns1YJ46CZacry0
87AwDsmORY9+S9hEi3Im4ZuiOGhrDV4zwjbRXJ5yaPEBHNq+dahZf99EE9nL
/r9kcZD3e28d7PcLDaSGupTq4aC/sEMpK34LZV0c4RbUx+HTrRK5c8+3kO0+
oxdkAxwc5eKn+fEtNJLCNXgM4WA/JWf4tmcLTSVHnm8+hcPn3NIFfslt9CMx
znDsTxzEwxYfS+LbaO+T9L2SzjjwvHl0iNeHA3rMN7BQFxxsu7M69odwQJrg
VYtvrjj0nY04pR7NAeLxEo/L3HFokhB+FZnFAbJxOVt6nuz5pvZFL99zgFpM
wfSl+zhUBSuQNcU5wSwSq09IwqEiVWxbq40TpvNdk/ak4DD7ILautp8TIhul
nVJScfC1PMF3YZwTGjYT+NLTcTjA6GrCtzhB29vz4vMcHHp6FL+raOyAQ/Ya
sxX/4GDG2rIyzNoBi3qV+1ubcBiMCd3Ggrkg7vKdabN3OJwLMDa9HccFR33l
a+nvcegf1iBrP+MCRyzlZtdHHLYpq6zd/3JBv6Jv+Wc6Dn/q3yoqYXJBw14d
o6kBHEY17aTGvLghaf71XZ4VHJzXLUrzQnigUMXe++4qDiWOmVM7nvJAlfuu
QAYTB9lA/5mg5zwwNG0TXriOw7s6Xg6sgQdUJjYyT2+x9d03lNLi3gnNX/U/
RuwmIChRdan76U5gtdfL7pQmYK3Hf0y8hhcEhV2V/pIhID2pZvMdnRekLURU
e2UJyH95riZ+lBeM2+z1ihQIUMtRcLvPtwtSPnJanlEmYGDkfnHftV2g2kgK
jtQhoJZbwaCbezc4EB/6dpoTkKnMqCz24gNBCcvadAsCJjuvG72M54OqR/25
Ry0JCNWp6N0o4QPhi9O3LawJ+IOxYKU8zAe1a3ybqVcIMCMvPDC24od9hhbS
Cs4EnG/sOqGHBKChsJf7tQsBxp3hAQq2AnBHyG7a5DYBshpvq22CBaBpyIO4
c4eAobxDURm1AvBXWOKpSi+C/T9JCmUnBOHjhx5X0gMCuDmXGnzPC4G32g1z
RigBUZbPNWL9hOBQxri6SxgBrpYfA7ieCYGv2w9WTAQBg2caohXHhUBW8Lf4
rlgCdFXdZsNDhIHmm+vpGE/A24TZ0K1CYQj4qnTpxxMCqk3TFVdowtBO1ZMS
TyLA2mV7iCIsAvfPX8NvphNA8VUruOEoAodrxlJXMgjQayjLMA4SgS65O0ER
z9h+uZUNwhNE4AjzPqk0hwCRtPx4njci0Pt3DmOpgIDh/z7hrpIUhYcch6vD
ighorE8wb9YRBZXb1ByxEgJIrStXrlmKQtjxBhfdMgI84gSFrCNE4Y98E7NP
5QQ88rSI+JorCv38nWrXKex9aas3m9aJQrjPlX0LVAJYLsyXHQOi8J/BkY0H
OAGxEiqjwixR+B9SU4nC
       "]]},
     Annotation[#, "Charting`Private`Tag$1132#2"]& ], 
    TagBox[
     {RGBColor[0.560181, 0.691569, 0.194885], AbsoluteThickness[1.6], Opacity[
      1.], LineBox[CompressedData["
1:eJwVjnk41PsfxUfGkpkpSfeilCU7UVSWup93KXFjoq7SVYohstxkC8mSkiU3
e3aTsSTLfL9f7iXNZMlSpKQskayJEGVffuP6/XGe85znvJ7nHFm7a6cdNpBI
pLh1/d97FWNjWlydwVwlmkN9yEFeB0ba8RRncJWfmR9cz1Qjg+1JDc7gp3ww
qzKZg/QdvuTayLgALp2i91cqByWz9J5NvXeBkP0RD0cyOOgPmcEvWwzcwK72
oeCOXA76cj9FupDsDiVji6qpZRyEiR8D6e3uQJ8SOZr5DwcFpH+3e7DXHaxD
NV6x/uUg8SLDx5427kD9+EOptIKDjjZP7jtU4Q6fOtJUhjgcxNx41LjF5TpU
1c/Sw+s56ELYN4+pNg/4dMGp9GsnB/WOi4wfG/eAw/4ljdofOcjGQo2Ryu8J
ry1yHoR0c9ClHa5njut4wvWqrQZyvRxkR0xqpyV6QhVPj993kIMcP0/PGll5
wY7yEQ/vSQ7yOLDgk9nrDYT2+4IKMhfdG+G/aTbqCwHDJxoddbnImaEoJ7jm
CxJqgW5kfS6i9xu/er7ND1pUfyZlG3CReHf0r1qGflByUsxr6Dcuym75tUws
0w/OhRB/3jjORZwytcnOM/7gtdA29/U0F03fPnPJruomxJ0aM538i4vO72IZ
3kgKBK98Dc2YQi66wjR8xy0KhHmV+id3irnIQ3bYhvwiEPrtJVj+bC6KlFfw
j50KhADFeyIupVxUqZSPFRoHQcTvzIqLz7hISqtQun85CL7umycim7nIjfnx
Zg47GPpslzXJk1ykgQ9Y0sRvw61DNeqBOs+RaURxhYnkXVh9K8ooGX6Oqsw4
K+Pf7oFIQXr0GKsKtb4Mm1xcjARL1tJrrnE1Wnq5hV27GA0TwoXHLrRUo2fd
15T97GLgdJuDtZBVDTIpKjNvGIkFicvlrdFrNYg0K3XVyCIenFbOuk0X1SLD
QNuxnLIE2GNeGlPn+AK1DwnU6B5IgguuASfCjOpQbFVz+FrGQ9ANZZx0U61H
GT5T+SvyKZB32FD6i24D0nGf2dORnQqSg/YNY9aNaK5PIm9IKB1UbUee9oe+
RAJq6V6KURmgv2r0anDgFUIToJ0+mwkVFzblSBk1I4e6RznyvzFhes5mVXv8
NZrWjmc48Ziw4qJhcpf9BsXYJKY24I+gX0HpUevhVtTpY3QLeWVDevl3UWX6
O6T9u3DXOxUWfCdvCdlj34baf6u88r6HBSbvfpxwiH6PdslufzSVngMmN+bV
A79/QIHP759RsMiFTRVW0RJWHeiOSTHuLZ4HcT2F0uZCXah1jmu4/0MejD9k
TIf1fER7L/aZbbqXDxKfbhQzuT1In2DzbzB5DFU+ENlE+oysamP2buY9BrSV
7/WIej96NmMbH8UugPK4lrQJzwHUIDbIlHN/Aio6PcoFpCE0MuGh+F6+ELaG
C7RfYQ6jkjLiXOZIIVwbbn7OyBhBUQf5+ULTi2DnvIWQ1tNRVNxj+LmDXgza
gxNZ+8TG0Q03bcfjoiXwU4UXas2bRGcrtZ1a6krgWO7JldG4adQ7OyjuHs6G
oC1k4bTYn8hSuWYZdmPQpWkoozs0i26He1jkK2AgErXAVzAxi7BR+cc0JQxW
m3gaMvOziFoQZtmtggHeKRYgu3EO1amYYp6aGHzzU4p315xD+9Q6HfL0MTDE
Gd3Hb82hTXsmWykWGNDzvrIXpOdRg7ZkfsctDMoU6qpHvRbQ2x+y2tZBGPiY
2geKhiygLrZqVV8wBgN6WSlHohfQN7VDnaOhGIS3Phatzl1Am3ZfElqJwEBN
1EpKs3MBnRPPddyVhIEL3hX/2WARjc5qKTuVrP8PaH/TQl1ClH+NC5Y+YxCT
HTfj/XYZRSYecbLvx6BU4c42895lJOKtr/R2AINg3+kZnfFlJKyjnpszjAHb
kyq+S2gFkYnNj059w+Cgo7+CD1pBvOKO5Nw5DLzI5wb1iBX0M9c+3IKKw1XW
VLMbcxV1JwY7FujhYJbT0nk9m4fSLQfINAMcqFkOLJFSHrLZdjTb/RAOB9Oa
yPgLHhpI4O89iHDI1xmsV/zCQ6PxYacbjuFws1Jnulx5Dc3H3j88dAoHHOxY
gcQa2vp38lbpKzi0CoxvCHciwQezJSzYcb3vc/qDe50ESbQ/6cNOOAgeZPot
+5NAKnp7RKELDruN3Rj3o0kgdz+Tp3cdh0h1ZtnmUhJoReaMnQ3AIeeam2/6
GglMw7CqmDgcvHb6xItl8cEYyylOLAGHxXKhizpP+CCsRsYhIRGHG3UVrZf/
4YPq1RiR5GQc5M4zJj8288F+z+uWWZk4JDCLgulLfLDz8r7x4ic4HHqL+5qd
3QBTemW/NNfiMJy7CCKS/HD/nOuYad06v8PfLUOBH1S9d3Pe1K/vq5DkD+/j
B3sswbbtJQ57DoQ4Z5/khy5F76KPb3CYzmrWEwzih+qtB46OduNQmzAfe2SU
H+Imy90EZnCQzCoQ46snQ676ZU+3WRziC2SbLNrIUOEi7Nc+t86X/H2G3UeG
vjGrO7mLOLCzRSRjl8mgPrKUepyHQ8AgTeqHlgA0fDZ4eXcjAa7zpq5aTAFY
flslJyhDgEFUCkUyQhDsiMZOQTMCdC7NfO+5Jwy07eacZDoBkRHl/e9ShaEi
tIupak5A1E/ajk/FwrDZcuwq/QwBe20nHmh8EAbOgshq4nkCDs3Rt9TLboRt
h+kyClcIaJLVEG2s3ggvGz84GQat982hIvk0CnhqXTRrDybA0+7azBZ5CuxM
+bLX8TYBm5WmgqN1KeDtPL8ceZeADLFf8ssYFJCjSUS3RRFwd4pOKa2kQMBp
a9w2mYDesN2lmDMVlJ4NJc6krPPbT5Y0BlOhTd7V/24aAQxKU1pvEhVU5gIM
CzIJ+J8l5cFKLRU6Hma2T+cQMPN04GiKFA1CSEqVt/PW/8qdw49o0UD9KjtT
/DEBmVSj3prjNLitX+2oW0iAqizDWtKdBhosY9OmIgLMNl872XKHBl2Ud1oX
SghoJvnuoqTQ4I7X+W3f2QRMtSpGpBXRQLN3YCkIJ+BBjvKTsec0+A8FF23k

       "]]},
     Annotation[#, "Charting`Private`Tag$1132#3"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{0., 4785.49990665102}, {0, 5}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {0, 0}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.7685447600793853`*^9, 3.7685447966506457`*^9, 3.768545791496478*^9, 
   3.7685467931681166`*^9, {3.768546874732049*^9, 3.768546887784131*^9}, {
   3.768548183759201*^9, 3.7685482229464273`*^9}, 3.7685604146813912`*^9},
 CellLabel->"Out[60]=",ExpressionUUID->"2b432197-b4be-4234-b069-ea0ceae093db"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"p", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{"1", "/", "2"}], "*", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"x", "^", "2"}], "-", 
     RowBox[{"4", "*", 
      RowBox[{"muonM", "^", "2"}]}]}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"delta", "[", 
   RowBox[{"x_", ",", "m_", ",", "w_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"p", "[", "x", "]"}], "/", 
    RowBox[{"p", "[", "m", "]"}]}], "*", 
   RowBox[{"m", "/", "x"}], "*", "w"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"resonance", "[", 
   RowBox[{"x_", ",", "m_", ",", "w_", ",", "p_", ",", "a_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"a", "*", "m", "*", "w"}], ")"}], "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"m", "^", "2"}], "-", 
        RowBox[{"x", "^", "2"}]}], ")"}], "-", 
      RowBox[{"I", "*", "m", "*", 
       RowBox[{"delta", "[", 
        RowBox[{"x", ",", "m", ",", "w"}], "]"}]}]}], ")"}]}], "*", 
   RowBox[{"E", "^", 
    RowBox[{"(", 
     RowBox[{"I", "*", "p"}], ")"}]}]}]}]}], "Input",
 CellChangeTimes->{{3.7682371323150387`*^9, 3.7682371468950562`*^9}, {
  3.768237196594181*^9, 3.768237231779109*^9}, {3.768237445770013*^9, 
  3.7682375674153843`*^9}, {3.768237656319708*^9, 3.768237776733774*^9}, {
  3.768237831092414*^9, 3.7682379940058627`*^9}, {3.7682392100197325`*^9, 
  3.7682392615389533`*^9}, {3.7682392963219557`*^9, 3.7682393604495034`*^9}, {
  3.768239612446779*^9, 3.7682396590901623`*^9}, {3.768239738097948*^9, 
  3.768239740696974*^9}, {3.7682398263849115`*^9, 3.7682398281471663`*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"6c45ba6a-a2e4-4dfa-bf7a-543610b4d110"],

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{
   RowBox[{"--", 
    RowBox[{"--", 
     RowBox[{"--", 
      RowBox[{"--", 
       RowBox[{"--", 
        RowBox[{"--", 
         RowBox[{"--", 
          RowBox[{"--", 
           RowBox[{"--", 
            RowBox[{"--", 
             RowBox[{"--", "--"}]}]}]}]}]}]}]}]}]}]}], ">"}]}]], "Input",
 CellChangeTimes->{{3.7683964786836996`*^9, 3.768396481869186*^9}, {
   3.768544862473649*^9, 3.7685448725247617`*^9}, 3.7685459358385644`*^9},
 EmphasizeSyntaxErrors->True,
 CellLabel->"In[64]:=",ExpressionUUID->"04e4efe9-ea24-435a-a90c-3bf1b8c0b336"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Abs", "[", 
    RowBox[{
     RowBox[{"resonance", "[", 
      RowBox[{
      "x", ",", " ", "jpsim", ",", " ", "jpsiw", ",", " ", "jpsip", ",", " ", 
       "jpsia"}], "]"}], "+", 
     RowBox[{"resonance", "[", 
      RowBox[{
      "x", ",", " ", "psi2sm", ",", "psi2sw", ",", " ", "psi2sp", ",", 
       "psi2sa"}], "]"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "3000", ",", "3800"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.768043080378785*^9, 3.768043114235282*^9}, {
   3.768043179283389*^9, 3.7680431802817307`*^9}, {3.768043605357559*^9, 
   3.7680436055181303`*^9}, 3.7680437130526094`*^9, {3.7682393713483644`*^9, 
   3.76823958358294*^9}, {3.7682396821684923`*^9, 3.7682397127018566`*^9}, {
   3.768239743711955*^9, 3.7682397922361774`*^9}, {3.768239843480197*^9, 
   3.7682398800733356`*^9}},
 CellLabel->"In[64]:=",ExpressionUUID->"44f1aff7-c138-4f4d-86fa-5ca7df87f6ff"],

Cell[BoxData[
 GraphicsBox[{{{{}, {}, 
     TagBox[
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], LineBox[CompressedData["
1:eJwV0ns01HkYBvBhlnJsKSypDlvhOMLG2ljJ+5ZaW7osqejIuEWtg9wVIXLp
Jpe2RXYSdr8jnSjGmHGbGbmMIsOw5bZpkWs71t2E/e0f73nO59/neXd4Bjpe
kKfRaE7U/Z+Szl3ytMsEo04G9etvYENhdPnja5TvFbkV397GhhC9wwdklF04
a7fohmxQjfCJmw8ieFajScPYjg0OW4tWJ0MI6vgE5/jFs6HN3UT2LoLg/YQb
0dzPbGiespwtiyUoyAqLFslxIOuBKNM0jqANC1cGVTngC+e+LaG87eoDPTVd
DiimRwU/uU5Q2rf/IP8HDhw2r5XmJRA086cV2N/lgPCq7dTdZIIbPiRLnHUq
oWbdiRHfewTDzgWo0Vy4sOo5pyGlHHoqw5IRyAWszbWLTCPI06Y3DCVy4WXo
BCslnaDFsFaedxkXXg+m+BVlEkwWiOIObuJBX1W9dPxXgl/Ej/4uEfNAW9Nv
R0gWwV52Xs6rcR64B6s6yiiDqcabf+hVMGTgUaacQ5DmVDAk2FsFE7+shhnl
EjT1zSvIYVbBcuD3Mv88gjeiN7MyI6rB+tX73fOUJ0uNbXMzqyFGP8U15jHB
Oa3K0p6SapDv765JzSf4203bINloNSgdDY0rKST4pjF808D5GtiiW6IwzSJY
ELvgWO1YCxZvdTeFPie40Sc1U9mCD3uOuFvNUFZbm0kcs+ODAe+hd/ALgv7r
i4RjznzYkqvKvVxGcCfPus4xig+LDJqnP5ug006l/A4BHyo/9pb5cAk26+RN
fjotgFJnzYFhyvzTxpH6lwTAEjmuv8AjmCL9UJ4cLYCsYpGrVxXBa6wEL1GB
AK4EchTcawhmD5opuswKwGoh46yLgOr7iJuFMVMIZhfbrr+lvBDjP3uoXAiG
75SenhUSHO9yVk9rEcLWqjjamXqCGf6pJ+oXhCCLCShybCAYpGSh98KpHqoV
j8rsRQSZYQ50D62XsF9D7pGNmCD7YovwvKQBDMV/HK2iLFcRUzU93QAad47N
W3QQjPdqTpCoNIJULuu4aSfB2kDeldvHGqFwyuTzri6CJ63Ee/8WNcKXDa4u
Su8Ibk/rMQrtaoKlWHmFJMrPNhzQyl5sghErVql8D0HFHhNt3e3NwC+dWfeZ
cqJrWHiIdzOE5t6q+NRH8Ot4v3bD5WboC6lUk7wnmFCf+KhmXwuITNzqfhqk
9kjy1jbxaYGKMbpfK+Xx8fX6+uktkMY4KWz8QNDjasZmzmgLHLIfucwdIiju
zhggzFfwdKd6G3OU4PCj5/YMg1bI7ude2TZG0PuGi2W+ayskZTH0sijn1zEj
mOmt4L6xODp9nOA3n/yTO1ZaQX35wO4bkwQzP0765f7VBlHiwJRLUoK5Sw4V
bj7t0GNj1CClPND28xP9W+1g9XRULnKa4KrJ/Ts3n7XDYpJHVNK/BNvMrb7a
s9QOodanAgpmCVYL1BQiH4ihs0il2HCO6rct4lxqnRjMNF9/fE75uKww3GtM
DNPThzz489R/dWfuYkAHBJC9Tv2LVB8M7eR9cx3Qqj6T5r1EUP2IivqIbicY
xZe0TlDmdARbxp7uhAlXA7vlZYK75R3asrmdcFF1q5XmCkGVxe+CfkyVQFNs
dziTMp2bsN+mUQL6UxlleqvUnqbmvV5rEhhuUjY2X6P2klrSfcK6wPPaio6D
HAsDbN2dQoK6QTDOdX1LWZlzxku1vBu+dg7PZsizUKrq8NB6sRsCQ2TyZ+gs
dNs+5Nsb+Sf8BzA1qkk=
        "]], LineBox[CompressedData["
1:eJw92Wk4Vd/bB3CUoZEyZiiUQlKGIuReJVJISYUMpQxFUpmSIlMJGSqFotoH
5+y9zfOUeZ5ChlD6RaQUlYzRs/5vnleuz3WOffZa9/291z5I2V4xseNgY2Ob
Ymdj+99P5UPqRl+tmMhafNih36sHDPV6utytmehPSOdmUrkHtA9KzD04y0Rn
zprlLuzuhs2IBRUXmMggP11rxustvCjzvnfIjom8CqbidfXegjgYtrdgFzyp
040SeAtC2hO2ffZMJMMZ6bgspxN4tHYH/3Zkos8SmW4xix0QWMzVduMiE1Vt
44ni6+gAds1eYbZLTBQnP+RqmtIB83tvslY7MVFEcNDcZtMOGFeraJa5zEQp
bn9mTua1g2N+tCCNHR4gERQZ2Q6f91ywVnZhomkjz/MuTu3wYTf3pPYVJqqf
/NuxbHM7vFEx4je7ykQrkAPX9tg3kLvznfl9dyYyjnNr+qHVBl+qtZ9xezBR
a+6qo+e52kDUgvEhEHvA+7/vsW9awS/wynlfTyb6trrn0T6HVjB4x3nZ7QYT
rR8e+zD1rAX+81f2s77FRMkCkkFbtzWDgMjTyn7sZ5k1OesXmkAvdXG5+W0m
unf7oUxYWxPQPfUhJ3yZqPx18Clx7ybwVLB5qH+HiditOI2duxphTXdosnIQ
E63zeHmh5VUDIKefX9KxBT1o0eu3G+Aa2+ntO4KZyF8gucHfogF65aUztt1l
ovi9Wyd6+RuA8C0oFA9hIjsfn4e7Quthr9znZq5wJpL1kN2/I7QO7Hy0f/c9
YqI5ydcrpmtrYGV1Ua3tY/x5LsTtd2QNZKxSi/uK7Ub57uSJqIH5uJ3752OY
6JpkTcIBsxp4UCgZsSGWifr0CM/VP6ohb5pju9lz/PudCedXb62GM9q+S4PY
1JvT7ZvXVgNb8EK7QwIT3R69aOkxXQVHBKe8PBNx/aJfXmqoq4IPKp9rY14y
0WrBWBbf5Srgulpn25XERBo+hmu8KyqBKtDZY5XMRAEC47Of6Uo4xla+4jM2
72v1nXeeVkJcREHGVAoTzRumzR1wrQTFNNYSP8lEVmnb+YekKuHUt9A4kzQm
0mu4bSgUXgHJdsc62vKYqEZB3KfetxzYq84mHs3HRmvZN10tB0vJq84t2Gf4
bm6Ksi0H/v4orqYCJsqtE5iu1S0Hv2OdGrVFTKSTVjvsuqYczDVPESWvmWj/
J312Np0yWL3O0o1Vy0QhB8Yux+eWgoOL837ZOia61VTZJp9QCpVNPmtTsH3I
ky3/BZfCjeDnTEY9EzFeOK2rMSuFkYUPA4mNTPTe65rWqaUSKB85pxvTykQW
nFVcJ46VwPUSB6GAbiY6Kmcz7LKmGKrPf7rwD5txRSSSOV8EAqussm/2MFHk
/p3eS6NFkGNucsytl4kmVUd5lyqLYGpa675dHxOZ9l2L0fUuAjel9UuHPjDR
aU7fX5XfC8EjueTz6hEmMrJ1zVsaKYBaIzXVEGxrDY74Ez0FIPQn059zFOel
j23L67oCyNNJ2bSEveX4KbdOZgHMfIwynxhjooYfnjbzTgXgJebQ0v4d92Oq
v5ngdD7ciFyXG/OHiRzUFT66ieSD1tdXjySmmQhZjuStW5kPSzqqbgzsaztD
fCsX8sB/5pRK1gzO4/cwXZvBPAizfpbZMsdET4+qnbyTkgcJCtvSli/h+dNn
eEZPMw+q67VSrnGx0LpIs5MfruRCsHRr8Bz2251cVVnnc+Gwj429HzcLhRc2
tseczoXWnf4yD3hYiLhw1ToccqEnpo5grWShIfW3Qkp8ufD1gsmLj2tZ6Muh
CwsaOTnAx+4Ye1SYhbILtD2TuHLgnmv4tWbsNZVzJeN/s4H9Y5bBEREWCmkK
/0/7dzb8er24qLuBhcK4L95Z9jEb3vo8tN0nxkIcT3Zrri7KhtjZ19sVNrGQ
115tnjWu2SA9Kfh6xTYWkvEyumc9lAVxNppP7mLP3hd0+tSXBfxtZ105ZVno
vig3XOnIguVplDS7HAt9sskOJCqyYMRp/91ZeRYyPaQk6P4iC8hR52Ojivh+
Ohh8mTZZoDJY9al6Dws5T75KUxvOhI3sbEZGaizEbjx5T2ogE3g2a+V3YXtl
x0gLvM2ED/Y5oSPqLKTR8cCavzoT7v9gqPJospBN//0QL0YmfFwMCjYAFjpS
+yKowz4TwsX15TsOsdDTDPNlV39mgJd24CMLfRZqIaTfWnzNANuz5UufsI/3
8owZDGWAOkO98/dhFppWT8vV6cqA4e1yPoKGLBT5R6/kXmEGaGiubDU7xkLe
ww0tNgEZMGrefHXQjIUOa5Zos4lnQG3NNSs7cxaSu64o8F4wA5KVNhz+hq2/
8WJzGW8G2K+wl5yxYKGRWO3amGX4+gX/Wtda4fvp0eUX6EmHT8IqCtrnWOhF
bXXcklc6vO+KHXl2kYX+6r7PGqlOg87jDmcsvFlI4WT1ezXvVHDsfIl6sB1V
35lLXk2FJdMBGdObLEQprT/F55gK8mbHJ418sL2007hOpYK/tWbQ/tsspGR0
esFAORWUnHjTZf1ZKLHkz0blcRoiAvM5ZkJYaPuD7gQ1OxoM8rnJh7G4f92j
xl5dosD0m5mkWBwL3TpK7I6ypcByExnzCjtq7TvjYAsKnO8aBmTGs9DD7Nb2
wCMUhJ2Osmx7zkLQIn6wQ56C1hlR3lWvWMj/r3Idc5yEY+qKbv4kC/1y9bgn
6E6CaaGp9tViFuIWk+qPecKC9cs00iawP5bZf3SPYsEbw00bXUrwfg7p+ZuH
ssDo49jipVIWEs96HKroywI9bt/SC2UslCbdqcPrwIK9J1laZlUsRA7s60lR
Y8HGn4sa0MRCzYmeM6IDTBiTTVFb089C9c/vchxUYoKoDKEegz193fJh+XYm
HJFK3LtpAOdzIvHe/q1MoDc80VR6z0LLe/2ijokxwWXlXTg5yEK3veRyX3Ey
4dc3+0PPh3B9E1xFBPpT4G/a1tM7xllINcT5x/qQFODdnexhtMhCvDot95sn
k0H15uq5Tuy8VLFqy/FksKi4dtNiiYWm5vndJ0eTgWEEfg7/cD9FfVWRHkwG
DfvekDvsJOoQOzCY1pIM9k9WPs/lJFFMXGPdJSoZXs+5VEmsJRFXgJK15cVk
cHmtxvdjI4neGSip0z+S4IiocpjTJhI17xfq0fuWBDIeCjxj2IE/LTiGR5Ng
YIcU27AkieSjxWdl/ksCg+crJ99Jk2hqX4xafWcSyPq8b6vZSiLDxz+rAwuT
YGhvQMRzRRL5Pgm8+DY4CcxyW9caaZOIWIx5EySbBMlVg2Jd2FpP3hT6yCTB
VPukrBWQ6EuB2gEP6SSI/LFe5zIikYCs8Par4klQv+20Z/gBElVyHlQL5ksC
tbjBwRY9EgVb/+ZTmGOAkO9kxlFjEl3huz5m1cyAt/rrTxw7S6LNJnpCG7wZ
0Oucb96G/e44+44eTwYMRFqePXqORFQvu9ETdwYM9yZfNrQlUe9Ec5LEVQb8
cdS8p3+BRHNGL5IOOOLrh1woRY4keh5Yq/H9FAPMGvK2KbmS6NeKhwNluxlg
+f2MYgZ2hwhHcakKA86uY9+98yqJrt8uDitVYoCjucGBHddIlMcTx121gwE3
xgYt5dxI5MQrceSzDAPiVqyIlvIi0Y2tYbeDhRjQf/jM33V+JKrP9CpXmyPg
bf6lwAfYw5lVzTdnCGiR8V696g6JPCSEmBV/CChjjxVf7k+ibS7LGk1+EUAU
dmtNB5CoTOdx2v1vBDjJmdzqu0siO5Xzccc+ELDAc3jpVQSJ/tzvUFCuIWDK
wyx4UySJ1giYGaZXEfB92GHtM+zznrqKipUEDJYHb3wcRaKSOxzqSmUEVHlV
a999SKKmfcPmhoUEhH4BX6cnJIq9AVHdqQSI1e9hU00k0c2vW1rVnxIQ/MVy
QzH2yYz4l+0xBPziCVA+8IJEezJe+Tg/JqD+cNuFYy/xfh7V0WdGE+DR5NDg
TOD1JB1QVA8n4E3r0+ikFBJdm574W+BPQGDX/BbhDBKlbdi7Zd8VAib+SGon
YPcbv/YQciHAQujQaZlMEm1wP9j405kApdMPQ1SySLRWp+t62iUCPvRu/2Gc
g/dntq8F7AlQH7DMv1dAoomEzfcqrQgY//T68Hw5iZCkYOqhowSs3XDD0quC
ROfc0064GhGw01jlyjR2cUBafpwhAdeKUx79qiRRkuC//KkjBMw+jBz8Wk2i
iw8+xJUeImCZrq17fz3uB/n0qVf7CdhyU/yeeQOJIt+t8BxFBOhldsf1YO9u
SM3fgX1PwqC8s5FE8zcV3Mq08edPq6xqbiaRuX7RwIwmAaIpnC9L3pBo0YYr
qGQPAVrvy7I12kkU75g5K45txe9dW4Dt2xWk4LubgBe+P77ldOD1bNXr0lcl
YKtZz57UtyT6V/hs508lApRXsJqf95JoHMoUY3YQYArnB8Xekehl2vHC9dju
7hK/YrE3Zv/kilIgoOC/KOGYPhL1dKiNR28nQLvI2/bBAIm6E7Z8puQIOOxk
OHv7I+7nSY1a5a0EDJyQeMb9H4lUpLnXVMkQ4Kr1AyKxn0QYXzqJ/XRN5N2X
n3B+X9Ve9NtCwEh6p1D1MIkcI0qCJ6QJ8H7KKDL8TKLM9W2hkdhr77hbd2Fv
F3UkVbB3mwgnj4xgl7bt9pPC/TFlsXvFGIm0U/ZGKUsSIPJ++7so7PTfe12G
NxFA1/z1Ef1Kosc/FRyfYnfGJFRv/0YiScbrHC5se98rjjnYr4QNeCo2EjDv
gFbvG8f1LknyvoUtvffTiaPfSZTsk1X6V4KAPKnsmW7sz+w2EWXYR1YGxtv8
wP0X2XonEPtqv8yQ6wSJUvYsNgpgL6+eDp7Dtgmvkv4ojtdP18n7T5LI1fnM
81Ts8lsXrz38ieslILl4FNvUXkNI/BfO74ZzQ9LYX46uKmJgX5d+9mNWjAAf
tQGrHb9J9NHBZlM7Nq9kKnsedmO9hgeF/YrndpL2FM5fc9nvu9h7fh49XIct
Rdx5ao9tXTkZ2fuHRDPDfBby2L/ICtVz0yTKWLPGixc7+GF07xj29F2ibFoU
95/PeZ9rMyR6PXZF6SN26gVVyQVs05u8rU3YyIizOmCWRN/VpWIKsd/u7nZY
PYfz/Fb3PgvbcWPKqsfYVO4mKh77L5dXusQ8ieglk7kI7MgJ/RPJ2BH3g9yD
sbf0bphRXCCR36S1pC92QfnXuHzsVUH35m9gG7CKtdFfvN4zqcs9sT9EhX2q
x+Y57A0e2Ne8rYKPL5Lo0tBz6n+vc51XlO/DTgyqP3wTO87gX4vtEolEj9Gi
d7AVVd9c/YZdLL+4MQS7QvyloNs/Ej2tdT39ENuU81rhX+z8xNfVidhfvh+w
CmKjUOdSokMatk83P/tadgptXGhDZdh8ZcOMGGxhn1XHO7AZKbn6mzgodFxz
Xewotlpk8HgKtkntyw3/sMsrK9YpLqNQ6hm5UyJ4//X//N2Tgx0rtqlABbt9
m7qlxnIKXXLQRsexLSyu3ynHXmtxcMIV+1NYWrIeJ4VqGiZqorGdysaamrE/
Hmavyftfv2w5K/SOi0LxgeeAE/fT8tPxmjbcFDqzTaNkJ3Z4SPfZz9iWaWEX
LLETfhhSP3koNJxeq1Hyv/4s2ItWrqJQy1NOgRjcz/rf3OwisaVFRJZ1YLdL
ZNwXWk2hvkqpLXw4L5/8t3ZJr6FQcs3Nf9HYnEb8lzR58f0jC5VsnLcHvkcj
KrC38ss1suG8CmWF5Bzio1B+0X+xx7Blhdn+nVhHodGY1NFp7CP/jT905qdQ
ubT9KWuc/w5+2cJf2MPsO1pLsc/onf/gJUChNwOHgzdtJsCZeicbJEgh1/iz
xWPYEe61pc+FKRQx5WIai+ePCJN9aLMIhczyD35diefViz4tHhJby8+i/jZ2
FmSb5G2gUNb9c+cub8P9vOLFaKsYhdRe8m2xx/POUqt/9UlxCv28tWrDV+zP
LkLK/dghy7KtXeUJmO4M8xmVwPVTsf8YgOelSMKNdWySFPos7jhQi+erlfIJ
TaUtFLordiKvVJmAZ97SFmnYuv4JlWdUCOir/Om1XYZCTb7+nxawzU5E5sls
pZBD1LpzB/B8N3VrURKVpZD2uXj3L2p4f/IOyS5TwOuXTRuO3UdAyJLQIT/s
1sb3Seb4PKnXG7FbxH6ALoWKAQG6PYHEzA4KXWheWZ+Ezx80W7FxfCeFal3G
t77VwfNCQ0uwS4VCHHrRt2Lx+ebuv0rVRBX3+8mGNT4GBOQ09pm0YVsXsned
xeeh0pkbkQ27cf247dYo4fNTwSdv1Ws1Ci3F/9AfP47n5+tdHCmaFNqyPH1F
mzmeVwdkJrx0KBS557CIsRPuL61ViWwH8f6wfgc9wuf3NbWfR+9hq6zqrBi4
jPdfoTTtiS6FunpGKXdXPD+FT7rkHcL1apHza3AnQPJ74PffBhTy5D3/dd6P
gF1Ph7+5mFLoXlBc2MNYAv6LboyfxvaMVS7fE09AdHiGwe2TFDq1MYLr/TP8
POTvQ4Wfwvm8JNul/ALPm8sCTrQZhX4z60uXpRCwX+fg1zFLCh25xPNVK5cA
kx/Elwt2FJIf8Gfb9JaAwiWJ6CHsOf64JK5ufD+8TzVt7Sl08ciOAz978PPK
zrAHNg4U4vrg7NTWT0DYVTdVi4u4/onplelDBNROHfQzvkwhUT3THI4pAjQW
PgtruFNI2Uf43zJhBrxcdbaiALtSoDHMeAMDeMT7Lql5UMgxktchQYwBPVqt
paqe/8uL1Bo9SQa43c47v/MGhRZuCV0vlWNAGvvd9C23KKT4Q9nUeB8DpHhk
D/EGUahfUvqboAMDbtuuvhaK/fRnr2nVRQb0lUw+4w7G/eJqd9TNmQHRVwt/
/cN+co+19BE/L3P0H0mYvEch+0vlOZ99GDCUevlPexiFvJ12RCU+ZEDSiWzi
0SMKedVEikxU4fenPWnle4zzaro9VLCOAdY8PnNh2Ms6XblQI37eLj1oHBSD
15PWsu3lGwYEb+2ed39KIbkA4fBH7xngMDd7/PQznGcT/RydWQbIJWqziTLw
PCoueTa0Kwlmqb0lL7Bf3C66OKyaBLUFql7bkihks9ib/FU9CWw75CdVk/H7
H0uFce9Pgmecwh+Nmbiez46VBxxPAj7nibIgmkKXfYc/h15Pgmn1F74/c/A8
d+IsWFuSBNW68ZpeuRRabfTlelN5Ejw0iZlZwjYzb4x4UJMEu5zDrqzOp9CY
VPgK2bYkuJjoZbOtkEKZOt3TzE/4+xTncbAqpdD4fLnf15XJUNnB/q++hkLP
9Vsvp51NBr0/kTl6tdjqHClf7JKhSVjyUjV28ak9XnJOyfDWUrurrI5COqDD
VeaeDCMj3lReA4V2/zZo1QhNhlULv08zWiiUztureis/GUw3j2T4duF+1w68
f00wBXp13e0XsSXV1JyixVLAynG5+M1ufD58+yVTJJUC9qnSdz16KMT58vxp
KcUUuKFmY3X5HYX+BNyfND6UAgkGvSvOvKeQZintp34zBUavN9ru+Uwhga8T
6ZFjKfBB35X1Gju08pGxzc8UeCshNKk3QiFz1SjBvXMpUFl37tbpUQoZcK88
t3YFExLE5p54jeH5udJD+IAcE0yrtrUUfafQwV9n6UEnJlSuD1KDaQqtaVF1
DZ1nQuGo/O06bLmVZXGjy1iQUfKm2ngG1zfaL+v4GhYk2IubnJ3F51fB4rpj
kiy4UZTtfGeeQoPWkt6zeizYZfvpZdUSzuuze9J8T1jwPBNW63PT6IAH234F
HRKOPLthWIe9XVvoVakRCbPB2WF6PNgbcm9Zm5FgYrltje4KGoULqYu1XSaB
i5t37YFVNOKL3wVXYklwsfzAq8VLI3mpibza3yRoc9/i3yVMI/UD/ZwVuRR8
+5l3Ih37y6u64DuVFDwdmHyoKEKjw21H/jvZRsGvzAsCOzbQaOTqmwaNMQqS
LQ0F5cVoFNF/aP28BA1rs8SEt2yiUdYzm4vuoTQMWhaJimyj0TmVbalZnqng
F2U98Aj7ZtOGy/PBqSBdy5GwTpZGp5flaJ+JSQX7nQZSq+RodHHIZ94xNxV+
sA9sW5Kn0fHwU6rFU6mwlLyoOqxIo9HiS+tcvdJA/CcyzthDo4bG5ScnH6fD
a5nPfDvUaCTsMT/XzUyHsxYhHSzsdW/OnxgsTgdGVftJQp1GFpsUPU99Soft
MbaWMRo0emo5wOkm8r+/lwVe9NGmEW8LcdztUgaYB9UF6OvRyC7S881DiUxw
cfbvbsO2f5diHSifCf4n9smdPkSjx91ZwZFqmUBKZbfa6dPIynnFzN/jmTBf
mrAh4Ah+3Zx7bOe9THg27Z5eepRGiuO2l43nMyHj/S4OXWMaHZOxNtqyIguq
q7+ZNmMfHH5XJiqSBePR5+b7jtGoRENA0GVPFuzbaaQ7a0KjRomCuSK3LPjg
sKVf+TSNhFZ2aJrNZsGvo4M7irAljA9I3l6RDVx74vz2m9HIdlna5RbRbNi5
nG/rMXMaDTpWh37dlw1+LxZcL5/BvlHouDM4G6TedXAxbWhEehBzVzflwIUj
d5QlHGlk4Dk0xu6dC4qH3tyOx/6tOV0nFJkLczqbmjZcpJFIi7eJXXIuPNhX
el7oEo0sV7wezerMhfxdc494nWlUf5/TVU4pD1YKX51hd6XR1glfJZvZPHjL
X67jh/38LP+ZO+vzIZGPN3IJW1B8/cXvCvmgujJVduEqjULTXVczbfPBavGL
+dR1nAc7i++N7fmQMXS2ZMSTRu1dv2YXygrA+2M6j50XjeJikYvVYAHovl8y
HcJ28WzP2/SvAPq6n30fvIGvr3enqAUKYVlj78Z3N2lE+eWUc9QWwqmMY3ea
fGnE6Rq8O3q0CCRTE1sO+9HI713O069riuEr68eGemwZ+QKTGtVi8CXCsqrv
0EjZnl9xdWAxMGPqh0oDaKT0OEpRQa4EFnyQXsZdGknl3NhhFVkK6grJm/fe
o1HQmejpg9Wl4Na/ir0Se6uIOy/HXCn82NtT3BlCo/7jif18dq/h4/RlpZlQ
Go1Pe6iK6ZeBeHLnWr8wGl3NYFTrBZSB+cm94zzhNDpjI5HhXlYGHdnLU0Qf
0OjUtxRNH4lyqHGNF9eOpNG2qgxmamI5cEiyz9di+xXxjeVUlIN2m32PcRSN
PrKxs1kOl0PBDuXoc9E0WpxPuxm/vQLIr/XcQY9olNbv0WhXWQERF6Z/NT2l
kZOEYwebUBU081u+MY3F+8mvOHnnYBXwVFWkvsdeCFcFDrcq8JcKd5yIo1E5
g9Us1V0Fbh82f1j/nEYPq+Nd+lOqISM8pDgeW3UxiDYZrIbvWhNPtyTQyP+a
yjkhkRqwjy86sSeRRkVByk1Z4TVgbna8wfwljaZVTE6k3a+FGO785E/YeQc3
PplqqYWOPPFAp1d4/0r2asysrwNDwS/atwg8n7rr5zxf1YF2x+2cF0k08pbX
rrj/rh687wxHySXTCKjD8+0KDbhfj1zJwh6ZkP12y78BdkUIylen4LxP/rvk
ursRNhvSiaMsGnWLLFgIlTaBqa7P0mWSRgqKp7L6tzRDsLah1RT26BG/lB8R
zTC2a1yUnaYRsyy9l3BtAVH5kht3sSXbAjRfjLSA4eaw3jWpNDrUp8mReLYV
0gUVYkTTaHRN9t23m+fa4OPav1MvsBMHFlrbvrTBep7mE9vSacSxf4dMiMYb
cJ93WqeagedpRNvd8tY3oDlIhRll0cimKy84QqwDnHtvfuvE3s014njNswMS
2g2OWGTTSOVfpjTq6QD26m/cDjl4vYK62SlEJyiXFtt/xz5r5xXazfcWLuSF
1lzPpVFZQRRTLOAt1DG3B/jl0SjS+Y+QvVcXzL5c+I8rn0bFon1T+/92gXx8
EwrHPuZ3gRgK7ga7ouXb7Qvo//9/+f8BqggZhg==
        "]], LineBox[CompressedData["
1:eJwV1Hs81NkbB/AhtxWixqhcopBKIiG78px+iqxQUmoq1bQVuXaztetSkpSd
CuMuI8b9e6YJIXeiRAoZbSpJQtNIZEohe35/nNd5vV/nj+d5Pue8jh4rwO2I
NI1Gu0jW//ewTgcr614KeWoNHHt59jmUircuU39DoQQz+0r9A92QPiX4O5b4
oEnm2K0WIUQqqnfN66OQ0RKzgAMbhLDT6HWE4lsKpa9T04637AIbq42vLxNv
o6/Gy6uegb59toVsP4XW8hbmGDs8g4nDfoOzxM0TQ+UrfDvh1ckO25B3FOIr
5LqFy3ZC4wWLxB/E9dabOpSyOiCOO+MgGaBQX7ZTbcWHdviLfyjj5HvS3+bR
rz6x7cCqbpocJTZvU7vnbdsOa1+y80WDFBootdU0gKewWDQm5TVEoUod5tFn
eU9A6vtO5nvid+kOHx9pPIF2hs7cvmFybmJwJVimDQLc+L5CEYWGcJbk6PxW
aF6wYtf+jxQ6EB6xRjuhBZYKs+A9sV+H+pJJ7Rbo9kieLxFTqO7cSk1X60dg
uog+HTxCoZGfO5t8HzfD1Z5r72U+UShQt6x3jNUMtvsjyumjFPKIYulxkh9C
ks7srTRi1zammpf1Qxh/cy5a/zOZP9wqW/XVA8hhBXiuG6OQSW0t7/CKB6Di
xZRx/0IhUXIRUnzbCF4rhCMvia+t6Qv25zRCg8j1+eEJCs1eKwt/7tgIQX6b
C05JKPQ1OjnUq+I+9J403Rb3jULiY3m+NwQNsH5dobXmJMnDoMZY+2QDxEoM
lmURpxVkuC23aAD7s5pfi75TSCu1Pyu0th74IXJpnVMUClINsRx8VwcKcOES
c5pCXmG2a+QEdcCiTfv3E1tfM09zCKkDxsWxjeMzFLpSaDf4i2YdhEa9GlKj
YTRInxFsaKqBfx09OpKJ5152n1fiVwNr53ZW6ElhZOBevFisUQOD7IdsM2mM
9AZ/YfqeqAbnuCLz7TIYjfzt5OJhVQU57qu1XxC7GvSPnBdXAo2RJ3dIFqOn
tuL7KlmVUJJ080WgHEaH+ek2pvRK0EqPCruhgNFKG3EMW7YCxHmerU+VMGrx
kQjva5fDdJgf30UZo9KPTdGzw2Wg5BEc00bc2iA7R+FuGRjLpni0qmBkf/fN
pK9rGfgdEr57oIrRXqvLxSrsUvi0aOt0NZ3MNzLhst3oLvz8zHxjo47RvDH5
yTM/S0Cl2buhkniH6eLWVcISMAm6HHWPgVFPtcJBiCiBgM4G+t2FGHkmHr9y
b7gYPl+1Xk1pYbSwaeD64+YioLG2qK7SxshpX1RtW04RqFp7fMkn3qvWwra+
VASmQ6fv5epgtK1XvNvRrghO2Ak2Z+li5GEZc4LReAfGpww9U/QxSvgtzHm0
WwDSzyw2LjIg+Q3wdrZUCGB+wSb9JOKf2XyJJF0A5rtZonhDjDT3Vt6x9hLA
6ZKbQTFGGCUVCmy407dB4ku/HmWMkZVyd4JnAB9+sNlfZFaT/uuyjb8gPtBu
y+2+QHxVx6dsdj4fFMcmdf8ywag4/sn1c+UYdM68LvIxxejxvLHTdQoY7ENy
ul3WYZTHnmJYVReCU/qS31qJW6U1RZHxhbCtNonrYIGRtHXPkIV/ITClor2Q
JUYVhgVGlG4h+F8K+GG2HqMx271B3VcKIIG9XkfdBiM34yXui33zIZV/JzyG
WOUuo77VMR8ynq4cUt6A0R2WaPTR8nwoUNMWyNmS/M728/IG8qAmXup/k4CR
Zef+rq5DeTB4s+XISztSX1EmWNc7F6z4+/AtJ4z+dEvSmMnOBv4FcabRVoyW
OnPkfa5kg+HO4KTbxJNBS5cz/LKBPp16sdoZo+OF2d90LbJh3PHlnheuGD0K
Dk2beMgD6v0eOTV3jD4sLnOb+pYFutoeB87vI/e9oEbeISQTEj8Pusvvx+jU
6IejmscyQaUx6PdrxIW0VZJF2zNh5niCRZonRhdjea1xBpnQUy6cW34QI9mJ
zV/92m9BvLt7+egfGOU0dS3wXH0LFP/ZrnbAD6OqpvKeS1IZsP7MGdpr4vm3
FXtV+rlwxDN5lOmP0TLZlKf29VyoM33btisAowmu/7Id57kQJAy86nKC5DUV
cLadxoV+nRtzbM9g1CwqqPqumA4Vd55ItEIwOsIrTTTamAbDKeMDKcSxuYsS
PxukgXoEo2thKEbeNmFhfYppELjLs4geRuYLnyi3FKaC4dSIv9IFjALiguuL
fVIhdpPy8FQERsOiKF4BNwW8nzv19PyD0RzLs/pX9ZKhamGP/zE2RkMrjjoz
FJJBlek1Z4KYl1Esbv6UBKWvIoyVr2N0s/EFZFQlgXR/dSjEYJQ7HGavwkyC
FPGaZbx4jOJao82zUhOhRYru65dO6q/qjSi1SgBtu0zaD+LHSlaxYUsTIDDC
ND6SS96nvtq6o8oJwJB3ruFmkPySPm6JfhcPLKVI1c5MjAbsolkmMfEwxZgs
scwl/9PNysTcMQ647I50vE+sUatxYk0vBzKT6b2ueRhx6PsDO1o4sEXLTN47
H6P6mOcah3kc4Ogd35NaiNGlvdy8st0cGGJNjiynMIo/eN9prj0HfuVFhpcQ
xzLeXv/TnAN9hlmFbRijH8ERG5NVOGDuZYaYfIxmVLx2uEzHQWR+bdcgsf3v
P100RXHwQuTsfeo2RuPON+Rp/8aBsfGrmVlieq5w23RTHPwHy6QAcg==
        "]]},
      Annotation[#, 
       "Charting`Private`Tag$2477#1"]& ]}, {{}, {}}, {{}, {}}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{3000., 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{3000, 3800}, {0., 0.3667515291749107}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.7682397534059963`*^9, 3.7682397927139344`*^9}, {
   3.7682398481955633`*^9, 3.768239880968974*^9}, 3.7683910661559143`*^9, 
   3.7683964843555346`*^9, 3.768542800650896*^9, 3.7685431297734756`*^9, 
   3.768543227688693*^9, 3.768544258591204*^9, 3.7685443333752265`*^9, 
   3.7685446624454374`*^9, 3.768545405955303*^9, 3.7685457964522715`*^9, 
   3.768548229600603*^9, 3.7685604172295766`*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"b94229de-f6a6-4f4b-96f7-f3a03ecef3ed"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"beta", "[", "x_", "]"}], ":=", " ", 
   RowBox[{"Sqrt", "[", 
    RowBox[{"Abs", "[", 
     RowBox[{"1", "-", 
      RowBox[{"4", "*", 
       RowBox[{
        RowBox[{"muonM", "^", "2"}], "/", 
        RowBox[{"x", "^", "2"}]}]}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"kabs", "[", "x_", "]"}], ":=", 
   RowBox[{"Abs", "[", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"mB", "^", "2"}], "+", 
      RowBox[{
       RowBox[{"x", "^", "4"}], "/", 
       RowBox[{"mB", "^", "2"}]}], "+", 
      RowBox[{
       RowBox[{"mK", "^", "4"}], "/", 
       RowBox[{"mB", "^", "2"}]}], "-", 
      RowBox[{"2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"mB", "^", "2"}], "*", 
           RowBox[{"mK", "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"mK", "^", "2"}], "*", 
           RowBox[{"x", "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"mB", "^", "2"}], "*", 
           RowBox[{"x", "^", "2"}]}]}], ")"}], "/", 
        RowBox[{"mB", "^", "2"}]}]}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"prefactor", "[", "x_", "]"}], " ", "=", " ", 
  RowBox[{
   RowBox[{"GF", "^", "2"}], "*", 
   RowBox[{"alphaew", "^", "2"}], "*", 
   RowBox[{
    RowBox[{
     RowBox[{"Abs", "[", 
      RowBox[{"Vtb", "*", "Vts"}], "]"}], "^", "2"}], "/", 
    RowBox[{"(", 
     RowBox[{"128", "*", 
      RowBox[{"Pi", "^", "5"}]}], ")"}]}], "*", 
   RowBox[{"kabs", "[", "x", "]"}], "*", 
   RowBox[{"beta", "[", "x", "]"}]}]}]}], "Input",
 CellChangeTimes->{
  3.7682396648567677`*^9, 3.768239765260335*^9, {3.7682416978732758`*^9, 
   3.7682417907539487`*^9}, {3.7682418554519696`*^9, 3.768241857059663*^9}, {
   3.7682419564249735`*^9, 3.768242000514103*^9}, {3.768242042934679*^9, 
   3.7682420544937744`*^9}, {3.7682421147277317`*^9, 3.768242146636447*^9}, {
   3.768242205581859*^9, 3.7682423308100376`*^9}, {3.7682423792465105`*^9, 
   3.7682423827242327`*^9}, {3.7685453661476965`*^9, 3.768545367009412*^9}, {
   3.768545620132639*^9, 3.768545648605514*^9}, 3.7685457311338654`*^9},
 CellLabel->"In[65]:=",ExpressionUUID->"a731069b-1b83-412b-b454-26db50b1cc27"],

Cell[BoxData[
 RowBox[{"2.9827682008970683`*^-22", " ", 
  SqrtBox[
   RowBox[{"Abs", "[", 
    RowBox[{"1", "-", 
     FractionBox["44689.96`", 
      SuperscriptBox["x", "2"]]}], "]"}]], " ", 
  SqrtBox[
   RowBox[{"Abs", "[", 
    RowBox[{"2.7873140061109193`*^7", "+", 
     RowBox[{"3.58795771727742`*^-8", " ", 
      SuperscriptBox["x", "4"]}], "-", 
     RowBox[{"7.17591543455484`*^-8", " ", 
      RowBox[{"(", 
       RowBox[{"6.793270969339969`*^12", "+", 
        RowBox[{"2.8114748180000003`*^7", " ", 
         SuperscriptBox["x", "2"]}]}], ")"}]}]}], "]"}]]}]], "Output",
 CellChangeTimes->{3.768545731553769*^9, 3.7685457986593285`*^9, 
  3.7685482313638887`*^9, 3.7685604172834597`*^9},
 CellLabel->"Out[67]=",ExpressionUUID->"58736143-4ff2-4ace-b27c-23f706b8596d"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"bracketleft", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{"2", "/", "3"}], "*", 
   RowBox[{
    RowBox[{"kabs", "[", "x", "]"}], "^", "2"}], "*", 
   RowBox[{
    RowBox[{"beta", "[", "x", "]"}], "^", "2"}], "*", 
   RowBox[{
    RowBox[{"Abs", "[", 
     RowBox[{"C10eff", "*", 
      RowBox[{"fplus", "[", "x", "]"}]}], "]"}], "^", "2"}]}]}]], "Input",
 CellChangeTimes->{{3.768242392305615*^9, 3.768242419546789*^9}, {
  3.7682424703938184`*^9, 3.7682425443890257`*^9}},
 CellLabel->"In[68]:=",ExpressionUUID->"6fc8f691-fb14-4e78-ab99-1dc2d62ef10b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"bracketmiddle", "[", "x_", "]"}], ":=", 
  RowBox[{"4", "*", 
   RowBox[{"muonM", "^", "2"}], "*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"mB", "^", "2"}], "-", 
       RowBox[{"mK", "^", "2"}]}], ")"}], "^", "2"}], "/", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"x", "^", "2"}], "*", 
      RowBox[{"mB", "^", "2"}]}], ")"}]}], "*", 
   RowBox[{
    RowBox[{"Abs", "[", 
     RowBox[{"C10eff", "*", 
      RowBox[{"f0", "[", "x", "]"}]}], "]"}], "^", "2"}]}]}]], "Input",
 CellChangeTimes->{{3.768239766453149*^9, 3.7682397666555705`*^9}, {
  3.7682425600979843`*^9, 3.768242569498862*^9}, {3.768242617189344*^9, 
  3.7682426786500216`*^9}, {3.768242747896909*^9, 3.768242771416026*^9}},
 CellLabel->"In[69]:=",ExpressionUUID->"cc0a0d90-49af-490d-8c27-e1b8c06612e1"],

Cell[BoxData[
 RowBox[{
  RowBox[{"bracketright", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"kabs", "[", "x", "]"}], "^", "2"}], "*", 
   RowBox[{"(", 
    RowBox[{"1", "-", 
     RowBox[{
      RowBox[{"1", "/", "3"}], "*", 
      RowBox[{
       RowBox[{"beta", "[", "x", "]"}], "^", "2"}]}]}], ")"}], "*", 
   RowBox[{
    RowBox[{"Abs", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"(", 
        RowBox[{"C9eff", "+", 
         RowBox[{"resonance", "[", 
          RowBox[{
          "x", ",", " ", "jpsim", ",", " ", "jpsiw", ",", " ", "jpsip", ",", 
           " ", "jpsia"}], "]"}], "+", 
         RowBox[{"resonance", "[", 
          RowBox[{
          "x", ",", " ", "psi2sm", ",", "psi2sw", ",", " ", "psi2sp", ",", 
           "psi2sa"}], "]"}]}], ")"}], "*", 
       RowBox[{"fplus", "[", "x", "]"}]}], "+", 
      RowBox[{"2", "*", "C7eff", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"bquarkM", "+", "squarkM"}], ")"}], "/", 
        RowBox[{"(", 
         RowBox[{"mB", "+", "mK"}], ")"}]}], "*", 
       RowBox[{"fT", "[", "x", "]"}]}]}], "]"}], "^", "2"}]}]}]], "Input",
 CellChangeTimes->{{3.7682427790974894`*^9, 3.7682429622488227`*^9}, 
   3.7685448390602207`*^9, {3.768545914397888*^9, 3.7685459390869055`*^9}},
 CellLabel->"In[70]:=",ExpressionUUID->"3223564e-c632-41a1-86b2-1263eda30155"],

Cell[BoxData[
 RowBox[{
  RowBox[{"dgamma", "[", "x_", "]"}], ":=", 
  RowBox[{
   RowBox[{"prefactor", "[", "x", "]"}], "*", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"bracketleft", "[", "x", "]"}], "+", 
     RowBox[{"bracketmiddle", "[", "x", "]"}], "+", 
     RowBox[{"bracketright", "[", "x", "]"}]}], ")"}], "*", "2", "*", 
   "x"}]}]], "Input",
 CellChangeTimes->{{3.768242980810162*^9, 3.768243023073167*^9}, {
   3.7682442257595415`*^9, 3.7682442795088024`*^9}, {3.7685453887732053`*^9, 
   3.7685453922469177`*^9}, {3.768545423906273*^9, 3.7685454277350364`*^9}, 
   3.7685456437873955`*^9, {3.7685457196475735`*^9, 3.7685457198251004`*^9}, {
   3.7685481204155617`*^9, 3.7685481246821136`*^9}, {3.768548206201191*^9, 
   3.7685482074268866`*^9}},
 CellLabel->"In[71]:=",ExpressionUUID->"5a733281-bb2b-4eac-bf46-072d6163fbcd"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"dgamma", "[", "500.0", "]"}]], "Input",
 CellChangeTimes->{{3.7685429465912046`*^9, 3.7685430409190083`*^9}},
 CellLabel->"In[72]:=",ExpressionUUID->"de858e80-3c44-4ea0-a374-ba372a5e43fa"],

Cell[BoxData["6.030442790985047`*^-8"], "Output",
 CellChangeTimes->{
  3.76854309618328*^9, 3.768543145064603*^9, 3.768543241926592*^9, 
   3.768544258752746*^9, 3.7685443334450407`*^9, 3.768544673532812*^9, {
   3.7685448196601267`*^9, 3.7685448444837556`*^9}, {3.7685454182085066`*^9, 
   3.7685454313813143`*^9}, 3.768545503597211*^9, {3.7685455481880245`*^9, 
   3.7685455732869415`*^9}, 3.768545805051238*^9, 3.768545884562702*^9, 
   3.7685459520143557`*^9, 3.768548129296776*^9, 3.7685482376221995`*^9, 
   3.7685604173732295`*^9},
 CellLabel->"Out[72]=",ExpressionUUID->"1893da70-0168-4162-8384-e0660be1306e"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.768545468760371*^9, 3.7685454914207835`*^9}, 
   3.76854552149939*^9, {3.768545556395049*^9, 3.7685455644196005`*^9}, 
   3.768545819343046*^9},
 CellLabel->"In[73]:=",ExpressionUUID->"d5c88e0e-f64c-4c7b-a72d-9d86409b15d6"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"dgamma", "[", "x", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "3000", ",", "4000"}], "}"}], ",", 
   RowBox[{"MaxRecursion", "\[Rule]", "50000"}], ",", " ", 
   RowBox[{"Mesh", " ", "\[Rule]", " ", "All"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7682430301472535`*^9, 3.768243057826251*^9}, {
   3.7682431094532557`*^9, 3.768243134577377*^9}, {3.7682435853432064`*^9, 
   3.7682436366859035`*^9}, {3.768243725722874*^9, 3.768243730282658*^9}, {
   3.76824379944876*^9, 3.768243800061116*^9}, {3.76824389602855*^9, 
   3.7682439722529855`*^9}, {3.7682440086985445`*^9, 3.768244048323577*^9}, {
   3.7683964521316905`*^9, 3.7683964531200514`*^9}, {3.7685429262236595`*^9, 
   3.7685429348515916`*^9}, 3.7685432550794444`*^9, {3.76854559178942*^9, 
   3.768545595049706*^9}, {3.768548251947874*^9, 3.7685482562224264`*^9}, {
   3.768548387532401*^9, 3.7685484618219843`*^9}, {3.768556212071046*^9, 
   3.768556214449185*^9}},
 CellLabel->"In[74]:=",ExpressionUUID->"2c1e3732-16d5-42ed-930d-8f72fb168b51"],

Cell[BoxData[
 TemplateBox[{
  "Plot","invmaxrec",
   "\"MaxRecursion must be a non-negative integer; the recursion value is \
limited to \\!\\(\\*RowBox[{\\\"15\\\"}]\\). Using MaxRecursion -> \
\\!\\(\\*RowBox[{\\\"15\\\"}]\\).\"",2,74,2,30453713367983442519,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{{3.7685484240998387`*^9, 3.768548462525111*^9}, 
   3.7685604174659452`*^9},
 CellLabel->
  "During evaluation of \
In[74]:=",ExpressionUUID->"2298e405-8520-4a21-b61c-4f7574b19f67"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJwU23k81NsbB3CRNUShbFkSshQJXXE6tqJEyU72hGyJrMlWyZKtss1gxjY0
yC5krWRJRKUsEYVCKUnW+/jr97qv3zXznXOe8/m8D3NF7D0NL9HS0NA0MtDQ
bP1v50sqLY1XAa5q1d9+sCcdFd8d/6SeUoB50wMmmtkzkDb7Ia6+lgLc8nS7
jY9ZBuq2bDl5bwb+mSOWTpySgT48j/kWwUHBIZ+nS0VUCYjw2/3LCxUKdnpv
YbjvOhHpf3s4qmxPwa2V6hqMsgRUlfv+0FIsBfvH+odZeRJQTpDBz9lqCj7w
uDpfrYKAbvKzP+YcpeDJ0heqIwsExHlz8kUEUyE+8VrQa68iESU92ja8X64Q
G66w37K+RkRpi8b2ayaFOP79RBrzYyIy2dZ6QjSsELsh7YTABSLa9t6H/zql
EJcrx5jJH8pE9l4n9Pf0F2Ke1V03Mj0zka64+q+fa4WY69XLaK/STCSUd/f+
8v4iPBL3xSLxWyYqmpqJO6tfhM9Y7XRWlcpCYr8TbvX6FmG3D4EJcpeykFGr
wUB2dhFWihuihuVmIZYS7fuFnUVY1iUodsfnLFRqH244tlCEg5nkXlxhy0bM
Q+8uOwo8wnHvDA8MSmcjQwHvIB6tR3g19afhp1PZSDRQi2bblUf44Fz9hSX7
bFSQcWj42INHOIxd686nm9mI/7ZNdWbDIzxg/8jHMj0baR2lTJ2feoRrbrM3
+VRno8UNiWg1Tioeo1xoM+3LRslJ81Kmx6h47sbKw5KZbDQ10y3aaE/FNwYV
7NjoSEiK5dfeC9FUPB3M127HS0K/T+qHCpdR8QW/+B3Jh0nI+1QrOj5ExR+U
cgzVHpDQGMeN1QjaYnzjslL9cU8Soss02MknW4wtm641+94kIe5FZ5FRo2Is
6Pw4xiaRhIqdnuR+DCrGhi+XLn8lk9AqQeSmWH4xfpFomXmoioSOEXpNCa+K
cVTEuSad5/A8Nil37JaL8ZMLX51vviUhMnogaSNagrXVtGq2T5IQ89SH9mjd
EtwpPdJy5CcJiV9QaqTzKcG6AeSEzlUSOhL1obwoowT7zb+pXN1ORs8+VOyN
bi3Bto1u3hasZLTjR97b6tkSTG+WYW23i4zUX40mcHKV4utrX2oZeMhISfnC
+uaJUvymTFCreA8ZqR48tZAdUYCT3VIiFOfTEVOVX3Hp2wJ8Xzal/PuNDPQm
IX8v83oB5sO8YeIjGahueO1y3gEKzu0Q9RdZzUR/PGjtb+pT8GjGvwgGEQIK
P1XDkOJHwcfyTvr+0COge5LZ3U+IVFypJSnRsR+e/1TUJ9kXVMwRMvth5A8J
HaQKET/NU7EM2+NIMQ/YD1Mn+td7irGSxShN8HUSos7dOfb+TgluvnjDNYWG
jFYfv7jEV1qCGW/KDo0xkZHzpIxf9vsSvPGUoZ60k4zErTUUHGlKsfrDed1r
XPD5aHnbmaRKcbrf44BD8Hnp/v0L6A0owA95cvaJjqaj/PE3ts1dBdhWEL/s
cs1ABksVpWd+F+DdxatGVs0Z6PR5rUIufgo+eu4K+XoTAenI1+7P0YLzLlj/
5I42AU3vLfCecqPgC9/vfdCHfLDi7fBmuE/FK/IfKzL2kdCuU98PF9RTMfuB
JtJPREI9/366s32hYr+eXj5DX/j8lVdDD7IV4/WZQ3dDvUnIRdgmaO5GCY4t
i3mrtkFCDwXOxdrnleBb4Q2hqwxkxGPz8xT3qxL8+oBhIzs7GaWK+0TGLZXg
kkan8YDdZNQ4tSrNJlSK78iVG5fCfrtGCaxkfirAu+S0DGWTMpCJdFrYCwYK
Fit7c6KBloDIVL1sSzMK7km9zakoSkDNyn7i7gVUXNWlfFr4IAkZnjoQX9lL
xa1SVhqRp3JQpLV+d9kyFbenPsv2dCehuhWLaB/RYszBEMMxHEhCttpqJ97U
luBS6yOGj1jIqM7FJV/6cwkeGmd9vYuTjCofcHxWYSnFTmvDu89yk9ETASX7
7YqlmCH8h20E7M9X2wto2qcA37ypGSf6Lh0t+YXWvXxRgCcOFPPG2GUg9kKV
Vbf5AuwjLxrwtiIDJe2Y1u3mpmCLP/1eRzwIyHtX8vaqExTMlPfg8fN7BJSx
LYsSf5mCq191SDsfJaC2xJHnMvFUnFYrfXObAAnRtgsMOlRTMd8DkUszR0ko
O9vcRGyMineJKB0wgvMeaf/36TvGYmxxLD6hyYuEXj7pJCz7w3ntm7JjXCch
xzf3s25nl2Cas1/GN+jJyF3A8PHrdvhnnme1x9jI6A0Lmyx5AfbHZk/QApzH
UAt363jeUnx56pSULezPk32dFwo+FuB3r3nejd3JQFyjQq+mtlGwr9lOD5f5
DOQtat4+J0XBZQrStLoFOag7TiLM/QIFe5X36mYIEZDE45iAZDIVcx9iMXUW
JyET+prKl7+pOMfHLs4Y9oencKxdRLAY47xD9WR/Ejq0Gma4VFGClQpOB3Yw
kxGnRBf71+ES3HK4J8aag4zmAgqjDehLMVsXp/QknJ/A+IL7tw+XYvZAHXVr
2J8ETVlHqeUCrPdBXJimJwPRpmQzXRSm4LfOyfkKF4lo37kRjURdCnYT/pfg
IUVAou1P7nxrpuI+71SZfZYktGYYdWtphoqX9xgoesB69tDPagTvKsbNceqj
7j4kZKMxNOS/WoJ9RmTCVmCe132i/0WLlWJ5iSXRRVivVP9dSnk7KNhQWEe+
kp+AJos5b4coUnBTteXB7GAS8u+62qR7kYL3H/J8ui5OQBZN9v0nqVTMvFOe
hUuGhOzbf3/9+5aKc/Wq3/u8J6Hymi8T6etUfJU42y4C+RNY5vef+M5SPDJ2
/XcxzOdy4bPD4SqleE0mu4kIn3+QK1TYx7sA32AM3XngTTrap2k6pfesAA+Y
lU+vWGUgd9qN7PXvBbjt10fNr8UZqHr6Pz7J3RQ8UL5T4KEhATFG0PxRU6Pg
N/4CWaOjBLT4J71twpGCKbevXyPLEdDSq8s/mmKp+NUpgUvlfCS0cltRprWC
inU9i1t+yEM/LMZ1lo5Qsff9mUjaKHj+/+qsZOmLMXtsdvEi9JXf7VDZkOsl
OCiC65nLGgmdkVUu3JtZgi+FaAiEwHzuXQ2flHlegscMIgid0BfzweZhMj9K
8E6TdPsimM+x3fl+vTylOGX729qjsN5jM9mh0oMF2KOvz2M4PANdU57w9t8s
wCULeZN+XzLQw82iFF1JCmb/GPFlNTkb2W7cv7p0joIdj1el58B8Wsxy54hn
UzFmtF/kOUBCKnc8Y16/pGJegUrT2j25qNnsu6TRAhUP7sATqTCvWqE5Tv58
xXjvaFlVrB8JcZmIWUWVlWBVpQZaE5jXBvo1Dc2PJTjszKTCNphXltk/ox9o
S3HGUEgtFeZ1SGesuF2mFEtenn18CvbrJ+Enc/2fAry+V0uevj0Difv8GqwV
pOCD7QoXbTYJqKL5tALfKQom3tWLr4V84LvzZkinkYo3/1jQduiQ0KX8hn9h
U1S8xvla7yrk89pwyfibncW4f/F+gPc1EgqulJCj+1eC//wQka6HebU+1Lar
W6QUJ88Rmfth/f72vR1SZIb9FeDMC2UnIDbVT6VvLChYdkLh+yUxAopbnV37
UEjF7+8Fa8RLkZDFHdbjJv1UvL6hzBN3lIz4d/CZia9SsekzZvUnsD4p+5/F
PWaFzzs3eisE5tPMYV2nTRnOh8i/mUT4vB6HHgZ47YX+lSMIoUQCGiEd4JLR
oOB7sWPhf20JiIVxeyT9Eyo+tf+d7e9jJKR89sDn1M9UbBN4+IFbAAkJkY5q
6wtAf5SpsdyF56+KtJIw3U7B6FZJrMByBlIWrEuYNAb/XnPMihUmoLmycKy2
RMWdlQrRZ+D5+N55BGkfKcUyBwslvOF5Rp/oH9u2n4LPEc/8NvxFRNGfPA2Z
9Cj4t+qOTSP4/FKMJIXYNir+Q3ZJd4ogoaTNC5xXZqmYX+huPyvMM6H+/Jiu
RClO66Zm7oLXOyZkkMDFTsFm++/G1EO+NPc2VYkeo+AJ9h4+r5FMNGGx/XWz
DQVveBUenjxIQOIRg4UVJVQcaGXJVCdLQq6MQ0lpg1T8X+4MayyVhHQdyzx3
0hTjKv629Rw4/yk/Y3jU1WD+DepH8+H9+jwecqpcLcAU8bBh8d50VE5pXWdv
K8DDH3iCqi0ykEDpp9/PvxXgaCFyDs2jDKTz3p72HScFf04s/tp3koAMjyu7
06pCv3oc5awXJCKRn3v74hwouPOK7pueQwTEMSsTeTeGihu7WHRY4LwvfVU3
uVoOz2fMHlsuR0KNgpRI32Eqjj+8wP0lDp43j586SFeMeT49b3sA69PjNvQH
+0K/dP6y4oPz3mjBVDJOKMEmMHFqcN49XDNChtpKsCcPW0cSnPe1qpuPFubA
F4RvgrfhvCd4aWy/yl2Kz/LcYNkH++1jwJoy964AS6bvr7MLzUCik8+e440C
HPJNeVl1PANJzvff/CtOwa/N+6qevM9CRdVdM3UGFJxiKb9bFObBh2VI+Wcm
Fc9Ek2s9xUhIaSqoKqSdiiv417u/34PXD6ef4P1JxW6qcjEDMC9D37dX7+ct
xh1R4nl0cN6jIv7Ynn9cgg+LD+byw3nf5cs1sjBYghmOuNf3ge/2HhBgztpW
ih8t3scJcN5j84JVnKXhvA0suxyH/crWfWMVsViA5dVYNAvaMlCtPX+LkwAF
/J6oQRomoA3/sEej2hTsnFMmKQBe45kapV9uoOKoBHqCnCbMb++7ilNfqfhp
t5p8C5zvJB0O+RD2Yvy7J/P4Mpz/6LmzFi//lmA397Rn6XDe2WNsstyES3FS
U0d4C6xffkjNoQ1GChZRcan2ZSIgs0VuyRRzWB8H+Q60n4DyecsTcyhUvM/H
6uFb8JXJ1N5L9G+oOENM/GJSMhn1Hj26a/YfFWsplSkkwPrkW32r8d1Rir9o
67I5wXlnak5766BUirPkS1aj4fPSCDbZ7dtDwd0X7BnTQwhI0sIq9iem4LnQ
Rb4hHwJiuMxRUllDxdcHTn5oUSKh2esnCNbjVPxuWdWQD/qyd723b56vFItn
kWd94Pkd5VGgEB0FH3548KLk7wwU/vYnU6ERBUvU1Pr/gzwnDxFWaP9Q8YuB
MonT8HydR3YwfZUrxdV7v911hufxovgJl4rA+0sILJplEtFASn5n92kKfh5M
w7sEfSx+ufji8VYqFrDO5512JaFTBz265L5T8V+hP/JUmGe9vfNu3w6U4rqS
6ko6eL3pLrmJYVYKJvN4dfwBnxo2TRN7lSiYZY9b4GGJbHRo9/FTt6wpmFHr
mL2zJAFhEhHdKKbil8z0S+PQ906Xpxp03sN+Fto7X20koY24roNvN+C+V/tr
0QbOu9gKx5FALgpWy1c6QLIhoIbOjuAABPfzBqLWymOYF/GDtLurwC9TNz0u
KJCQUec1jm+jVMxpSq1VhbzazhlemnGQgumTijXoC8FbJ63/yRtScDafCnsQ
rNeOw30e935RcReB7r8tL0vNaopIHCrFChyS+hfg892/1fGTRYiCA/USrV0O
EpEzE1e9sQ4Fm9+0TfkB+bBP4kVZWhP8/LLsecdzJPQ4mHyiYZqKi9wlLc2u
kpDOgcUodxYKVn9Tt/SVi4CuE1Q8ZI5S8JU5Sf5is1xkLUHDWTAA51El2XgN
7o9O45tE2zXIkyoP2m/wPBQGuwbR/0rxTMAV6VR4nuET7omtvBSs/eoGc04+
ASXNzGdc1qRg/Zls5nvnCYi94K+5SR3s/10WnszjJBS+f4dk/wQVN7wVJ8jB
+Z11jNgrbErB19+9PHEW7oP6BfJPHF7DvH+sWX87mYuORTnf9vtLxY9vFTub
wfvrdR04d0GMgtPDNAM77TIRlvjvo+ZZCh566FfkCPt9qOA+88gzKvbgDp5h
hrweffufad4c/PzUiegY2L9Dnx++8vmPglWL6Ey/b8tEHgU9zt8+ULETM2v1
2SzIxwwRFqJXAU6SlP4i8TodBVsuM3m1FuC9NQZiD8wzkM+tc4JSkN97ee9t
JBZmoAR9vyUZyO+hlhIOC3UC8ksxLHQ8TsHGeH8F7UkiOjb3/v1newr2cJDw
mJElIMGgeP35aCq21lPgquYloUUene8vy6iYsGJ3VwjyO/VIld3gEBWXGnmR
LBPBM5bdFw5BfjsxOwk4w7yvfGKyIMP9/u6kW8Qv2B//BS+5k5Df2gEtvgKQ
37zkgq4jkN+jgYIl1yG/I98kYSPIb/7RpFfXIL+nMxNOt8D9/vW13dd3wflF
tBpkfcjvPbZWYuI3M1BgneKP3K37OyHOTvZTBrp6afqaPuS3tEyoXYt1Furd
NXVAEPL7K3eGy23I77A0NVpjyO9bLU0XBSC/j+4JP/0O7usZ4XJ/jyiQkX7K
rubgH3B/p6kf+LeVT6SqiwF7i3Hd2dGcJriva/BZT1TB/Vxop//ZDbifP6UK
S+lCfl9kVeh+AvmdpX44cxvkN5LkfhsG+T3Z3xlaC/fzmJ3/GBRh/paFNG9P
wf3b9T+N8M2WDLTt4+Gqp3D/rjl1LNfzFQHtf0uYwpDfGhPlPVGYgLj2fEk3
gPymFbe1msOw/pMTCUVwv16T038VCvcJ643mwX64Xzu797M1QH5zstj9kIL8
XkaJ9tGQ305Pir80wH2Z14AmrBbWb8NI0swe8vvnjlKreHoCunFO3vEf3I97
7zSvVcE8VlC4rJggv6nFSd6BkN9Ge4RWrPpgvkuaHj8fIaMiXfUf+pDfz1jj
X4XA+sgGrOsNwn3X6klhtjnkt+79fmNfHgq2zSdTDlwnoLPVfi8NIa/T5jVS
f4YS0PYUHk1WyOuPV3ZxdSuS0IsYvYAWuI/SODF25NwgIUkZob0nIK83DzNm
u8LzNjCW5kTQQr/wuCW2/sxABx9p9T+UpuAd31JZD1nloQJ5Vu0dkN8sx4JU
2yCPngjFcjFAPt/S4wvuCyKif79bGA5DPttxcWbdh/wUS6n4ldRCxTsNvLcp
2pOQKAqiJH+j4ptP7eOY4L6We1K1QgXy+dWDDpo1eH9MJrGqQz6L25pEasPr
ZyitJ0pAPns5s8xEVWWjC54mEVNwH2tTDH7xWAI+rxwN82e4j+Xc2jYSAPns
mDJ3PvsdFVt5tb4d7CChQfIgwzHIZ7nQVD4M59v96gHxPrg/vdNVO/DDjIA+
i3+pHYL701ztkqncCwJaMzmc4lpJxUfXHosQjkAefAqd09vKZ/GQgNBbJHQ2
6mHOL7ifSPjbFnTtICFVU5clvw4qPkiq9T8zn4e+ao6eWYD7x4NCqZ9hsF8j
ayOTdvsomEdEdt1sNxHdczz8twbuBxU70rzy5QlIwLT75ALcD8JmiWw6euCB
xFIBIchjlwWx8FHI4+xVx/I+8L44b0lFAScBWZky5G+A53+cypf/A/OOl1iD
2sDzd6jV1b3wfgdqwm7xQt4eKj3SfC+DgGq1WQY7wOsfs5Ql2+Hzsr5fny8G
rwt2qPMoq5CQ1bG99EqQt13H49rU4T6/qyXTKNyEggd/uz2agvOr915Zv76H
iknRPP0l+floF9HX4SP4nElgg3AO3u+BhopWEXi8mDlF7uuBTHQpwognHzzu
+k3utQ74ROZ9euZRyFuHz956K/dJqPbbyaNMkLcxuReKX8J++Jk/aukCb1/P
fvw9xT4Tub0fWFeHvKWTfdBzLY+EZD7xy2eDd23rXDwaVwloufp4Cy/cX6d2
u/7Wg/vrJFNGwS64v+a/iLXuiiahT9ebYrIl4HluLW+6K2WjWxaPi83g/tgg
oSTzHubpWI2vsATcF5fcrinvEM5B6e7L9DXgR4+9V40q4fN4fmd4tRPuc5dX
eYcU5whoUMDJM/gk+KXZc07uGAHtrRn0yn5KxbEiV+k6tEkI13e8mwTfCYwk
dHOC78xqdvUWg7+mJbzEHS+S0ZKQhOLdFfj/7zHU5sDr1+QVubeDr0ouj+4k
3iEg2g187p46BRdU2KRpuRIQkzIL78VaKn4YLdP4Shn2w/5I0yr46v4DnuaK
QNivgLfdJ+G+dGdyZVYV9kdC3kHGTBT64+Syy93nRFR6POae6xkKjs1KwHEH
wG/KBS8nwEfcZ5PIdb4k5MZlN/8KfCTJuTvYDPqCWWegMEAZ1usSUrUKzkLO
yuMSC+CbBxa6ep2VcH/xKMF7TlBwgiTLJCWTgAKrnBTufKLickMZOrkwEqpP
c6/XkqLgKXe6+2E3yEiO71mi9G+4vz1xDbWFz3uGvWqtAjxyUeeh1Rd1ItKk
Lof8Bo9cqnbziJYhIOG/TQIazVT8z3SIS9uYhN5dfmBvPgPzLK40PgL5IHtA
PPAt+KM7bYelaCYZ2f79F2TwlorPcN4Y+DwF5zmu34Zei4Jf3OQLlzkD9x+/
hWcb4IlLOL78kBp4hH3jtNckFbMOKbUoQZ+0MzkPPwVPaNLGERrAE+f8pMVY
e2G+7T8FiHzNQdbXLAu5lqnYM0Sp0x6en86I02oNPKFxP+T2ekEm4usSzBwH
Tyh+p/+XCz8v53fP59ZzKmZhifNn6AK/tcgEasxT8SH/mj9GMN9fXMM/CKpQ
8MN0t6WiSiLy8nvZn/SRir/ML9ItpIFXjly+rQ+eCDdYuErpSUdjR5PvCYMn
2NkopdvBE9k/mY4MzxTgwkHB4iFKBjoaHeY1yEHBo1ZXNKoRAUVr27ZzgyfW
RJXY/lkRUdsVZcs48ETIxZo5P/BEbmtzRTZ4Ql+s/2fMlicOdVX5gCeMuym3
Bw+TEP3T4713wRPK2/u5UDL4oSpX4gNtMb6WLtR6HuYjIU2I9QJ4Iv7V4bhR
8MRb1byryxkl+FxD9F5W8MT8U1nv0dYS7Fxkdd8ZPPHSViCGHjzhnfQtzRU8
kTM0LLH5tgAblru2R4VkoPwD83Sm4IdtL6ljkqMZiIf0a+7fAQrmY7xNIEpl
oc/Ni5o9+hS8fuj2u1aYb6Zxj8uM4Afu06caf+0nIYXD9Nzh4Ier1HMWeZzQ
VyJphw6BH05yuAayw3rHrsXJi4EfEjhtXxBhvzMvJYk5gx+CckpTfoAfkl7m
9y6+L8H34j+2F4MfPr+unS6nKcWxR64nB0KeisWspqeBF7oXPljKghfyuVUO
u4AXmvd/uKYP/aAUUUT/C+ZtWZluYFGDgHYfdr60Wk/FFhY6/76eICE7T9dN
S/BCVmSc2BR4Qcs8jRoCXrBc/K2YCV4wSv323/BSCX74/NOBcPCCVPOYGxf4
wNo5zW//dgJK2EDhReCDC7N/wnzAByo/tVvqCqj49YtUfBZ8cOEw+RkT+CAq
svRnKVsO2jPQfpEGfOAXeHz+OsyrF+OdlLvgg8gAHesL4IP51yOvhcEHtL7u
pQLeBGTreOsrHfggn0U4PxPyh640ja22moqDZe1vmYEPFHj6p7zBB4buOQPX
Q0ho2vtyrSJ4gCsizef+jwzU81X/rAZ4AP0ZuBNyJRfVfeJ1aLhAwbWVQ5eL
IV+tLR/eKxOm4PDq3JadbkTESn4tNqZLwbqp7u0BB8FbpUpPEXhATGBHUbA1
CcW/3LZHGzxQfb5ORgXOuxubhcjvHZBfHTX8AYJw//JmZO1XpGCVFxdFjDng
viHoSEqF/t/F1VVkAv0/4CfvcR/6v4eDXHIS+t+h1EtAD/pf1/JAsmkvCWnL
2vjOrMP9rVuhSQ7mY1fx4exg6P9SIsNGjTEBLQ78ehcL/R80rvOZ7Q0BrXqU
2/FA/xuxb3ukBf1PcGG02oB+acqhVXpym4R+tWULZkL/X2I5k8bQnY2Oe842
74f+Nzm7De/+mIsu9DjWkaH/t40fGY2F/Qg/Nm3NBv0fw6i/Yx8zEWXvUuN2
gf6/1VritV+BgPg9DTiJ0P+93ft+XT9NQsv8SjG9U1QsEZ+qIQ7z0u0k2B4B
/Z+tkqzKyUFAlp53vYqg/08TXnzo2E9GJaHJMz7Q/zc2Gpqebd2ff8VsPttL
wfxps3oGDwmo8/eepUDof3WfdtVGK7g//gh9Yg79L3UhWzrxP/j3WzOvTX2m
4swehZ5+6H8xT6q3AvT/R+Ga9lo4fyohWo3G0O93qseQ7I5M5C9vf9QU+v2M
Ew+lA/pd+oegy1gbFdeEOO3PjAF/+L4iP5mF1+e2GJyB9RZ4UHXcD/p9JDHj
pjIhE135YfpybpCKf+37XTNKIaETcQcNzkO/s3VIczKxEtHfnsaY5+Vwnzz/
uW0D7jvnk0V4ng1TceJ465WzsSQUWd5Brwf9Tla6ZCS0MxulMjV2MUO/9z4o
vb/1+wTl18vf3rRDvlZccczpISPG5RAfZ+h3qxDPN22wPlwSnKfqBSh4jJuo
0TBJQN87r8XJQr93FMlf2KNCQHteWx7Xg373K2wWvKpFQsXa57QfQr+PitHc
uwz9bvr6Qbo59Dvrff/nAdfJyKGoOEcF+l3ocE1FGry+1XjllDf0+yJTUjdT
BAFxU/U9TkC/v+wavufgSUCMBnETTNDvJiN1fVrQ7zdna2WLod97OE9M2QeB
PyQmvv0FXzu2rWfwwvp3u7nn00G/J2mdtdQvI6LWmVQqP/S7uba43nFx6H8D
qeR46Hfev2Y2dHB+PmZOWd2Efg90PagXBvldEjG6/4Dy1v57H2obyUKXDRx/
EaDf225R3195QkJFAg3WLxEFnw1V2QguIKCpE3OkY9Dv+U0eXw6Hk5DtWz/F
nwcpOKJl8vYQCxll5p6mGflFxfMu92wvw+f9qdt03Qb6XYw1ZemCIhGZWDAQ
c6Hf5w2bzSugf4R2ak7NNVExua1K7YkhCWn+4u9jgX5fPERXwwu+fWMrtxgK
/b6wvbI63D4H2ey8cXJ1gIoTHqkHuiyQUOjx8RvVmhS8TyiKT0ufgHbG+8QX
Qb/XHi1stFUlIe60yCci0O8FAvunH0Pef6z2XfGAfic+DWyOgH42iC+Zr3lN
xRvHF4q1jXMRJZ7/XQH0eeAlhxifiEy4T0k3JEGfp4pM+tNA3h6O5zwhD33O
8MNVw7SOhK6eqZj+DV51V+P09YR5js6PvdDxHwU/Hkn7cOQDEXnGn8lE0OcB
bl/toghw/5jRiqmAPu42GnHdlCOi7oOj/8hwP76d/5agkJqFnjp+YHSD+3Hj
v9zZc7C/ipsbN6VgXhVaWAfkrMhonDX2fBf02bxTdNRXWF8OifiOeLifXud+
YqAOfc+9adebA/dT8SvyxvQaJGS8mXG4DPpgSmmdMb2KjFSYXItIkO/s/3Gq
lAcSUKSnivl3yGvxxUFPvqh8dD9iQ9wC5uvqfzw5w3BeuCbD0i3g/rbrW+qO
uFgi+qpdOX0D7m+X3+frN0K+HrgpJzkFeX2QVu3jBScSWlr5tScY7mO+v/nD
qp2z0aWbrncWIW/JX6QIX9tI6LX/+3QpmKfXr75qydUTEEUv9uApmJ+n0iJZ
fUYk5P/EdOQEzM/uDyyc/vD5gqpFPjDAvCQ5GKliOQISPKEbrA3zwtagcQHr
kxDfoNb0ewUKVghcpS9JykMXT0QMGcJ86CZWTz9jJKO8lLUne2E+zJ0du9gg
z882VUbshP0u2//ahCYnD/1kkOvcgPyysJ9OKtPIRIZvTy4vQn6d22nxtw/2
W7Zpz0AU3E/qmsL2lWfBPAqSTgnD/jL0F69eP5yJ3JvOKzyE+8jlsnYNbTIJ
qRyv4vwD+TPvKlbpbZ+NvKnk+ALIH8W/Z1LrYD3/y/plHAT3C32PF0fFonOQ
7MTVFHrog++7mxjz4fOaZ8XT00J+O7wpCEo/RUYHhcQ2lCE/jbRf/hMPhvlh
7TzcA+fV7cU5mnSJLMRj7aCeBr5e1PzUIUKfg9ri1fYHgq+7FtITLeD1RD6p
Gj4EP9PfmaxnNyOhqBl+Ginoy/6l4DW+TfDIp+tlOeBllttNpdqfIE9tKOa2
4N1eWsPT6vfAu3sZP58Df+mO+0o2vc5ETIIcoXLgL3a5uRfn4fzI2+wo+gDn
QSjZQZw8Cvedw/423uDZv0rbn3cmE1HM+l0tCfCrf4J7nTP4NYlk27fcUoCb
fYbCG8wykJL4JdUq8Ov8i/C2UPArA+W2bAn4tdzNjN9VjYCSBS6bzMPr8Wk5
Naq7EBGT9BTBHvz6dOZb8AnIj8tO5l8DwK/+lq3D17f8WseacRr8ivYWplWB
X1Vpb96xB78ye1L+/gf3S/eJhMQS8OujRJHKk5B/lktSl6TBr68EKuIHVrf+
/n4iqQ/8emcmZIoO/Pqhy0u7Evy685ZrpB341USnyPU9eFVz346jKuBVVcvT
7HLg1f5If/eCkQw03fi36DWsV9S51ysXdmWh2ZTVN/mwXuKxK00TcJ67IhnT
x4hU3NBQszCy9f2U+pMrpuDVbVfL8vbRkdEvp+EvDHC+t+/TMuSFPNk8GJjx
b08x/jOWOhYL+aUXI6OtDl7dfKKtMg1evfFuJqcTvEq3cbO1ALxa12ya5w0+
vX/kVeWz5gzkr3BiFoFP5c4R8z+0EtAJi7L7neDTfg9bzf+0CGhXfeHpPvAp
x8TOCRfwaWplcfVR8KmFYhCTANwHq2/f4jECn9a8EbaIB280/1H1m2Og4B1r
vDqedASU9labMwI86n7z1B0jOC/ryaO/k8GjTrslGpW2vs9RP1E4CvOk+P3o
yJGDOah57nPsEjcFS0e69O73IiCX47G+w3B/LI3PcxeNgfvu8aFfceDNlQth
bsLgTVdLnqwz4M0zdxikJ8GbIZS2Vnbw5galKX5jHu4/241Y9kJ+9dOt8iZs
z0VyARln7oAvGYPDDt12JCLugHcVteBLxQ8rUr/Al6LH/clckFcVahtfW61I
6AWbdo8g+PK2udvnC9CPtJrV0t3gyzrd4LoHArA+CawdFDgvk3kPtjFB33su
uspdBV+6dJWzSUH+pfB3ZruDLw1fZJ2R3fp7//HybaKQdxd5aiT6+kko2Jss
0wq+nL89rXQQ9nOy8+hFQ/Bldv3Ba2JGBLTx6N+lS+DLDv41R5NBAlpZ+X1y
toKKXbD5FC/48k3J6d3vwZf8BfR1lndI6Nl43sdr4MtQK7c9QqXZSGVVNv8f
5EndYJToB69cZGQyxz8hSMF2UuPVh7cTUdGMZq46eNLRgOfrN/Ak32rCkjd4
0uFy7usMXRI6ZBegUQiejB42tjoD+2tPCOizAk9KmWcp6+wkIIvVTuNQyCOb
l7N209JkROBIbEgHPxZya3K3JIOvO/fdugB+tDcW0n9jTUAsQa2Zh8GPIzsi
9dTBjzRMn6WawY+Ka81GQgHwficrtrGCH6/MB1uQ4TywEnSDpCF/+XLenLiz
QURhTC5thyF/Y+KrV9m2/l4bZKNfDflXdn2w2xY+v1Udh0US+NFpc5ZvFdZz
vrxuWQ/8GNgt9bK7PhO5BqUVt4EfW6ojqdQiEtJxlBs9CH4kOrRsL+cioiUL
mqsZ4EfJ0zRXBsGPUYOOrwjgR/vLYr48cSRk0tvTKgL5Xb2jcWHHShYifXpm
8Rn6OG7PEisbPK+SpUovBfrYqI4cOkUlI6TWKq0OfvziO8zQDXn7VYtneyL4
cfjL5vmcMbg/kAIFGMCPI7mlb7iPExCPJVFYFPzIIeTH1aMJ8xCdvOYBfhSt
aBhLBj+aWL6/eRj8eCbluJ1/JBkd2nySqwNeLGPLHuwOJaB9Ecbre8CLlclH
3u2G+xrDUM/L0RoqNnMTyaQFL1bLr03eAi9e/agW9xS86L8xE9ALfc6nVCLK
DM+fGe3t8AH6XG933ktqARF1NWcdX4Q+b2k6xJGy9feyoSv+l8CLVomuZIct
X0xKcZmBF2e6xYMeQl6GkibGVqDfG3bdJa8zZyOnobz33lt/70IqnuP1JFTx
u3UuG/o9R1Qw7gAVnveWnDkneFFG0Tu3Erx45Hfwr+fQ962iRMQ+SEKtl58p
KIIPW/Qq0k8fJiIbeTaZG9D3uo0O6/pbf4/KiOxug77/G/lZP+o8CQUq20RN
TFPxzozo4mPgw9wQxTIT8GFiraTm5XfQVxmN0n3Q/3+SnJbuLZHQ3V3TrPHQ
/46uComfDQiInb/aJxR8+KORh2PjOOTpaLvXygT4UEGP9Sfk60Rc8N1T4EOr
VxcUnLf+nsT/sz8OvHDgjtIrzdRcFFgsnHwDfHh9n3ObjncmUjp35MwV8GFz
KY2v/Nbfj/gviDFCH1rdP1T/vZyEKAlLWd3gw0uLPHo3YF6tniVyZIEf/isW
vL05S0Qe/HERXOBDdbYfO74QScij+QxfLPjw5mAo78njRJR5aqLaH3wowq8h
9jwwC704PvZQC+YxL4nmjSvs59EYfIkG5vHAb0lyzkky2mv4mikX+oP/r9D2
n1vfj9h4YOYMPjw+dYZsewI8F5NTFQA+fPa7YoefOgkZxYxw3QEf6hcTQgde
kpHWJ6//AsGHD8oFNttuEJBJ3AmeFsjX++/rqh7O5KFMC/rqIzBPuftbb74G
z3zddsdUDubpck1tKXskEf0UqPczgXmqOrNUILH1+/0d3mxPIW8T3wo0KjmS
0Mu9ywPnYX54jwdmxcdlI8cdVLdOyMvanfsPvX5BQoMGowfpYH5mSr5+ftBM
QIGuiV8FYF6Y7Oaus2z1eZcEJxnu0516GybTefkoMk7S+RPk27GHWok/wYcC
XdF2/DAvO59/pz9xloRm/+jmFYMP/TU/C3Cez0dWXc+aJWE+5HdEJsuww/3Y
g053AfIr61Alw4QJAemdWXL+0kPFsUntFGW+fNQqpnhnAPIp0K6Ak+FoJrKi
nr3zCvLp1YViKdat31+fMX9iDT707Qt6K54O83g8f/sS5NF7fR2VUzqZyO1M
8l438GHcSpXJ2xwSYjWv6+mEfDnp01bw6kw2Co4o3BcKPtzTkif9cuv3z9JH
j5+HPG8MTHpy9XwOcpSs0h+FfLmHwpmKYT/NpPu+vYd8ENsmlCpkSEaSjy5N
sG3dr8ciXQfhfq4k9Xo2B3yoaCu2fvdcFpoOct7wAB9S6tMdHRrJSPhKiLIb
eFBYtL+KZEJCFowiTzfh/JQOuY8kniYj2yu1DwPAg8N6X2kDJkjo3PcZpAT9
HTqp/Le2JgfJKbFmiYNvBotKha5C3nIu8uxkBt+wNehaBW79vvO7XnwpzL8F
4VCu2wAJ7Tp7Q0AH/FZ4RDX1VQ4RfZiZq/wD/76xqejTnK35zWk7TAX/xOw4
8ilDhIyMcmw/HYF5fJCrop+7SEZrL3UT3GD+7K8XRwWO5CIPI1qDdejjN/fu
mq+PZSPh/Azbl9CHSxN7Qoa+ZSPViaKhszAfTEpuEZ6meUjGYTiCBfZnWJcj
QCyBhO7J16jKw3oru46QprbWW7OGfh3OT/vaffdWeD8zTaOqCFhfYbIbT5AL
rG/bdPETyMNXA6d2dEI/Dl/64GQA63skwDnL+FEWulfj7c8H63vXZN/woDoZ
0WhKKz+C9Rz+6Ty+skJG5+jcmT9DfkRqdZ4y4AOf0PWaW8L6PNfW6bBohf26
9jTdCN6/2G8o7aAPGenUnZCmhZ8PrNXd+eJ7LtJXuJb5H/w8/8KS/ZxqHnJ5
m0HzGp4/422gawk8v8ruMi9IcmxBb0ogdOUgLWwd2QbPI9qlJRtomIOqu3g0
jWC9Yn1uVz2F+6phv5NPNOznr2wR+h0uOehIf8czOlh/C8+o7sYZEpJyyMjZ
vfX75jT5Tzrg9VnT9l+T4HXXnW88U8Hr29STW0jg9SEmM+sL4PWfSsysaeD1
QkULz+eqBLQ36hvxI+z3jcoT24ieRPSIpUntLHhdzK7NbD94Xac276YTeL3n
ooih25bXQ4vjFMHrwgrzJ3PhfL/PFZ7QB6+LO6WV1YHXSzc4YtvA2wNu3Au7
wNuM1/2TBcDbTsO3vjqBt/vOBt6tg3kUDnWZS2XIQl50N/mSYb429Da+LMJ8
Ud9dUHgF3k6wlI/t3/J22NSsBnzeT4F9k6ZrJFSLvc5Nbvla9p9nJOR/5NAt
KRvwc4gdY3AC+NmutaXzIPi5xGHb4oEWAnr2/OhwNfgZSTw49mnLz2Gn1RrA
z7V0UZ0y4GcFU4NdIuDngIkr35XAz17EoC4EfhZhPv3zDsxPDMeVtQ/g5zGm
c1nc4Od96h1mnuDnUfWcei3oj9EvUnGh4Gcl1Uvh0lt+DoskdMB+cb88Xreq
nIMKgy9dGQc/6+glKEh5EtB8/s6nL8HPtpZsO37Ggp+3BX/xBz8rS6i1bwc/
02/UrCiDn/Mve7/hv0lCum0Lp1a3wc+fp3ytAD/fTX43TQvnK7Dhe2PXQA7i
Zdon5Q1+PmOiLPvOlojSJ81Xc8HPzCvrQ4e3vl+7jffBJuRH0xvmjTlLuD8+
/nqQBfxMn8y0bAt+fnVx8nEt+Pm4j2m9GPh5dV98QTLMX4kX50qxDey3zbzf
RfBzlODnod3gZ393HhUL8PNJoWf5wlt+3nZhkRX6QGt+JpH4Fl7PvOa4GniZ
sXGBbHWBgML33Cg+D15uvEpvkTQEXm6+/9978LLVU56W7eBlj/OzHm3gZXZ2
h/sl4MViL8kn1pAP4dNfzU1J4OWWntRJOC+Su01LgjRykdoXz5VX4OVrD1yG
2rYR0SHj9u+y0Cf952wWzm19P7NFfuYieLl97uwlefByJx1/5gPw8styn/+s
YT9lagpMdMHLqoL2eJIdvNxyRecKePmLOgc6KEdGt7J2Z9wCL//l3aYckkRA
Zkrk7hPQN/l/tOs+24CX1R0S+cDLD9wCefnAy4SiiNtU8PIfL5an2uDl8jIl
3WVj8Pz0dacUmOc/Ut3mPNBHZuw+ajbLRMTgM8siAH3UXmcrLr/lZXUadRL0
QW/Hz3zLWyT016q8KgS8LGRuPk0HfnuvrvZeGfpJ3LJuxbwbvKyukl0CXnao
8xVJfwTr/briAA94+SWPh6kUL3hZmORwG7ycKXjH7Tl4eX8Qw8G74OXDzsE7
b4CXT9icJ7FCn81U3Z6nTGchJZGh1mTIsxtiKeRAIhm97k/9HAQe7rOjCcwa
JaAAtg2HX+CR8sG9mylbHhZBu9nAw1VhV4pCwMM3ZRXNLMHDHElzzkVbHhbx
9+aDPG76nLBnLIaM9txTjVAAD7t7WBmU3ySgnPlBre3g4Ri/4Wifa+Bhkmd9
B3i4LKWG/bMSCZ0fIGJv8PDvv3+KvoOH7/D+UX8GXom86WH9j0xE52T2R4yB
Vwr0SqRebPmXxHrlHOQ9T2VvqIgbCVFvd3tqgX/PfHDTJsH6ucg5NX0Bv3gw
WSSe2wv+JZ3svAj+RceOM4k/JaFjNicV48AvV2ojtaJLCOiJxJvK9VEqZrt8
RmMZ/Msfvtn7GDxj0qTsvKeZhB4tndwlDP4NVR3kMJYhoklCwjUX8O+hedfC
zC3/2h1oKgHPjK8+ebp4joS4P9JMvgb/tl1Xc9cB/8ax1carQ1/Q6L17lsad
i6ztbPY1gG9mS+lMVP+REHvIFaMg8O+E2ZH86nPg33GTy1fAv2vxBo2vwL9T
oS7dX8G/5pvHxDn8SMj532avAvh3uVk232LLv+MJL/yhfxJs6GnR01zkoJpz
1Rn861pZGbXokonq66uSTMG/IxsN7ae3/Dv+m+cX9K1gn7OnaykJaSqkrtaC
f2nypa0jwb/aF7h+3AX/VhgXa5z7C/4dl72+CT6SiXcdTMkkoV/Sr+38wb88
75nxA0xEUe8uPXAA/5aYqV1NcctCZwocRuXAv7VZs07+W344P276HfI8qyxd
hqpKRpkuvNlGMF9ppROpv7e8e16L4gje/Tcimyu/5d3zIYze4At/1h8Mu/vJ
6OXpxVAn8G5djNKYCcwTFmhdfAR5WGU8xTOdnYd6E5iV+GFe7sVECU7fJCKf
KNSkAfMinatfcnHLt327aChb319p/y/xiAOc14LAClWYj9/bWQdFIW8c+/Rs
qiDPIh7tbDZ4SUIX6+2u/oD80re8Z2nVRkCO/3Y/Z4R5MNXsaBcUBN8mvqSL
A7/ceZl8z1c7HzH/e1TaBXm0uCb+6vTW90ESpY0YYB7izkjotOqR0LvI7ohU
8K3xV/Nt1+rAt4mXKnbD/lsHrehGcpKRwsfbCSOQN9PLbkf6TcG3nCmWveBb
qtGT9hsX8xG1ouZSE+SJeuymyBeZTDS2p6/nCeTJ6/3DslJbvuVceaQL+yl5
awTdSgE/dUpOjEF+CB/LK1EzBt9yHt1hCvtnoWlZT84loaUp1eIqyIO3HYsN
SzgbHVsYkFWFvL14oMd5BvrLbOHqSCucZ6NbD15TTMFXVxneLcN54zLfWagJ
8yf85HzXPfCV0pW/FbucstAblu/D5uCZhZ7IUnXIE2F5ISlT6J8o8nXnBWMS
0nqYk/EN5l/ChHhU+xr4Vd78riP41WQPz3j+F/BrWbS8yJZfwxt7NHpyEN/z
ezc4wQu3TuW8IZRlooc26SZ/YH5vsW44x2/5tWwuLB181u0dvo30moTWhjZX
joBnJjnDWYWpRGRtfHVtHHxhxOHHUbE1fy6Ooikwf5Y919s8ecGvLtveCMJ8
bfwpEh6jyUFfTLs9TGGeHLYLL2/LyUWmE5Gy09CHnN7N9QeWsxGz4r4TFdBP
Mybbs+WGwK+FZ3qOwf6/T/DL+8kNfj144/ofyPOV7WGd7yFvjzGbLU9Bvnas
arLwj4NXmf8UeMB6Ijdr+mFPWM+7d7Nzt9ZzTKfQDvLh+aKltgqs5+nTsTV7
n2Uhv5N/TOhhPZW4Gz51HyCjmb5i7oewflNqyqMGsjnoXAfbWjecd1W1QJZP
SuDVDi+9k7Aegiecdmo3glf/s447Ae+/rvlN5+NlMsoYqwmkkYF5ETZAhaP5
SF63lXkO5vHdwbmd0eBV/X+7E8Xg9SiKi4+ZPfKQyofzDgkwDyOCFiHUwhwk
+3rMvgSeh/k3b/OcVA7KskgTOQHro3h9KmnPMPgnjfGyL+zfcfK+Gokg8Gqa
a808PM/uYZW/nJMkpKIYV8EO6/XkYJ3Ef7C+Sg/SRp7C8zGM+q9oPiUjfWOD
cTK8f1btR7Zi31zkvtGg/A72Y6G9hiNaKR/RXW5/0QrPq+S/HxfP5yG9W8bX
V+B87OiyN1UdzkPHLiXrOMHzXpqfHe5yy0Gqgiwz9vB+ontWP5/MykMr3GL9
zODlwR1VC8fAy3srAumGwMtmK/sv3wYvO+lLbU8FL1sN5peqgpcdQ04H3QMv
n9VS9L4GXv6oRm3phfly+7zv3aA3EQWvSN3TAC93rnSd4QYvL06sM0iDj713
mudngI/Fri9waIOP2yUKnq2Cj4fNnxGegI9lHp46v3kjA7lF7Hy3C3ysX08q
PQk+PvmgfrkU5t1pd9T4AE0WOjKp7HQM5pXgxZfACf3wbVrA3Ri8Kxp3xuk6
ePcdM52KEHj3/XOx5rNw/74dtLRJBe/28G6XUNAG706mV5eBd+3MwlLZwbuh
ksJ4D3j3HU14gBZ496cg7Yle8O31Xs+v07QEZDg5wN0E+1ewr/6/ghM5aN+z
/I/vwK8GSrp7SjwI6LFV31gj+PVqlNTHD3HgV8duaw/wa8SVEba/R2H/X0dJ
HgK/ltppfFUGv5b0RkUvgF8nms/l3gO/bnzcY7YM8+M/Eh/+7GkOumdgk3YZ
/HpyahnR2RCRVTq/UAb4VXPXC26DLb86urD/gfwoUln0uQV+LQmrMaUBvxqm
P068DH49ciXodwn49WBj/eVZfgLycurluQvzOBOeorMG9+9q7scZF8CvgTYq
nUzgV5eHDbfPgV/LNK+Z8m351THHnRbyviJsvP/OOxLSIC0XHQG/0i3eU/pn
SEBar1DfKcj/dQ1NN6UR8Ov+L+Xd4NftIomqK/Ik1CD6lPAE/Mo2kOoyCH4V
rXWRNIJ8wB6lNoR08KuY0K6PMI96YfuaqTK5qFymzKUN/MqxvaI8gYaIfrLc
3SkGfYFXc+Mit/wqFm5/AfyKrjb/XtMhob1vv3fFgF+9751vuAx+TXqG0k6A
Xyv7Tl4u3/KrWEOLLfjVrzyg3OQIGa18NWIJAr++/snSP5IIrzfr8FMR+oQk
+TzRbeu/78mtYOEAvyYdpS1hAL9+CzKqJINfn26w8JuDX1292Tp2QL9oj5+j
Gv8hopcHSw7tgn75cX5gTWPLq7nn6h5Cvk1M0kjaRJKQ3t2Qk7LQL5/2hSb1
D4BPc6N5c8Gnm1rpt3dRSSjhlq8uK/j0iCEP0ZMffBqxMBwEPmXRSVquB596
6b03uQk+jRR8NfkY8rIyP4mDFvqott6ef/8Y+DRSWvUu5APT0+abex+SkRb/
x/NXwaflrx89yoL9UG1rujMNfti4+qz/w5ZPI+/F0oJPdS/p+puCT7s36G+f
B5/WNTI8r9ryaeTL7zshf12mogpL48ko9u2NX5Lg08SmookD4Amn3HTnFfAF
t9Fjbomt/z5ie/PRJvDpnuWP033g08w8Dw9n8KkePj1HGwx9Y6tU9gS88fgF
T29+NhHt+jyZ/w68QXfjCOPolk+3X5zUhny/x/elLuIKCY36FRwaAm9IHHmw
t0QIPLo9RfMCeHSN8fJwJ3h0hYvJPAI8qnI/8H1zKbz/P/Ln3+DRpktK9nwR
sJ4dGOeDP7R3WQZ6V0L/oqgwbvBo7JeazREpIso3NEyx0dn6/tH43xdbHn3m
eywX/KH2yECMCh61f9TH+QI8mlDuaH0O+obW+N+GEvTJYPI7ZzL0h/Wz0pQy
8Mj86NxS2goJ9SrI+l8FjzJWK8+JbH2/X6tgzBY8Gh2s+LYKPHrUUGVlGDwa
q5lpLAwe0NcaP+EB+b3u6nWxvDcXDUZMjl0Ef4rtfqD/8VImCrnqV30W+pur
cjPfdMufWloJU+CVXahfprqYhOrSja1vgjddbGxM76+BN7VCfiyCV0QL1x1+
gTcrxxNueYA3e5zDRXU0iWhNsGC7OXjz8b6eHNKlLHS0Ur5vDPKy+tqfs2PK
4AuhobbTMB++U8bf1ba+D155X8QCfMlwrc5q55YvK3tCL0P/O1067XJxkIyi
aapzrWD/F6WOfhINJaBKBzpHEvRP042XSpwheUhr8HQeB+x3h47MwMcbRKSE
aCaOwX6zhx1hDtjy5VGHACL4MuGw2r9yexIS/1EvegT2W/itR4/KI/DlUeJg
EeTNSL2L0u4O2L+r+x98gXypIp/JffaMgAax8dF1yI+M11yDdlzgy8W9N8Kh
v7y/ttd17s1H3SXu/S2QF5f8iVf9t3y5GPRqFfKi0CZ5KBJ8abVYJcMM+/db
r/bOf7vJaClep7Ifzv/5YdoKPTPwpM/UwHPoS6fxaAb9sHwkvLo8UAnn/a5p
pZ/dwUxEfpX0qxTOO0N3jut/W570OX3gBOwPMz/XNrWHMC+7XA3fwXmnD57k
GbUET/pERurB/uzb2faKJQ/y43CIwCM4v8c2TnZQjoMnjQ4UyEP+mSc/s2qX
Bk8atZrUwvnj76HhDLUA/wyYGMzD+fDVTqq6cZ2E7tMkq0Ru+WfWvX/f1Sx0
0lBGzwD6Ime0pEwpDjx53zNHD/rgoqB5Szp4snZkkmUM5jcm9p+MWxTs9/1C
WgvwpJROOEfbV/Ak11ARD/SZ1sWZEfqP0DesvbMM0K/j+wzyXlAzkRm3VeB3
mMejd24Kp295kgutxm/56SlP29lXkE/71V0loe+j628ZWT8moqPdlRmxMF/q
0kO1W/99pVH3+TO7YX5IOX6/jzLmIIsstlE9mJeNGaaJydu5qFvqBWEE+ii5
VozIv52EQoJtSihb+/upcjimH/zISjh1CPY3dX3gkNnfXCQT++rHN8jX9aUp
cm0s+NGNcnUE8s+hu+d53Xvwo5u2iCOsXwRnLt9+8LPk7AfedFi/snC1U2XQ
hyqqxMrDsH4M+0PMe3uz0Gyc0osVyAu12wNv5PeSkTXvj8gYWK8nO+QPRZ8C
P8pY+7bCeVX5rrw+oA1+lGl5fhw+f2eI3hndBvAjoYRBEd6fPDgZs2BHRhyS
/2b+wOczDj8RX1+ej/IT6AI+g8dC18zyp53Aj1aOLHzwegOeoZZ0UeBHVfLw
Ldj/A9wF27PSchCZR/RNDuxn6O0rLn17chAX+WPCUVgfG497NvZz4Ic10zFX
2K+6N02NyZHgx7U6uUl4nn9SsrmG4/D5Ukdktm+dh4fJDWF3wY/LMyYV8Hyf
U39kLVWCH2uzLFLh/Xmv7XMohfX/rL2R3wXP++Dtm8DwnflI1fl9ChV+npWk
MUf/MR/5lzAp1cLzv7h7zZVZKB/pTef9+AHnI4F77nXaCnjy5USLJTx/8pfO
6XFr8GSohb0Z/HzB8+9pEYF5iPet74Fk8OJC/UlzefDiPtHvR+6CF13yau+e
Bi/W743r74L5abXfcTHQh4i6Rj9IYPDgZNnRYIEHJCTEYvmXDfx308kwRQn8
R/Of1eMimE9/NqrAv/VMVM49JqoA88ZN/9HF9y8JETyd750D/xFO8Vx0Af8l
FOcs8IL/fqnnsX9uIqDLJu85KeC/u49bLC+C/5K4Na9RwX9ntVofbCISkkhk
Mt0F/mu+33zqDPgvJ/esbRf4L2zkQ8Ib8N8b7htvnsD6L3CyeMdr5qDBy6eS
+sB/3Jk2il7gv0Tl8l9PwH+/2UQiDe4RUPgDH34X8N9Re3fln+C/f8pXTxwE
/42PG9KcBP/ZmQ48mgX/iY8WMwaD/2pDanf+hv13MPyd5Fueg3R+b56zB//Z
Bj1k4rcmIjXfNbmH4D8NG+kua/Af+QFn1w8472oHaQrPg/++LLzgKgLvFTeO
6vSC96RYrvZEwPz0cXNd1ITP8/nBGal1OP95jSNr6+C5Na00e1nwXPQSb1Y7
eE6cIjShobb1fbZeGpZRArKPyauoBL81Vvp4zIHfhi1ZxvThfDpV3xrSSclG
dcbtPQOw/4Gp4T/3iOYit97zIo3gt03XV6/+bsJ6fXEREYI8ZvncHJYJfks3
lhLShzy2st///TP4TetMdJ0KeM0n/VxEGnjto7FjiCV4zfmGSZuXAhlVRv9p
8QWv0d8wnogDr6290KKVh7y23Vz6WgZei35r2c4CXgtwL6BbPwb9P6naRQSv
Kc88GLAHr4nxlIQyQH6vPRa73v+biIreLl9PhPww6Dd0/AxeYJAV2pCAvFaM
onUX/JCJvr098i4TfLZ02OGAAfjsqq2JJSP4jHGBIeG+ABHtqWpyDwSPqQ6+
X+gEj3mcVni5Bv36KYp94ttQFqqu6l+OgPN21EwyNSiJjDYaA5ndwGOpBXmn
hcFjgpnZaRPQt6wX98+vgMceVinprkPfmhwK3tQEj72r8qKwQH4pNLKkZCaS
0cmAfUX7wV/zgwUTn0II6NTNoMA/0L+zNzqnQ8FftxWdF56Av/ZJnxntAH9p
s5vecgB/7RZwy2AHf1Vvf+9cCX0cPZ6rdzOLiH43PK/tgz6OnXJl/gb+ylek
J2HIRw42vXxF8Fcqv87Xt9DHuefZsev+bPRV8QSNPvgrf53puyDcl4d+TLuG
gL/Wp7WOXHxMQMV2Mf/mwV8Zdmm3pWE9T7uP/82GPH1/2G+jsoSEPo5L/scB
/rrGG2G3H/x1R1ah0AL8lbN5XW0I/JX5R/BPJviL9o+oaBT4q2E5rVoe8pdx
o2vo9qlcNPrH1IgKfX1DXjzu9CoJlbGzR7uBt8K4csZMwVv3fA2IluCtDRGX
n4XgraE6Uc5B8Na1+VVfSfBWj+/dNWfIO42hcy3FH3NRksStVDPwlhphnbfP
IROV+M7qfYb+vlahz/4J9tsHLXEHgq/ah5x91DaJaN5XnPoD+rtugnNRKIuE
HjZcTXEBX11hML2Wq0VE1S2nGi7A/nP40r4/ZJ+FKr98iBva+v5TWeT4LUUy
0lipf6sN+12/9+1pB/DV/S+qo0aw3wcNai5sYBIa+HJdzR76UTxDqpVtiIzc
hh5Wm8L+djooRN4FX3kw5vITIJ/JL5vqt7nnoY2g71Y7YD/PGOcXvQ0mor28
44sKsJ83nBq6b4Ovci/tUEoFXzm/a+64Br4aS7L6KAP7SeuXv2nxOBtNXtJ+
kLf1/cXz15p9/2/rzKOpfN81bkoSEiUkShqoJJRv0kNkCJk1CGkgCUUlCYWU
ilSUoWz2YGdKUZEypYgMSUIiFalMSUiSzvX81lnne9ZZ58+Wtfd+9/3c9+dz
vbZ9h3zFay5w8z3mveXGAvGAcvCwa6T3J+Z73t3K2RLiOI9Fj9cGgecm12Z3
G/Nzya1zFp8KMc+ix0tdo5GvbixSOj+Kef7p8NnFDfmqbZHzV36c1+xPClFS
s1ik2UO5oo5+PnDgc9055Kua1OjLpfBHREnJ5+cxXNLKSIi8jfk0izhZxVrC
IFmpQx/W4jyaCy79TAT/C+8LC7/EfFou6qnW2skgfanLNxjjPDb3rVLTRZ5S
eiXfmIo8FTBt45Rpa1PIg/DaPcvBo3vSHa/aFoHP4R6SdzFPzNX2gYM7WCRt
yt/pX9HvvrM/1s9FnrLgaAwFIw886Bswz/dLJjyjNfym4O+CqcNR1uEsklI+
x8UYfBWy5as5iDzlGxJe1or+tB6zFlsdwyLvy61L7MDT+NUjwTafcT9gGLZv
JvyQ9j1gt9B78PvOIQ4vfMX/XemYSwaD3DHszj8PX6fM2s5ZWsMkKs8+LFgI
/21M2NMUmptE7v1x3HYW/SP4MU8hGvm08c+vqWLoj2N3vBO8hJEHDLKvGqMf
xjbLGcw9ziFn6tzsWsB3R5Wm6pXCyEPzedzYOL+eZ3O3lL9IIYXrNwoo4/z6
bs05GfmFQ7KKjmZ9Av8GTPJfqCA/PUiyUm0Br/Ty9ixZ/IpFWpL6251RL0PJ
ecdqjqBe80KbYlEv/uOlAg7IT8Md6z1VtOj3O+2N9ZuTCXdVc9CIMs0H9urB
M1lEuiha5wzqY36oX3vAjk3qnQTVizB/kvy304M2c8gdp/2ha/D+/2hvl7R8
BD7VbnmyEq9/5fQVQ5GdLPK8OoE5iPfnmOBsoJ/EJS6aHLV25I+MeIXwpMBU
Uhcl8mwWnm8LpzLO5WoqeehlmnAK5/2nd13ei8ts4viw7ByDfl4821HLU4xN
6gxPmKxEfThn1z1SG2WShhU8Sa44H1/15hW7z7FJ7oo9/e9wPcMaKpfHOnB/
VRX+dRLPp97SECISxCF5ylcks3F9Q9IGBz/fYZG6r8ZzYvD6K2pf97lFc8iN
rwznClzv1xKf2v0TqaQw4aTdTdRbXN3RJrmES9Rt0vtzcf25FtqMJi3cn8hY
ZPWg/1mR3DhdYS558OtC8Ba83psplwZlt7FJYQ6fgi0er3vNx7zcPZVMHH50
8Qzy0rXC37JrkJeWnwlyrEK/CHgnTTM/mkTkJS3urkdecjgR1bkc89J+frlI
GvrN5+77isgJBvFV1QxRq6C/71YWT8f79/p7ePMc5KHq8k7LLOShlifPElOR
h3abx4b7Ig/ZqMa9TEceyj7TIzOGPDSx78dgJfLPwYtP9z1D/rmi+lItH/W7
PV5hEGHIJvcL5gzUIf8M+B60tEb+4UtlnchH/nn7juf8eeQfc89Pzxcj7yhE
7p+5GXnnan3EgUGcT3L4xalvs9hkkWF91k7km/4dsutvOiWR3dmuyv2YN/+6
rvdayDPkvmfqTeSZiHYVnTLkmamOG1RCcJ6qj4cFv/kzCSs7OWIc86f066ft
2mYmeXNtf7EK8ozNTqG32cgz12aLe+qBb6ONHNUjyDN5JTK6ucgzTh0BP38i
zxTcal1njnnZuLI1mnk1hQR4z13RgPNgnNkU80GOQzYpKgY9Qp7Zqfr3TCby
jL6tQ4Yc+PfWeGt4JvKMg3dwkBn41/b52dBr5BmGmJrZP8gzs1YGSEYhzyR6
P2jbhjzz7qQiZxJ55nJb+bzDyDN3M37/dkWe8dRcfUEVfNxTPCOrEXmG1z24
jR88bCnyn5RAXvEYMH99EfNao2fhwoFfO4O/bV8EHnKihNIM3zJIxsAZzRvI
K7+Cp17ygr8+1hp1CSCvFNZnMpbMw+O/avw6hrzSP5OXtwV5xewGv9I4fLVw
2wzPl2+SiV/dkq0htL8VeiRyL7HIW3Ez1/3IK60qYY/et8Gv0y8rfoC/gti7
z05H/22pO39zHP46Iy1tpo28cq3uqaAQeMHrO9keh7yzsHpgygLklTNaM74X
Iq886DvA/wM+8yWjI7eQV2K23HyaA38FScn/7mUkETeL7X/Xgy9/54i9EEP+
OOolrf8KvtJ/IFVgtiiFcC1inM2QPy5HNFw3Rv4wd3rz4wTyxyGtg5xM5A+e
a39P9yF/HBTqm6eJ+iwpybFlgEcjfI4H3mYySb7Gr1hR5A+/gFvmFcpJxFnM
1/YG8saluWc1/JA3Ekw9JFaCV7L8A1kzrTg4v6zcdPjLZ/NXUzbyhmLgpIQH
8sagpOeAL/KGxUX2723IGyUTfbwc5I2LF9sd3MCDooLLGoUdHHLwuP2wPfJF
Z17tsdrdDOJ9US+zAz7rKZmx/gjOx/pShc8x5IsyefuaGB7kj4sBwv3w2anS
2ro1yBeCKrvn70O+iO9fbG9rmERiJaWlrXFeUXvEfI+5JJMjkyv6WjDP355N
D96hySK5+7KcDXA+v76uFvJDvrCbvKRjg/MJaG9+PYJ8ETv5/PpO+OPUssGn
7m0s0mZyVteefl/exoqnCPnCfPuRYwng16/AA4GX9qaSttpHD4RwPowgqSO2
yBd7g3fGXkWeCGE93maHPFHctXyNCs5n/Vt77cHcFMIJTvjBQp44epu59iby
hG/B8Mp3mLdTvcekXyNPHFpTYTyC+apaXX+/X5RJAsns+ADwbTtP7ZKCsVQi
/5N4P8Q8rXxu9yURecKR+H/5gXmKSHBQ34o8cYPkGvPiPHq2r88MR54waZGx
rMG82Na2H05Dnogs7vxWBJ4uD9hQsxK+eDC4/0sW5ifeNEH76WIG8Sw20tNC
/WXSff1MwcfEqW9d6zA/Eq/9RT+7IG8Un0reiPq3b+84th15gm/Rt5Vs5Im9
28N6H2qlEH/GgjJl8GD90QeDg4psEs8o9rmDfn91tPdLuiOLuMvbenSjfx9X
TVvYi/ur5aP8FoHwY88BhkbE8WTyzjhppzF4d//LlUG3EBZxeXegeCP4ttPi
jrsD8oTVK3v5FvRf0leNeTfjkTfepc6zAc+MFudHJiJPXHJuLhcDbw2EcmdJ
dLLJYr4NPH/RX5JWoXGq6QxyyFlbKgL+MhW8qxZXzSRCMrlB8+n3FXwW+167
m0SOSOXknUa/LFw/HpYgwSKxUhZu09EP7gqHIkVF2OR7bPD3jTj/yzIb7Y19
OWTPfK07r8FX1vXazXo4r9m+9Y9TcH6eim/+dtWmkMAt8bsX4/ymXbTPN/yA
+72mKuFO8GnX/S67rgtM4p/HiXwNnlQuVr30sZ5F4vM26DiiXlMbDJJZR1Gv
NU2aV1CvMzNCGKrIEy9XiVUuQb3+Of615VRrMgkLu/lmCPMrsDPlsYoYi4yI
ucSHoT5LD8c75u5gk0t+DpceYt4eekrkOVlzyCG/ovcaeP/XHQLfNTxkEr/u
jAUr6OdlnEG7JCcWuTnP408/3p/uy9yNnCtcsv7MkchW+NhsV8yqxLOp5CLX
ZZEEnu/FfaPpS2+kkoAzjJ9BOO8eEZdzI1FsslYk5tN1nN+7qtDsJ9PYJDPO
jLMc9fkc930ra5xJrhjZ/d6N8xFUPsd8cYFNfI3yzdpwPbcioj6bI08EfHxj
PIHn41XwutvjxyF++t0+mbg+a4G9sZ+yWeQi342jl/D6JjOecepjOOT4uoMP
n+B644RsOT+HUklg7qtcDs1vYVobsnO5ZAbn+KY7uH6j0zkibhu5JFKdJfwF
/T+X062+SoJL/CXet9ni9ZTvcW6WIX8FPt8SZIXHG/srnxd1SSVdhSPzK9Ef
DWFesveRH0xWCB3XQX44cnC+7GrMR83dWQ+5yA/Nvncmm34zSIO2s9pK9I/p
7mIrA+QHx6WOv2cjPzz9Nul3GflhSPt3x328/yvThdrDjdiE8/Hv9VrkAVMz
rY4NyAMP38rq5SEPmJtpyKQhD5z+dUjmG/39xr7VI+IZyLPeDxyc4f8JV+0D
XvB/2cPoll7MRwuPkcAy+H/0tHLLSdRb4JXnl4oAzOejKq1fmA/LGl3xrfB9
ZYm1pzJ8/23q1wAWfG9hHqGmC/44CJ9zTofv05rKvpqif6v9bp6pi00hrQGG
bfWoz9/jvHZmshzyj5nwsofw+xTFAO9w+P1hQNpyU/DnxNTY4y/g8/GAjxe2
wt9Vx34+NwRvT03JrvaBvyP6PUVs4G+lA7csV4BHc0L0ZnfD34OX3M7zUX/f
2tqbO5REno/7nY7CfBziSToaDj813GyeqgT+eBlvTHRoYxDh33c6E+HvJwMr
awLhh0SBkjR++DvfwdvMB/6O5R8o9oO/7cRLz7+Dv7XLe1+OwQeCtmPbyluS
SVOrteAp9JNjy+mIumjk57Uaku7w9zlhw4lS+Hu0deC+IOZvzz2P0Vm4X5EY
brqrAD+HnZ17IhN+PrSw69l3+EBicUD4U/g59GS07x3wf0qLpTcLfq5wCMvW
wbyGti449NuDSfZe4xl8Cf5nKwdbrV+cQgR3FImaws962inBjvCzQKz1gwD4
2enLk4EG+FnqSwJvEuY5Xk2852AG3cf1QV8EPrZLr6qPg4+LZW/wXYePlbob
9Lzg44jDNqUrMO/uco83O9lxyKRs06408P9k5fWKF/Bxa65fszt8PO/StoBT
8PFA4lRhV8xPSe/pTOOPHOLMWs+yAx/HfgZ2mcG/dYkeDu/A/+pTN/h1UV+9
gmwFP/i3YT45oMvLIGLXOQW94P/A1DY9ffi32KW1yxX+zUyMH86Af0N1eCos
F9PvF1rr3tiZTBqFt11vxjz0njD5W4E8Nyw8YuoMnp7YH/hHqh18PfmX3w71
zKmvrWyCX7VDHRfFY553MJ94OmH+no+y3aaivlbXTQ5JwK9Pz58ziIVfee99
St8Iv2aJzu5aivqeWvJCUedeChG48ITFhF+Nc6z5S+FXlUHP4Tb098d7ljl8
FTeI887skR/o764vSXvbRHD/baprfBzzHpjfsSb4B+5HTFlxQ+jn6VVlDEv4
dMK0bfQv+rlCxmVXEXzKM+OKZDX6N6dkXk0+fNr3XDS5kPo0xHxKH5NLuDI2
1zLp34tf5VtkD5/WPD/0bTXq6ca6tkURvIhQeyJZi3526OgwKdnFICLVGRYG
qKdjyDamG3w6uLm5nQmfxsSIWc2CT1syzHyWYv5eWfdYBi1gk7GML/Nvoz+N
JAL7BcHrSpUgmU/ot0e7FrU3wKeyCn2/A+CH13PdX30PSCY1h8JFjcAPoZdV
iYdPskhpz1UvA/CiQr650BI+JRPra5rQP1KqE7F7ruP99b6osQIvujLt04vh
00EP3qOi4NUqB37xhC7wZ5ly9h/0y+30f0aupjFIvcfe8jPgN+umQZEFfDq6
IXGZAni5KVQ5bv69JNKo+NY1DOdv4y++QQE+HVb0ny2M80/LkD70Dj59X+SW
bIDzNr81NfKDN4fYbFJwbgS/Wvd3Ln8uhvm4/uBgMs4rR+bZxZGaFNK2u2zm
IpzXw2Sj0fR2Dqn5+KPgA3ihqcur8Rg+bSmrW9eI+Z6/PMar7AWLjJUd6HVA
vfZ9C5UQ80O99Hm6LqFeqaMi49LwaanD6NHFqFcvs/H8grZk4psVffY75m+Z
4s+kUBEW6dIyMQpFfV7frOSR2ckmg6GndQswP9trVt+1suWQ+tDP0ep4/7v9
9/+JhU+bhhrrl+H1+Y60V3Qj71wzssnsw/t7WWbU4hjFJctvO2q/gY+KfpWr
L4hKJQO55xvF8XxPLc5KXElOJa1XKrmBOO93S++XT49kk6WrT8Qk4vxGwo99
8ZnKJomlGjbLUB+XS7aV7RNMMmQTnLEL5+NvfGtgKvzbYPNhvBXXs6mG3+Pn
O8z/AP/PcTyfv9ir7o7DHNJkIb4gA9d3cL+mRPctFhkQq1CKxut7B+hMm3aN
Q9z3bnUvw/VmKjoJGQ4g3xb/2cVGvUuZkj8Ms7hksmbXz2xc/5dpMXMfmHJJ
H6kp6Eb/8y3pr10lxSUt8sKRNni9GWkxRjK2bNLWdGq5JR7/aDDO/tGOVBI1
wrTmUUN/uF/05XzgEr2WVWHP6OcZzidzvsKvBpVCM6lPHbdoaQTCp3wbfpxV
Rf8sCTaKFIdPX+lY2VCftk9bJnIaPlXccEmL+nSkZcGogDGbSP3cwqL7mI9Z
vuvMgz9/i7oeGkD92p2/5melscmF0Kwcuj/YQlTvrBX8OVSqpEb9+WLtyzgl
+NMrYd5K6k8DnveOd08wifpj5+gxzMPattghV/gzY9/TKwQ8sT5cOPEEvlTo
zdOjvkxTklZ2hS+FQ1o1qC/HCl968cpwyO2dkyEF8OUF8m2vL3ypEnKsk+5b
3XBfyHM//Lgv1eED3Veadedbbyj8N84j0kb3k0avmSXoB/9tKqx2XghePA8t
MPkE/2nzGq2j/nsYNn7gNPg8qvBbmPovaJPasmvwX86bjqXUd7OaPI7rwneC
H7qdTtJ8xxwWd4LvNlgs8dgH36kueCifA98t/hAmSn3XuXM82Ra+qyHzVajf
Rvf5v3gBv/2KCaui+xXPSXdmB8Fvoy7SgtRvr6yannyH395mDBtRv4nNmWeb
B7+t2WXjugl+S09jDbnBby63o2Sp37L6dOO74beoiahtNzBfVvXLTzxMZ5J5
Pq8TpsNvDdr+107Bb3/PGs2hPtOYEZV9eguHqC44UED3CcqqZk/vhs/kmTd2
7UW/isTMPH6xk0Ne5636Rff52R0lKtPhr8mUvzl0f5/TQbavEuqTX8vyOwp/
pcaEbr0Kf8kGShlTX1kN/nrRBF+NWw3LUV+589eVq8JXAjP7h+j+Ob9a775A
+Epp5nmWE3jFai0x9oavBhPXudmiPttN0m93wle516wC4zAvZoECxobOqUR/
enwR3XfW3bR++PuJJDJ8Sf56DHwlvaSjUQe+UlKcto76arPvUo+P8JXm5W3j
dB9ZWZ7PgRr46sz09C3UVy//fAmVga+aDrPMqa/EXc/+sIWvRKwbGf6Yp4VH
Qho6B1PJcmsfC+on84L1Rm/gJ7mXrFG6jyv36Ot4cw6XyKsa9dP9W8HGadx5
8NFE/RRj6qOrt+osxOAjHsN8jxr0l9CMrVHLUD+dl7pcffho1DJYyws+ct9d
rUl9pJmSKP94TQoRynlfSfdj7W4ul3g6H3zKCfLPRv9MU90Xtho+GlOT8Omi
93eyvsUV8NEVjfc21EeJjLQB7olksjH8+F5DzOMG6wSvgGAWyRNXX0j9o6Am
kV12g0XUvu9Vov7RejqnqP4z/XuGazV0/9JhUf45Tp/YJFJ33hTqn3rVA/mb
4B8en19y1D8luZJXVeAf7y0XQ+XBE650SIkr/COgfKiI7gNaNK6m/3EmzlNZ
1GsazjPAvlDOWJRNtjY4jNJ9NR+7g3JWenFIlZPkfbpP5pT43tPRM5jkxO2s
CrrPZdnj6wLB8I/I/q37lejn6wXO4+1vOWTiS8RM6p/qqN602/CPUNWemFeY
t29HPperwz9Lq3gMqH8OF/F3GcA/YyZX11H/qFueWCYG/2gd7Kmh+zciGneR
Ifjna2nYO7pPI+fU18Hs6Syyw1wnKQT1eZ5vfrzWhU3mn5OJe4D+35x1b9AU
eY/nXOBnuo/CIlp5rh/8I/hr/1LqH7eLcx5egH9EdxjxUf8UfBrlf3SeS9hP
rC63gN+VuskfuqNTiXyBwnLqH8e1CVkmzFQinOA0SfcT7FozN38O7s+SNx3s
ofsF0qO3cDYJssnM10vS6f4AhmSD8sQfJlHcLslH/eNZdT2AA//wbfezof75
lf2rtAb+ER6Jt6D+MV02S8AR98uC9qn+6bg+STuPTDf4R372jsCLeH2GoYTC
4jgO6fA3K6HfN3QdWb2lvzeViJRfKaDf91voTupnp3PJ0Q/2FtQ/dX/Mijo3
c4mc4a6Z1D85wjUSmtJcIrQouZN+f05X+5GhoQ2biLTPPm2BxzPCf6fEbk8l
h4Wds/9WZ/2PfwY+WWpS/wSU3eUR8EsiwQ2hhfT/m+guTnllAf/IPi3SoP4Z
N25S+zmC+4mnDp30/ydYr/nm6kb4Zo10x1zql9DPcjXjN9mkInKxE/WLlOgp
KX34JY0d10r3xZfLyT5QDmQS043K66hPfk9U7/WCT7Yd/WhAfZL82m3zS/jE
YUStj+4371CecVoNPlEUSO2g+8sD3KsrHaU5RF9AP5r6Y8R1QvcE/MG63RhJ
/dFLygc84A+J001nqT9maT+x9IY/TlcYClN/9C6wvKXSziCcFTKFdD/0HB6P
pDT4oq9zeyPdr3x5zvzYieZkIr83Zxr1hXLuWXcl+ILstSig+4L5tOUmfOEH
TcafI9QH53/0WB2AD6Tuew/R/bg7b634OAT+VwoKCVC++4optuwE37edsDSk
fN96XGv8CPiufaW0jO6LTbbc8iR+K4cYcVa7Ur7zbRJU+w6+r1HaJ0r5PuWD
v5J6F4eIP87kUL7fvTVW34O8LKX03InyXf3R7BVzwPehZjlFyvej6m/tSsH3
ynOaY3R/qlrQ1dI+8H3NDu9Kun80toJn5phzMpGLzWdQvs9SnzewFXzXibWz
oHyP3LxvLAF8H0h5voTynPvTM3KtUyoJni2xn/JcbvXt/A/gucvysG66b3JB
I5e1+n4KMZ+1kEv5/WDYNqUB/NaUeTb1P/zW5+1XAr8lguXG6D7IiddNlTPA
b6VGhinl9/DPA3vjvqWSjY3rxum+Rv+Q+es+g9+adl4syu/PToIhQalc4vBP
aTzl9ybWtR5e8HuWXf0Q3Ze4dtcPcwHwW8dylRTl91wR44Jd4DfH0/A93W84
6/63yDjwer5m5hHK65T47vh94LWepokS5fWiZONX5uC1+L02Ocprt4Yld0rB
67r1DpN039+RvCXRvIHJ5FT0jxmU14meTqN/glhkSDarju7bW+og6uubxCIm
R1bW0317LsP2lW/A63+GXfwpr5fwDemUgNeVm+Lu0P14b0e2JS0Er6WHn1TS
fXjfQ3eES4PXGS5CqpTXF+YpbY8Ar+WytPdTXk8wd1U/Aa91spqlKa9jC5+H
nACvr7Q1Mimvz++LEdzqySEC7lEulNejVStNh8DrvILFvpTXj4ulzfTBa6W+
BbMprzsWzDw61sohszxti+g+qlvCcZ+44PV8c1Vdyuvi3W4ivXUsomdePUD3
N/XrlTQ4g9e+x/f6Uz6rmbZkO7QnE6XqP+con9VfntUwBJ+vbckzoXz+mOH4
9+YuNvlHpEuf8nl2ib1uJfgsLWIcQ/ncaJN4yw18lg/UbKT7fA7K6eg7g88m
rqW36P4dr6CByx7nuORd7XMdyue8OhFD08upZI1OTzPdV+MaxNypxEolinJL
Myifv5V7dyWdZ5NW269XKZ/ndH1Jvz+FTcw7mHaUz+ZWk6cvTIJXbe236D6X
80d+vxkCn2Xb9Cbpvpb0M+998sFnRf8943R/xKZpaw4n+nCIfMtBJcrnyfAI
sSXg85r4RUspn3MLJTvvgc+yYeUelM9J+2s7H/ekEiUTJ1fKZ+N9icXBN7nk
Tu+LX/T71jc2Ga8OsOQSzcplRfT71GuyYu6uluGS+Wn7oymfXypJ+by1ZhMl
l/crKZ9ruzw6Qrb99/1Bzb98/g+vV/57v/D//vx/8fz/3l/838cLV4mr8C7/
998fmI2PM+v+ffzj+YfELmv8+/OtDC/R7Kp/f/5ffaW59Q==
    "], {{{}, {}, 
      TagBox[
       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
        Opacity[1.], LineBox[CompressedData["
1:eJwl1HtsTncAx+FTr8sal5o1NJmKakNiiFJdNkWHjshCdRoyl+rYUqxrTNKy
KWJdu6ITVcTcpYpYNotWYk2mtigz2yqoS1ZxWaTZjGXYXJc9J/543ibn+8nv
pHnPexLeLsh6PyoIgiof4d+GSBB0JqtNENS4kMBx10uZwWiGYg5kwTEfXXjT
hb3aPjS6XsZMxpDCSU13cnS1msGcc30DC5hMW97x8Z0uhinafbpETtg+ZRZj
eU93WhPPPN1RTRpXbbsoZBjzdD9oejBbV6dJ5rxtY3gORZoL9pco0jTZJ3HP
VkcWUe2CYKHujCaRAl2jZiyt9v3Easo0t+zpVGha7XnEsEwTrVmqabEnU6y5
YJvKU/u3SIKHmnf5XtOVbN1+TRIn7eWsFnZiue66JpWPdS2amUTI0ZVoIizS
ndX0Y5HutH0CdzTb7Uls1kQxVVOjifCJ5qCm0t6dcvtt+2jW2v+wLyCWDN2/
mpX8p8vnJ10v5msbNCO4puvYPgi26nqzQfckfO7YpXuqWc7u8FnTrdU8T4mu
VZNGue43zVw68oEulhrdQPZoOzNXW6vt1ubZc17kzFzd4/D7ogNLdJd0A1mi
O6vLCr8zXQKHNel8o0tima5JN4gHzivQdGGnpi/bNO2YrjugeY5VzjuizdO1
pUoXR4Xub904qnR/hf8LpZr+HNO8wQldMmt0V3RpvKgr1MTxhSaFLzU9KNAd
1cSzU1epeZlfNG9xXpPOFt0dzURSw9+t5r59BU/s8zmliSNXc9hZQ/hQ15Ov
da9wSBdPoa5Rl0Rz+BvTRPO5phfrNY80mWzTPHTWR5RrBnFck8mPulTW6W7o
XmNT+I7QBFRoYsJnTXNT8yplmmvO2mwfSbM9l18146nWPNDkalba+1Jvz6DB
3p8STbMmhQ66hZoX2KMZQLWmE3N0hzRd2WHPoMWexw37JA5Enr0oK53zmX0o
p+zZNNmHs9H+uzPGh+87+0Ruhu8bbtmnUauJJt9xs+2PbKW0Z7H9om0Ai51x
xr02adI4Z8/hsuZ1dmv+0WRqiu29qbOP4og9kWLNz5qvbNn8Gd6Du7Yc6sN7
2attE7huy6fVNoWD9nrX/wf+wcdn
         "]], LineBox[CompressedData["
1:eJwV1nmcjfUewPEjhmlC9oxtUGO50mQpTOFa7lgmaowspRLSLUJFZSnj5VpC
6DKSKYQaDbJkmSxXSBERNxmUXRoqRiMhpvfvj8+cc76fz/N7nnPmeZ5zavQZ
3GVQoUgkUr9wJLKpaCQy0ePhWyKRlZ53x4WoSGQELheJRJ7GRr4p/yGXjFPc
IORyj2IlH4VSFn2ZL4/FfAIW86Xxbz7bGhUwn0/CUf45nOYfwVJNIcy2zjS+
MXbz3bGfb47Z/C/WSMbH/MM4y7+EX/meWKOJwWjvsQ9/g3sTMXidP8o1xBvW
yLGv2ZrmOMD3xhFNOyzUXNH00Izm70Q23xqb+NoYo/m/pgFu2t8QTSks0NTB
fE00ntIt1xTHXL4tfuCfxUm+M7I0BfxU+3uLb4idfFfs5RORrsnVtMNm+3te
UxSzNJUwTXNZkxw+R02+tV7BRM3d2KbpjJ26xpimO6FriTjda5pK+ERzP1Zo
KuFl3TZNDSzSzdQ0wz5NLxzUtMFc3SVNCuhIL3/y+TdwjX8WX2hKo5smy1rx
GK6rjBW6plilq4Jhuu26O7HTegM00ZitqYK3NVc0nTBH80d4D5ik7acr0ExB
CaRpT+ma4D/ao7onMUYXj/Wattisq4uxugO6RiiiK6NZxP8DC/iYcH1oVvIl
MV3T236v8uNQBEM132nqYKjmG00y4qy1mm+BbE1NjNJ8o6mHi9YpoZnHx2MO
fwt6aDLD8WC8dVbqZvB3YBJ/gW+D6fwv/EBs5jpiO5+ASfwPfCLK81lcA2Rx
ZfE8v56riDl8kn0UmA/CHi4OA7gtXAs8YNs1XByGc19zdXDSdjW4WeY3zFOx
IFwntknDQr5M+Kz4XK45JvFnuGfQjtvK3Y3xXA53f7iGuI/MSuIZ8zVmZZEe
/vfWHGF+2PwejDD/zjwVG8zjMdp8n3kCrurnmhVFL/OlZrdisn4997tZe6Sb
5YX7Gxpq3zI7ZtYcVcwqmg0x+9zravjArJXZe2YXvX4YTawX63Ufz9fx8zz/
PZxfnr/q+Veex+NgOCaP18yvexzl8bTH1igWrms+2myQx/0eU7CXq4/LZlmh
w5vma71Oss0Mry94PTgcK14I936zyng/fD66u8yGeP2V1/9CblgLaTjGNQr3
VP4Q1xMF5p+jGLZzZdCdX8LVwtfmk/E02qFpOG7cihjchuLhPEdJ3I5SKB3O
C5RFOZRHBdyBiuEzRCVURhVURTXEoXo458K2jqWHY1nqOGpjl9mU8H2A9mgW
1tD01WTzjZFj9i4Goyvi+Rf5HeEawDmzJRiD41xjpPGHucfCFyq2eKiJL7ly
6MkvM6+D3eZ/mnVGhvmVcN8L93rz0+E7FuPMj5k9hSj05YY5jw5wdTGM32P+
EPK4DK5w+N7jFnNRmMCtCp+9+UXztnjb/Nfwf0eFcA67tpbw5TCA38DHIoPr
YNsIP5jby1XHQG4r1xIPcmu56hjB7eLq4lR479w75jfNu2KheUG4x2IRP50v
i/Gac3wLTNb8xPfHJq49tvH1MYE/FK4dlOAzuXuQyd2O/vxarhxmhf8z1w37
uAfxDn+e64jh9v8XNxHRGMkf4RMwkj9g+66o6b67jm+FjZpaSNPs19wbrstw
rmo+4GtjHl8MT2iW8TFowx/h+uME1wkf8ze5KfYxhW+AHXwq9vDNMFPzsyYJ
G+znOU0U0jWxmKrJ13TALM0la2Xw/0QO3xdHNR3xkeaaZqhmAl8PW/lO2ME3
wlTNcU0LVA33IU0slmnuw3JNLF7UbQm/EbCAb4/j/AD8xHfBJ5rC4R5vnRl8
U3zLP47v+dZ4n8+zxiPhPsGnIDd8T+IC/zjWaYoj0ft/kv+DG4Mb4VzFLq4y
+lnjM/tK1yRiv+YJ5GjaYr4mX3OfZiRfFav4RKzmq+E1zQ5NLSzkO+AEPxBn
+VQs1xQJ56HjeYGPQQYfh5n8dS4F88K1zqXiPPcK8sJ5gc/4647jPa4VDnH9
cIxLRib/l+1fD79d+AR8yadgN98E/9Wc0XzK9UQeNwpXuL74X/ju5zO5TjjD
DcF5rjs+5aPD7yTvoz9fCFP5UhjLn+UewMRw3+S64Bw3DBfD9yOy+ZOOcQ7X
Ege5PviR64AP+au2X809hkvhNzf+5Pphc/j+s/1YrjY2cknYwtXDOP6g7Vdw
3fAbNxz5XO/wHs3+BvmbmAw=
         "]], LineBox[CompressedData["
1:eJwl1X9wz3UcB/CPYZ1Wim52IdcQNqylG9Zifp6wNbnIr+TIz6tZI+6alFGn
iI2utfIjp4lifnWJtl3HylF+ZM5dTFSW8mNXHOnox+N9/fHwfr9fz+f38/n6
/lrixJnDcxtEUXTNPxWxUZRnfc2axJeNo2goXzeKolSWNIyi0zFRtEM2it9l
BVyXTaRSniHfIMvinGwmF2Qj2SGP5V73eFGewCfyh9kib0GuvMo1ymVPckk2
lyuy8eyRt5avkw3ijGwGdbIn2CKPYY17FMu7c1g+muPyTN6T17tGNh/Lh/Gr
fBb18jF8phPHCK/NM/IbsoXEkC8/JutAvmsccK93dDKokY/ne52BrNO5pjNI
Z568DTvlvdglT6RA51udZC66X67OHazWace7OlF4DfXKdBqwWt6Pk/LJnJVn
sVHnb/lC91siT2W/fDiH5D1ZqXNeZyCbXXqqTkOKdeJZrHNZpy/LdH5zreks
0ulElc5gqvVSeEPvlF46zfRm67Rgk04qm3SaM11vt04CpXpFOmkc0hlFjU5v
SvUu62RR6PlN0vlH/iZxzNf5USeNBTqnXGssr+q1Z7dOf6r0OlGod1yvW3jt
9JrprJcnsU7ehAk62+R3slxnr/veplei04oines6WZTG/P99mUtXnWp5Dgd1
0ijS+UmnD/frbJX1YLusFbPl1bK2lMnbudcr8puyaeyXxzMmfA/kSRzSKdW5
jxU6f8lzWKVzQ/4Sy3SayhfI62TpvC4/K5tALJ3l6+Vx4Tsr3yG/i2LZsx7f
WD5HdkKWzBzZEVk2V+SrZI3Cd0S2SRbLYtlO2VLzP8wHssK83jz8uCSwWRbP
c7IKWUtWmQ/xuDzzo+aJPG++z7wP52T/hu8gH5pH4fNgXmZ70SyTpWbnzabS
lLtda6r5LvN4SsJro19rlso8sxNmI7hpPt6s3CyOt8JvoNkQsxLnq+F3ijbk
m+0N7yfrnTP0Jtvvse/OAfuOnDT/wHrLbD6/2A/gffOfnSeR5twkfA6tK62D
rXny6+ziBbOWrGWk80HrY1wKL4HuIvst9hesM2jOQudCWa39OGLYZ9aeseZb
nZM57Lw8PHeG0tb/K1deLetLndkGCniU1vIp8i9kPTgVPgPhOfIUD9CBjnQi
iWQ604WujPP4beGzxxHnIqaQRS9SeFpnu7wLR52Lw+8S2fQOe3mFrCe1zqvJ
ZxQPhvcx/F0K30e+c17BNB4nM+zllbJ0TjuvYRajSWWneQrH7FcynRy+MusX
3kf7j5hHH6rMHuEH+7XMZgzfhPeTy/blLOKh8BzNaqxvMyM8L+er1k8ZRn/n
89aNvMxI51vWSvqS4XwmfLa4PXwurX/yeXh/qWcr94S/39ZquvEfm8T5BA==

         "]]},
       Annotation[#, "Charting`Private`Tag$2527#1"]& ]}, {
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.], PointBox[CompressedData["
1:eJwl2dMSINgOBdC2bdu2bdu2bdu2bdu2Ztq27Z427+q6D6tSdT4gydlJ2Kh9
pXYBAwQIsDFYgACB1MAEISieAgQnBCEJRWjCEJZwhCcCEYlEZKIQlWhEJwYx
iUVs4hCXeMQnAQlJRGKSkJRkJCcFKUlFatKQlnSkJwMZyURmspCVbGQnBznJ
RW7ykJd85KcABSlEYYpQlGIUpwQlKUVpylCWcpSnAhWpRGWqUJVqVKcGNalF
bepQl3rUpwENaURjmtCUZjSnBS1pRWva0JZ2tKcDHelEZ7rQlW50pwc96UVv
+tCXfvRnAAMZxGCGMJRhDGcEIxnFaMYwlnGMZwITmcRkpjCVaUxnBjOZxWzm
MJd5zGcBC1nEYpawlGUsZwUrWcVq1rCWdaxnAxvZxGa2sJVtbGcHO9nFbvaw
l33s5wAHOcRhjnCUYxznBCc5xT/8y2nOcJZznOcCF7nEZa5wlWtc5wY3ucVt
7nCXe9znAQ95xGOe8JRnPOcFL3nFa97wlne85wP/8ZFPfOYLX/nGd37wk1/8
5g8BNICABCIwQQhKMIITgpCEIjRhCEs4whOBiEQiMlGISjSiE4OYxCI2cYhL
POKTgIQkIjFJSEoykpOClKQiNWlISzrSk4GMZCIzWchKNrKTg5zkIjd5yEs+
8lOAghSiMEUoSjGKU4KSlKI0ZShLOcpTgYpUojJVqEo1qlODmtSiNnWoSz3q
04CGNKIxTWhKM5rTgpa0ojVtaMvfxt2eDnSkE53pQle60Z0e9KQXvelDX/rR
nwEMZBCDGcJQhjGcEYxkFKMZw1jGMZ4JTGQSk5nCVKYxnRnMZBazmcNc5jGf
BSxkEYtZwlKWsZwVrGQVq1nDWtaxng1sZBOb2cJWtrGdHexkF7vZw172sZ8D
HOQQhznCUY5xnBOc5BT/8C+nOcNZznGeC1zkEpe5wlWucZ0b3OQWt7nDXe5x
nwc85BGPecJTnvGcF7zkFa95w1ve8Z4P/MdHPvGZL3zlG9/5wU9+8Zs//B38
AQlEYIIQlGAEJwQhCUVowhCWcIQnAhGJRGSiEJVoRCcGMYlFbOIQl3jEJwEJ
SURikpCUZCQnBSlJRWrSkJZ0pCcDGclEZrKQlWxkJwc5yUVu8pCXfOSnAAUp
RGGKUJRiFKcEJSlFacpQlnKUpwIVqURlqlCValSnBjWpRW3qUJd61KcBDWlE
Y5rQlGY0pwUtaUVr2tCWdrSnAx3pRGe60JVudKcHPelFb/rQl370ZwADGcRg
hjCUYQxnBCMZxWjGMJZxjGcCE5nEZKYwlWlMZwYzmcVs5jCXecxnAQtZxGKW
sJRlLGcFK1nFatawlnWsZwMb2cRmtrCVbWxnBzvZxW72sJd97OcABznEYY5w
lGMc5wQnOcU//MtpznCWc5znAhe5xGWucJVrXOcGN7nFbe5wl3vc5wEPecRj
nvCUZzznBS95xWve8JZ3vOcD//GRT3zmC1/5xnd+8JNf/OYPf5f+gAQiMEEI
SjCCE4KQhCI0YQhLOMITgYhEIjJRiEo0ohODmMQiNnGISzzik4CEJCIxSUhK
MpKTgpSkIjVpSEs60pOBjGQiM1nISjayk4Oc5CI3echLPvJTgIIUojBFKEox
ilOCkpSiNGUoSznKU4GKVKIyVahKNapTg5rUojZ1qEs96tOAhjSiMU1oSjOa
04KWtKI1bWhLO9rTgY50ojNd6Eo3utODnvSiN33oSz/6M4CBDGIwQxjKMIYz
gpGMYjRjGMs4xjOBiUxiMlOYyjSmM4OZzGI2c5jLPOazgIUsYjFLWMoylrOC
laxiNWtYyzrWs4GNbGIzW9jKNrazg53sYjd72Ms+9nOAgxziMEc4yjGOc4KT
nOIf/uU0ZzjLOc5zgYtc4jJXuMo1rnODm9ziNne4yz3u84CHPOIxT3jKM57z
gpe84jVveMs73vOB//jIJz7zha984zs/+MkvfvOHvx/+gAQiMEEISjCCE4KQ
hCI0YQhLOMITgYhEIjJRiEo0ohODmMQiNnGISzzik4CEJCIxSUhKMpKTgpSk
IjVpSEs60pOBjGQiM1nISjayk4Oc5CI3echLPvJTgIIUojBFKEoxilOCkpSi
NGUoSznKU4GKVKIyVahKNapTg5rUojZ1qEs96tOAhjSiMU1oSjOa04KWtKI1
bWhLO9rTgY50ojNd6Eo3utODnvSiN33oSz/6M4CBDGIwQxjKMIYzgpGMYjRj
GMs4xjOBiUxiMlOYyjSmM4OZzGI2c5jLPOazgIUsYjFLWMoylrOClaxiNWtY
yzrWs4GNbGIzW9jKNrazg53sYjd72Ms+9nOAgxziMEc4yjGOc4KTnOIf/uU0
ZzjLOc5zgYtc4jJXuMo1rnODm9ziNne4yz3u84CHPOIxT3jKM57zgpe84jVv
eMs73vOB//jIJz7zha984zs/+MkvfvOHv2FfQAIRmCAEJRjBCUFIQhGaMIQl
HOGJQEQiEZkoRCUa0YlBTGIRmzjEJR7xSUBCEpGYJCQlGclJQUpSkZo0pCUd
6clARjKRmSxkJRvZyUFOcpGbPOQlH/kpQEEKUZgiFKUYxSlBSUpRmjKUpRzl
qUBFKlGZKlSlGtWpQU1qUZs61KUe9WlAQxrRmCY0pRnNaUFLWtGaNrSlHe3p
QEc60ZkudKUb3elBT3rRmz70pR/9GcBABjGYIQxlGMMZwUhGMZoxjGUc45nA
RCYxmSlMZRrTmcFMZjGbOcxlHvNZwEIWsZglLGUZy1nBSlaxmjWsZR3r2cBG
NrGZLWxlG9vZwU52sZs97GUf+znAQQ5xmCMc5RjHOcFJTvEP/3KaM5zlHOe5
wEUucZkrXOUa17nBTW5xmzvc5R73ecBDHvGYJzzlGc95wUte8Zo3vOUd7/nA
Rz7xmS985Rvf+cFPfvGbP/wN+QP+zf4JTBCCEozghCAkoQhNGMISjvBEICKR
iEwUohKN6MQgJrGITRziEo/4JCAhiUhMEpKSjOSkICWpSE0a0pKO9GQgI5nI
TBayko3s5CAnuchNHvKSj/wUoCCFKEwRilKM4pSgJKUoTRnKUo7yVKAilahM
FapSjerUoCa1qENd6lGfBjSkEY1pQlOa0ZwWtKQVrWlDW9rRng50pBOd6UJX
utGdHvSkF73pQ1/60Z8BDGQQgxnCUIYxnBGMZBSjGcNYxjGeCUxkEpOZwlSm
MZ0ZzGQWs5nDXOYxnwUsZBGLWcJSlrGcFaxkFatZw1rWsZ4Nwf5/t9rEZrbw
P6OVUR4=
        "]]}, {}}, {
      {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
       Opacity[1.]}, {}}}], {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{3000., 3.9943029533487583`*^-7},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
     Charting`ScaledFrameTicks[{Identity, Identity}]}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{3000, 4000}, {3.9943029533487583`*^-7, 6.819193219987319*^-7}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.7682430588455515`*^9, {3.7682431166410356`*^9, 3.768243135486945*^9}, 
   3.768243529546382*^9, 3.7682435606681747`*^9, {3.768243600825775*^9, 
   3.7682436376044683`*^9}, {3.768243716202337*^9, 3.768243730741431*^9}, 
   3.768243801253951*^9, {3.768243911465245*^9, 3.768243916832893*^9}, {
   3.7682439581324754`*^9, 3.7682440494006968`*^9}, 3.768244086632202*^9, 
   3.7683910898293076`*^9, 3.768396455565508*^9, 3.7683964962287903`*^9, 
   3.768542936127215*^9, 3.7685432555342445`*^9, 3.7685442275064507`*^9, 
   3.7685442588754454`*^9, 3.768544333540786*^9, {3.7685446748403177`*^9, 
   3.7685446829695807`*^9}, {3.7685447153480153`*^9, 
   3.7685447230703335`*^9}, {3.768544820728234*^9, 3.768544847998358*^9}, 
   3.7685454336392765`*^9, {3.7685454945943155`*^9, 3.7685455079067216`*^9}, {
   3.7685455509376764`*^9, 3.7685455963662157`*^9}, 3.7685458074518476`*^9, 
   3.7685458868156705`*^9, 3.7685459540898037`*^9, 3.768548131293475*^9, {
   3.7685482401863365`*^9, 3.7685482571400075`*^9}, {3.76854842462544*^9, 
   3.768548463396747*^9}, 3.7685604185101895`*^9},
 CellLabel->"Out[74]=",ExpressionUUID->"863e5267-4792-4c73-88f7-cf403592ee7c"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   RowBox[{"dgamma", "[", "x", "]"}], ",", "x", ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{"xmin", "<", "x", "<", "xmax"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.768243742674552*^9, 3.7682437731480513`*^9}, {
   3.768244104933549*^9, 3.7682441630531597`*^9}, {3.7682443392505217`*^9, 
   3.7682443620356054`*^9}, {3.76854672978654*^9, 3.7685467321013517`*^9}, {
   3.7685484997306385`*^9, 3.768548500964332*^9}, {3.768548546767845*^9, 
   3.768548575816181*^9}, 3.7685551691781044`*^9},
 CellLabel->"In[75]:=",ExpressionUUID->"1f401330-d01a-41d9-92dc-375db5b59cd7"],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.7685442239689226`*^9, 3.7685446205664234`*^9, 
  3.768546593251582*^9, 3.7685467905790057`*^9, 3.7685485642550945`*^9, 
  3.7685669445683055`*^9},
 CellLabel->"Out[75]=",ExpressionUUID->"56317672-98f8-4319-a837-3565218e385c"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"FortranForm", "[", "%", "]"}]], "Input",
 CellChangeTimes->{{3.7685564379611263`*^9, 
  3.768556440645036*^9}},ExpressionUUID->"97837c7e-6026-4a3f-9cfa-\
273fc5caf3c8"]
},
WindowSize->{1536, 781},
WindowMargins->{{-8, Automatic}, {Automatic, 0}},
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 5635, 164, 1150, "Input",ExpressionUUID->"037a651f-7204-4dd6-865b-a90153e54791"],
Cell[6196, 186, 4010, 129, 291, "Input",ExpressionUUID->"8d2adb61-5946-4313-bd6e-7ea714537656"],
Cell[CellGroupData[{
Cell[10231, 319, 775, 17, 28, "Input",ExpressionUUID->"1bb2a22e-ee7e-457c-9fa5-3ad0259e7bb1"],
Cell[11009, 338, 13970, 252, 261, "Output",ExpressionUUID->"2b432197-b4be-4234-b069-ea0ceae093db"]
}, Open  ]],
Cell[24994, 593, 1734, 44, 67, "Input",ExpressionUUID->"6c45ba6a-a2e4-4dfa-bf7a-543610b4d110"],
Cell[26731, 639, 591, 17, 28, "Input",ExpressionUUID->"04e4efe9-ea24-435a-a90c-3bf1b8c0b336"],
Cell[CellGroupData[{
Cell[27347, 660, 967, 21, 28, "Input",ExpressionUUID->"44f1aff7-c138-4f4d-86fa-5ca7df87f6ff"],
Cell[28317, 683, 16902, 297, 236, "Output",ExpressionUUID->"b94229de-f6a6-4f4b-96f7-f3a03ecef3ed"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45256, 985, 2266, 62, 67, "Input",ExpressionUUID->"a731069b-1b83-412b-b454-26db50b1cc27"],
Cell[47525, 1049, 783, 19, 59, "Output",ExpressionUUID->"58736143-4ff2-4ace-b27c-23f706b8596d"]
}, Open  ]],
Cell[48323, 1071, 591, 15, 28, "Input",ExpressionUUID->"6fc8f691-fb14-4e78-ab99-1dc2d62ef10b"],
Cell[48917, 1088, 843, 22, 28, "Input",ExpressionUUID->"cc0a0d90-49af-490d-8c27-e1b8c06612e1"],
Cell[49763, 1112, 1352, 36, 48, "Input",ExpressionUUID->"3223564e-c632-41a1-86b2-1263eda30155"],
Cell[51118, 1150, 838, 17, 28, "Input",ExpressionUUID->"5a733281-bb2b-4eac-bf46-072d6163fbcd"],
Cell[CellGroupData[{
Cell[51981, 1171, 212, 3, 28, "Input",ExpressionUUID->"de858e80-3c44-4ea0-a374-ba372a5e43fa"],
Cell[52196, 1176, 618, 9, 32, "Output",ExpressionUUID->"1893da70-0168-4162-8384-e0660be1306e"]
}, Open  ]],
Cell[52829, 1188, 272, 4, 28, "Input",ExpressionUUID->"d5c88e0e-f64c-4c7b-a72d-9d86409b15d6"],
Cell[CellGroupData[{
Cell[53126, 1196, 1072, 18, 28, "Input",ExpressionUUID->"2c1e3732-16d5-42ed-930d-8f72fb168b51"],
Cell[54201, 1216, 512, 11, 21, "Message",ExpressionUUID->"2298e405-8520-4a21-b61c-4f7574b19f67"],
Cell[54716, 1229, 45465, 766, 216, "Output",ExpressionUUID->"863e5267-4792-4c73-88f7-cf403592ee7c"]
}, Open  ]],
Cell[CellGroupData[{
Cell[100218, 2000, 636, 11, 28, "Input",ExpressionUUID->"1f401330-d01a-41d9-92dc-375db5b59cd7"],
Cell[100857, 2013, 282, 4, 32, "Output",ExpressionUUID->"56317672-98f8-4319-a837-3565218e385c"]
}, Open  ]],
Cell[101154, 2020, 192, 4, 28, "Input",ExpressionUUID->"97837c7e-6026-4a3f-9cfa-273fc5caf3c8"]
}
]
*)

