DarkRadiation.nb

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

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

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    142898,       2951]
NotebookOptionsPosition[    138592,       2872]
NotebookOutlinePosition[    138966,       2889]
CellTagsIndexPosition[    138923,       2886]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"Quit", "[", "]"}]], "Input",
 CellChangeTimes->{{3.745660166781345*^9, 3.745660167757299*^9}, {
  3.7488467175602493`*^9, 3.748846719023501*^9}, {3.784016038017746*^9, 
  3.78401603864161*^9}, {3.7845525560743*^9, 3.784552556592206*^9}},
 CellLabel->
  "In[336]:=",ExpressionUUID->"aacafbb2-6c83-4b7a-ac50-85d6ebd12eb8"],

Cell["\<\
This notebook can be used to find the thermal history of the Universe in the \
presence of massless dark radiation. 
This code is based on Escudero 1812.05605 and 2001.04466. \
\>", "Text",
 CellChangeTimes->{{3.784289428040432*^9, 3.7842894580615473`*^9}, {
   3.7845529786904488`*^9, 3.7845530561184473`*^9}, {3.784553518573971*^9, 
   3.784553520496688*^9}, {3.7845538147186527`*^9, 3.784553817206806*^9}, {
   3.78455452636663*^9, 3.784554559216297*^9}, {3.7866894918462477`*^9, 
   3.786689492491478*^9}, {3.786690651041564*^9, 3.786690658697277*^9}, {
   3.7866906888350143`*^9, 3.786690690147058*^9}, {3.786691273273799*^9, 
   3.786691277159902*^9}, 3.78669425558433*^9, {3.786694316759899*^9, 
   3.786694316893906*^9}, {3.786694651231285*^9, 3.786694664693595*^9}, {
   3.7870330006946917`*^9, 3.7870330171571093`*^9}, {3.788074134235181*^9, 
   3.7880741353144617`*^9}},ExpressionUUID->"59c482db-ea09-4e94-81d9-\
312f917c36af"],

Cell[CellGroupData[{

Cell["Basic Instructions", "Subsection",
 CellChangeTimes->{{3.74574668207099*^9, 3.745746689244767*^9}, {
   3.7463515306408567`*^9, 3.746351533698724*^9}, 3.7583737804255657`*^9, {
   3.76225827696426*^9, 3.762258285063748*^9}, {3.7624475759609547`*^9, 
   3.762447591042226*^9}, {3.7624517841319838`*^9, 3.762451791183298*^9}, {
   3.762455016080265*^9, 3.762455017214542*^9}, {3.762457724739826*^9, 
   3.76245772567489*^9}, {3.762685205023032*^9, 3.762685206310626*^9}, {
   3.762690970431069*^9, 3.7626909751148443`*^9}, {3.7653058744048*^9, 
   3.765305884212072*^9}, {3.76562752201742*^9, 3.765627526285842*^9}, {
   3.765696720125304*^9, 3.765696723059465*^9}, {3.765712801041895*^9, 
   3.765712813918517*^9}, {3.765720461911504*^9, 3.765720464523703*^9}, {
   3.766236221765518*^9, 3.766236224466456*^9}, {3.766236277151506*^9, 
   3.766236277154964*^9}, {3.766395423860496*^9, 3.76639544075397*^9}, {
   3.768278278583282*^9, 3.76827828351609*^9}, {3.781424955335102*^9, 
   3.7814249616364813`*^9}, {3.7842761259936438`*^9, 3.784276126750905*^9}, {
   3.786689169478221*^9, 3.786689171484433*^9}, 3.7866892107530823`*^9, 
   3.7866940903434153`*^9},ExpressionUUID->"a910b207-0afe-443b-a530-\
630ec823168a"],

Cell[TextData[{
 "The code works as follows. \n\n1) Evaluate the Basic Modules Section. \n2) \
Evaluate the Dark Radiation Section. It takes ~\.1d20 s per example \
evaluation.\n\nNote that everything is in MeV Units but the rates, the Hubble \
parameter and time which are in seconds.\n\nIn this notebook \
\[CapitalDelta]Neff is an input parameter. We describe this dark radiation as \
a massless fermion with g=2. \nEntropy conservation tells us that the initial \
temperature of the species describing the dark radiation fluid \[Alpha] is:\n",
 Cell[BoxData[
  RowBox[{
   RowBox[{"T\[Alpha]0", "/", "T\[Gamma]0"}], " ", "=", " ", 
   SuperscriptBox["\[CapitalDelta]\[CapitalDelta]Neff", 
    RowBox[{"1", "/", "4"}]]}]],ExpressionUUID->
  "97238487-4a0f-4700-a207-dd66f93b46df"],
 " where T\[Gamma]0 = 10 MeV."
}], "Text",
 CellChangeTimes->{{3.784289428040432*^9, 3.7842894580615473`*^9}, {
   3.7845529786904488`*^9, 3.7845530561184473`*^9}, {3.784553518573971*^9, 
   3.784553520496688*^9}, {3.7845538147186527`*^9, 3.784553817206806*^9}, {
   3.78455452636663*^9, 3.784554559216297*^9}, {3.7866894918462477`*^9, 
   3.786689492491478*^9}, {3.786690651041564*^9, 3.786690658697277*^9}, {
   3.7866906888350143`*^9, 3.786690690147058*^9}, {3.786691273273799*^9, 
   3.786691277159902*^9}, 3.78669425558433*^9, {3.786694316759899*^9, 
   3.786694316893906*^9}, {3.786694651231285*^9, 3.786694664693595*^9}, {
   3.7870330006946917`*^9, 3.787033012162702*^9}, {3.7875843153027067`*^9, 
   3.7875844582953043`*^9}, {3.787644849507729*^9, 3.787644852336698*^9}, 
   3.799555293273086*^9},ExpressionUUID->"e68671e4-6b35-4720-aee7-\
68c7917205e3"]
}, Closed]],

Cell[CellGroupData[{

Cell["Basic Modules (need to be loaded first)", "Subsection",
 CellChangeTimes->{{3.74574668207099*^9, 3.745746689244767*^9}, {
   3.7463515306408567`*^9, 3.746351533698724*^9}, 3.7583737804255657`*^9, {
   3.76225827696426*^9, 3.762258285063748*^9}, {3.7624475759609547`*^9, 
   3.762447591042226*^9}, {3.7624517841319838`*^9, 3.762451791183298*^9}, {
   3.762455016080265*^9, 3.762455017214542*^9}, {3.762457724739826*^9, 
   3.76245772567489*^9}, {3.762685205023032*^9, 3.762685206310626*^9}, {
   3.762690970431069*^9, 3.7626909751148443`*^9}, {3.7653058744048*^9, 
   3.765305884212072*^9}, {3.76562752201742*^9, 3.765627526285842*^9}, {
   3.765696720125304*^9, 3.765696723059465*^9}, {3.765712801041895*^9, 
   3.765712813918517*^9}, {3.765720461911504*^9, 3.765720464523703*^9}, {
   3.766236221765518*^9, 3.766236224466456*^9}, {3.766236277151506*^9, 
   3.766236277154964*^9}, {3.766395423860496*^9, 3.76639544075397*^9}, {
   3.768278278583282*^9, 3.76827828351609*^9}, {3.781424955335102*^9, 
   3.7814249616364813`*^9}, {3.7842761259936438`*^9, 3.784276126750905*^9}, {
   3.786689169478221*^9, 3.786689171484433*^9}, 
   3.7866892107530823`*^9},ExpressionUUID->"44eab734-09c8-49d6-a694-\
0efa772da8af"],

Cell["\<\
This Section loads the main common Modules: 
\t1) Thermodynamic Formulae 
\t2) Finite Temperature QED corrections
\t3) SM \[Nu]-\[Nu] and \[Nu]-e Interaction Rates
\t4) Constants and Parameters
\t
These modules are loaded from BasicModules.m. One can alternatively copy the \
content of BasicModules.nb here and load it locally.\
\>", "Text",
 CellChangeTimes->{{3.784289428040432*^9, 3.7842894580615473`*^9}, {
  3.7845529786904488`*^9, 3.7845530561184473`*^9}, {3.784553518573971*^9, 
  3.784553520496688*^9}, {3.7845538147186527`*^9, 3.784553817206806*^9}, {
  3.78455452636663*^9, 3.784554559216297*^9}, {3.7866894918462477`*^9, 
  3.786689492491478*^9}, {3.786690651041564*^9, 3.786690658697277*^9}, {
  3.7866906888350143`*^9, 3.786690690147058*^9}, {3.7866941159667883`*^9, 
  3.7866942301119633`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Load", " ", "Basic", " ", "Modules"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"SetDirectory", "[", 
     RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Quiet", "[", 
     RowBox[{"<<", "BasicModules.m"}], "]"}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.786689174046639*^9, 3.786689202901915*^9}, {
  3.786689509741823*^9, 3.7866895382585373`*^9}, {3.786690944976067*^9, 
  3.786690961052801*^9}, {3.786694212884294*^9, 3.7866942277537947`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"9494557c-ee11-41dc-82c2-e9af64174292"]
}, Closed]],

Cell[CellGroupData[{

Cell["Dark Radiation", "Subsection",
 CellChangeTimes->{{3.74574668207099*^9, 3.745746689244767*^9}, {
   3.7463515306408567`*^9, 3.746351533698724*^9}, 3.7583737804255657`*^9, {
   3.76225827696426*^9, 3.762258285063748*^9}, {3.7624475759609547`*^9, 
   3.762447591042226*^9}, {3.7624517841319838`*^9, 3.762451791183298*^9}, {
   3.762455016080265*^9, 3.762455017214542*^9}, {3.762457724739826*^9, 
   3.76245772567489*^9}, {3.762685205023032*^9, 3.762685206310626*^9}, {
   3.762690970431069*^9, 3.7626909751148443`*^9}, {3.7653058744048*^9, 
   3.765305884212072*^9}, {3.76562752201742*^9, 3.765627526285842*^9}, {
   3.765696720125304*^9, 3.765696723059465*^9}, {3.765712801041895*^9, 
   3.765712813918517*^9}, {3.765720461911504*^9, 3.765720464523703*^9}, {
   3.766236221765518*^9, 3.766236224466456*^9}, {3.766236277151506*^9, 
   3.766236277154964*^9}, {3.766395423860496*^9, 3.76639544075397*^9}, {
   3.768278278583282*^9, 3.76827828351609*^9}, {3.781424955335102*^9, 
   3.7814249616364813`*^9}, {3.7842761259936438`*^9, 3.784276126750905*^9}, {
   3.784379517871928*^9, 3.784379523470381*^9}, {3.786691290920047*^9, 
   3.786691291093912*^9}},ExpressionUUID->"6001c685-a11b-4e3b-a1b3-\
1976b70c9246"],

Cell[CellGroupData[{

Cell["\<\
Hubble Parameter and total Energy and Pressure of the Universe\
\>", "Subsubsection",
 CellChangeTimes->{{3.738664316543005*^9, 3.738664318701496*^9}, {
  3.747637483529047*^9, 3.7476374840383997`*^9}, {3.747639114123144*^9, 
  3.747639141305916*^9}, {3.74764575876602*^9, 3.7476457666042557`*^9}, {
  3.781427521782346*^9, 
  3.781427546699238*^9}},ExpressionUUID->"4b7b5908-48c7-4d67-a338-\
91b888461aa5"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"\[Rho]tot", ",", "ptot", ",", "Hubble", ",", "RHSec"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Rho]tot", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", 
      RowBox[{"\[Rho]BE", "[", 
       RowBox[{"T\[Gamma]", ",", "0"}], "]"}]}], "+", 
     RowBox[{"4", " ", 
      RowBox[{"\[Rho]FDM", "[", 
       RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}]}], "+", 
     RowBox[{"2", 
      RowBox[{"\[Rho]FD", "[", 
       RowBox[{"T\[Alpha]", ",", "0"}], "]"}]}], "+", 
     RowBox[{"2", 
      RowBox[{"\[Rho]FD", "[", 
       RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}], "+", 
     RowBox[{"4", 
      RowBox[{"\[Rho]FD", "[", 
       RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}], "-", 
     RowBox[{"Pint", "[", "T\[Gamma]", "]"}], "+", 
     RowBox[{"T\[Gamma]", " ", 
      RowBox[{"dPintdT", "[", "T\[Gamma]", "]"}]}]}], ")"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ptot", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", 
      RowBox[{"pBE", "[", 
       RowBox[{"T\[Gamma]", ",", "0"}], "]"}]}], "+", 
     RowBox[{"4", " ", 
      RowBox[{"pFDM", "[", 
       RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}]}], "+", 
     RowBox[{"2", 
      RowBox[{"pFD", "[", 
       RowBox[{"T\[Alpha]", ",", "0"}], "]"}]}], "+", 
     RowBox[{"2", 
      RowBox[{"pFD", "[", 
       RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}], "+", 
     RowBox[{"4", 
      RowBox[{"pFD", "[", 
       RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}], "+", 
     RowBox[{"Pint", "[", "T\[Gamma]", "]"}]}], ")"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Hubble", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"\[Rho]tot", "[", 
         RowBox[{
         "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
          ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}], ")"}], " ", 
       FractionBox[
        RowBox[{"8", " ", "\[Pi]"}], 
        RowBox[{"3", " ", 
         SuperscriptBox[
          RowBox[{"(", "Mpl", ")"}], "2"]}]]}], ")"}], 
     RowBox[{"1", "/", "2"}]], "FacMeVtos1"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"RHSec", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"-", "3"}], " ", 
    RowBox[{"Hubble", "[", 
     RowBox[{
     "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", ",",
       "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"\[Rho]tot", "[", 
       RowBox[{
       "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
        ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}], "+", " ", 
      RowBox[{"ptot", "[", 
       RowBox[{
       "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
        ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}]}], ")"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.7456552410298*^9, 3.745655363540448*^9}, {
   3.745658524575444*^9, 3.745658535755413*^9}, {3.74566016027192*^9, 
   3.7456602162364492`*^9}, {3.745660439526217*^9, 3.745660446592559*^9}, {
   3.745662075102223*^9, 3.745662087809806*^9}, {3.7456624686920557`*^9, 
   3.7456625215131683`*^9}, {3.745662571534659*^9, 3.745662608545594*^9}, {
   3.745665694072517*^9, 3.745665695815156*^9}, {3.745665984429*^9, 
   3.7456659846020193`*^9}, {3.7456662846772947`*^9, 
   3.7456662850025167`*^9}, {3.745667968982484*^9, 3.7456679743021193`*^9}, {
   3.745669758425774*^9, 3.745669763456752*^9}, {3.746264588735017*^9, 
   3.746264612749408*^9}, {3.74626467695746*^9, 3.7462648062644243`*^9}, {
   3.746350280341352*^9, 3.746350313135069*^9}, 3.746350358089541*^9, {
   3.7463515469036293`*^9, 3.746351652901125*^9}, {3.746351702876404*^9, 
   3.746351880349824*^9}, {3.7463519284020863`*^9, 3.7463520300678463`*^9}, {
   3.746352079156046*^9, 3.746352185719494*^9}, {3.746352292161129*^9, 
   3.746352423699407*^9}, {3.7463525149485607`*^9, 3.7463526968150597`*^9}, {
   3.7463527306860723`*^9, 3.746352891499804*^9}, {3.74635297250379*^9, 
   3.74635300686102*^9}, {3.7463530377542133`*^9, 3.746353040616948*^9}, {
   3.746353489150091*^9, 3.746353495850409*^9}, {3.746353575508936*^9, 
   3.7463536343674593`*^9}, {3.746353673545128*^9, 3.746353773029422*^9}, {
   3.746353813896914*^9, 3.746353944471158*^9}, {3.7463540767092943`*^9, 
   3.746354136413411*^9}, {3.746354424345664*^9, 3.746354445053331*^9}, {
   3.7463545397485647`*^9, 3.7463545462827187`*^9}, {3.746354709734251*^9, 
   3.746354711411407*^9}, 3.74635475472985*^9, 3.746354795025867*^9, {
   3.747560593576756*^9, 3.7475606042421713`*^9}, {3.7475731974419928`*^9, 
   3.747573201256133*^9}, {3.747645803094768*^9, 3.747645840512867*^9}, {
   3.7476458936640263`*^9, 3.747645926861782*^9}, {3.7622589537282476`*^9, 
   3.7622589539504023`*^9}, {3.762450479570212*^9, 3.762450481538858*^9}, {
   3.762455226789645*^9, 3.762455235552943*^9}, {3.762522491135914*^9, 
   3.762522515940456*^9}, {3.765301355715255*^9, 3.765301415057746*^9}, {
   3.765301469233863*^9, 3.7653014998908577`*^9}, {3.7653015383761053`*^9, 
   3.76530156814037*^9}, {3.765301610968997*^9, 3.765301670872017*^9}, {
   3.765302729323125*^9, 3.765302748483686*^9}, {3.7653032633812933`*^9, 
   3.7653033002096357`*^9}, {3.765303353951025*^9, 3.765303365857939*^9}, {
   3.765305486046583*^9, 3.765305486219564*^9}, {3.765696828167473*^9, 
   3.7656968324297543`*^9}, {3.765696865391028*^9, 3.765696947991127*^9}, {
   3.765697031780899*^9, 3.76569713707465*^9}, {3.7656971723531733`*^9, 
   3.765697183864621*^9}, {3.765697223085705*^9, 3.765697256688649*^9}, {
   3.765697297025337*^9, 3.7656973695222797`*^9}, {3.765697786360318*^9, 
   3.765697826624411*^9}, {3.7657062697704077`*^9, 3.7657062949300327`*^9}, {
   3.7657063708583107`*^9, 3.765706421119624*^9}, {3.7657183907335253`*^9, 
   3.765718392477724*^9}, {3.766237393624736*^9, 3.7662374052923403`*^9}, {
   3.7662374427731667`*^9, 3.766237592004656*^9}, {3.766395492567287*^9, 
   3.766395505436418*^9}, {3.766395557656231*^9, 3.766395580243944*^9}, {
   3.7663964117791862`*^9, 3.76639642838614*^9}, {3.7663967504865313`*^9, 
   3.7663967679573593`*^9}, {3.766749544908289*^9, 3.766749555645052*^9}, {
   3.7668127507597218`*^9, 3.7668127858334713`*^9}, {3.766812820810549*^9, 
   3.7668128707291327`*^9}, {3.766813211904985*^9, 3.766813220406299*^9}, {
   3.768278547727191*^9, 3.768278586064493*^9}, {3.76827868191238*^9, 
   3.768278720455064*^9}, {3.7682810391110783`*^9, 3.768281076280697*^9}, {
   3.768305455351981*^9, 3.768305463208398*^9}, {3.76830550395992*^9, 
   3.768305509839654*^9}, {3.76830579031933*^9, 3.768305790841123*^9}, {
   3.7814275145537233`*^9, 3.781427516164678*^9}, {3.7814275788949757`*^9, 
   3.781427592422401*^9}, {3.7814276522094927`*^9, 3.7814278021485147`*^9}, {
   3.781427880321557*^9, 3.781427887386751*^9}, {3.7814279318537397`*^9, 
   3.781428032327965*^9}, 3.784276108320692*^9, {3.7842761558434343`*^9, 
   3.784276190223151*^9}, {3.784276712762424*^9, 3.784276713354579*^9}, 
   3.784277991001051*^9, {3.784379658084126*^9, 3.784379705217134*^9}},
 CellLabel->"In[3]:=",ExpressionUUID->"8ceb9456-0881-4f2f-bed9-20971dfdd32d"]
}, Closed]],

Cell[CellGroupData[{

Cell["Temperature ODEs", "Subsubsection",
 CellChangeTimes->{{3.738664316543005*^9, 3.738664318701496*^9}, {
  3.747637483529047*^9, 3.7476374840383997`*^9}, {3.747639114123144*^9, 
  3.747639141305916*^9}, {3.74764575876602*^9, 3.7476457666042557`*^9}, {
  3.781427521782346*^9, 3.781427546699238*^9}, {3.784276254206637*^9, 
  3.784276276872074*^9}, {3.784276734216667*^9, 
  3.784276736342863*^9}},ExpressionUUID->"493c6bcd-d8bd-42af-8c28-\
b11f430f58ba"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{
    "dT\[Gamma]dt", ",", "dT\[Nu]dt", ",", "dT\[Nu]edt", ",", 
     "dT\[Nu]\[Mu]dt", ",", "dT\[Alpha]dt"}], "]"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"dT\[Gamma]dt", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"-", 
      FractionBox["1", 
       RowBox[{
        RowBox[{"2", " ", 
         FractionBox[
          RowBox[{"2", " ", 
           SuperscriptBox["\[Pi]", "2"], " ", 
           SuperscriptBox["T\[Gamma]", "3"]}], "15"]}], "+", 
        RowBox[{"4", 
         RowBox[{"d\[Rho]dTFDM", "[", 
          RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}]}], "+", 
        RowBox[{"T\[Gamma]", " ", 
         RowBox[{"d2PintdT2", "[", "T\[Gamma]", "]"}]}]}]]}], 
     RowBox[{"(", " ", 
      RowBox[{
       RowBox[{
        RowBox[{"Hubble", "[", 
         RowBox[{
         "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
          ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}], " ", 
        RowBox[{"(", " ", 
         RowBox[{
          RowBox[{"4", " ", "2", " ", 
           RowBox[{"\[Rho]BE", "[", 
            RowBox[{"T\[Gamma]", ",", "0"}], "]"}]}], "+", 
          RowBox[{"3", " ", "4", "*", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"\[Rho]FDM", "[", 
              RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}], "+", 
             RowBox[{"pFDM", "[", 
              RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}]}], ")"}]}], 
          "+", 
          RowBox[{"3", " ", "T\[Gamma]", " ", 
           RowBox[{"dPintdT", "[", "T\[Gamma]", "]"}]}]}], ")"}]}], " ", "+", 
       " ", 
       RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]e", "[", 
        RowBox[{
        "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
         ",", "\[Mu]\[Nu]\[Mu]"}], "]"}], "+", 
       RowBox[{"2", " ", 
        RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]\[Mu]", "[", 
         RowBox[{
         "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
          ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}]}], ")"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"dT\[Nu]dt", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Hubble", "[", 
        RowBox[{
        "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
         ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}]}], " ", 
      "T\[Nu]e"}], " ", "+", 
     RowBox[{
      FractionBox["1", "3"], " ", 
      FractionBox[
       RowBox[{
        RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]e", "[", 
         RowBox[{
         "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
          ",", "\[Mu]\[Nu]\[Mu]"}], "]"}], "+", 
        RowBox[{"2", 
         RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]\[Mu]", "[", 
          RowBox[{
          "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e",
            ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}]}], 
       RowBox[{"2", 
        RowBox[{"d\[Rho]FDdT", "[", 
         RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}]]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dT\[Nu]edt", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"Hubble", "[", 
       RowBox[{
       "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
        ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}]}], " ", "T\[Nu]e"}],
     " ", "+", " ", 
    FractionBox[
     RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]e", "[", 
      RowBox[{
      "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
       ",", "\[Mu]\[Nu]\[Mu]"}], "]"}], 
     RowBox[{"2", 
      RowBox[{"d\[Rho]FDdT", "[", 
       RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}]]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"dT\[Nu]\[Mu]dt", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       RowBox[{"Hubble", "[", 
        RowBox[{
        "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
         ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}]}], " ", 
      "T\[Nu]\[Mu]"}], " ", "+", " ", 
     FractionBox[
      RowBox[{"\[CapitalDelta]\[Rho]SM\[Nu]\[Mu]", "[", 
       RowBox[{
       "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
        ",", "\[Mu]\[Nu]\[Mu]"}], "]"}], 
      RowBox[{"2", 
       RowBox[{"d\[Rho]FDdT", "[", 
        RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}]]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dT\[Alpha]dt", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"-", 
     RowBox[{"Hubble", "[", 
      RowBox[{
      "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
       ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}]}], " ", 
    "T\[Alpha]"}]}], " ", ";"}]}], "Input",
 CellChangeTimes->{{3.784276745200658*^9, 3.7842767909186487`*^9}, {
  3.784277960766592*^9, 3.7842779644217978`*^9}, {3.7843797108113737`*^9, 
  3.784379758015731*^9}},
 CellLabel->"In[8]:=",ExpressionUUID->"e6b93132-914c-4e2f-ac34-041af44745a6"]
}, Closed]],

Cell[CellGroupData[{

Cell["Definition of observables", "Subsubsection",
 CellChangeTimes->{{3.738664316543005*^9, 3.738664318701496*^9}, {
  3.747637483529047*^9, 3.7476374840383997`*^9}, {3.747639114123144*^9, 
  3.747639141305916*^9}, {3.74764575876602*^9, 3.7476457666042557`*^9}, {
  3.781427521782346*^9, 3.781427546699238*^9}, {3.7842768094290943`*^9, 
  3.784276817034499*^9}},ExpressionUUID->"80a01eff-b49d-44fa-9421-\
4af4ea950cef"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{"Nefffunc", ",", "gstars", ",", "gstar", ",", "OmeganuMdiv"}], 
    "]"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Nefffunc", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     FractionBox["8", "7"], 
     SuperscriptBox[
      RowBox[{"(", 
       FractionBox["11", "4"], ")"}], 
      RowBox[{"4", "/", "3"}]], 
     FractionBox[
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"\[Rho]FD", "[", 
         RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}], "+", 
       RowBox[{"4", " ", 
        RowBox[{"\[Rho]FD", "[", 
         RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}], "+", 
       RowBox[{"2", " ", 
        RowBox[{"\[Rho]FD", "[", 
         RowBox[{"T\[Alpha]", ",", "0"}], "]"}]}]}], 
      RowBox[{"2", 
       RowBox[{"\[Rho]BE", "[", 
        RowBox[{"T\[Gamma]", ",", "0"}], "]"}]}]]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"gstars", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     FractionBox["45", 
      RowBox[{"2", 
       SuperscriptBox["\[Pi]", "2"], 
       SuperscriptBox["T\[Gamma]", "3"]}]], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", " ", 
        RowBox[{"sFD", "[", 
         RowBox[{"T\[Alpha]", ",", "0"}], "]"}]}], "+", 
       RowBox[{"2", " ", 
        RowBox[{"sFD", "[", 
         RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}]}], "+", 
       RowBox[{"4", 
        RowBox[{"sFD", "[", 
         RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}], "+", 
       RowBox[{"2", 
        RowBox[{"sBE", "[", 
         RowBox[{"T\[Gamma]", ",", "0"}], "]"}]}], "+", 
       RowBox[{"4", 
        RowBox[{"sFDM", "[", 
         RowBox[{"T\[Gamma]", ",", "0", ",", "me"}], "]"}]}]}], ")"}]}]}], 
   ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"gstar", "[", 
     RowBox[{
     "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
      ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
    RowBox[{
     FractionBox["30", 
      SuperscriptBox["\[Pi]", "2"]], 
     FractionBox[
      RowBox[{"\[Rho]tot", "[", 
       RowBox[{
       "T\[Gamma]", ",", "T\[Nu]e", ",", "T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]e", 
        ",", "\[Mu]\[Nu]\[Mu]", ",", "T\[Alpha]"}], "]"}], 
      SuperscriptBox["T\[Gamma]", "4"]]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"OmeganuMdiv", "[", 
    RowBox[{
    "T\[Gamma]_", ",", "T\[Nu]e_", ",", "T\[Nu]\[Mu]_", ",", "\[Mu]\[Nu]e_", 
     ",", "\[Mu]\[Nu]\[Mu]_", ",", "T\[Alpha]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        FractionBox["3", "2"], 
        FractionBox["rhocrit", 
         SuperscriptBox["T\[Gamma]00", "3"]], 
        FractionBox[
         SuperscriptBox["T\[Gamma]", "3"], 
         RowBox[{
          RowBox[{"nFD", "[", 
           RowBox[{"T\[Nu]e", ",", "\[Mu]\[Nu]e"}], "]"}], "+", 
          RowBox[{"2", 
           RowBox[{"nFD", "[", 
            RowBox[{"T\[Nu]\[Mu]", ",", "\[Mu]\[Nu]\[Mu]"}], "]"}]}]}]]}], 
       ")"}], "/.", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"rhocrit", "\[Rule]", " ", 
         RowBox[{"8.0992", " ", 
          SuperscriptBox["10", 
           RowBox[{"-", "11"}]], 
          SuperscriptBox["eV", "4"]}]}], ",", 
        RowBox[{"T\[Gamma]00", "\[Rule]", " ", 
         RowBox[{"2.7255", " ", "Kel"}]}]}], "}"}]}], "/.", 
     RowBox[{"{", 
      RowBox[{"Kel", "\[Rule]", " ", 
       RowBox[{
        RowBox[{"1", "/", 
         RowBox[{"(", 
          RowBox[{"1.16045", " ", 
           SuperscriptBox["10", "4"]}], ")"}]}], "eV"}]}], "}"}]}], "/.", 
    RowBox[{"eV", "\[Rule]", " ", "1"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.784276859569652*^9, 3.784276922151506*^9}, 
   3.784276964645434*^9, {3.784277064675457*^9, 3.7842771781045322`*^9}, {
   3.784277250208219*^9, 3.784277329435068*^9}, 3.784277390559286*^9, {
   3.784277501753249*^9, 3.7842775188550053`*^9}, {3.784282995223583*^9, 
   3.784283000183496*^9}, {3.784379570231998*^9, 3.784379609715898*^9}, {
   3.784379648635922*^9, 3.784379650627304*^9}, {3.7843808578638983`*^9, 
   3.7843808591264067`*^9}},
 CellLabel->"In[14]:=",ExpressionUUID->"99e71276-99cd-4934-a58a-86806b6fab9a"]
}, Closed]],

Cell[CellGroupData[{

Cell["SM+Dark Radiation Key Functions", "Subsubsection",
 CellChangeTimes->{{3.738664316543005*^9, 3.738664318701496*^9}, {
  3.747637483529047*^9, 3.7476374840383997`*^9}, {3.747639114123144*^9, 
  3.747639141305916*^9}, {3.74764575876602*^9, 3.7476457666042557`*^9}, {
  3.781427521782346*^9, 3.781427546699238*^9}, {3.784276254206637*^9, 
  3.784276276872074*^9}, {3.7843797829032907`*^9, 3.784379785493688*^9}, {
  3.786691371788769*^9, 
  3.786691373827161*^9}},ExpressionUUID->"77d4f13b-2bce-4415-a123-\
b5e9cd28b555"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Function", " ", "of", " ", "the", " ", "case", " ", "of", " ", "T\[Nu]"}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"FuncDRnudec", "[", "\[CapitalDelta]Neff_", "]"}], ":=", 
    RowBox[{"(", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{
       RowBox[{
        RowBox[{
        "\[CapitalDelta]Neff", " ", "is", " ", "an", " ", "input", " ", 
         RowBox[{"parameter", ".", " ", "We"}], " ", "describre", " ", "this",
          " ", "dark", " ", "radiation", " ", "as", " ", "a", " ", "massless",
          " ", "fermion", " ", "with", " ", "g"}], "=", 
        RowBox[{
         RowBox[{
         "2.", "\n", "\t", "Entropy", " ", "conservation", " ", "tells", " ", 
          "us", " ", "that", " ", "the", " ", "initial", " ", "temperature", 
          " ", "of", " ", "the", " ", "species", " ", "describing", " ", 
          "the", " ", "dark", " ", "radiation", " ", "fluid", " ", "\[Alpha]",
           " ", 
          RowBox[{"is", ":", "  ", "\[IndentingNewLine]", 
           RowBox[{"T\[Alpha]0", "/", "T\[Gamma]0"}]}]}], " ", "=", " ", 
         RowBox[{
          SuperscriptBox["\[CapitalDelta]\[CapitalDelta]Neff", 
           RowBox[{"1", "/", "4"}]], "\[IndentingNewLine]", "So", " ", 
          "that"}]}]}], ",", " ", 
       RowBox[{
        RowBox[{"when", " ", "\[CapitalDelta]Neff"}], " ", "=", " ", "1"}], 
       ",", " ", 
       RowBox[{"T\[Alpha]0", " ", "=", " ", 
        RowBox[{"T\[Gamma]0", " ", "=", " ", 
         RowBox[{"T\[Nu]0", "."}]}]}]}], "  ", "\[IndentingNewLine]", "*)"}], 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Tmax", "=", "20.1"}], ";", "        ", 
      RowBox[{"(*", " ", 
       RowBox[{
       "Start", " ", "integration", " ", "at", " ", "10", " ", "MeV"}], " ", 
       "*)"}], "\[IndentingNewLine]", 
      RowBox[{"Tmin", "=", "0.009"}], ";", " ", 
      RowBox[{"(*", " ", 
       RowBox[{"Finish", " ", "integration", " ", "at", " ", 
        RowBox[{
         RowBox[{"T\[Gamma]", "\[NonBreakingSpace]", "~", " ", "0.009"}], ".",
          " ", "By"}], " ", "then", " ", "all", " ", "electrons", " ", "have",
         " ", "annihilated", " ", "away"}], " ", "*)"}], 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"T\[Alpha]0", "=", 
       RowBox[{
        SuperscriptBox["\[CapitalDelta]Neff", 
         RowBox[{"1", "/", "4"}]], "Tmax"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"t0", "=", 
       RowBox[{"1", "/", 
        RowBox[{"(", 
         RowBox[{"2", " ", 
          RowBox[{"Hubble", "[", 
           RowBox[{
           "Tmax", ",", "Tmax", ",", "Tmax", ",", "0", ",", "0", ",", 
            "T\[Alpha]0"}], "]"}]}], ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"tmax", "=", 
       RowBox[{"1", "/", 
        RowBox[{"(", 
         RowBox[{"2", " ", 
          RowBox[{"Hubble", "[", 
           RowBox[{"Tmin", ",", 
            RowBox[{"Tmin", "/", "1.4"}], ",", 
            RowBox[{"Tmin", "/", "1.4"}], ",", "0", ",", "0", ",", 
            RowBox[{"Tmin", "/", "1.4"}]}], "]"}]}], ")"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"(*", " ", 
       RowBox[{"Solve", " ", "the", " ", "ODE", " ", "system"}], " ", "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"tCPU1", "=", 
       RowBox[{"TimeUsed", "[", "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"solEarlyUniverset", "=", 
       RowBox[{"NDSolve", "[", 
        RowBox[{
         RowBox[{"{", " ", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"T\[Gamma]", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{"dT\[Gamma]dt", "[", 
             RowBox[{
              RowBox[{"T\[Gamma]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", "0", ",", "0", ",", 
              RowBox[{"T\[Alpha]", "[", "t", "]"}]}], "]"}]}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{
             RowBox[{"T\[Nu]", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{"dT\[Nu]dt", "[", 
             RowBox[{
              RowBox[{"T\[Gamma]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", "0", ",", "0", ",", 
              RowBox[{"T\[Alpha]", "[", "t", "]"}]}], "]"}]}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{
             RowBox[{"T\[Alpha]", "'"}], "[", "t", "]"}], "\[Equal]", 
            RowBox[{"dT\[Alpha]dt", "[", 
             RowBox[{
              RowBox[{"T\[Gamma]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", 
              RowBox[{"T\[Nu]", "[", "t", "]"}], ",", "0", ",", "0", ",", 
              RowBox[{"T\[Alpha]", "[", "t", "]"}]}], "]"}]}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"T\[Gamma]", "[", "t0", "]"}], "\[Equal]", "Tmax"}], ",", 
           
           RowBox[{
            RowBox[{"T\[Nu]", "[", "t0", "]"}], "\[Equal]", "Tmax"}], ",", 
           RowBox[{
            RowBox[{"T\[Alpha]", "[", "t0", "]"}], "\[Equal]", 
            "T\[Alpha]0"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"T\[Gamma]", ",", "T\[Nu]", ",", "T\[Alpha]"}], "}"}], ",", 
         
         RowBox[{"{", 
          RowBox[{"t", ",", "t0", ",", "tmax"}], "}"}], ",", 
         RowBox[{"PrecisionGoal", "\[Rule]", "8"}], ",", 
         RowBox[{"AccuracyGoal", "\[Rule]", "8"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"tCPU2", "=", 
       RowBox[{"TimeUsed", "[", "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Print", "[", 
       RowBox[{"\"\<tCPU = \>\"", ",", " ", 
        RowBox[{"tCPU2", "-", "tCPU1"}]}], "]"}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"T\[Gamma]oftfunc", "[", "t1_", "]"}], ":=", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"T\[Gamma]", "[", "t", "]"}], "/.", "solEarlyUniverset"}], 
          ")"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"t", "\[Rule]", " ", "t1"}], "}"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"T\[Nu]oftfunc", "[", "t1_", "]"}], ":=", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"T\[Nu]", "[", "t", "]"}], "/.", "solEarlyUniverset"}], 
          ")"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"t", "\[Rule]", " ", "t1"}], "}"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"T\[Alpha]oftfunc", "[", "t1_", "]"}], ":=", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"T\[Alpha]", "[", "t", "]"}], "/.", "solEarlyUniverset"}], 
          ")"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/.", 
        RowBox[{"{", 
         RowBox[{"t", "\[Rule]", " ", "t1"}], "}"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"T\[Gamma]oftfunc", "[", "t", "]"}], ",", 
        RowBox[{"T\[Nu]oftfunc", "[", "t", "]"}], ",", 
        RowBox[{"T\[Alpha]oftfunc", "[", "t", "]"}], ",", "t0", ",", "tmax"}],
        "}"}]}], ")"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.784276388638754*^9, 3.784276436000514*^9}, {
   3.784276482741776*^9, 3.784276517077971*^9}, {3.7842765572169037`*^9, 
   3.7842766679618473`*^9}, {3.78427671094425*^9, 3.7842767429809628`*^9}, {
   3.784276796863751*^9, 3.784276798983423*^9}, {3.784277529869104*^9, 
   3.784277531128708*^9}, {3.784277593280299*^9, 3.784277748764703*^9}, {
   3.784277799984528*^9, 3.7842778248160353`*^9}, {3.784278154076406*^9, 
   3.784278192096439*^9}, 3.784278408469638*^9, {3.784278572480801*^9, 
   3.784278578367305*^9}, {3.784278948081503*^9, 3.7842789641685677`*^9}, {
   3.784286914475436*^9, 3.7842869164830637`*^9}, {3.784287128508923*^9, 
   3.784287128616454*^9}, {3.784287250028783*^9, 3.7842872750907993`*^9}, {
   3.7842874895272617`*^9, 3.784287489641776*^9}, {3.784287694734706*^9, 
   3.784287728874187*^9}, {3.7842877904001827`*^9, 3.7842877916779127`*^9}, {
   3.784287846513122*^9, 3.7842878476295557`*^9}, {3.784287879853305*^9, 
   3.784287910884501*^9}, {3.78428848579189*^9, 3.7842884877107773`*^9}, {
   3.784379797890789*^9, 3.784379934317215*^9}, {3.784379966300375*^9, 
   3.784379973232546*^9}, {3.784380180372258*^9, 3.784380217407034*^9}, {
   3.784380255817383*^9, 3.784380371826427*^9}, {3.784380564811451*^9, 
   3.784380573642219*^9}, {3.784380721106922*^9, 3.784380722897764*^9}, {
   3.784380757823461*^9, 3.784380792111672*^9}, {3.784380951568275*^9, 
   3.784380953797742*^9}, {3.7843809961365967`*^9, 3.784380999135358*^9}, {
   3.784381097407529*^9, 3.784381098813191*^9}, {3.784381137454171*^9, 
   3.784381155667902*^9}, {3.784381330942957*^9, 3.784381334277706*^9}, {
   3.784381477775881*^9, 3.784381477989942*^9}, 3.784382369480582*^9, 
   3.7843825295640163`*^9, {3.784382711353929*^9, 3.784382726822091*^9}, {
   3.7875842288537083`*^9, 3.787584229298276*^9}, {3.7875844673345957`*^9, 
   3.787584484945245*^9}},
 CellLabel->"In[19]:=",ExpressionUUID->"d0349130-a313-4165-8202-eb6e99433ed6"]
}, Closed]],

Cell[CellGroupData[{

Cell["SM+Dark Radiation Examples", "Subsubsection",
 CellChangeTimes->{{3.738664316543005*^9, 3.738664318701496*^9}, {
   3.747637483529047*^9, 3.7476374840383997`*^9}, {3.747639114123144*^9, 
   3.747639141305916*^9}, {3.74764575876602*^9, 3.7476457666042557`*^9}, {
   3.781427521782346*^9, 3.781427546699238*^9}, {3.784276254206637*^9, 
   3.7842763087390413`*^9}, {3.784288171272852*^9, 3.784288176661298*^9}, {
   3.784379786925622*^9, 3.784379792319029*^9}, 3.784382470625646*^9, {
   3.786691374867038*^9, 
   3.786691378139057*^9}},ExpressionUUID->"40d4ace1-368e-4d46-a7db-\
e71149dadc8f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Run", " ", "case", " ", "with", " ", 
     SubscriptBox["\[CapitalDelta]N", "eff"]}], " ", "=", " ", "1"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"resvec", "=", 
    RowBox[{"Quiet", "[", 
     RowBox[{"FuncDRnudec", "[", "1", "]"}], "]"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.784277772039658*^9, 3.784277775263968*^9}, {
  3.784278221518071*^9, 3.784278240498653*^9}, {3.784278416059729*^9, 
  3.7842784189191713`*^9}, {3.784278836699937*^9, 3.7842788719260607`*^9}, {
  3.784279190280746*^9, 3.7842791982520227`*^9}, {3.7842792624573383`*^9, 
  3.784279276919652*^9}, {3.78428208196564*^9, 3.7842821052513123`*^9}, {
  3.784380379117023*^9, 3.784380380674103*^9}, {3.787584158352742*^9, 
  3.7875841841358643`*^9}},
 CellLabel->"In[20]:=",ExpressionUUID->"4319d5e2-33d4-4dde-b6ec-0667c929a531"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"tCPU = \"\>", "\[InvisibleSpace]", "15.954581`"}],
  SequenceForm["tCPU = ", 15.954581],
  Editable->False]], "Print",
 CellChangeTimes->{3.78428215908705*^9, 3.78428770273249*^9, 
  3.784287759792061*^9, 3.7842878199359922`*^9, 3.784287888460031*^9, 
  3.7842879492977057`*^9, 3.784289504220636*^9, 3.784380493322295*^9, 
  3.784380580069489*^9, 3.7843807868252573`*^9, 3.7843809726634398`*^9, 
  3.784381011964217*^9, 3.784381120668913*^9, 3.784381167177724*^9, 
  3.784381357701438*^9, 3.7843824001675463`*^9, 3.784382487211404*^9, 
  3.784382556232923*^9, 3.786691323599337*^9, 3.787584148082653*^9, 
  3.79673378613068*^9, 3.799555344980586*^9},
 CellLabel->
  "During evaluation of \
In[20]:=",ExpressionUUID->"09480279-923b-427c-9917-b2602fde0576"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "T\[Gamma]oft", ",", "T\[Nu]oft", ",", "T\[Alpha]oft", ",", "tofT\[Gamma]",
     ",", "T\[Gamma]dT\[Nu]ofT\[Gamma]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Gamma]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Nu]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Alpha]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "3", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tofT\[Gamma]", "[", "Tgam_", "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"InverseFunction", "[", "T\[Gamma]oft", "]"}], "[", "Tgam", 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Gamma]dT\[Nu]ofT\[Gamma]", "[", "Tgam_", "]"}], ":=", 
   RowBox[{"Tgam", "/", 
    RowBox[{"T\[Nu]oft", "[", 
     RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T\[Nu]dT\[Alpha]ofT\[Gamma]", "[", "Tgam_", "]"}], ":=", 
    RowBox[{
     RowBox[{"T\[Nu]oft", "[", 
      RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}], "/", 
     RowBox[{"T\[Alpha]oft", "[", 
      RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Neffval", "=", 
   RowBox[{"Nefffunc", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"gstarval", "=", 
   RowBox[{"gstar", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"gstarsval", "=", 
   RowBox[{"gstars", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"omeganuval", "=", 
   RowBox[{"OmeganuMdiv", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Tratio", "=", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], "/", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(N\), \(eff\)]\)     = \>\"", ",", " ", 
    "Neffval"}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(g\), \(\(*\)\(s\)\)]\)      = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"gstarsval", ",", "4"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{"\"\<\!\(\*SubscriptBox[\(g\), \(*\)]\)      = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"gstarval", ",", "4"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\[Sum]\!\(\*SubscriptBox[\(m\), \(\[Nu]\)]\)/\!\(\*SubscriptBox[\(\
\[CapitalOmega]\), \(\[Nu]\)]\)  = \>\"", ",", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"omeganuval", ",", "4"}], "]"}], " ", ",", " ", "\"\< eV\>\""}], 
   "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]e\)]\)   = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"Tratio", ",", "6"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Print", "[", 
    RowBox[{
    "\"\<\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]\[Mu]\)]\)   = \>\"", ",", " ", 
     RowBox[{"SetPrecision", "[", 
      RowBox[{"Tratio", ",", "6"}], "]"}]}], "]"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"LogLogPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"T\[Gamma]dT\[Nu]ofT\[Gamma]", "[", "Tgam", "]"}], ",", 
     RowBox[{"T\[Nu]dT\[Alpha]ofT\[Gamma]", "[", "Tgam", "]"}]}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"Tgam", ",", "10", ",", "0.01"}], "}"}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<\!\(\*FractionBox[\(T\[Gamma]\), \(T\[Nu]\)]\)\>\"", ",", 
      "\"\<\!\(\*FractionBox[\(T\[Nu]\), \(T\[Alpha]\)]\)\>\""}], "}"}]}], 
   ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", "\"\<T\[Gamma](MeV)\>\"", "}"}]}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", 
      RowBox[{"{", 
       RowBox[{"Blue", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.784282164027906*^9, 3.784282279509446*^9}, {
   3.784282320628632*^9, 3.784282345259815*^9}, {3.78428410458879*^9, 
   3.78428410487086*^9}, {3.7842842030978107`*^9, 3.7842842729937563`*^9}, {
   3.784284336300796*^9, 3.784284361155418*^9}, {3.78428451854011*^9, 
   3.7842845192634478`*^9}, 3.7842868102056828`*^9, {3.784287687102829*^9, 
   3.784287687995302*^9}, {3.78428933398415*^9, 3.784289341230422*^9}, {
   3.7843803849139967`*^9, 3.784380460841353*^9}, {3.784380888328226*^9, 
   3.784380888424686*^9}},
 CellLabel->"In[21]:=",ExpressionUUID->"e9047461-364e-443c-99ae-825d269dffed"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(N\\), \\(eff\\)]\\)     = \"\>", 
   "\[InvisibleSpace]", "4.0523758740798455`"}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(N\), \(eff\)]\)     = ", 4.0523758740798455`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.799555345149469*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"400b2419-6969-4730-b85e-b56b5cfd935b"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(g\\), \\(\\(*\\)\\(s\\)\\)]\\)      = \
\"\>", "\[InvisibleSpace]", "4.5704099460833331747`4."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(g\), \(\(*\)\(s\)\)]\)      = ", 
   4.5704099460833331747`4.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.7995553451603193`*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"db7577a4-dcb7-45cc-bb54-160f7c44d740"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(g\\), \\(*\\)]\\)      = \"\>", 
   "\[InvisibleSpace]", "3.8406484298266816424`4."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(g\), \(*\)]\)      = ", 3.8406484298266816424`4.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.7995553451707077`*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"d53ea4d1-2aae-48b7-ace5-5c25eb65e77f"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[Sum]\\!\\(\\*SubscriptBox[\\(m\\), \
\\(\[Nu]\\)]\\)/\\!\\(\\*SubscriptBox[\\(\[CapitalOmega]\\), \\(\[Nu]\\)]\\)  \
= \"\>", "\[InvisibleSpace]", "93.0884682899143029999`4.", 
   "\[InvisibleSpace]", "\<\" eV\"\>"}],
  SequenceForm[
  "\[Sum]\!\(\*SubscriptBox[\(m\), \(\[Nu]\)]\)/\!\(\*SubscriptBox[\(\
\[CapitalOmega]\), \(\[Nu]\)]\)  = ", 93.0884682899143029999`4., " eV"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.799555345178672*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"c5861b73-6470-4ae1-9890-578ae0499f77"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(T\\), \
\\(\[Gamma]\\)]\\)/\\!\\(\\*SubscriptBox[\\(T\\), \\(\[Nu]e\\)]\\)   = \"\>", 
   "\[InvisibleSpace]", "1.3959713995655487206`6."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]e\)]\)   = ", 1.3959713995655487206`6.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.799555345185912*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"419323ae-504a-402c-8916-360b7a00e31f"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(T\\), \
\\(\[Gamma]\\)]\\)/\\!\\(\\*SubscriptBox[\\(T\\), \\(\[Nu]\[Mu]\\)]\\)   = \"\
\>", "\[InvisibleSpace]", "1.3959713995655487206`6."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]\[Mu]\)]\)   = ", 1.3959713995655487206`6.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.784382400310816*^9, 3.784382487353376*^9, 3.78438255638686*^9, 
   3.786691323754562*^9, 3.787584148274643*^9, 3.7967337862526007`*^9, 
   3.799555345193894*^9},
 CellLabel->
  "During evaluation of \
In[21]:=",ExpressionUUID->"bb25756a-610e-45ea-bb13-a370900e03a1"]
}, Open  ]],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          AbsoluteThickness[1.6], 
          RGBColor[1, 0, 0]], 
         LineBox[CompressedData["
1:eJxF13k8VF/YAPAZ+zJjbbOkLClaSSS6z5BKWiQUEUIikqwVlaUQCRUh2ffs
e6F7FLKmshbJvif7Gt47v/f9fN5/Zj7fzz1z7vM8c+5zzhU1vXXhGgOJRCok
PujfYAm5OdZ8KMqpl/2pfguYN02lHpf8fz/BEmJ/dPOisP9zdrxcYWgIL3L4
P7ezVdZqa/Ciaof/dYZMSAqViRdRjf/XnpeveX8u40E+fT3/2TlT51TvHR7U
rNb1n60Y1CirB3nQqa9V//lDjFi7x3duNOQq95/fKpe9Yj/KjVbDV5jpdqs8
Ni4lyI1exAsy0X3uXC1NY4EL8YMTme7tbedDbjRzoct3clf9CU8bt4345XAh
U3+9RboNEgQ3CMVwoTc+A7N0Vwxfwd4+40LJN+om6V70e6g985gLaTK7TNC9
Z2+spdJ94veif8foDrnd/7zOmgv5/j44SHctP0sKvxkXCtSf6aV7rWBnmcFl
LhQ7bvWbblm9U9/jtbjQOtuLDrotlm8MjalzIXe/G210u9etTAWqcSHvO/It
dEdEPv13kMaFEv6Nfqe7EcvkdVPgQmrPWb/QfdhrSokiyYWs5Wc//pffwfPt
F0S5UC5HHE53Rn+mY7gwF3rEtLuUbrETNhmS/FxIsbE2n+7O+dpTN7m4kN7Y
72y6Q5OlBvPYudCZoJp0ujnZh0RUSER+y5QEuqveHS/xWaEi0SHT6P/ivZFw
6cs8FUUPeEXQPVt3NcjgDxXdstoeRHfXs06y808q+ugefofuMFCKKm2hIuPT
W+3p1p4MP8L4jYpu+hpa012tdckh6DMVFe2fMaTbk6GQu+0jFf0uUdOlWzlv
Q/rWD1Q0Yatylu6cDd/70/KpyEZcVplum8oDHlNZVMR5XfQg3ZLOgVsPv6Wi
Eu1cqf/q2XbmYmUsFSlTczbS/Sji8+ffz6iIfKBt2I8wdlrSXNKPilCWWyfd
SyuPSDcfU9GhF7GNdNteUVVcdqWizG0X8+nexRXbQnOhIq/MzkS6ez+s3/ax
p6I5y6+hdF/aXpa2wYqKvnVPONPN903opIE5YUWqBd31Hvf6Yo2p6Guxuw7d
tD4F4f0XqajircZ+uqWTcgNOHaMi/46CjieE+y/ySgdhVHRXvaqC7mhWu6pW
RSoC5qV0ujdY7V03O0BFLaJh9+j+IvD0ddpuKtqkO2lM95PaUYUpSSrqSlVR
o3ttd4rdg61UVBCZzk53cScLtXILFW2vLRj1JewQcC2VcwMVfT7wqpbu4Qmx
3jAOIh7PCm+6v+W80c5boCB9ec2/PvT66Z5b456hoNSrr6voTltcS7GZoKDf
a9aRdHuCydqOAeJ6id0xuid6eVI9uymo5lzwBroNvMu1f3dQkJQsd783YbkG
sdSw7xT05kHvfbpj7Jq05xooaNtDijrdlA2P1rRqKMj9/CkeugcuD2hzIgq6
FVYc8Zi+noeT1+5nUFDgx88hjwgzPNVL/ZlCQdGbO8/RbbufXUchgYJsLgUz
033S+UbqZAQFOQpftfYinC8gpHM2lIJ4LqoK0r29rG4tNZiCGgK1qzzp/YFx
r46ZLwVde562kW7z5F9ruBcF9XtXFnvQ66PxLFX4IQWNovd6dKcF/11rcaKg
bt3RZ+70fLflpqqbUZCWg/D3B/T1/9FUJ9GIgnaceqxPt5wF/zr5MgU9iu3p
vE/PN8NRp+Q8Bf11m211I3zn/I71zWcoyG3U9izdAzMtqY4nKWhrqDnuSu+f
RxTW92IU9NThcsg9+nrqGkr1U6Sg4GTppbv0+niE6QzJURCwbNKj+1b1UmrM
biKfhnKWO/R8L5bq8AtQ0JKJnKgT4adLNuu3NlBQJ16v50jPP3JrWj03kZ8X
t58DPd++B+uPWSho4/Jq+236evDZn9ZLJuqlajxvR89fulsHVjmRVVkZhW7P
2yppizOcyHvitrgtfT1smNbR/cuJBBw1xW/S61Mct54zyonCfRWFbOj1WGfU
te7hREZvzaat6Pk+rVr/9YUTbffNOXyNHt8eThe/Wk704dgldnP6+HrNCfkq
TnS7yuabKf1+nD86A8s4kavwwgkTwiJ+Y+9U0jkRb0RY+WXClVIHZCaSORHc
KjypT+8nNY4pEfGcqKLp9adLhN+zrYXORHCiw0KGr3Xoz7cPj2PiE04U2fet
/Cz9eZLUHdN6zIkGgg79OU04qSrcdM2dE/2MXuXUoPdHZnGtS3c4Ualig9gJ
wkGPDu1jv86J/p3VyDhKWEHiXmKBKSf6LRNgrUTvn58+CJsacSJO1zuCivT9
h1Gds0SXE8mntKnLEa7xuDxso8aJTugndkkRthONNhagcSK1K7NbdhLeUt7X
WqnEiRSH1VQk6PsN6WalyEFOlJfxRF+E/jw8fBj3VZQTKZzgiOAlnCpSIXB/
Kyc6qmZ1iIvw+Q9swVICnGjRseU9B+Go1eCHnjycaPW3licjYSW3BEO5dQ5k
IROzOKNH9D+hkabuZQ5UwSjc/pewX8ne0wHzHOi8hP6bMcLty4WKg3840FRy
QGsvYce7NZvCOjkQf+jvmUbCQgJcz9TaORCz/q6pWsIfiy8wTzVxoM+5I42V
hHkWO2ZO1XEgkSKOAyWEM5wnvq6840C98fqmCYRn50g/NQo40I6Gge9vCCs7
8feFZ3Oga8aaW14RrnU4PKeQTFwf1hd5QnjYzlPA8SUHOnkCrd0gvH/ypdin
QA4kPKj2z5Sw863k3Xz+HMjp1Lmqy4SZbeuPZntwoA9TgSWnCItZbzIdt+VA
OtdluHcQthrdZa10gwOp1thpCxHOtlJy9LvGgY5ear7MSxgsTbx3GXKgI6+V
M/5dagHDa2lp5qc4ULfaQ+WvhBP6S/Py1DgQw6U9ThWEx8waSxloHKhhU59F
MeF7prNfYuQ50FuXjJvRhF8ZYzOd4hzoV41/mxXhrq7z/3Zv40ASUY1JBoQl
jcyYXQU5kJkARews4XxD380CvBwoV9Fl837C3/W/K+musaPrUVqfJy62gMCP
frWEJXbkHsZR+Yuwid7C2ZlZdnR5KMS4nvDfi8ImwWPs6OKshnkKYYqOxaOG
dnb08ZWjuwFh7aY7AcLN7OiY58FDJwlHXPAPtW5kR2vVo86yhKW0slPYq9jR
vTtKr1gJnzi3XH88jx2ZyShT03VbwEP92cYPAeyIa/bVQodOCzBVyI0ZPmFH
9S38F8oJe0MHWnnEjowDeSWTCPvJ77JRdGNHOxc5DtoSfi7x8WP+DSKfFydy
FrSJ/TDGMkzbgh0Z3E+IbyccKsRtO32VHWn1964WE47gMxQ4oM+ORD+1d7kQ
jifN33p7kh2VHjpMnbhA1M8t8rjGMXaEXBYXaginLKgKjWDsyLHMSSGRcPpE
YJWkPDs6kzqUpE84v1Nqa5wEOwo0bfv7TovoJ8VXasPI7GjvoQ01uudbQF2O
KUZhlQ3JnNQM3UG4NivNqXWRDb0Z2+Yzq0ns10kL2zdMsqFg2RPVgYRbXwa7
BHaxIU4f8fsl51pAj+fw2X0/2ZD98aMbHxH+6d8l1tDChr5WixSeJtzlufsL
RwMbahzfyth+tgWG7ColvEvYUKxdn03/mRa4MW69LFHEhtqfao0nEB6/zvf1
Uy4betwWcdmc8KSRsSs5jQ2R1G3+/j5N9PczS9/cwtgQf9qxuzUaLcAutfeh
gxMbGtIllZuqE/391Q6zXbfZ0E0D3iZewvbMIid/2bChknK2SPwkEV8PF/dJ
cyL+hhSVzYTTIiajBLXZkG2DlUPhcSJethHPxnNsKIfV5PtlwnwuPRaPNNiQ
QPTNulU1Yr1rf983ocKGziXRHICwBiX/Q/kBNnRqZH4lW7UFHtxLj3Pew4YC
tcgBpwjnDid4797Fhq7jbs+6VVpAsDLkXMg2NkS2L2hiJzx636XLkosN+X8a
IZ2CFvD/e2SNe5wVnVI+u99JqQXQlYO9FUOs6NA7ldjpI/Tz7O6qu32sqDe5
VN6WsGGK8LO+n6zIXsI+11iR6MdX17YW1bCisRNS53YrEOe57+VHjZJZ0QWT
L29NZFuApPJelD+eFZ05W1hUJkPsR1m5zNVRrOizBmvAFsKR/vENMqGsyHV/
pMfn/cR+o/b4CvNjVvTocC2ZY28LUAvV76ebsiJF2fAgoZ1Ef2XknmE0YkXP
+p6dU5FsgUKtFksDfVZEUt7BYL6jBYz+XNXhOM+KKhd2iMWIt0CmhOtuy6Os
aIvdG53pbUQ/f5HxQ3wLK9oY3Mi9uLkFJHocNF35WVEAf4HiyCbi/953pPI7
FytitL7zrXVjC8TWVmV7MLOilWLv/cn8xPsFudvn9zQL8g1s/iPE3QIv7fjk
XzewoOlat+4zzC1g+aH97XQ1C2peLeMiMxH9lxItqlHBgoIYXHflMBD7dfJu
6tJ7FmTkePAPE6kF5H+r9V9MYUGnS97H3V1uhvazLs/5vFjQinGCadBEM6S/
Psp24wELStTeWrk23gzuI4wPyu+yoIAZfXfLsWaQehxkZWfHgsS8Tg7uHW6G
e6WptC9XWNAVdFjtYU8zCO/unHhymAWFHncyk2tqht1O8jIycizIveeYqeO3
ZjiCBzm072dBbdVyzVmNzaCnfXxRcicLWuilLm6pb4aXrpnkio0sqDXSIi2v
ohmo9e4b1qaZEbq6vImpgJh/U8fFxAlmdCKpcMdibjPsMTkUfmaUGe3cmUAb
ym4GjdmRra97mFHRxb2OhenN4C2svfPwV2a0cHK2fXNCM6zZ7Dhin8GMuII8
pXuCifmLHroJpDKj60VedzwDifnJPz+gBGZ0U0vWSCSAiDfk2THuSGYUnZLa
ety3GZzLFs6k+zEjK9bF6QsPmmGCWms8eJ0ZOfy8yXPVqhl+Z9o+vizGjLZ+
Zv8srtQMz9pD9RO2MiP+i24PahWaAWPA905sYUZfvH5Rbsg1Q6Qud6sHNzPi
Pfg97sVeoh7/siRT/zGhuWY19qBtzfD91NTnhVYmZF7F47uLoRmq+h04QvyZ
0B7hn8HF5U3gxBXZ1eXNhMRK9JFdWRNIHK7I3eXJhNy2AIi+awIv/w0GpXeY
0Mouhw3W2U1Aky1423eNCU3ssyuMjWqC9w/nz8jSmBBHgKV0xN0myBC8G9g4
y4jc10Q4c6SagEth5vb+SUYkFmNQ5yvRBHbaN3UCxxiRa640VX9bE8g+vSqg
2cuIPM5wi/dtaILC1VPxXxoZUefGTb5v1r8D6hYsakhjRMFNxyUzmr7D96TS
rjoTRvSZopJwxvk7LMiS9lXXM6BIEwXTkKRv8Pw01y+XagYUmL3JSyzmG+wz
F366s4IBuXBXyqaFfwOLkMOj3iUMqDfiyMXUp9+gdcEu6XgaA3ous8/9uv03
KCztEfnkw4DOkI6OjSh/A+cTn7g/qDKgY1RRW+UvX2FO33sqr5iMfn4q9e76
3Ajrlo92qQWQ0FMfI24Z9zpwpTKyB46t4YemuvbMXfsMXAJ/krsC/uH2qQKZ
9pKfIPz+11uHlJZx8ta9755Xf4Cn6zke9jyLuJjgh1auvGJoVgsJp7Au4lIW
UjIoohiE/e7kJK0u4I3z77W9PIshfQP0/BxdwEMN0v5iF4qhQaqedqxyAfdM
zNTyny4CLp3Bdf57C/iqpguzi3wRBKcIPijom8dLbD0/+NQVQKiW173Fwjmc
8sd+Y6ZUHugEPtstmzmHC/3qucWxOQ/4GsI7rRPncA85B3MnpjwIUM/Gfr+Y
wzOF3zC6/M6FR/CLocpuDmcaKnfcFZoL9nsV/F9Iz+ERvu/uAksuaLKPvd4X
NYu3bfKt2TCVDWzlF8rMH8/gzU807P7+zgDpfrJ/rdsMvlWgWHLvxww4zZqj
d8BxBmfQlOrwSsiAZ2e5Z1dMZ/CKPRzsblYZsPFnnfQL2gwucd6mM2Q2HcRm
1MLKV6bxtbdmFmbc6aC8Q8FBxH4an9rv8oPpUhoYqQ/SHllN489ln+ueVE6D
h9YhXKMm03jhXmVqrGgafMyZSS3QnMbNb15dDR9PhVNYdveZvdO4z8FCsxiv
VLh4Seqc6/AUznn0h9t0QQrcfiIk1W40hdf/2ffp5O5kaLK2EPTSncKlqjr3
avImw6FzOZz7zkzh9wu2NN1aSIIlvpMTXopTOLVl0PPfpyTwiLTP279hCo9l
8klWMUqC4Kyaoz41k3h3kl3ijVeJkNPickFebhJfYQmGIqEE4C/+eKxHehIn
GTegj0wJ4BRBPfRUdBIXdjA7OPAnHhSN4zf3ck3iQw1P0mxRPHh+e/NOqOkv
/uTkrtVgi3iIWzj9pfboX/yqi+4ejqI46D2WsrCLfwIHXn6ymmUsBC8FUD6Q
J/A/DA4PzHVjgZblIKo9+QdXk7ZifaUaC9GCcPp+wx+cM5iHU3ZrLBhOtUR9
8/mD/3hYfh1rioHWKMYT9/6N4wfOvu6JV42B2iWjl3UDY3iOHdc1T8MouJul
lmrSPIY3Rd3axqsUBbuuSX+Y+ziGL21Kc0cCUeDzdW5oWwwx3iPzm3vbG1BL
fqrkeHkMN1DynXmk8wY+6LzvFW4cxQMYNOO8dCIhL2ujjG3xCK7GKyeVYB0B
vX2yAQGJI7i4qKFd5OkI4NtyfiT9+Qje/u37iaLdEWDv7hc7ZjOCkzpzBHTG
w0FWa53XSnQED7t25p3krXDInR6eNvMbxgXPLnZscA2D3EOl+fqGQzgo2F9S
fhsKvVY/eO6eGsLdr2yvOBUYCnxR8zZh8kN48E1Kn5NDKNizyOxo4xnCxaUU
yqWUQ0G2LTFEp3IQP949wbj6JQRy7wQ6a+4bxHmjPBgOL7+EnBLTw8dJA3hQ
o/K6osULSHn/uaH0Tz/eEyUhoHfuBcS822Mm19GP+x7zsgiSfwFBRfMB4oX9
ODVH75Am6wuwzfPrI9v04+4sqayMac9hz9vcZx/a+nDJhyJue2eDISWCYUAh
qxf/cZN2ui8iCGLCr9/LiuzF0zz9Dbd4B0FYWD33Tr9ePN36Ybj57SDwCQ09
sulaLx5s2XBMVT0ILJ5LB80K9uIyu3cObJsIBHG/C0o53j34J9Prh8pznkHU
3bhgaaNu3FLpyvA9/afgYXVJv+dMNx5ww+z+Z4WnYK5PEQ1T6sa5u0qTFTY9
BWlF52xmgW488cVOFc8mf7CZehFVnvwbd90w+bdU0x9e/+0V2PeiCy9ZUZ1M
O+cHy2PuPKw3OvG3gXF9vtd8AdWNRjue7MSP+6R7sqv7gvdbnf29Ep34CIi8
T5X2BT7rXefKfnfgrToGt9QnfUBqtNHfQbcDP3/te8eKmw/oDYuwdtN+4rqO
fUvbI72hqL90tXhzO14fHzQSNvkIHpSHCClNt+FsFAnhlOZHcCLKVrGsvg0n
3WBQbi9+BM2XtjuVe7ThZt80siI8HsF0rdd4zVgrvrvGLy9hwyPYk6PR8RNv
wclJRRlZNC+Iu99e/O96E753b/2u91keYHU5p/WBahPuEbRFKOmVBxxQ8Jsl
bW3C64kXwsyHHoBPKh1g/v4dX3vJYbfpvAd0mkUnU5W/4wVHCsWPTLvDZg2L
UBGeb7ik8Z7Zd0rusGXv8JC7XQO+mf9NXmTlfdBJ+zHuId2Avy90oM5H3Yeg
nXVTnv31+KbkM+dc79wHVrHMlUd69fhLtc7fbHvuw/wmR+4ntDp8F9f3/N0h
btBEIikE89TgAi4mOVa3XOFZ82bv2OwKPPjX4GC06l2o1WL3j7tRgZ9cz5KS
EL0LzI3LgfESFbiRaIxI4/odcKv5FZ4Y9gl3iylxCSu7A9Yf4jNS3T/iY2HL
ObNH7oB6yv6WnPMI3/3gYXKOsguwLD5Tf9RThO92P5Qlb+AEJerZwup6Rfhq
ZtJsO+YEduHfJjkbC3G3PRcnX4s5QceRDeEvSgvwoEOTrl6jjpB9P3wk7lUe
nicnp9ni6gj5LOobvV5n4g0h/oqWCQ7w7ki/7MpQBG6j3uRrS7GHZ/etvDg+
heFfpbYWjMzcBjM00SQQFYq35p3M8ei4DdQTS46HdZ7jFd3UONa3t+FjnLDt
mLsXvpyewmd5+jb8NBbhh8j7cFHV0MQ02A60CxSNqB0J8OdrU+CUqS107X6o
2sOUBIlf8G2bD9mCZVyFZMG+ZOhaKta9w2oLboGafw09U6HU+kLmePpNSLC6
5p4ulQldPQP+Ics2sPJM/duO0jxYlkybyYyxhjPHHGcqOMrgbftE+/FNVhC8
FjqQGV0G739cNjD7awmt7961hcl9gHc7ptLyqy3B5MB6iY0RDtd3jpZPuVrC
e67jqhw/EPipaR9t7L8OX1YaOXL3fAIeBoHQAdwCFloGXpObP4PWv6FTS+Hm
MH6rqO6uQDWMqy3s1XE1h172JytTRtVAcoiRajU0hy9H9xj0jlRD18/1d+e2
m0NCor3gp/Ua4NQNt6OmmsF5p/WwR7vrodos8cObclPQ8aPUeal8hQS7jtEp
6lXQvDTXC+ZfgdW8PKt7xgQ0JLqWV7y/ghOsvJj6YQLwIUvasf4rvHa/vu1u
kglIT1/wM79EnNOeyG8UUjEB0uXwU8dvfgeFdY7pjUtGkC4lWc0S3gxKbWmq
qpsMIXEpWivzQzNMcwwtXB80gKgagQ7d/mYQIg25NBUaQLAVdSKBeC97h4/4
39MzAOfUWf5jVS2Q4qIevh51GVSkPhk/nGmFYH6z028O60PrLpOFhbM/wIix
dGr+9UX4uvjDPdrhB4h5vrQ67HgRaqq1OU6G/4DUWfaSgjMXodTy5NbQ/h8g
eWs61XRNF2JT9h2Tc/0JKmuFbffNdcFm1+ozu5QOYN8lyumqpAPkXRGSIwxd
gMkt+8TxXQCJTY6n+SS7oIJd+tT6jBacYDpnp3SqC7gkYsJ9W7TgaTf5fUBg
F7TztgpuD9eCLeFWZ2WFf8Psy2aBK2JasPmx8PsywW4YVGfh7VE+D/vmx7pZ
5HrgNkkmyFD0LND+6Dm8Pt4DPobGto4TZ0Crv5L5wKUesOOv/tVecgYcv0dL
69/rgYyb5iOPL52Bd5najm9RDxxZ1JK8/vw0qF0vYT1/thdeLkz6zW3UAIN2
v33hFn2QUVlFqdU6CQlM2T4tLn2w2X5iIHTXSZg40NLN+6QPatndEuPWT4Dn
k20v/N72gfTcz7dOmScg5Uj+gutkH5iuHXlF5j4B85G/kNG9frDr3hpv06EG
QVcPaEs8GwCKbqzCo2hV+Bmgm24SPQBrbXuqX99XBYn395jfZA+AhmH5yrCB
KhTxVRZtbBqAOttiMXEBVej6pC/MumUQhh2TV7lCVWC3pNfASOwgXE4tqX0S
ToOq0RaXrIIhOBF5kMfN7ygEhn16crV6CG40PnI3PXcU9E7kvObvGALDiwFf
8/iOwmjsU9yFNAxzD468eBmpDFz6x1hpp4fBMeuRKFagBBc/Z4d+7R6GQGnJ
hZQlRRiI9y+Y4hwFJ1MIaauVh0ytu5/jRUaBY7S0eD5CHpzXLX7oyoyCnZCF
tq21PLAZqK6+uzgKJ+fhjzlVHvbwLal5xo7CPf5W7au6h8DR3aKZT2EMXn+n
MYnPHQSmKyozB03HQTWH7cJpOxnoFfj9PclhHCby7zkynJUB1OqWK/B4HCZv
cE4sSsmA6/ni2+vJ4zDSP+0Q138Apo/tn6z9Mw7p2YOXfxkegG7prX+u3v0D
eVf2PDqWsw9KFhcHA4MmYDnlWEvXnDSEFYRWMcRNwMvwLxnlZdLgbC+X5JQ3
AQkVs0sUb2mQHb91zaB1AsrnFk1/bpaGlN+DfZLCf2Hru1E9jCYFoVXN3WXJ
f2EoyDMmOn4nKGro/XEynoTGvgM7nmZLgKWVv+Ko5SQM3ToxE+kuAaG+Hx4b
2U+Cq3nHGK+WBMxUSYicfDwJd22EXu+YEYeMY1PntrydhG0Xnh/tVhIH0aN+
2e/nJ6HYHf3m+C0K5w3L/u0jTcG4Xp0SpUAUHrhOqsezT8HjxRzG6/6i0PHu
Yo+/8BT8DmJQkVQUhRB5cb4rqlNgaJLy5F3YdmA/UOqwHjAFmrsa9p5vFIHD
mn9xh1dTMLrJ8M3IYxGwsBWjDMdMAR4qKyp2VAQ+pfsmfM2bgqKXKmHWGVvh
vpRuS+wPYj5xIyGVEGGYEpuQV9sxDcWKEtQVd0H4sXH7sm/pNKgFO/BsLd8E
CZpZUt8qp+HpOHblnt8msH2C6Qs0ToNP3LDAeZ1NwLRqWJzaMw3+TxbNr45s
hH394c61LDPgamftxiCwEbxy+Wc4tWZAO65b/G4AP5wbjxPTuTwDbm/f181c
4YctkrIXIs1mYOLnwEz7Pn7IDD+fs8d5Bv5khr/5/Y0P2t0Dbp19PQNx+iY9
n4X5YK8m6/izgRnQVXubpVXJA4u+oUJtEzPQbZbUuSGSBz5+3HF62yIx3kdg
+ZADD+jJq6VlccxCyMuuzYFiPOC51cPy6/5ZEGp/pc7DyQ1tY8v9vPdmweLR
9uU/RRTw8J349ZJrDpTEBOZYnrJBa42u49HNc5A6qczscoENdnOWcQxum4Pd
kc85rwmwQUuAv4KCzBxk7x84XJTCClIvdgX/0J4D2xRnIcEGFvj2xlRNJJwY
7/9HmXMPM0j+rvlZFTsH9srCi2NLTOC6Xeb2rbQ5oKxUCSlXM8GO+PU3qGQO
okzst++0YIK7KZELpl1zEBfC5MGdzAiiea1pyeLzcDyupuOGMgM4zx5VOb9n
HkyrtXZ18DBA3aHEtkW5ecjyG7TNGCCDU7EDk8aJeWjxfVjhGESGmjIeo3HL
eYiwDii4M0qC29UavDKZ8/D+YN7bK8brWPmvDy7vDy/AsVeZs1kzy1irsd/o
VtUFMOBTki1JX8ZGe3SveJxegLKEgqcuFsvYhsFxVXWjBVALuqe92rGEWUwI
crV6LQB3ucMVrsZFjLLunDD1ZQGM9akHKc3z2HYP1c267QtwWrv/nv2reewQ
I5dfcc8CkHnd828bzGNGrIm3Hs4uQPrJ5bmDA3NYDnfTEargIpROWCYmr89i
F7fv/7br2iIE2sp7XtGcwW7ELR97arsIHTm9UplCM9hDiarCvy6LUHP2b7Dn
8DSWsutKZKHfImyu+1Wc6zWNrRzwv66WvQgaAZevT+FTWIzK0KrxyiJEv2SS
2nR+EhszjZZ+FbwEFhvmW0J9xjC/cVUmg4glqPA8WtSsMoZJOw/+EolfAuVM
3tq3/0Yxyyd7g5Pyl8C8yWP4oOMo1p9VslDQtgQJNa3P1m+MYJ3/2iqaty5D
mOmlXcdchjBXb9eoMMlleNZSLvlKaQgT5Nl2x3D/Mlym3q5yJQ1heuIWu/tp
y+AolHOw038Qaz41Gzxjtgwa89RzhSkDWH0ojzFf2jJMMISI7//Xh1lvzz/c
mrsMF7Qi4taq+zCOtEt8ESXLUH1yXF47tA9TL4uqFG1YBjvVN8EPZfuwir49
e2Qml8H01dSMo10vVnrg1JKm/ArAikKiNbkHS6/1eBHwcQXOCPauVY51YkO2
Hud31K0AU/w1uBzXiYnxe1DLmlagQFomPE6/EwszcPcZ718By/ArzOdrO7DH
4w9cT7P+g8hLLUcW835iRlxuZuxn/sH1O3lcgsntWHiu6/ZYnX8wL96fp2TZ
jjVfdP11+Mo/0PkqqTkt1Y6djr53ydL2H5j1D3QIZbVhhw/cPf056B/cbR+v
s0GtGO8F54OPWwh/pkL2WjN2Zt5pUrjrH5RpmvYtVjdjvhFOGfmD/+Dil0am
ihfN2Fqv487+hX9wsC1OekK6GRt1cBA6JrgKoaf5lGSMm7CKEDtGsvEqeORE
ZIj8/oY5/7jR5Da8Chxue+o9b33Bbn/4YtM4tQqLJodV1kS+YDbxsixiK6tw
7Ubp55n6BuzqzeXDn6lrUNG27erijgbsDNkvivfgGsyeLZlRbqnDxKVSrZLu
r8GmNnThvmcNJsJFZVzyXoNd3n1nzI/WYFtm7CJPB61B+OIv7cuL1RhX2eHG
ybg10DeiLerYVmPL56vllKrXoDk6y/yP/mfs652h9Ua+dbjzTTjZ4EolVnfl
dLiY8DqUqBW/MBCrxKpUs2SddqzDgoCyN/twBVZCcbkmeHgdfIpSvZmcKrCk
GJY68yvrkP5LoYrjxSfsfvWO0KXkdTDzNk1fHCzHNlJkQ4NFSLS/DCky15bL
MD2G85vGxEg0zWUZW9m3ZVjE4s1QtZ0k2rjN+9ofhmXY9oHU0IX9JNr0m5sZ
IR9LMekPoq+uqJBo+tZux6tDSzDMjjdMypxEi+M4X7f36jvM3WL/Fi9LEq2D
dT71usg77JPh2bBOGxLNZ8vRNsGuYkz91JOwQCcSzcG+w8TeuBjzh+QtI3dJ
NFJexJEdIsXYl0OVYaoPSDSIs7XT+FWEaYuRw+e8SbTTb04N8lsVYaFbtglo
+pNoNOH78skiRVg719HwlEASLaw66Nm15UJMiNlAgOEliXbsh/MblblCzGjl
TrhBGIl24XaItev7QixmKlSgIJKYj3lb2JVLhVjfUH44VyyJZnLcJMi2qwCT
7PouYJlIojH9GnMNu1qAWTZPhpenkmiDk5sqTi7kY2m1XIJCmSTa0aejGRLJ
+dgftCfCMZeoj3j5qrlNPnagSEPwSyGJ5umwcOmaaj5mn2EZsbOERMPTVv71
SORjBfHegh44ifaA5j1P5cvHFsMTIn5+IuLJq+6dY87HjgR9FJSrJtHs0mIa
51fysPve3REB9STazvqTsiF/8zDktiY49JVE0wo6uCenIw9jcBB+TWsh0W7S
jNRiS/Kw41ZHhCJ+kGgWrlM/tgXnYb7Geq9nfpFok8xdQXpGeVitrrPQ2V6i
XiWFp6/tyMOoZ16+Thok0X5x3cH/DuVizw9/fa3/l0QbDqJU3zTPxZr3TQjl
zZBocppjKSOiudjmHZRIyiKJNisQzrW7Jwe7LCQtbPGPRCvz4XEzT8jBInnV
I3ESmVbV2LCn1yYH62K1EBZgJtPC2sXYGxRzMNE1r0h7djLt4degag2OHCxp
FI/cwUem+RdP+nLkZWPD3b+EH24i0w4l4k/S/bKx3W0rke2CZNqW44+WHptl
Y9mfFN74i5NpWVq/6oQ3ZmMlkVlvpuTINManNw9Nd2Ri8kydb+oVybSrfRpD
Ym6ZWLY1W1QyRqaNr/OUZghnYklHrkZdUSfTjrorhm0xy8CC2vmiawzItIiH
IoCxp2OcNIhOMCHTJChxzTqFbzHvZOvoh9fItE3O+wa3m7/FXJ0rog/dItO8
VZ9YW5SnYRYbnGNiPck04fzhyA1eqViPa1yMmw8x3y/fuWa5VMyw70vMpadk
Grvx9dvOf1KwVo2VGNlgIp5ABj7DmBRMK3dnLDWUTLvYUQ0KailYvYBO7HAE
mWa24O147W8ydsLDPfZTNJlWkmXEtTctGSsfTo+NSiDTFH57iz9TTMaUz/+I
vZtK5N/Xg02kJGFFRcxxOplEfV9c1Kr2S8Jkt8nG7c8j00Z6dzoWuydh6d5G
cRzFRP60pYfGlknYzgm/uIFSMm3K6IJz2bEkLE63KA6Vk2mjp17XoA1JmHBZ
X9zrKjLtRWK3wHhnIvZKgifeuY5M83LRUxmKTMT4nirHa30l034lm+ceuZSI
BcxYxu9pIdPmVqpcA7gSMTaDkHjWn2RaSpZZhk1lArYqNRFf1kem2VW5XPsj
n4ANf45OkFog09gya/N3KMVjZvvrE5j+kWld51iY8xnisa7QxYTfJAaaHvfe
oG91cVjTNa3EEA4GWv1G6vYNV+OwD4yMSWQRBtpnWbu/L5/FYorW+5M6xRho
1gOHY2mXY7H87wZJRTsZaKNff3TKScZiqbH5STdlGGgbDm9N/InHYBLsPUnq
8gw0je2dAczPYrAoO2qyuBIDbQ/puV6eYQz2EiySf6gx0D582Z7XNxONcSc/
T84/Rdx/fJPtzZfRmB8Xnhx4joFWeFn3wkfZaIzJeSz5hjYDbTE8ly3oUxT2
8NfmlON6DLTuu71tzPZR2LKaWsr2Kwy0NdmytEG2KMwx3S5l5Srhpl9HvcXf
YH/536S0WjDQdrAFON1riMRuuNak5Fgz0A70FkbsXHyN/Q9ZeUOI
          "]]}, 
        Annotation[#, "Charting`Private`Tag$138464#1"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          AbsoluteThickness[1.6], 
          RGBColor[0, 0, 1], 
          Dashing[{Small, Small}]], 
         LineBox[CompressedData["
1:eJxF1nk41NsbAPDZlT1LJElKthAKFXOmULjZGSK7pEhKqCxZUihaLGVtsXep
bEmW71FKtoiSlK1QJNkHY8zvzL3P87t/zfN5nu85c877vuc9Z4vbactjBBwO
54rH4Ti/wAuUlngLwV+EyfMBwbagM1P8MdPwP3t0TRcabP/PsdScB58H18Ex
yr/uWfO62cr4P0faH7vSWCv4/+8DH1sbfTsvCCdWf//jEwR9XpaGIJxc/tej
Nlok0y3/2aNQcSWL/z/X3ZfpiegU+P/4v3Vq76zV/c/aUdP7eLfz///7y2mN
jQMJfHB2/F+/L8m0KmPwQnzjv3bYXFpo6M4LjdkT/3jk+ht23zseGJ0z9o+b
Iux/+ujzwMfm/f+4OHCyg1nFDdM8n//jO87U2a9bueELBeV/HGGYIFoXvxbO
Z1sGcvz6uWPzXfxaeNzgqx/HaxWUL/kHrIGrP0JdODbmLa+r37kGqud6HOb4
2p+9qwITXJBPJ3QHx62d9bpO+Vzw5s8A5jlkvmeGoUVuXDChbLmSY6ffrtbc
5lywf8IkiePH24KVvHS5YKDY3WMcsx2SCG8UuSDtzEE1js0Tiz9vFeeCBrxP
l/2RHzS/eRpB5oKHWi2rOZ7BD14dmKHAicyc4xynZfr4J4xR4IdgpjrHenuW
nHQHKfBGWcLyWeQkPyHN9DYK3LV7MIZjHd57W4wbKNBp+rYlxyP5SnxLLyhQ
eXu9GMeaA/rD9AIKBIKv75xB7r/4vp18jwL7XZNsOL663qm6PJkCH82PiHC8
s3Q8z/06Bfb6Rnf5cerLJOi2UBQF3l+ai+e4KF13zckwND4lyYvj8DFiWP0F
CtTbOk/jWCH65gk/P2Qun7HTyKud9IFGbwosy7pbxXGn9CabzccpMFCw8CrH
F2sKae8cKXDKW2czx5JKXydjtSlwXy3r3ClkpQBNNbVdaP/ir/Q43ovd9O9R
pcA4Bzkhju2sDBa3y1Hg33fli3048Ql+jG8QpcA1fOxP3sjZr9foe6+jwE2H
Xt7juFTA/YoQHwU6Wmh4ctyRLcbjSqJA3d6+yZOcfLeGi6zOkKHa5tLfJzjr
Wf+FnjtJhpK+OoUc73DZnXp4nAwV9QjuHBvPjW1KHyJDk1Tyey/kK5JWctod
ZMg8fTz5OHKyZ/GJgRYyTMgy2s9xzlOuoiuNZOhzjD3uifxSv2bnxzoyJP2p
0eJ41Ud279liMlzoFqnw4Kyn8lLIhkIy5PmgZMqxJL63DuaQYRMxe9idE4/k
BD2BDDK8sRi0luPAWsbhojgyvH7KVM0VOZrL8obVFTIc3f2nzIUTH4ui98sR
ZCgdflGD49IRF7rRBTIMCGDKOyNP8jU7jx4nw8ap+F9HkVm22x7Gu5PhSR1X
K455H4YN73Imw0jh+EoHZEVNjRNRdDKU0bl+zh7Z0ynjrLQBGa71iKyzRQ4o
WChvpJGh1WAYhePLM+YMXx0yHEuWNqYjP7xCDqnVIMNcP9Yra+SBx77R9jLo
/8IXwiyQE3pSjuRsIsPxPu8Ec2QqAVOeFCdD/8m2O2bIGTYC3RECZPjLb/yO
CaceVp5sL1whQb+jrkZGyFzbe5ZnGCRYnOewzRD5mRmuXWeWBE2thBcPIotm
mwd1jJHgudKdsfqc+jSabmR0k+DLXWKmgNNf/Ddk7O8kwYMicFyXcz4y9/td
ayPBZIuWMB3kG1O3xDc3kGA4we76HuTDd9RPHCwhwcyD7doayExor3uziAQ3
1dRcU0N+NB65rjefBH88iv6gyulP1K6qU1kkeJ/WfnAH8pthf+7kayR40YFV
sI0TP/6M/v4rJJg64Jgrg7xNu6FUPpIESx3YydLIUddEHGrOk+Bdk3ZLSWSa
esXf34+R4DazTRrCyH8c+i4pu5KgdmR/kyByVjTZOugoCVbwfrPk5+S3x2aF
24oEpSor965FfnFp4bA6jQTLznjbrV5E90Wh1JaQfSTIzea+wEQW7zo4/1qT
BLO+349dRA7afifziDIJyr29cWEGeVeb1mS4BAnG1pfGjSAXS1y40T5HhJei
3iq9RebXmj2jOkWE4aWC2g3IflanrG/8IkL3A52qEFn9uusGs29EWF8bP1iJ
/IxllP2unQjL0vn25SGLSTREq7QQ4aRx6t0HyOc1qV4Jb4iwKPFBfwbyvtPq
yqa1RNiVNiuTiAwHJSrbHhFh1HG/8DDkLaykVOU8Ijzb2Lr+AnLkBoGQ+AdE
GNd85LY/sr4lcb/JXSJMFMjV8kJuevWruTWaCL+t5qmaI3fm1fS3uBDhu423
Tkgga7zSrFc6SoS3hsQ2CyMnDTzNvmZLhOlhx0p5kOniuV5/mRJhh8NOv5UL
tuBLbPxM8z4i/GxStu8rsk4e10dFLSKkiVI6u5CzXkZUxqkTocurzQdakF2Y
ASHGCkS4ll+huAp52MeJ0ryeCC/DrHtJyAaxPT8VhNB+u8Lm4pDzci1bYvmI
8JNjilQE8on+gzeNSERYsX70jw/ypJmqRNM0AV6lrwYdQGao41TethKgRXtX
8vB5W3D7L/6+oLcEuIrb1tmDrOIheV2ugQDFtuUNtiJ7JmuPX6kmwCit3jPl
yN0MvzyDRwRI42qlRyA/qxmSenWVAPVeNpSsQ7b8+KftbBQBjlqGthOQJ3+z
QmQuESA5tOvZbJAt2L5Z4ktEIAHy3tEa+ICcEmF5h3aMAF0/MSKTkAMPvhKo
O0CAi+FDyWuQhZzf152iEuDWTzdX5wPR/Rs0cGrTXgI8AM9v/448WsBsCVEj
QKkTuNIaZDrPrph90gQot7Mm3wdZsyMHV8XCQ2GT475YAMrnz9InXkt4OGFe
FVaA7IuvdxKfx8NU+/t6t5Dz1PtqgibwMOPSj5+uyOuTRS9ofsHDWaGMm+xz
tmD+yJXpsud4aIqbyFdCVrd0U6eU4+Hzh0MmfMi+xlR/uyd4qCrQHD/pj9a3
d2FuJRcP54vLep4gd2/0XDRIxENr8ZorysjP+wzYn07hYXyM/RGRs7Zg7uMW
muIJPPz16ofGnzOo/7xjhYd4oPmS1qc2IRfUVRC2OOAhtJ5khyCn3pOlnDTE
Q+XIQ0r9fuh+dSXzrcjgoUpdGG/saZSfI0MmplJ4aKcbsPUo8oxFbcL9DWg/
Og8rVJBPHggQNBDEQ+f+nEOdvui9tnVEOIGFg1MlZjzCyDojDRLSPTiYPVK1
46IPeu94XZbXj8dBNz6rhU9eKF+Oc7tHY3BwaT+xIwN5wsJDL+YyDr46cFvC
Fbl3r75TazAOVj90wf08jtbDS0qyPomD90vTGZOeKF5Po9geh3AQxyjpG/RA
75mcWV4uPRz8E0T3yUA+muouUUjFwaapOSdbZONIvd2/d+Ng6JM8nxZ3VE/W
RO+AbTi4UHJPvMgNvXcYkd3ReDR+58Z1Zi6o//aMNJotszE746JYHPLd+tff
EqfZmM31bPsSZ1twLjFaXHKQjZVWPyeuQ96hRb68o5aNVZV3erx1RPdNGMnB
JIiNVd4suLvewRYE8xHX3vi1iqmtU3RUskHzV4QdEv+2it1WbcXVWaP3k+Ny
9IOeVUyZevmVGfKux9OE8ter2H5/6Yd+Vug8mA6u9NxbxXanU0vzLdB5v1k7
JWO9itk8anAYN0H9SPh8z7MaFvZbt+BBjYEtOFs9u55WxsIuNjllqiHfdz9t
01TIwiSX32fk6qN+XubZ+SWFhbm98c6J00PxtbZpxZ1hYUPMIsxwP3qfpahD
Y1kW9irjYW+SDurHG37n98evYN97N6YcV0fz5fESky6vYM0pSTN1aqied+1w
MgpewXjuPD4kimxr5i1S7rWC2ffFfalTRd9Hj4XH6q1gXaJV8WRlW6A2M3Jk
1xITa2CO33aRswXWLf08cR5MLLyL291oI6p/u1VP4MDEzphrloVJoP49uunl
nAUT23xLb6x0A4ov3vG8M2Bi7KKrJHFxlH+tL8O7JJhYhH3DwgcRVF85n2oH
2pext4ETXhv4Ub2HdpzevW8ZO+vfFvCWTQf3ieIecerLGJFvHTa2SgcFsc52
AwrLmBEljcKN/Dx5khYrvoy1HuO9fmiFDnqKeYX65pewF3dzHZ8t0sH6fsPy
y0+XsJ4vEZ/cpungNrV+sUt2CYto4nk8P0QHaQ1rfstvWsIqQodX5gfp4KGx
+VCoyBK2WHSFOj9AByX0gSY54hIm2fk9daqPDjp8WWnBg4tY4M+U6t7PdCBw
b4/u1vRFbEbKLiLkPR1cZ5dEnBVcxKSm1+nMYHTwQT85lZdrEXs29ffvyjo6
kIw7X5LHYmDMO762IbV0UCQChnrHGdh8+FZzfDUdtCm00vReMzCjM1c7mRV0
wG89yha+yMB2XmSPlv5NB7cKJMIqvi9gy5694nIpdJBiEXVx8dk8lsrfsCh8
nA6sbyQoqT+ex44Y739w7xgdCLWlfvXOncdUfY5LKHrQQbzhU+pA4jw2IekQ
vs+VDi6DPsIbv3lsdtljvbkDHZxV1rqWqDiP/XabE6OZ0YHZ2l/pKllz2HTV
ky4tTTpYU29Z6xE9ix0Q+84MxNGB4jD+WnPILJZ5wcGwdtUG/MVVYrfz3Cwm
wCC4EFg2IMFEYI7pNovJTLdOxSzZANHeFsVE2iwWeVs99+qMDZCZ1b9bz5zB
juY++Wo8bAN0ZLX8pc7OYKnPfC2G39iAM7EbFXqcpjGXn/ZW9jE24JteAUNe
eBI75jbwaHjSGpQ9EVXzfT6GhbzQFyk0tAYl1W7aBrgRzKxKRQt/0QpkXXh4
S9FpEFN5IfdUNtESVA7XsJ6L9WBcA8JOug8tgLjyzx/hfm3YEdtvp3QizAFl
McHw8lAlpqQ/UaleYQqsKvY48X3JAePc7i3PjA4DxseRdPyHRhChFqEaEGUE
ihS2v6WkfgCZO254P/lwEODl07aPEfrBVFLkQtE3PeDQE6eS6vkdzPXRFnKV
9oM34x+DnlT8AB1c9V1T2rqA5Lh/VsNtAkxuNnMo4tEGe4ztfgc4T4FwMFj9
tlgDfBaVXo6pmQHb9twUHBBSBTlmTxTev54By7tKE8bmVIBvLPXIhvYZ8PSE
h9KaTyqAxDr6vHBoBixqvHgTka4CVIZTA5sps2CTvUeMp6wKiCoVnuWxmAWG
/M+lFnSUgbIZ10TCyCy4sFhYpBSsBCJiJvuS+OeB23CQS4SGPOhusjmnKzYP
gvDpPzQl5YESTy336OZ5EMlwDGWR5MHH+GtaWmrz4GTI1obabjmgkCh/67PV
PCiUOnaAO1gOvM9005dKnQekrOcx795sB1vKuh/lb10A4UbBqt9OyoL6vrqg
F9oMIBzS4cw9LAO6nePGNx1gAJeXQid122XA+JCNY8RfDJAk0tsf+UIGiIxO
HDB0YgC9UKdB/VsywHNSgr87igGsitsCcoAM4GUH5ky/YwCJ29Yub+9tAXRp
1ffyxxbBDpXTPO6npMEvt3uKd24tAa/OpodXgyRB3MQBkkPaEpjwXm1/YScJ
FANH+6Syl4BeDPjIvVcSeMUq38orXwJnXnR8Ya9sBMNPqhkVn5bAO/mD0xqX
N4KvK58aPmxaBpGyPBY5dyRAa4qgs9CjZWAuvhMGdYoDb+ly7e7SZTAon66I
rxQH3I9shdKql0ErQeRldro4MKzNer2lbRkc5aKX7DomDhq+79ihNrUMvgof
fhOzJAZqdhotmWkyQbg6TnR6uxgoao5IjH/JBOXhTZ5jiaLgh2+EuWwLE7w/
AfhWg0WBjHAEX20XE+DWyyoreIiCuw7hVyeGmUBSdYPty12iIHoiLPgvrhWQ
envg4Gq3CHDiD3Ffe3gFVHns9UzYLALWWQZqRH9cAT3Ssq+fQyFweCFgSrJ/
BZw+pK4Z/LcQiEkLKC4fXQGsZDkFyxQhsPrtnNwwYwUYaKm/0fERAuP+/hv1
JFjgpO5GAfIGIdCQ7EfEO7NAV3dTh0ngOhD4+WRXyE8WULegp2oDQXCm7p1P
+zQL7M0NXK1VFAQ+2eoUGSYL7KsZZrutFwSup5a1G/lWgWWMvcTarwLgMD4u
a53GKpBW83y1oiMAtioUnsgLXQXUmXs8cWR+0HH+B7tdiA00rrGEiip4QIvj
X6kykmywYag3cfoWD3hz4Il6gCwbjDUZuDv68oBq3qBjEtpsUPxDLv6RHA/I
u09p8XBkg6ZCXSmFDG4Q+lY2ZSmfDUytRnlNr68ForzqKbekcLTsj+w+n6tc
wI5gvv6XDI5ml0drcTjJBdIWT6Xoy+FoVQPNOl4mXEB6pDCFoYqjmfFbFH0W
5gKKdVvuOO7H0ezFN41LPaQAqt+6uwoeONq58JRpsQYy8PowlVpfiKPJsks3
sqRI4FEzv8TGxzgaQerv621kEvgNd6SdK8XR9gUsHXWfIIKzxV5pctU4Wtya
qj/yVUQQemUwLb4VR1to0ylOtSSC29od6Uf+4GhNylzZJjcIoDrjSeb0Ljzt
q2VNo64kHmiSvma27sHTpPIElXWIePDUe01WPhVP25+TL3dkHAfy9rpmORri
ab1reRNGnuPAzR6he00OeNqGFw6j+2xxwFMk8P6DSDxtUpE05KzOpgpd18m2
6MDTNrnlXmx1XaHGz3pl7/iIpwU+MO8r0VyhrnFIzubqxdNGDtl/ecezQmUp
TGbXfsfT2sr6XqY+Y1J/Nt7LUWDgaeInfXlk+JnUOiIxDy9FoLlYqp8oeb1E
3eOtmvdVhkCTjyoRi8xcopZ3OuRVyhFo4Xb2p6MClqiFD8rzTqkRaFwf5a4I
b1+iJgHP/M/6BBqr/uBH99hFqkD+7fxyIwKNmGn6odBtkRrHj+XfMCXQAvbP
/bVNZ5F6qU+swMCOQMuQlufL/sOgLuvrF0g7EmhFYhFe1U0M6rkivwKmK4Hm
Hxr8EpfLoP4Rzizo9iTQYmQO/jG7xKCeDG4qKPEm0IRVdSfe2zOo/wMstmkw

          "]]}, 
        Annotation[#, "Charting`Private`Tag$138464#2"]& ]}}, {}, {}}, {Ticks -> {
       Charting`ScaledTicks[{Log, Exp}], 
       Charting`ScaledTicks[{Log, Exp}]}, GridLines -> {None, None}, 
     FrameTicks -> {{
        Charting`ScaledTicks[{Log, Exp}], 
        Charting`ScaledFrameTicks[{Log, Exp}]}, {
        Charting`ScaledTicks[{Log, Exp}], 
        Charting`ScaledFrameTicks[{Log, Exp}]}}, DisplayFunction -> Identity, 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None}, DisplayFunction -> 
     Identity, DisplayFunction -> Identity, Ticks -> {Quiet[
        Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , Quiet[
        Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& }, 
     AxesOrigin -> {-4.605170185988091, 0.}, FrameTicks -> {{Quiet[
         Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , 
        Charting`ScaledFrameTicks[{Log, Exp}]}, {Quiet[
         Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , 
        Charting`ScaledFrameTicks[{Log, Exp}]}}, GridLines -> {None, None}, 
     DisplayFunction -> Identity, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, PlotRangeClipping -> True, ImagePadding -> All, 
     DisplayFunction -> Identity, 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None, 
       "ClippingRange" -> {{{-4.605170045013494, 2.3025849520194486`}, {
        4.167758571277503*^-7, 0.33359051671308204`}}, {{-4.605170045013494, 
        2.3025849520194486`}, {4.167758571277503*^-7, 
        0.33359051671308204`}}}}, DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {
       FormBox["\"T\[Gamma](MeV)\"", TraditionalForm], None}, 
     AxesOrigin -> {0, 0}, CoordinatesToolOptions -> {"DisplayFunction" -> ({
         Exp[
          Part[#, 1]], 
         Exp[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         Exp[
          Part[#, 1]], 
         Exp[
          Part[#, 2]]}& )}, DisplayFunction :> Identity, 
     Frame -> {{False, False}, {False, False}}, 
     FrameLabel -> {{None, None}, {None, None}}, 
     FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, 
     GridLines -> {None, None}, GridLinesStyle -> Directive[
       GrayLevel[0.5, 0.4]], 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> 
     NCache[{{-4.605170185988091, 
         Log[10]}, {0., 0.33359051671308204`}}, {{-4.605170185988091, 
       2.302585092994046}, {0., 0.33359051671308204`}}], PlotRangeClipping -> 
     True, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.02]}}, Ticks -> {Automatic, Automatic}}],FormBox[
    FormBox[
     TemplateBox[{
      "\"\\!\\(\\*FractionBox[\\(T\[Gamma]\\), \\(T\[Nu]\\)]\\)\"", 
       "\"\\!\\(\\*FractionBox[\\(T\[Nu]\\), \\(T\[Alpha]\\)]\\)\""}, 
      "LineLegend", DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0]], {
                    LineBox[{{0, 10}, {40, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {40, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{Small, Small}]], {
                    LineBox[{{0, 10}, {40, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{Small, Small}]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {40, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[1, 0, 0], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> RGBColor[0.6666666666666666, 0., 0.], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"1", ",", "0", ",", "0"}], "]"}], NumberMarks -> 
                    False]], Appearance -> None, BaseStyle -> {}, 
                    BaselinePosition -> Baseline, DefaultBaseStyle -> {}, 
                    ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[1, 0, 0]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[1, 0, 0], Editable -> False, Selectable -> 
                    False]}], "]"}], ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0, 0, 1], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> RGBColor[0., 0., 0.6666666666666666], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0", ",", "0", ",", "1"}], "]"}], NumberMarks -> 
                    False]], Appearance -> None, BaseStyle -> {}, 
                    BaselinePosition -> Baseline, DefaultBaseStyle -> {}, 
                    ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0, 0, 1]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0, 0, 1], Editable -> False, Selectable -> 
                    False], ",", 
                   RowBox[{"Dashing", "[", 
                    RowBox[{"{", 
                    RowBox[{"Small", ",", "Small"}], "}"}], "]"}]}], "]"}]}], 
             "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{
  3.7842842820875063`*^9, {3.784284350636279*^9, 3.784284371386569*^9}, 
   3.784284536464419*^9, 3.784287709460148*^9, 3.784287773013843*^9, 
   3.784287831315805*^9, {3.784287897289638*^9, 3.7842879162470284`*^9}, 
   3.784287958085731*^9, 3.7842895215752287`*^9, 3.784380497576807*^9, 
   3.784380586797553*^9, 3.784380838162087*^9, {3.7843808691763687`*^9, 
   3.784380890925535*^9}, 3.784380976156371*^9, 3.784381016557641*^9, 
   3.784381124065312*^9, 3.784381172580337*^9, 3.784381362865377*^9, 
   3.784382405463233*^9, 3.784382492333268*^9, 3.784382561975569*^9, 
   3.786691328732959*^9, 3.7875841533804913`*^9, 3.796733791558237*^9, 
   3.799555351977353*^9},
 CellLabel->"Out[39]=",ExpressionUUID->"22183634-7a38-4e4f-bdfb-5a0f0091f41d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Run", " ", "case", " ", "with", " ", 
     SubscriptBox["\[CapitalDelta]N", "eff"]}], " ", "=", " ", "0.5"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"resvec", "=", 
    RowBox[{"Quiet", "[", 
     RowBox[{"FuncDRnudec", "[", "0.5", "]"}], "]"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.784277772039658*^9, 3.784277775263968*^9}, {
  3.784278221518071*^9, 3.784278240498653*^9}, {3.784278416059729*^9, 
  3.7842784189191713`*^9}, {3.784278836699937*^9, 3.7842788719260607`*^9}, {
  3.784279190280746*^9, 3.7842791982520227`*^9}, {3.7842792624573383`*^9, 
  3.784279276919652*^9}, {3.78428208196564*^9, 3.7842821052513123`*^9}, {
  3.784380379117023*^9, 3.784380380674103*^9}, {3.784382391724833*^9, 
  3.784382392715899*^9}, {3.787584196592946*^9, 3.7875841993592453`*^9}},
 CellLabel->"In[40]:=",ExpressionUUID->"da9b7e77-b9ff-46f3-b07d-ff45106f0d03"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"tCPU = \"\>", "\[InvisibleSpace]", "16.021849999999997`"}],
  SequenceForm["tCPU = ", 16.021849999999997`],
  Editable->False]], "Print",
 CellChangeTimes->{3.784382512258757*^9, 3.784382578541851*^9, 
  3.78669134494081*^9, 3.7875841743007507`*^9, 3.796733812684442*^9, 
  3.799555368138257*^9},
 CellLabel->
  "During evaluation of \
In[40]:=",ExpressionUUID->"588fe484-f155-4b22-96e1-9ab66d686dfb"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "T\[Gamma]oft", ",", "T\[Nu]oft", ",", "T\[Alpha]oft", ",", "tofT\[Gamma]",
     ",", "T\[Gamma]dT\[Nu]ofT\[Gamma]"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Gamma]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Nu]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Alpha]oft", "[", "tt_", "]"}], ":=", 
   RowBox[{
    RowBox[{"resvec", "[", 
     RowBox[{"[", "3", "]"}], "]"}], "/.", 
    RowBox[{"{", 
     RowBox[{"t", "\[Rule]", " ", "tt"}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"tofT\[Gamma]", "[", "Tgam_", "]"}], ":=", " ", 
   RowBox[{
    RowBox[{"InverseFunction", "[", "T\[Gamma]oft", "]"}], "[", "Tgam", 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T\[Gamma]dT\[Nu]ofT\[Gamma]", "[", "Tgam_", "]"}], ":=", 
   RowBox[{"Tgam", "/", 
    RowBox[{"T\[Nu]oft", "[", 
     RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"T\[Nu]dT\[Alpha]ofT\[Gamma]", "[", "Tgam_", "]"}], ":=", 
    RowBox[{
     RowBox[{"T\[Nu]oft", "[", 
      RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}], "/", 
     RowBox[{"T\[Alpha]oft", "[", 
      RowBox[{"tofT\[Gamma]", "[", "Tgam", "]"}], "]"}]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Neffval", "=", 
   RowBox[{"Nefffunc", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"gstarval", "=", 
   RowBox[{"gstar", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"gstarsval", "=", 
   RowBox[{"gstars", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"omeganuval", "=", 
   RowBox[{"OmeganuMdiv", "[", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}], ",", "0", ",", "0", ",", 
     RowBox[{"T\[Alpha]oft", "[", "tmax", "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Tratio", "=", 
    RowBox[{
     RowBox[{"T\[Gamma]oft", "[", "tmax", "]"}], "/", 
     RowBox[{"T\[Nu]oft", "[", "tmax", "]"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(N\), \(eff\)]\)     = \>\"", ",", " ", 
    "Neffval"}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(g\), \(\(*\)\(s\)\)]\)      = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"gstarsval", ",", "4"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{"\"\<\!\(\*SubscriptBox[\(g\), \(*\)]\)      = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"gstarval", ",", "4"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\[Sum]\!\(\*SubscriptBox[\(m\), \(\[Nu]\)]\)/\!\(\*SubscriptBox[\(\
\[CapitalOmega]\), \(\[Nu]\)]\)  = \>\"", ",", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"omeganuval", ",", "4"}], "]"}], " ", ",", " ", "\"\< eV\>\""}], 
   "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
   "\"\<\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]e\)]\)   = \>\"", ",", " ", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"Tratio", ",", "6"}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Print", "[", 
    RowBox[{
    "\"\<\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]\[Mu]\)]\)   = \>\"", ",", " ", 
     RowBox[{"SetPrecision", "[", 
      RowBox[{"Tratio", ",", "6"}], "]"}]}], "]"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"LogLogPlot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"T\[Gamma]dT\[Nu]ofT\[Gamma]", "[", "Tgam", "]"}], ",", 
     RowBox[{"T\[Nu]dT\[Alpha]ofT\[Gamma]", "[", "Tgam", "]"}]}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"Tgam", ",", "10", ",", "0.01"}], "}"}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
     "\"\<\!\(\*FractionBox[\(T\[Gamma]\), \(T\[Nu]\)]\)\>\"", ",", 
      "\"\<\!\(\*FractionBox[\(T\[Nu]\), \(T\[Alpha]\)]\)\>\""}], "}"}]}], 
   ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", "\"\<T\[Gamma](MeV)\>\"", "}"}]}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Red", ",", 
      RowBox[{"{", 
       RowBox[{"Blue", ",", "Dashed"}], "}"}]}], "}"}]}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.784282164027906*^9, 3.784282279509446*^9}, {
   3.784282320628632*^9, 3.784282345259815*^9}, {3.78428410458879*^9, 
   3.78428410487086*^9}, {3.7842842030978107`*^9, 3.7842842729937563`*^9}, {
   3.784284336300796*^9, 3.784284361155418*^9}, {3.78428451854011*^9, 
   3.7842845192634478`*^9}, 3.7842868102056828`*^9, {3.784287687102829*^9, 
   3.784287687995302*^9}, {3.78428933398415*^9, 3.784289341230422*^9}, {
   3.7843803849139967`*^9, 3.784380460841353*^9}, {3.784380888328226*^9, 
   3.784380888424686*^9}},
 CellLabel->"In[41]:=",ExpressionUUID->"2473a44e-6671-4310-855b-bde5549eb92c"],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(N\\), \\(eff\\)]\\)     = \"\>", 
   "\[InvisibleSpace]", "3.548844619679207`"}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(N\), \(eff\)]\)     = ", 3.548844619679207],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.799555368261116*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"c6ab6768-bdc3-44db-a55d-f31debb9122f"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(g\\), \\(\\(*\\)\\(s\\)\\)]\\)      = \
\"\>", "\[InvisibleSpace]", "4.3111459158600835551`4."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(g\), \(\(*\)\(s\)\)]\)      = ", 
   4.3111459158600835551`4.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.799555368270293*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"01c1257e-c5db-45e7-9baf-7fe50b650c15"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(g\\), \\(*\\)]\\)      = \"\>", 
   "\[InvisibleSpace]", "3.6119371647366316225`4."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(g\), \(*\)]\)      = ", 3.6119371647366316225`4.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.7995553682781067`*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"6c42f006-0226-4cc0-a20d-086796ddc4c1"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\[Sum]\\!\\(\\*SubscriptBox[\\(m\\), \
\\(\[Nu]\\)]\\)/\\!\\(\\*SubscriptBox[\\(\[CapitalOmega]\\), \\(\[Nu]\\)]\\)  \
= \"\>", "\[InvisibleSpace]", "93.0706240224313745557`4.", 
   "\[InvisibleSpace]", "\<\" eV\"\>"}],
  SequenceForm[
  "\[Sum]\!\(\*SubscriptBox[\(m\), \(\[Nu]\)]\)/\!\(\*SubscriptBox[\(\
\[CapitalOmega]\), \(\[Nu]\)]\)  = ", 93.0706240224313745557`4., " eV"],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.799555368281281*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"f69e3031-bb5b-4ede-a901-7017e78bc47c"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(T\\), \
\\(\[Gamma]\\)]\\)/\\!\\(\\*SubscriptBox[\\(T\\), \\(\[Nu]e\\)]\\)   = \"\>", 
   "\[InvisibleSpace]", "1.395882195251403779`6."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]e\)]\)   = ", 1.395882195251403779`6.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.799555368284388*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"b591cad5-5e0f-4d5c-9b7c-782562183a32"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"\\!\\(\\*SubscriptBox[\\(T\\), \
\\(\[Gamma]\\)]\\)/\\!\\(\\*SubscriptBox[\\(T\\), \\(\[Nu]\[Mu]\\)]\\)   = \"\
\>", "\[InvisibleSpace]", "1.395882195251403779`6."}],
  SequenceForm[
  "\!\(\*SubscriptBox[\(T\), \(\[Gamma]\)]\)/\!\(\*SubscriptBox[\(T\), \
\(\[Nu]\[Mu]\)]\)   = ", 1.395882195251403779`6.],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.784284273653009*^9, {3.784284341962997*^9, 3.784284363435701*^9}, 
   3.7842845199678593`*^9, 3.784287705357251*^9, 3.784287763946207*^9, 
   3.784287820240755*^9, {3.784287888782918*^9, 3.784287913572201*^9}, 
   3.784287949447773*^9, 3.7842895045322247`*^9, 3.784380494856699*^9, 
   3.784380584312818*^9, 3.784380836839382*^9, {3.784380867889707*^9, 
   3.784380888807678*^9}, 3.784380974219545*^9, 3.7843810120862713`*^9, 
   3.784381120920039*^9, 3.7843811673559217`*^9, 3.784381357818151*^9, 
   3.78438242582051*^9, 3.784382512373486*^9, 3.784382578670459*^9, 
   3.786691345047328*^9, 3.787584174448915*^9, 3.7967338128349524`*^9, 
   3.799555368287431*^9},
 CellLabel->
  "During evaluation of \
In[41]:=",ExpressionUUID->"3a902286-aa49-454b-87bc-2a858a4c151f"]
}, Open  ]],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          AbsoluteThickness[1.6], 
          RGBColor[1, 0, 0]], 
         LineBox[CompressedData["
1:eJxF13c8VX/8OPB77XGNNM0ie6dBqfO6ShIqMhp2JNIgImmQbKFBSPaO7FHG
eRvJKFnXSNlkZyQr/O79fL+P3/efex/Px7nnfV6v13mf1+tc4Su3z1+lIxAI
OPWD9g22kJdrz4PK6ch0HhcpwEIvm2N0+f/c+nbX+zXN/7N121z6SXEehP7X
/lhSfHf/lv9/PCfxQFF42Bbk/r/uYvnUoK+1Ba0R/sdZ+8LSOBi2oI+z8J+f
XL7q87mcGwUY/49d3hucHrzHjRwHVf+zHZ06aX0/N9LtFP3PFXEiXZ6tXAjV
hRBpfne0/DXrMS7U5G2+8ZjqB59OTEnxcSF/XZ9/NJ8920DWWuJE2KG+FZqL
X3DYfx/hROcv7VmmeU+nbtj1dk7ExrblL83+/K/w1SpOtE8vYoHmefPO8YBc
TtQuHTtLs3ES3zb+OE5UqCAwTXPNmCn2LpgT+QvRj9O8HPBYf8GbEz0dyxml
WVYu3lb1IScqqZIapjnMcfhFoz0n2hV9s5fmhq1MaVutOJGIsWIPzRuFEuXG
lzmR4nBGJ81KF0+3JupxokTGzjaabVav/5rU5ESjO7K/0ezRuDYXos6JTi+7
faU5Kjro334yJ3owK9dI8zfs/ZYHytT1WM7W0qziNadKEudEdF80S//Lb79u
13lhTjQEq8U0Zw2/d44U4EQfZH0KaBbRuJElvpUTCeRKZtH842/D6ZucnMiT
61g6zeGpUqP5rJzot45EMs3srL+E1AicSPegYzTNtR9OlvqucaCqf42v/4v3
etKFpr8cSMhz/gXNfxotQ42nOZD/ZJQfzb3BP4gu3znQ+UoJR5ojQDWmjMKB
TMqyrtOsPxt5hL6FA3UYTFvRXKd3wSn0MwdSNQgwpPkJXRFXZxUH+hM/cIbm
o/nbMgUrOJBP6shJmnO3tQ5nFHCgfeTxAzTf+KToOZfNgW79HJKhWdwlRFDl
HQcqVPcU+a+enTpGn+I50OD0G06an0Z9/twXzIGauj8MPKIa0xa3Fg/gQJ7h
uhSaV9aeEm56U/P7YV1H8y3T44dX3TlQa/jfTJolOeMpZFcONNp5I4bmwYpN
R987HOjFofMhNF/YU56xzY4DFa2Z3qaZp4X/lLE1B3oZ7GZG8xfP+0Px5hzI
cGBEh2bykLKAghEHmnnwUoJm6ZS8Z6dPcCB7mdC2h1QPG22RDsU4kMSR+2U0
xzI71HYc5kBejU+TaN5mJ7dppciBFrFxR5qbeIPeZMhwoICcIxdo9m+YUJ4T
50DD8q9Uad6QSXN4JMiBKBG8RJpLfjBxfNrFgc4UyQ8+oNrp2dV09m0cKEd+
ZxXNYzMigxFsHOjh6yOPaG7Jfaufv0RCfSZ1o+60+hme3eBaIKE7h8UraM5Y
3ki7MUNCrN9YX9H8BCw2xEZI6DyfrCrNM4Pc6U/6SUjLUY9Es7FPpX5fDwlV
fajvuU/1ga8i6RGtJJR42PcuzXEObfqLX0nouUcR0Eza9nRDr56E7G3oWWge
uTyiz45IqPpX0gs32n4eS914mEVCZ61Yg+9RTRd0Mf17GgkRmXafovmWAquB
chIJ4ct1665Un3K5nj4bRUI1SlFWNBfw8hucCSch/fAgHpr3lDdupD8noTyl
ogoXWn+glzOw8iOhZwVanDRbp/7cwL1IaKuHT+5dWn20gtMFHpNQx7E3ejRn
PP+9QblLQubSG77OtHx356VrWpGQQblR4x3a/q+6YpBsRkJ8is3naT5gs3WT
eJmEai+f6HCk5ZvlbFCqS0Ifp8KbHai+pyu2uVOHhI46qJ6ieWSBku58ioRG
Giw/3Kb11yPKm3IYCbHs/Bpyi7afen+lBxwmoXnxwYWbtPp4Rhj8OkBCRbbc
+jTfrltJj5Mhobt6tcQbtHyNygy28pLQmdOm/HZUB63c2Ly9jYR8xdLO29Ly
jxbM+MJFQj+j33hdo+U79GjTm4mEGKIWWq7S9oOvQsYgkYQeRe/6bU3LX7rf
ANbZkfNpc0aanziqZSwvsKOw2WX+K7T9sG3ewPA3OzJQxfktafUpSdjMnWBH
f1mfbLWg1WOT3tB+gB15bvBMmNLyDard/NnEjiYPJ8pdosUny+4a0MCO1vOy
1y/Qfv/l3MyhWnYk5Hz7kxHteuzdP0LK2ZGR/SFVA6qFAiY/qGWyI/2QgaKz
VH+SUtw3k8qOLCfnjp6h9ZN657SoRHY03vC6RJvqjywb4QtR7KhLqz5Ek/Z8
+3I7J/tTXb6tWI32PIkbTup5s6PU772DQHVKbeSVDQ92JMv9j4jR+iPjXr0L
99gR/8n1XUeoDn16UJ71GjuyyPBJ2Ee1suj95MIr7Ohk07qFAq1/VlcIXDFj
R7WfT3HL0eYPvSZ7qSE7av78FiSprve8PHZDnR3FkV07BKl2EI415yWzo4LZ
J9z8VO+qHOr4pMqOiHSuKrto84Zw85PQfnYUoB2ny0N7Hh4/TmgWZkcLY6kv
GKlOF6rhfShIjVfMS56Oat0KludSvOzU9wHLgo0LFIhZf/74CTc1/mK3B0tU
qz5IMjmwyYb8Rhrmf1E9yD/e1r/Kht4vWrUOUR1QKqf97C8bosNaX/dR3bVa
dHh0mg1ZX09q6aDa2a1+R8QPNnTnBv9cNdX8vJzB6l1syE5F8HcF1VUl5xnn
2thQiGFG0wequZd7Fk43Us9X6FLIpjrLZaZ57QN1/QBLywiq/ywSvmsVsqEg
0+Mtz6k+enfrUGQOG8oz79geSHWDk8qiciobihH5zPuQ6jGHJ7zOr9jQp+SL
S2ZUK8y+EqkOYUPBlmqLRlS73E6V4QlkQ67FVfhZqhlvfTmW48mGOFO1C49R
LWK/48rULTZUlBJN5KXabkLSXvU6GzqhMXyci+ocO1XngKtsiGU7syYj1WBr
4SNpwob2h4pHzRpRwORqRob1aTZ0TVlcvobqpOGy/Hx1NmQ4aGdTQvWk1bcy
OjIb0t9aa5RJ9f0rf5riDrEhckfr5ZdUvzbHFn7sZUPSFEK9KdW9vbr/ZHaz
obbyhshzVIubWTG687EhHrnR7WpUF5j47eTdwoZ4HUtZ91LdeqlV1XCDFQln
MhYPGlKAt3tYPWmFFd3zNS5qodri4tKZhT+s6GF4y3lE9W8jAYvnk6zoS2Xr
5WiqSQY2T792saIw7ianc1Trt917JtDOikxupUqoUh11PjDc/hsrigv0thOn
WkovJ421lhVJsO0PWjOggMbZ1S8n81lRr7MgSxzVnprB2yueUeOxlJn/pk8B
hpoDkyb+rIh/T8GZQqp9oAetPWVFX49PikRRHXBI8sbhB6woeXuh0hWqX4hW
VRVcZ0X1C7pFk+ep8zDONkLfhhUJbn+S/oXqcH6uW/OWrKhjSJ8hi+ooHhNe
xUusaNLTYMSe6kTC39vvTrEiwgsGjkE9av0eRJ/UOsGKHKTIf8qpTls6zj+O
sSJ9C3bFSKozZ0JqxQ+xIlX1R2/PUF3wQ0owQZQV7ck5NJulS+0nJaYNEURW
pAzDGafPUUDzAEOc8joL+vG72Y2X6obsjLsdyyyI5e72W2NnqfM6ZWnPtlkW
ZOSdWORFdcer564hvSxI7gjZOvsMBS5yq5yR/86C/DwSll2p/h7YK/KVwoK2
nfv9GqjufSLTxPaVBZmm1ow36lDgl8MnUZ9SFhTnEGHRpU2B61P2q6LF1PV4
uXsjqJ66xtNcnceC9umUnbpI9ayZuTsxgwVxx71va9Oi9nedlZYHESzIJND1
WvlpCrBKyT12usuCVrwl0o1OUfv7azErSUcWVFOhXsZI9R1GoVM/b7AgGadu
33wNanwDnFynrFmQd+/8AXaqM6JmY/j0WVBsbMqtDHVqvCzjT76dZUFir7c0
nqWax3XA5qkWC/q3ex+aO0Hd7/qt8jNqLOgJ0w7r/VRrkQoqKhVZEFdg/VyS
GgUe3c9McJGl1lPUzPsY1XljST4ykizIY8rFv51MAb5PYWfDdrMgZKPXtg4U
mHjo2mvLyYJyOr2ZjmEUCPx9ZINrihklpWoo2B2hADLdP1jzixntjpGJHD1M
e5+VqXUbYkbCchaSV6g2SRMIHvrOjK4cJKboq1D7seWGYHE9Mxpk+W68+xD1
fa618phZKjOixNT4GuyjAEHto/DWRGZ07KV3ZK4idR5l5zHWxTCj6GvKLiSq
owMTv+4LZ0an6ZqcS+Wp80bd25TRmxl9uCfVuy5DAY4izYeZV5iRvpzqUy5x
an+l51qgN2NGaj/+qBwUo0CRHsXW+BIz2mfqNHhBlAJm05YGbLrM6HzD+NpL
EQq8F3WXsT3GjJgIG9qjQtR+/jKre+8uZiS4dvzv5A4KiA44nXPfyoxEz5RK
/NxOvd/yRz61cjKj+SNfahu2USC+oTbHk5EZJTS9lX3DQ/1/Qez37ZtnQhV/
a9u4OCnwyoHn0JuvTEi0L7qXzEAB24qud/N1TGiJ7MG+SEftv6RYYa0aJlRU
UyaaTKTO61QZjpWPTIgotG1ueaMdDvWpDxulMaHIiCP5N1faoeuM6wseLyb0
hCHR/ul0O2S+OcZy/RETUhgK6ZqfbAePcfpHlW5MyJ/RPsJkoh2kvEPtHByY
kHfYjynhX+1wvyyd3GTKhE5Eqgk49beDgMyPGX8VJrTfYlZFsrUdZO4e2rfv
ABMi3f+kbNvcDkfwUKcuBSYksfQkJ7GpHS7qn1wWl2BCibnHxkmN7fDK/T2x
ZjsTamMueZFa3Q4cXzy2bcwzIvWR5G3L+dT1d/QYJc8wooiLAsJTue0ga3Ew
UmeCETF95FTqyW4HrT/jgm8GGBF9surZjHft4COgL6HSzIhuboVq9sR22Lgh
duROFiNq2vuKhRJKXb/48QPedEbkc8jgokswdX3i9wqUxIjcjAuO8wRR4w0L
PsEVzYjOJd2rO+zbDi7lSzqZAYzoaq/kb42H7TDD0WA+eo0RnRVr2mNo2w59
7295XxZhRDqjtxJ2HmmH4K7wS0mCjKiowMiy4lA7YHS43Mwu6vH791ZN97dD
tCFXhycXI3p35XCkryy1Hv+yxdP/MaDhkOMDXkLt0Hp67vNSBwOy+/TcUoDY
DrXDTmxhgQxI/nhQSCZqg7uc0b29PgxINsQ/z7qsDURVavIknzCg4gJGse0l
beAVuM247B4D4u29tWiW3QZkpcJ3Q1cZUI1GUEHY2zb4+PivjhKZAXV9a1AJ
vtcGWXxuId/+0KOSizF/kyTbgFN5wVFhlh7dCv+GP9jbBg76Nw1CJunRoxf1
LGeE2kApyJL33CA96l8y4u3c2gZF66cTm77Ro9BLyXeeb7QC6ucr/ppBj8IC
KzjjW1uhNaWst9GCev6ju6/Jd1thSYkgX/eFDlVe13T2T26BF9qcP13r6FDs
xarQHbEtIG8tECRRQ4duegWpxkS0gE2YyoRPKR0qvv7A5G1gC3QsOaSczKBD
RGONaGPHFigqGxCq9qVD72ZqOHtVW8BFo5qr4jgdehoREar4tRkWL/nM5ZcQ
qfPs3J7Pn77Bpu1TSfVnBPTs7Zl5vkeN4M5BzxoyuYF3xFdcGbD6DJy806m9
z/7hfjfNt1iJVkPkw+bbB1VX8c0O/eXHNRUQtJnreYd7GSdLdh8kZZdAu3pY
JIl5Gbe+d7KvPLwEBALu5aasL+ErHIHnPR+VQOY2GPg+sYQLhcu8PHa2BL5K
fSGf+LSER6aKqgXOFAOnwejm1vtLuOve1DA3pWJ4nsb3qHDoL16UkN0e8rkQ
wvW87i8XLeKhN0tNq8TywSAkWEbp/SK+I6x71+6t+cDzNfKHffIi3jgjfCuA
kA/PNHOwvpeL+HB9WlNQTx48hZ90tQ6LuPnVX43wIg/uyCkHvpRexFPzTv25
RJcH51gn38jH/MFV5SjVB6dzgKXyfLm19wK+Zl+ScKYrC6SHiYENDxbwzNXj
J4I+ZoE2c+5FRecFXPa90+BQdBYEn+H6s3ZlAQ9avTTy0zILtn9vlH5JXsAN
hN+Yr01lgsiCekTl2jz+YtvLq83MmXBUTNlJ6M48Li5dauF5LgPMNEfJT+3m
cfqak5SaAxnw2D6Mc8JiHk9O1PTfzZcBVbkL6YXn5vHh8aPufCPpcBrL6deR
m8d/baGbEX2QDkYXpM66j83hUzkpRx6/TwNHf36pLrM5fC2zxH5kbyq02dvw
eRnO4dnMzZvTrKlw8Gwuu7zOHM5VWtbENZsCKzynZrwOz+G84ys9MWUp4Bl9
J19h2xwesTP005hRCjzPrj/mWz+LJ2uy3+YPTYZciuv5QwdmcX03wqz7tiTY
WlJ1YkB6Fhc6xVjst54Id6M4DgYJz+Jy2rtvZ44mwmHzxJ2DnLO4ZSm7s8qH
RHjS8vYDf9tvvBhMFy3MEyFhSbup4dhvXNfqS+pQdgIMnkhbktw6g1cK5gvs
t4yH5yvPSBXEGVzt7US/wdl4IGc7CevPTuMll5wK/FTjIZYPtB9+ncZvXGo1
EN8RDyZzlJgW32n8vIPEpPKXOOiIode4/28KbzlsuJyqGgcNK2avGkcm8Rdm
i62FZ2PALVs93aJ9Ejd/vMvglFwMSF6VrlismsRjc856rrDHgG/z4q/dcZO4
58sqO7z+LainBqk6X57EeVTOvW/SeAsVBh8HBb5N4HGxksRhjWjIz96+71bJ
OG7B8XdGySwKBoeUnj1LHsdr/gVcl8eigGeX7njmi3Ecv/HSSUcoCu54BMRP
3hjHDWo+JP7ujQQlvc0tdsLjuCzf2P5vlpGQNz82bxUwhjOZqvp23oyAvINl
BZdMfuE25ovoQEw4DNp1c7ud/oX7MmWFn3wSDjwxf29EHPqFT7E+33vXJhzu
MO0T6+T+hfctcYwoKISDUmdymMGnUTxoe/w8d1UY5N0LcTknP4pTKNZj16Zf
QW7pFZWThBFcgrthe+Kll5D28fPXsulhXOHsIqpRewlxH2StDvQM44ninLEE
6ZcQWvz32d6iYdyxdmtY7eoLuJUfMES8MYznuMu0+r19AbLv8oIrOofwSZcx
q+ZfzyEtim5EOXsQ950W/VgZHApxkdfuZ0cP4iemNF6suoVCRMQXLomAQZz8
WCJC62oo+IaHH9lxdRDXEUjsVDgaCjYvpEP/8A3iEsKDD3g7QmBvwHnVXJ8B
/Pek2OXlyGCIcUt4Lm3Wj8/qHnNQOB4EnnYXLg3o9OMG6JST154gsL5EEo5Q
7ccvup1hX9sIBOnDLjmMvP14eXGtu0FZINyYexlTmdqHv45//i9TJRDe/B7k
lX/Zi2vorJzoVA6A1UkPbubrP3BbLfc8hXN+gBonYp1P/cAHbY9zlOzzA593
BgqDoj/wQyVjdBbb/IDHXvJseV8P/sa/GuP87gtSE98CnQx78BcOdAP9Nr5w
cUyIuZ/8HT8xfbNW3ccHiofL1kt2duHnPOIOH+15Co8qw/hV5ztxrZseIgb4
U9CIuXW4/Esn3tRGvxGU+BTaL+y5W+nZiQuF8W9a3HgK8w1eU/WTHfi+pgKl
4E0vkM3V6vmOU/BRKV1Cv4wXJDzsKvl3rQ1PcNLee/+NJ9hdzu14dLwNr44v
Ervt5QmKygF/CIJteNrvwQ0Pe0/AZ1UVGVtb8aNvpXJXVT3hh1VsKsfRVpx/
EnqO93rATi2bcCHuFnxd6BczSdIDdsmN/fJw+Iq/XWTqs4t9CAYZ3VOe0l/x
c2GD2e33HkKoROPck+Ev+LfKg0725x8Cs8j7tacXv+AjaQFbBZgewt8dzlz+
5EZcbPcReu/bD6CNQFB+zl2PP7vmsnFOyx2C23f6xOfU4L7wLjFByA0a9FgD
E67X4FKBOg0q6/eA8dtqSKJoDX6HOzJlsecePKj/GZkcUY333lLrHIy8B/YV
iVnpHlX43UPPm3123gPNNAVKri7CmxNf/2zjdQWm5WDNpwPFeN7zAUb3w3eh
VDNHQPNiMR4bdnW3guBdcIhsmWX/VoSHyeXM0xPvQs+RbZEvywpxv+RCjKve
GXIeRo4nvM7HLdgfZ6VddoYCJs3tXm/e42WeHTyNXk7w4ciw0tqvKHx4Y7/f
3JQjBD+082KrjsDZueLxR+2OYIVm2nhjwnG7b78+S5Y5AofGirOKwQtcJ/hY
V2ugI1QlCNya9PDCn1nwvnOXdYTv5kJbIfohdL3W+8rk5AD6hYfNOHqSIPjP
w6RHkregV+bx8QGGFPj9uDSqcfUm2CbUiBfKp0Km9Pw+g6834UHIud8mT9LB
9d2Mxk6nm5Bkd9UjU+o9mGfPFZtW3oC1YM0WsbJ8uHtfdG+erT3onHBeqGEr
h1wn0YCKflt4vhE+8j62HDRF3PToP9pCx4cPnREHKkC7Oazf/aUtWChult4w
w4HuTaoL8ylb+Mh58jhbNwKJuUmmJznXoGntG1uebDWkur/MbAmygSXKyBti
+2cQ3zd3zNDOGqZuFze68dYB/UqLWbuWNQyy+q/NmdVBk6z4i/uy1tB0TNZ4
cLwOpF1svFRnrSAp+Q5f9WY9nOQNrA5wswLdu5sRT2W+QFXdLk3Vl1fAIIDU
6KXWDEzJwtGNAxZw7sLiIFg3g+X3womfny1AS7R3dc2nGa7HfDvL894CoCJb
2vlLMwSFLNT137cA6fnzAdYXWsDbw8KWe4cFEC5Hnj55sxXexZsJHUg2g0wp
8TqmyHZ4sMvndmmtMSSvxOq9r2iH5OlsDtnXxhBTz9tjONwOb5p2Wg5eM4bn
dhwzSQoU0MgT7RNiMwaX9D9bT9RSwEe0q+ar3mVQk6o2f7zQAfzS+tUfpy9C
h6TF0tKZbjhoq3k2ztAImpe7PWKdukHhZc1rKXkjqK/TZzsV2Q30MwsJM0xG
UGZ7SjB8uBtUugubpT8YQnya/IkD7t/BLUtZwmqPIdyQXA92SOsBM0s/4bAV
fSBKRomP0/WCr1uhKhdFD0R3OGvziPeChUPc63dFeqDBcNZB9XQvxD6OuOUc
oQdB/cSPz0J6oVBUIDPARA92RdqdURLoA3V+c+/qUV3Y6S3wsZyvH+q+hP0+
TK8L8n8n+5kODADbvb6s9EIdIE9fdHpzcgAKX/vKHPbSAb3hT4yKFwZAXK+J
brueDji3xkpfuj8AeaYz0jnT2vDhvb7zOzQAQirlCteltUH9Wimz7plBWIw7
k/429zQYdwXIR9oMAZ+zy1b/PxqQxJDjS3EdAlV2hp9VjRowo0jp3+I/BHeM
lwxVEzXgif/ulwHvhmDnK1vuX+c1IO1IwZL77BAc1jOs2VN4Ev5G/0Rm94fB
toYt+PtTdQi1VNQXDR4BTSza9fK54/D9mWGmRewIyLZKi3TIHQfRj/cZ3+aM
gGWVipQv6TgU83wq3t42AqVytVJ+DWrQW31JgHnXKKyvhGAbmmogI+41Mh4/
Cr8r33U+PkeG2gmKa3bhL8jMZ+/fSncMQiKq/S3rfsE1vf1X8NqjcFEj983W
nl9A/hnKSwk8ChPxQbgrYQyuLrsXHt15FDgvnWAma4+Bs50TNqmkCkafc8Kb
+8fg9OS7CHh4GEYSAwvn2CdALVvR8eOlQ/Bez+1zotAE3Hi738lP+hC4bNp0
G+6bgBrS4sn3aweBxfj4+gejCVCSktm7GXMQZHlW1J/ET4Cgx+7Q1PED4Oxh
086jPAnZegooNng/MJiqLey/MgXbpOOO5WzdB4O8fa0pTlTjP2YWZxUBdTzI
4/WegouWR8/fb1IEd90Sx83UKUjUtWs8FqAI8ycUZhump6DvkJZ/MoMi9EsL
Tlu6TcOLxHXCWVZ5KF1eHg0JnQE/1cwd/0ylIaIwvJYuYQb2ehwx6BCXBpc7
B1Lu5s+A98nXjTtnpUBp6vZV444ZkHObnfnqJQVpfaND4gK/YcdnYdLfbEkI
r23vL0/9DewTL2rS+STgsNbF6bvms+BhLlxsfVAUbO0CD0/YzkL99mHKZSZR
CPer8Da7Qz0uf78lv3MvLNSKCp3ynoUVt8C8QPe9kHVi7uyud7NgWXCuB/8s
AsLHAnI+/p2Fn7KW/fGOwqBrUv5PnjAHKtjGe2dNYXjkPquZyDoHJgo92bG7
haHng9FAoMAc7Lk7/XmtaQ+EHdrLY3p8DsTlrtp+3LcHWBXLnDafzcG3iLnD
zziEQOXcb9zp9Rxsxx7dSfwkCDa3REhjcXOQ2ffkvtEjQajO9Etqzp+D02Fa
14PnBOChlCElvnsOJL47O4v18cOcyMwhdbF5oITUF/+l8EL39j2rfmXz4NvE
2ZAivQOSzmVLtXyaB/o3OvYHF7fDLX/sEu+3eVhcrg4/grYDw7pJSfrAPFyM
0JhLurAd5IcjXRqYFqD20rr4QOA28MrbusCutwAenhZX9rFuhbNTCSIGlxfA
8P4yO/6dB3aJK52PtloAtX9X7QszeeB9pG6urMsCiNwNnejS44Euj2e3z7xZ
gKSETZ/huC0gd455KnhkAdpvPbhmrscNy37h/J0zC8D1rJGJV4obqqrEtHcv
L0C+flOrPJEbLh5Sz8hm+wMWBsZqXdZc8ETQ07ZZ4Q/EBoSI7FfihM7J1eEt
9/8AZXt0/YlpdvD0m/n5inMRmm7ffGTTwwwd9YbOx3YuQlfqt5QTicwgw17O
Nrp7EST3eF+LtmcGyrNAZeV9ixCx3/lb2ToTSL2UfN6tvwgMStmf30owQcvb
K+pCkYuAX7Pr2vKcAcT76r/Xxi+CoPfr8yoWDOC+Z5/j7YxFMKJH440KDCCW
uPkWlS6CeK0Xn3IzPbilRS9d6V0E5vNaN5h20INwfkdG6t6/UFlP+fqpgAgu
f46p6cr+pfZHBrYLAURoPJjcuXzgL3j6EJdPWRDhbokTg5bGX9iSq+enwUGE
+nJusynbv8DxeWo9yZ4AjnVaW/a9/wvGwl96nO9sYJU/K1w/qiyBfV+mCb/Y
KtZhHjAheHwJji4RBHb0r2ATA4amntpLsO5r+g5/s4JtG506rmm2BD0iDgd2
bF/BbGb4ODu8lsBfJFo4gmsZI226JM01UX0wW99d5C+2x/P4TsOuJchija3y
/7WIHaTnDCgZWIKdsWXnFbMWMTPm5NuP/yyB6rrbjpNHFrFcrrYjHHzLcNM1
0ELX+A9mtEehRfLqMsye+07aLJrHriesngi6tQyxOn7Tn57OY49Fa4t+uy5D
yC6tapL+PJYmaRpdFLAMaYudOrvm5rA1xcBr6jnLMBods7Fz/xwWp/Zr3Xxt
GSblShNvvvmNTV6JlX79fAWsmYYEvdMmsICp4wzGUSugVHKZoc16ApN2Gf0p
lLgCFH3F8grhCczWX+55SsEKuAqNOdx9O44NZ5cuFXauwM2jvrMPosewH/86
a9oFV8GxTYXxV/Eo5u7jHhMhvgoeR/zFZh+OYnzcu++ZKKwCeCtp+KiPYhf3
2sgMk1dB3Dwp4HzbCNZ++s/zBatV+Ll0b4B/aRj7Es5tzpOxCrGVuikBF4Yw
+z0FKh15qxC+elRoWGwIY8u4wBNVugokg5Cp0j+DmGZ5zCfhr6vA5W5yZeDl
IFYzJCu7b3YVIu8K74jsHMDKFE+vnDu0Bkax65xLTv1YZoPny2dVa3DjclbH
7HIP9uuWp65Y4xrkrC+/bMrtwUS2enKUt63BO7XH/Fb2PViEsYfv1PAavFTt
ufqu/zvmPfXIXZv5H0j4EjnM2roxM84HVqw6/2AyuDe5ldKJRea574k3+AfT
lpVZLJGdWLuR+08V03+w8uSNSr9JJ6Yde/+C7a1/sPPSnz670Q5MRdFN+3Po
P9iyelAmf5OCbTnvst+b8g8K7JT1I7TbMZ2/d2cFev9BRs8uh/id7Zhf1N2s
gtF/oBsf4XxouA3bGHSWGF76B5FifOZzj9qwCScn/hN86/BiqOL69IdWrCbM
gZ5ovg6ev0rG955swVy6r7c9GFuHK0oukSZWXzHHiqYb3+bWgWPLuy6mLV+x
G4lKTCJr68B8eZwUjH/BLG+uqnzm2IC9pO9f9XZ/wXSIATFb9m+A/XwFaWy8
AdsrlW6X8nADIvm3uSon1WFCnBz0Kz4bEEEserNsVoftWnCI1g7dAGW5w+Ov
+OowznKVb7MJGxCm7Vqt+uoztqpbd0C1bgNsk54eqwqtxZrv/dr8xrMJH4Od
hhkya7BGU+1IEYFN+GdtxF/rWIPVHs9Wuiu2Cex18srvVWuwUpLrVT6VTVDO
j3k81F2NpcQxNVqbbkL1NcHmj3NV2MM6sfCV1E241JJlPK9Vid3L8lfQyd2E
9kvCanPsldidF9N1MR83wSaOccSrGWE3XfWuzFVvgsjJlIqCCITZmhSunfi6
CUmfNDyeWSPMSo03LLxjExpdWwV3HkCYmfhD+fG+TSixYU9IYEaYway6ZfD8
JtRkPMM6HHHsHCVttX9tEzrk1JTiJXFMzJSpIZKOQC5l7iHs6q/A9hy/bjXF
TiDXFipf29CswLaTlMKfCxHIUQKGb4Tpy7GLdLo7JkUI5KmWakViYRkWtXwz
XF2CQDYV4s0ysynD9oykhy8pEMjWOm5b8PpSTLpC+LWpGoG8NhHqPhT7EcMc
tkRIWRPIFqFtotoBJZiHjcIuL1sCuVz79Zra6RKs2uRMxI8bBPLRbsZb11lK
MM3T/hEhdwnkg93jKRW+xZi+CDFy0YdAhpgsj8hnRVj4rt285wIJZEsTTUZl
vSKsi/NYZFoIgWzUnrLPdXsRZrZ2L9I4gkCmk0XBkrGFmG37bGRlOoHcYqca
UZhXgGU0cPLxvyeQTzk+nMRdCrBpJBvlnEcgD01e+T51uABTLNbiayoikJlF
xi+Wr+Zjd7JsoyRKCeSY03Mv7IvzscJEHz5PnED+2m3p9dsxH1uOTIr6Xk2t
T8zVxkWJfOxIaBXfgToC2X/hzi2J1jzsoU9/1LMvBHJjr0+dmXMehh5s8P1q
JpDbjPzq7eTzMDongTdkCoGssK9s7iAhDztpd4Q/qptA/lE9bfl6KBfzM7/4
ZuEngezmNH09rDUXazB04T8zSCCrXkrasP2ci3HovHqTMkogF3vnHCouz8V0
j+fxEyYJ5Ed2kzoy+bnYC5XmN5d+U+O/OvEmOTkXa5ef4c9fIJCDtujN73qZ
i+0UI0WTlglkj/RnW5rv5WKX+aUFbP5R71+NxvebBrlY9BbNaJxAJA+21b4b
lMjFepltBHgZiWR5z8ID5xZzMOENr+g7rERySeUPv6qKHCxlAo8W4yGSj+my
ruTp5GBj/T8FHu8gkjlUlw0ctuZgMp1r0V18RLKGdvNYSlE2dusrr6DSbiJ5
wdLF5LVWNpZTrfw2cC+RXEcYejY5/B5b+GAoOCJBJG8r67M29n2PHcpxeovJ
EsmEj+jlD8X3WGl09tu5A0Ryb/PPgI+vsrBDDD/efjlMJJfz+eWUnszCcuxZ
YlIxIpnf5FWl3GomlnLEMsZUk5rPVvSo1yoTC+3iia03JpIj6Zocefa8w2y2
ucTFPyGSz1+uhNTkNGzAPSHugS+RLJYtNiptkoaZDDXFXQgikl9dq267zJOG
6eVJxHOEE8ll5gbecY9TsaO63fFu6USywYBkgo1ZClZczJhg8J5IPilgdkdW
OQVT2q2UoJBPJO8x+tsWwZqCZfqYJbCVEMlmbfJG7nnJmMRMQMJIGZEsvFKy
pnYkGUswLE5AlUTyzZjzg0nFSZhA+VDCm1oiOV7s6adW7STstSh3oksjkbzX
gUD3iDcJ4wk6mqjXTK23/m5vF6Yk7NmCbaIshUi+onvnd9VSIsZiHJbI/J1I
HmUOMWAbS8SeVFUmDvYSyX4rj2b42hKxdamZxPIhan3ktoU15ydi957zJUWM
EcmVfa25Af6J2PyyRpLTNJH841tl21ujROymhVPS2Xki2UYp3eS4YCI29jk2
SWqJSObyu+Cwqy8Bs1L4ksTwj0hmGI5a/vE2AesNX07qI9CRiacHtl01TsDa
ruolh7HRkZ2tKf1j3fFYBT19ClGIjnxzZfoL20IcdtheIeWHCB3ZyYFhYr0i
DitoNU4plqAjT31jEXkYFIelxxek3NxHR74RyPKgUDoOE2UdSNE8REf+vfd8
7ofFWCzGgSN1ryod2TuOUS0mPBZ7BTap3ep0ZA3Xz5Z13TEYV+qL1ILTdGSP
21pfXz2MwQI48dSQs3TkbO0C4+i9MRiDy2TqdX06suqCPmGh/i32+OfOtJMX
6cgFeWMz+Q5vsVV19bQ9pnTkixlWuz6IvMWcMx3S1izpyPJtKn3Wn6Ox31vf
pnXYUPPZrzoQcjMau+5en5ZrT0d+jTSkZ52jsf8HWN2CaQ==
          "]]}, 
        Annotation[#, "Charting`Private`Tag$291448#1"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          AbsoluteThickness[1.6], 
          RGBColor[0, 0, 1], 
          Dashing[{Small, Small}]], 
         LineBox[CompressedData["
1:eJxF13k4VVsbAPAzmTKmlCgZUghF5mSfDN0jyZAGcyEpkSRzynhDlAxFoSJT
SRlKob3q05FZgyKRhEoSDg7O9K3VvU/3r/P8nrXPu9693rX3u7aCxwn7wyQC
gaBPJBDQL+aDVT7wlQR3pYukrd3o2FUW2y7hqCS4/a+nU59MOR35z4U1uloU
z/8sSFa/v8/pP7/Olb7Hov0Xz+vNVKnF+v/GE00Kb/YOLv0z3iP4omXPzqWg
6l+Xa2WWiFL+c4zT4YSmBok//w++52A5FCrxZ/woyVyEs0XiT7w7xg1XhLaJ
/xk3iJ3aKrJe7I/jcpqaPqWKgvx//epB7p4qpgj46187r60spXmKgOwV/3jk
Ap3X3yEMmsT/cXO007fj5sLAe/72P/kGT3SxHi8BpmX/+Iq7CeOj0hLQlFT4
29G0VKmnKUJg9vGN335R69pylSgE7E5n/raQqsbZU6cFgXTnqd9O/mXEFR8X
ABGf5H5b9CHtzF0PAUA3JK1Etk0v71WSFgDbkqKkkDMCJPWutfMDQ/yp5O/1
tA65LBnLDySeNIghr974cSLRgB+clTcU/H19xD1ioxQ/aDQfoPyO33ZuOXea
DxjsmCYgJ6zes8Ggiw/obpNl74LmHlc2CiznA6biGnPIRpmpZuLX+QB/U900
cnADc9fdJD6Qsn9iArly5OA+yzA+EBzbMIY8IdriPnqEDxTbrh1FVtPbcjR2
Hx9wTRX6jOztdj1Q3oIPGJ8P7EO+lcAX2bCFDzynWXcjf7rnH++kyAfYtqkd
yCYkXGNCmg+Q191qQb6+V/xdtDgfKL0VQ0c+wK5YX8qmAKb9jQbkhzaETmMG
BVy27a1FliqwDen6TgG0gO4q5NeWU03MdxSwNka4FHlz7vaA5HYK8HalFiBf
nEyTXttIAbErNHKRd13RPrrjAQVUjctfQmYBp22X7lKAJ+POBeSysZilH4op
4KWDYiKykMmbx355FNAv2xmNTB8+tSQzmQJEibKByKfFrg8MJFDAFwuGH/I6
g8ZKlRgYz7rqKHJs8nLn+lAK+DXy7SAyVbvmzpfDFHBdb9Vu5F/O/Wc1DlEA
z1/eEjkvns8hxIUCZBWEzZE5PXvZS/ZQwDOjREPkJ2fndmlTKeBE7Uol5KOl
cgqRWyng4oaNa5Cl3+yYfaFHAcZ/rVmJHLL+Sq6jBgXcCLi2BFmnXX/inAwF
6PtO/rJC+18m7GLnDBlI9zlXIYvpM05umiSDD926d5AD9vg5XPxBBpEuI7eQ
tS8cWmUzRAYRHc/SkB9yLAs6Oslg3TEXP+SVMo3xmq1ksFKuwxM5VM/EJ5VO
Bh5W4k7IW09oa+xuIIPvUWwLZDAo86i9jAyexh6XRVbgZGRrFJEB94CZBHLM
KvHIlJtk0M4ZoiCb25O3W18lg8cvQn/uRM/z/360tMXD+U+O1iG/LqofaD0I
x09wdyFv+Z/es40uMJ7rzDbkjE/3C5L3k4GGRaEm8j7p2z5Wu8kgqE5dHLkv
MWW6ZSsZmAaDNkto4yKBbjV9MnhddrgOOe959KMkbTJoCq8vRT7IOh25U5UM
3oeExSMPH3fjb1lBBmpHhA2RLRJ7vqlKkkH9kLQyctFt+9ZEUTKoCuqQQD46
sOOSJYUMXFQ1v9LQ82azSaZ5igRiPbIuIjO1CZov20hgteKP9r+gL1uJ9Ye8
JIG++ds1yJpeqy9saCSBH1tJucjemQZjCXUkUEW9dgz5HTOgyKKMBHZlaxGR
H9Z/lvvf3yRAOzUkswPavvtXe2AsCfDLabIt0Pw/OZGKZ0nANV6wH3n9Wpm+
6GASiNh3OBc5K9r+CvUwCZw5cFMaOXjH/8SfmpIAn8AM0wxa0v3VUz8TEiiI
zXmFfC/kk98aIxIQNKorQx4tYbVGapGALuOwM/I+YZ3zW+Xh/URuf2gKrddV
SHjMIYLg272221E9v1VW+CwQgfe+vxWQ/YnP3KRnicDnUc4UFa2ndn99yDgR
mFqPXURekSkVptdHBGE67EYMetYxYaqqlgi2WtaKm6D9a++hzV9NBItWK19t
Q/F2mpw6UEEEk82FacijRnMz7NtEwJQ0l0B+J+s9b5FOBJaKDmRj6Np+C957
PyIwXliBG0LPdCtQ1Y4SQeZX/0DkzR2cc5FeRGC27t065JKnNSQFZyIoaWpK
MIDOzlfmP0YjAqEviyb60OGH+ETZikQgG8WN0kH1cfxsvVuOCNaFWysiT9s1
pN5YRQR8ohcat0AfMz0tYSFBBJW6bXzIzkojy1I5BMB/SCpCC+3XkUYZ+R4C
MAlQ26QJzfOJUzFPIYBKDcGK9aherjO6o+cJQNX1lw7yuJ2X2fk4ArDIrnis
DP3ByNytLYIAeh0q8XUoHxFKhsMxAjgcnt6giNbrfizP6y8C2FgfdFYO9ddC
hoiAGQHEsZWm10C7ZHvKlJoQgGlJuifyzhgz3Z+6BNDlVmq6Gu0nB7Lv6XUE
ELM3aGoV9AAz5l08kQCGFYwkpND7t2ekyWaRh1ebL/NfDn312Yuh9Cke7jSQ
0bIMOig9Xnr1IA/fnH3hrCS0uj5fnHoDD6/9QukVR/0kiuJsHcLDy8pv0ZZA
R4iShS7+4OKxeWniXFcYvybqL+khLr4rrUKDA93juhh/s4eLfz1I2MmG1rk3
Rap+wcUHaHpnFqEndg+ye/K5eER9Rvsc9MFLDZOKDlx84Y3t8l/Q5stCex7W
c/DPR8o3DUAH1jFWUKs4uOfMOYl+6BueJ/Y2l3Lw8oxdv/qgOVXer/uyOHjI
A1DUC/3QYW8b4SQHz45isN9Cq2Zpg53KHHz4Zr9FK7TYqp/FAyls3C2Rhj1E
8YpEyBlxbPytzcBQNbS2jrqbZQQbb+i3i6uC3m/ju7zah42/Ot4G7qPr47+f
SzRj4+9ectfegdaaHnHUWWDhPUKbg/OhHVoHhJO8WLi/LSE8AXr0ANcbc2bh
N3atfR8HHTq65vmMHQu/+2BOKxb6OtE11B1j4WKU9sGz0CP6fcM6MixcqF1c
KQw6uPB9w6fORTws7YmBL3T2ma4TulsX8YLvP25ao3zI0l5J2ou4f8rlbivo
kkT3A59UF/FrDHmBndC1mRPUROlFnPZz0nMHqk+5iGT/7ALulbqFgkGvGKBV
x91fwKsciZ82QV82eTb/RnkBX4h6tSgOndMo+FNlzQKeJ7JxTBT61k7bz2eW
L+DKLw+8F4Z+sO9T8wbyAn5+Qu+OAHSXPycnYnAeL+sXMeS50DHxfMNtStfm
8a/7f/74CX2B9yA6UGIeH6wpLG6GfmuemS0iMI9nReBrmqBXJ4U+KOIwcZt6
rbRG6LvLsc8fxpj4c1fpABy6XbWNavaCiXcXz4jWQIs5jPKWhTNxtYwzzfnQ
aSUyUTVf5vACQ9LWU9BZdrHh8w9ncdcq709LoR0upm7UvjeLP6FY9otBS7Zn
f/S9PYu/HHnUIwydQrtv8il9Fs/OxV7yQcdh/SR6wCy+9f7XxAVnuB819JPT
1WZxYZe+rEFoG6Ef1zTzZnDuwY7Ou9CCz+wbvOIZOGhUyt8GrTZMTG6JZOCO
ikdFjaCtBB4c2BzEwMcPbgrTg061Fp9heTDwycemlpugpT60qqVTGXj6SGGr
PLQiw/zqM9Y07uzYtZ8Mbaysf0oucBrfOFAbT3eiYycTZVV73KbwIjF3ZQvo
IbMSpsqyCTzbb9UDAUc6VlUhpeVf+x0Py7T6WLEf1qvOw8CCMIJ73tmbRNkH
+2/YrTQ1t0Hc18PIa4MDHXs0XM+pXdmDn7/64YWuPTz/aHz7ei6gHScJh1t+
tqFj/POptLjPj3CrgYTWaGs6tqfG0E20rxA7YpU3TraC/bV75BrxbRNWbul8
V5UG66e6/iV/9lsslHK93ceCjhFVctZ/Jw1gMyUBUV6m8H3bk6SZ7f0FOyty
xzAPg+fBse6QipqvmNFzjgVpGx2juG5nbPEYx2RfjYQcNKJjhjsP/DztPol5
DAVdUzWgY71S8ovn66cxXZ1jdkv14PeYTYXqqxfT2IxCKhCC9k80cVzVOY05
5vTqkKApHJfa0s/TmGXYsDpDF/b74ezgFn4GFjj3yqkbOrZyGUPYjoG5Vyne
y4HWsBEYTx1hYHLJNZtVoaPPT/RniM1ivo1lMXt1YD9r3hu0beUspkPwKrSB
3ijcsGR07Sy25GNQiyV0d0qyvr7WLCYWk6NkAq2arpLWu2cWc7XrFdoA/SrX
w1wuexbTjxe6vbAFnteq3pUVK81hXipEkVvQz/qfhjwxYGInhR5b87ThfO5J
Y2tMmdiUdLrkIvTY572u0VZM7BJL4uMM9PLRcVOaGxP7IciKHIP2npARexfL
xNqqwue6oUV4wYVTHUzM8OZ9ejn0PvlNr1QOz2N/8x2hH4T+4ZGvdiVtAXN2
bHj0RouOJY2bUpxzFrAeq/HeTmi14NF+uYIFzLPvLq8V2idRI62oegEzbn2+
txF6uKKOWfN+AVOPP7WxBvoj+33j2zWLmL5/cvgV6LYsCXfJskUsJ88rzw3a
V77a4F3lImYj4afuDL2kbL9kTt0iVjFV0rAfmtaQ90KhfRGzu1P13Ra68Yu6
utbkIrZ1JsDHDLp+s+WCjR4Ly1ynwVWBvtsSnZ7ynIXtkDYNmttMx776R9sq
t7KwUt9zQQxoxWXRog1vWNhSpeSQSeirzuf+Hh9mYZ8rWPFj0PHjURFWAmxM
tnRb4wC0m1ikp9AuNqaNX6U3QS+1D94S383G2reLtl+D3jV3enL1ABvb23hJ
KRv6fM7p8upRNqbZkROZBc0dCtowzGRj8zM0vTTosVOnZM1kOJjgYuDLBOjG
zAAy0Z2DjedTvAOhg3uPvYn8xsHCMkVOWEGffNpxvHOKg8WXFw3SoI8XaPMr
sjhY/Uq2ww7oQ36LBk2iXCxQ+YnFdpQPMSlv6RYuZq/M3aEPraRaerToDBcL
OaTarATdFfqV1ynJw8Sn7Ku5m+hYq6tVtuJqHiYaTgtlQ9NNK7RPK/Mw3pGP
xovQdSIhh2UMeNiSN4MvZ6GLbvC3ernyMGcFxrKf0GdeKmctFPOwEJsLyz5C
S4loZ6XJEaisLxPVddAHSLYrfigSqCZBV748hs6Z98sy30CgekScXF4LLT9S
msXcRKDet1eKrIZWe6pwxXU7gXqFPu9eDm0SsPSqqheBKljgEpIP7fN2MvtZ
KYFamWlJiIcuaxGTkb1HoJLmVXxioX8C9ZygSgKVs17sdTR0YLlPzoY6AlXl
dXh5FMo3YTAnpY1A/ZZ2KjEU+rJB1zXHXwSqgUrKuuPofq9X5E7pEKlHpgrG
HaD1KB9z2wyJ1CdPvJz2QN/3FcwrNiFS3y4pa7ZD62F0KM+VRqRKpP8q3w19
qUcyv9mZSFU9HnuJBu29PPjGzRgiVc/U5Z4xtOQF4wK7LiJ1zc4Qg/XQKQyf
AvVuItUi2ahnHbSgc2aBwAci1f1qe5gSNEd1oqDhC5F69XZuozz0t6b8QlUm
kWosNRQhC/2UTC4iypGovk7h1KXQhr6bij4qkqjnxua44tDVr52LHm0gUbvO
dj4Vgy69WV3kp0WiqpvVmYpAZ2Dexb3mJGr3iy0uAtDixZeLqy1JVIymuZ4f
OkkML764m0TNDImbokCf7V9ZYnGARL3JunKBBL1obl4i70qidlTmOROhg+4G
lLAOkagZQ0x1AvSvZbkl77xJVNUzFTyuJjxvRzSXPPAlUT1Gi9o50P8HWhtT
Ng==
          "]]}, 
        Annotation[#, "Charting`Private`Tag$291448#2"]& ]}}, {}, {}}, {Ticks -> {
       Charting`ScaledTicks[{Log, Exp}], 
       Charting`ScaledTicks[{Log, Exp}]}, GridLines -> {None, None}, 
     FrameTicks -> {{
        Charting`ScaledTicks[{Log, Exp}], 
        Charting`ScaledFrameTicks[{Log, Exp}]}, {
        Charting`ScaledTicks[{Log, Exp}], 
        Charting`ScaledFrameTicks[{Log, Exp}]}}, DisplayFunction -> Identity, 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None}, DisplayFunction -> 
     Identity, DisplayFunction -> Identity, Ticks -> {Quiet[
        Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , Quiet[
        Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& }, 
     AxesOrigin -> {-4.605170185988091, 0.}, FrameTicks -> {{Quiet[
         Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , 
        Charting`ScaledFrameTicks[{Log, Exp}]}, {Quiet[
         Charting`ScaledTicks[{Log, Exp}][#, #2, {6, 6}]]& , 
        Charting`ScaledFrameTicks[{Log, Exp}]}}, GridLines -> {None, None}, 
     DisplayFunction -> Identity, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, PlotRangeClipping -> True, ImagePadding -> All, 
     DisplayFunction -> Identity, 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None, 
       "ClippingRange" -> {{{-4.605170045013494, 2.3025849520194486`}, {
        4.2360986294842946`*^-7, 0.3335266134236399}}, {{-4.605170045013494, 
        2.3025849520194486`}, {4.2360986294842946`*^-7, 
        0.3335266134236399}}}}, DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {
       FormBox["\"T\[Gamma](MeV)\"", TraditionalForm], None}, 
     AxesOrigin -> {0, 0}, CoordinatesToolOptions -> {"DisplayFunction" -> ({
         Exp[
          Part[#, 1]], 
         Exp[
          Part[#, 2]]}& ), "CopiedValueFunction" -> ({
         Exp[
          Part[#, 1]], 
         Exp[
          Part[#, 2]]}& )}, DisplayFunction :> Identity, 
     Frame -> {{False, False}, {False, False}}, 
     FrameLabel -> {{None, None}, {None, None}}, 
     FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, 
     GridLines -> {None, None}, GridLinesStyle -> Directive[
       GrayLevel[0.5, 0.4]], 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
       AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange -> 
     NCache[{{-4.605170185988091, 
         Log[10]}, {0., 0.3335266134236399}}, {{-4.605170185988091, 
       2.302585092994046}, {0., 0.3335266134236399}}], PlotRangeClipping -> 
     True, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.02]}}, Ticks -> {Automatic, Automatic}}],FormBox[
    FormBox[
     TemplateBox[{
      "\"\\!\\(\\*FractionBox[\\(T\[Gamma]\\), \\(T\[Nu]\\)]\\)\"", 
       "\"\\!\\(\\*FractionBox[\\(T\[Nu]\\), \\(T\[Alpha]\\)]\\)\""}, 
      "LineLegend", DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0]], {
                    LineBox[{{0, 10}, {40, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[1, 0, 0]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {40, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{Small, Small}]], {
                    LineBox[{{0, 10}, {40, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    AbsoluteThickness[1.6], 
                    RGBColor[0, 0, 1], 
                    Dashing[{Small, Small}]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {40, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[1, 0, 0], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> RGBColor[0.6666666666666666, 0., 0.], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"1", ",", "0", ",", "0"}], "]"}], NumberMarks -> 
                    False]], Appearance -> None, BaseStyle -> {}, 
                    BaselinePosition -> Baseline, DefaultBaseStyle -> {}, 
                    ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[1, 0, 0]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[1, 0, 0], Editable -> False, Selectable -> 
                    False]}], "]"}], ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0, 0, 1], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> RGBColor[0., 0., 0.6666666666666666], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0", ",", "0", ",", "1"}], "]"}], NumberMarks -> 
                    False]], Appearance -> None, BaseStyle -> {}, 
                    BaselinePosition -> Baseline, DefaultBaseStyle -> {}, 
                    ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0, 0, 1]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0, 0, 1], Editable -> False, Selectable -> 
                    False], ",", 
                   RowBox[{"Dashing", "[", 
                    RowBox[{"{", 
                    RowBox[{"Small", ",", "Small"}], "}"}], "]"}]}], "]"}]}], 
             "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{
  3.7842842820875063`*^9, {3.784284350636279*^9, 3.784284371386569*^9}, 
   3.784284536464419*^9, 3.784287709460148*^9, 3.784287773013843*^9, 
   3.784287831315805*^9, {3.784287897289638*^9, 3.7842879162470284`*^9}, 
   3.784287958085731*^9, 3.7842895215752287`*^9, 3.784380497576807*^9, 
   3.784380586797553*^9, 3.784380838162087*^9, {3.7843808691763687`*^9, 
   3.784380890925535*^9}, 3.784380976156371*^9, 3.784381016557641*^9, 
   3.784381124065312*^9, 3.784381172580337*^9, 3.784381362865377*^9, 
   3.7843824307274237`*^9, 3.784382517346518*^9, 3.7843825831170177`*^9, 
   3.786691349373405*^9, 3.7875841795320063`*^9, 3.7967338176166897`*^9, 
   3.799555375021298*^9},
 CellLabel->"Out[59]=",ExpressionUUID->"29c4df7b-9537-47f4-b606-71139576e023"]
}, Open  ]]
}, Closed]]
}, Open  ]]
},
Evaluator->"Local",
WindowSize->{945, 749},
WindowMargins->{{92, Automatic}, {Automatic, 39}},
FrontEndVersion->"11.3 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 5, \
2018)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 344, 6, 30, "Input",ExpressionUUID->"aacafbb2-6c83-4b7a-ac50-85d6ebd12eb8"],
Cell[905, 28, 948, 15, 58, "Text",ExpressionUUID->"59c482db-ea09-4e94-81d9-312f917c36af"],
Cell[CellGroupData[{
Cell[1878, 47, 1219, 17, 54, "Subsection",ExpressionUUID->"a910b207-0afe-443b-a530-630ec823168a"],
Cell[3100, 66, 1648, 27, 242, "Text",ExpressionUUID->"e68671e4-6b35-4720-aee7-68c7917205e3"]
}, Closed]],
Cell[CellGroupData[{
Cell[4785, 98, 1216, 17, 38, "Subsection",ExpressionUUID->"44eab734-09c8-49d6-a694-0efa772da8af"],
Cell[6004, 117, 819, 16, 173, "Text",ExpressionUUID->"bd9d2478-5942-4b26-8b7e-ef42b186c61e"],
Cell[6826, 135, 679, 16, 73, "Input",ExpressionUUID->"9494557c-ee11-41dc-82c2-e9af64174292"]
}, Closed]],
Cell[CellGroupData[{
Cell[7542, 156, 1213, 17, 38, "Subsection",ExpressionUUID->"6001c685-a11b-4e3b-a1b3-1976b70c9246"],
Cell[CellGroupData[{
Cell[8780, 177, 417, 8, 45, "Subsubsection",ExpressionUUID->"4b7b5908-48c7-4d67-a338-91b888461aa5"],
Cell[9200, 187, 7909, 158, 205, "Input",ExpressionUUID->"8ceb9456-0881-4f2f-bed9-20971dfdd32d"]
}, Closed]],
Cell[CellGroupData[{
Cell[17146, 350, 458, 7, 37, "Subsubsection",ExpressionUUID->"493c6bcd-d8bd-42af-8c28-b11f430f58ba"],
Cell[17607, 359, 5993, 157, 489, "Input",ExpressionUUID->"e6b93132-914c-4e2f-ac34-041af44745a6"]
}, Closed]],
Cell[CellGroupData[{
Cell[23637, 521, 420, 6, 37, "Subsubsection",ExpressionUUID->"80a01eff-b49d-44fa-9421-4af4ea950cef"],
Cell[24060, 529, 4698, 128, 392, "Input",ExpressionUUID->"99e71276-99cd-4934-a58a-86806b6fab9a"]
}, Closed]],
Cell[CellGroupData[{
Cell[28795, 662, 524, 8, 37, "Subsubsection",ExpressionUUID->"77d4f13b-2bce-4415-a123-b5e9cd28b555"],
Cell[29322, 672, 9626, 206, 846, "Input",ExpressionUUID->"d0349130-a313-4165-8202-eb6e99433ed6"]
}, Closed]],
Cell[CellGroupData[{
Cell[38985, 883, 597, 9, 37, "Subsubsection",ExpressionUUID->"40d4ace1-368e-4d46-a7db-e71149dadc8f"],
Cell[CellGroupData[{
Cell[39607, 896, 899, 18, 52, "Input",ExpressionUUID->"4319d5e2-33d4-4dde-b6ec-0667c929a531"],
Cell[40509, 916, 804, 15, 24, "Print",ExpressionUUID->"09480279-923b-427c-9917-b2602fde0576"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41350, 936, 6477, 173, 542, "Input",ExpressionUUID->"e9047461-364e-443c-99ae-825d269dffed"],
Cell[CellGroupData[{
Cell[47852, 1113, 1068, 20, 24, "Print",ExpressionUUID->"400b2419-6969-4730-b85e-b56b5cfd935b"],
Cell[48923, 1135, 1102, 21, 24, "Print",ExpressionUUID->"db7577a4-dcb7-45cc-bb54-160f7c44d740"],
Cell[50028, 1158, 1078, 20, 24, "Print",ExpressionUUID->"d53ea4d1-2aae-48b7-ace5-5c25eb65e77f"],
Cell[51109, 1180, 1253, 23, 24, "Print",ExpressionUUID->"c5861b73-6470-4ae1-9890-578ae0499f77"],
Cell[52365, 1205, 1176, 22, 27, "Print",ExpressionUUID->"419323ae-504a-402c-8916-360b7a00e31f"],
Cell[53544, 1229, 1182, 22, 27, "Print",ExpressionUUID->"bb25756a-610e-45ea-bb13-a370900e03a1"]
}, Open  ]],
Cell[54741, 1254, 34994, 636, 213, "Output",ExpressionUUID->"22183634-7a38-4e4f-bdfb-5a0f0091f41d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[89772, 1895, 949, 18, 52, "Input",ExpressionUUID->"da9b7e77-b9ff-46f3-b07d-ff45106f0d03"],
Cell[90724, 1915, 450, 10, 24, "Print",ExpressionUUID->"588fe484-f155-4b22-96e1-9ab66d686dfb"]
}, Open  ]],
Cell[CellGroupData[{
Cell[91211, 1930, 6477, 173, 542, "Input",ExpressionUUID->"2473a44e-6671-4310-855b-bde5549eb92c"],
Cell[CellGroupData[{
Cell[97713, 2107, 1065, 20, 24, "Print",ExpressionUUID->"c6ab6768-bdc3-44db-a55d-f31debb9122f"],
Cell[98781, 2129, 1100, 21, 24, "Print",ExpressionUUID->"01c1257e-c5db-45e7-9baf-7fe50b650c15"],
Cell[99884, 2152, 1078, 20, 24, "Print",ExpressionUUID->"6c42f006-0226-4cc0-a20d-086796ddc4c1"],
Cell[100965, 2174, 1253, 23, 24, "Print",ExpressionUUID->"f69e3031-bb5b-4ede-a901-7017e78bc47c"],
Cell[102221, 2199, 1174, 22, 27, "Print",ExpressionUUID->"b591cad5-5e0f-4d5c-9b7c-782562183a32"],
Cell[103398, 2223, 1180, 22, 27, "Print",ExpressionUUID->"3a902286-aa49-454b-87bc-2a858a4c151f"]
}, Open  ]],
Cell[104593, 2248, 33959, 619, 213, "Output",ExpressionUUID->"29c4df7b-9537-47f4-b606-71139576e023"]
}, Open  ]]
}, Closed]]
}, Open  ]]
}
]
*)