(* 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
[encoded data, 4575 characters]
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
[encoded data, 9943 characters]
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
[encoded data, 33977 characters]
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"]
}
]
*)

