(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 58447, 1595] NotebookOptionsPosition[ 56915, 1537] NotebookOutlinePosition[ 57273, 1553] CellTagsIndexPosition[ 57230, 1550] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"Needs", "[", "\"\\"", "]"}]], "Input", CellChangeTimes->{{3.580451243442741*^9, 3.5804512494024143`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ SubscriptBox["En", "\[Alpha]"], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"22", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"4.6", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"26", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"6.5", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"28", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"7.6", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"29", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"8.1", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"30", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"8.7", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"40", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"15.8", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"42", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"17.4", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"47", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"21.8", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.580449460866584*^9, 3.580449564485923*^9}, { 3.580449826421934*^9, 3.580449836807396*^9}, {3.580450483567335*^9, 3.5804504936734123`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"22", ",", RowBox[{"{", RowBox[{ TemplateBox[{"4.6`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", RowBox[{"{", RowBox[{ TemplateBox[{"6.5`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", RowBox[{"{", RowBox[{ TemplateBox[{"7.6`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", RowBox[{"{", RowBox[{ TemplateBox[{"8.1`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", RowBox[{"{", RowBox[{ TemplateBox[{"8.7`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"40", ",", RowBox[{"{", RowBox[{ TemplateBox[{"15.8`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"42", ",", RowBox[{"{", RowBox[{ TemplateBox[{"17.4`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"47", ",", RowBox[{"{", RowBox[{ TemplateBox[{"21.8`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.58044956692302*^9, 3.580449837706252*^9, 3.580450495131833*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ SubscriptBox["En", "\[Beta]"], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"26", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"7.1", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"28", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"8.4", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"29", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"9.0", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"30", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"9.7", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"40", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"17.6", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"42", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"19.4", ",", ".2"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"47", ",", RowBox[{"Quantity", "[", RowBox[{ RowBox[{"{", RowBox[{"24.5", ",", ".3"}], "}"}], ",", "\"\\""}], "]"}]}], "}"}]}], "}"}]}]], "Input", CellChangeTimes->{{3.580449574168458*^9, 3.580449601908291*^9}, { 3.580449841191415*^9, 3.580449850022829*^9}, {3.580450497265316*^9, 3.580450499720324*^9}, {3.580451332832963*^9, 3.580451333664521*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"26", ",", RowBox[{"{", RowBox[{ TemplateBox[{"7.1`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", RowBox[{"{", RowBox[{ TemplateBox[{"8.4`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", RowBox[{"{", RowBox[{ TemplateBox[{"9.`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", RowBox[{"{", RowBox[{ TemplateBox[{"9.7`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"40", ",", RowBox[{"{", RowBox[{ TemplateBox[{"17.6`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"42", ",", RowBox[{"{", RowBox[{ TemplateBox[{"19.4`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.2`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"47", ",", RowBox[{"{", RowBox[{ TemplateBox[{"24.5`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )], ",", TemplateBox[{"0.3`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"keV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: kiloelectronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Kiloelectronvolts\""}], "]"}]& )]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.580449602808928*^9, 3.580449850625502*^9, 3.580450500484*^9, 3.580451334380588*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ry", "=", RowBox[{"Quantity", "[", RowBox[{"13.6", ",", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.580450121755941*^9, 3.58045012743661*^9}}], Cell[BoxData[ TemplateBox[{"13.6`"}, "QuantityUnit", DisplayFunction->(TooltipBox[ StyleBox[ RowBox[{#, StyleBox["\"eV\"", "QuantityUnitTraditionalLabel"]}], ShowStringCharacters -> False], "Unit: electronvolts"]& ), InterpretationFunction->(RowBox[{"Quantity", "[", RowBox[{#, ",", "\"Electronvolts\""}], "]"}]& )]], "Output", CellChangeTimes->{3.580450127837018*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"a", "=", RowBox[{ RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"En", "=", RowBox[{"First", "[", "#", "]"}]}], ",", RowBox[{"\[CapitalDelta]En", "=", RowBox[{"Last", "[", "#", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"En", "/", "Ry"}], "]"}], ",", FractionBox["\[CapitalDelta]En", RowBox[{"2", RowBox[{"Sqrt", "[", RowBox[{"En", " ", "Ry"}], "]"}]}]]}], "}"}]}], "]"}], "&"}], "/@", RowBox[{"Last", "/@", SubscriptBox["En", "\[Alpha]"]}]}]}]], "Input", CellChangeTimes->{{3.580449740289617*^9, 3.580449786470715*^9}, { 3.580449948755124*^9, 3.580449973073613*^9}, {3.580450016206522*^9, 3.5804500362624617`*^9}, {3.5804501255952387`*^9, 3.580450133407486*^9}, { 3.580450515174827*^9, 3.58045052447425*^9}, {3.580450772677656*^9, 3.5804508912213173`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"18.391174353957037`", ",", "0.3998081381295009`"}], "}"}], ",", RowBox[{"{", RowBox[{"21.861865804880157`", ",", "0.33633639699815626`"}], "}"}], ",", RowBox[{"{", RowBox[{"23.639448585188376`", ",", "0.3110453761208997`"}], "}"}], ",", RowBox[{"{", RowBox[{"24.40467650459882`", ",", "0.3012923025259114`"}], "}"}], ",", RowBox[{"{", RowBox[{"25.292407602933753`", ",", "0.2907173287693535`"}], "}"}], ",", RowBox[{"{", RowBox[{"34.08466966074855`", ",", "0.21572575734650987`"}], "}"}], ",", RowBox[{"{", RowBox[{"35.7688658571373`", ",", "0.20556819458124886`"}], "}"}], ",", RowBox[{"{", RowBox[{"40.03674782584854`", ",", "0.1836548065405897`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.580450527957347*^9, 3.580450894768358*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"b", "=", RowBox[{ RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"En", "=", RowBox[{"First", "[", "#", "]"}]}], ",", RowBox[{"\[CapitalDelta]En", "=", RowBox[{"Last", "[", "#", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"En", "/", "Ry"}], "]"}], ",", FractionBox["\[CapitalDelta]En", RowBox[{"2", RowBox[{"Sqrt", "[", RowBox[{"En", " ", "Ry"}], "]"}]}]]}], "}"}]}], "]"}], "&"}], "/@", RowBox[{"Last", "/@", SubscriptBox["En", "\[Beta]"]}]}]}]], "Input", CellChangeTimes->{{3.580450045620962*^9, 3.580450065270432*^9}, { 3.5804501386073914`*^9, 3.580450143223666*^9}, {3.580450530118785*^9, 3.580450533728911*^9}, {3.5804509014482203`*^9, 3.5804509064660788`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"22.848606599296417`", ",", "0.3218113605534707`"}], "}"}], ",", RowBox[{"{", RowBox[{"24.85250608738542`", ",", "0.2958631677069693`"}], "}"}], ",", RowBox[{"{", RowBox[{"25.724787771376324`", ",", "0.2858309752375147`"}], "}"}], ",", RowBox[{"{", RowBox[{"26.706465399180907`", ",", "0.2753243855585661`"}], "}"}], ",", RowBox[{"{", RowBox[{"35.97384670922507`", ",", "0.2043968563024152`"}], "}"}], ",", RowBox[{"{", RowBox[{"37.76864557056943`", ",", "0.19468374005448164`"}], "}"}], ",", RowBox[{"{", RowBox[{"42.44373438135827`", ",", "0.25985959825321386`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.580450068627583*^9, 3.580450146808926*^9, 3.580450537112287*^9, 3.580450910033532*^9, 3.580451340646832*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"sa", "=", RowBox[{ RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n1", "=", "1"}], ",", RowBox[{"n2", "=", "2"}], ",", RowBox[{"En", "=", RowBox[{"First", "[", RowBox[{"Last", "[", "#", "]"}], "]"}]}], ",", RowBox[{"\[CapitalDelta]En", "=", RowBox[{"Last", "[", RowBox[{"Last", "[", "#", "]"}], "]"}]}], ",", RowBox[{"Z", "=", RowBox[{"First", "[", "#", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Z", "-", RowBox[{"Sqrt", "[", FractionBox["En", RowBox[{"Ry", " ", RowBox[{"(", RowBox[{ FractionBox["1", SuperscriptBox["n1", "2"]], "-", FractionBox["1", SuperscriptBox["n2", "2"]]}], ")"}]}]], "]"}]}], ",", FractionBox["\[CapitalDelta]En", RowBox[{"2", RowBox[{"Sqrt", "[", RowBox[{"En", " ", "Ry", " ", RowBox[{"(", RowBox[{ FractionBox["1", SuperscriptBox["n1", "2"]], "-", FractionBox["1", SuperscriptBox["n2", "2"]]}], ")"}]}], "]"}]}]]}], "}"}]}], "]"}], "&"}], "/@", SubscriptBox["En", "\[Alpha]"]}]}]], "Input", CellChangeTimes->{{3.580450098054573*^9, 3.5804501145388327`*^9}, { 3.580450152373061*^9, 3.580450196707265*^9}, {3.5804505413255053`*^9, 3.580450586404819*^9}, {3.580450621500594*^9, 3.5804507588227463`*^9}, { 3.580450936523467*^9, 3.580450988367432*^9}, {3.580451045868162*^9, 3.580451050371841*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.7637010720591277`", ",", "0.46165867234654084`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.7560917851299358`", ",", "0.38836781869030873`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.7035159916943492`", ",", "0.35916426326717965`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.8199069011682703`", ",", "0.3479023839361942`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.7948433239849209`", ",", "0.3356914560461504`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.6424135922546981`", ",", "0.2490986481502867`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.6976713375483357`", ",", "0.2373697049566188`"}], "}"}], ",", RowBox[{"{", RowBox[{"0.7695457305383542`", ",", "0.21206630398835613`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5804501983173113`*^9, 3.5804505454723883`*^9, {3.5804507366275167`*^9, 3.580450763100288*^9}, {3.5804509658432207`*^9, 3.5804509922907352`*^9}, 3.580451055185425*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"sb", "=", RowBox[{ RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"n1", "=", "1"}], ",", RowBox[{"n2", "=", "3"}], ",", RowBox[{"En", "=", RowBox[{"First", "[", RowBox[{"Last", "[", "#", "]"}], "]"}]}], ",", RowBox[{"\[CapitalDelta]En", "=", RowBox[{"Last", "[", RowBox[{"Last", "[", "#", "]"}], "]"}]}], ",", RowBox[{"Z", "=", RowBox[{"First", "[", "#", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Z", "-", RowBox[{"Sqrt", "[", FractionBox["En", RowBox[{"Ry", " ", RowBox[{"(", RowBox[{ FractionBox["1", SuperscriptBox["n1", "2"]], "-", FractionBox["1", SuperscriptBox["n2", "2"]]}], ")"}]}]], "]"}]}], ",", FractionBox["\[CapitalDelta]En", RowBox[{"2", RowBox[{"Sqrt", "[", RowBox[{"En", " ", "Ry", " ", RowBox[{"(", RowBox[{ FractionBox["1", SuperscriptBox["n1", "2"]], "-", FractionBox["1", SuperscriptBox["n2", "2"]]}], ")"}]}], "]"}]}]]}], "}"}]}], "]"}], "&"}], "/@", SubscriptBox["En", "\[Beta]"]}]}]], "Input", CellChangeTimes->{{3.5804509206442013`*^9, 3.580450931403729*^9}, { 3.580450998171288*^9, 3.5804510031431103`*^9}, {3.580451053963378*^9, 3.580451054308941*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1.7653929994607047`", ",", "0.34133249296534224`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.6399366241947284`", ",", "0.31381027828339614`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.7147421834125431`", ",", "0.30316953129541624`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.673515822072929`", ",", "0.29202561008172245`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.843973569812377`", ",", "0.21679560471697512`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.9403019012286435`", ",", "0.2064932891689246`"}], "}"}], ",", RowBox[{"{", RowBox[{"1.9816214000914343`", ",", "0.27562272612188915`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.580450924319539*^9, 3.580451006857436*^9, 3.580451058590249*^9, 3.5804513470831614`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ SubscriptBox["\[Sigma]", "\[Alpha]"], "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Z", "=", RowBox[{"First", "[", RowBox[{ SubscriptBox["En", "\[Alpha]"], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[Sigma]", "=", RowBox[{"First", "[", RowBox[{"sa", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[CapitalDelta]\[Sigma]", "=", RowBox[{"Last", "[", RowBox[{"sa", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Z", ",", "\[Sigma]"}], "}"}], ",", RowBox[{"ErrorBar", "[", "\[CapitalDelta]\[Sigma]", "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "sa", "]"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5804510301403017`*^9, 3.5804510367165747`*^9}, { 3.580451084704589*^9, 3.580451089270938*^9}, {3.580451133219759*^9, 3.5804512155666733`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"22", ",", "0.7637010720591277`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.46165867234654084`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"26", ",", "0.7560917851299358`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.38836781869030873`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"28", ",", "0.7035159916943492`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.35916426326717965`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"29", ",", "0.8199069011682703`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.3479023839361942`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"30", ",", "0.7948433239849209`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.3356914560461504`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"40", ",", "0.6424135922546981`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.2490986481502867`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"42", ",", "0.6976713375483357`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.2373697049566188`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"47", ",", "0.7695457305383542`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.21206630398835613`", "]"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.58045121635207*^9, 3.5804512558526707`*^9, 3.5804513507805862`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ SubscriptBox["\[Sigma]", "\[Beta]"], "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Z", "=", RowBox[{"First", "[", RowBox[{ SubscriptBox["En", "\[Beta]"], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[Sigma]", "=", RowBox[{"First", "[", RowBox[{"sb", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[CapitalDelta]\[Sigma]", "=", RowBox[{"Last", "[", RowBox[{"sb", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"Z", ",", "\[Sigma]"}], "}"}], ",", RowBox[{"ErrorBar", "[", "\[CapitalDelta]\[Sigma]", "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "sb", "]"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5804512278549633`*^9, 3.580451237534871*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"26", ",", "1.7653929994607047`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.34133249296534224`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"28", ",", "1.6399366241947284`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.31381027828339614`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"29", ",", "1.7147421834125431`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.30316953129541624`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"30", ",", "1.673515822072929`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.29202561008172245`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"40", ",", "1.843973569812377`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.21679560471697512`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"42", ",", "1.9403019012286435`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.2064932891689246`", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"47", ",", "1.9816214000914343`"}], "}"}], ",", RowBox[{"ErrorBar", "[", "0.27562272612188915`", "]"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.5804512380375566`*^9, 3.5804512573344517`*^9}, 3.580451352216655*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ErrorListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ SubscriptBox["\[Sigma]", "\[Alpha]"], ",", SubscriptBox["\[Sigma]", "\[Beta]"]}], "}"}], ",", RowBox[{"Joined", "\[Rule]", "True"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"20", ",", "50"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "2.5"}], "}"}]}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\<\[Sigma]\>\""}], "}"}]}], ",", RowBox[{"PlotLegends", "\[Rule]", RowBox[{"{", RowBox[{ "\"\<\!\(\*SubscriptBox[\(\[Sigma]\), \(\[Alpha]\)]\)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(\[Sigma]\), \(\[Beta]\)]\)\>\""}], "}"}]}], ",", RowBox[{"ImageSize", "\[Rule]", "Full"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"GridLinesStyle", "\[Rule]", RowBox[{"Directive", "[", RowBox[{"Lighter", "[", "LightGray", "]"}], "]"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.5804512602527227`*^9, 3.5804512990171843`*^9}, { 3.5804513650374403`*^9, 3.580451370348692*^9}, {3.580451439656364*^9, 3.580451597295989*^9}, {3.580451627325652*^9, 3.580451639574376*^9}}], Cell[BoxData[ TemplateBox[{GraphicsBox[{{{}, {{}, {}, { RGBColor[0.24720000000000014`, 0.24, 0.6], LineBox[{{22., 0.7637010720591277}, {26., 0.7560917851299358}, {28., 0.7035159916943492}, {29., 0.8199069011682703}, {30., 0.7948433239849209}, {40., 0.6424135922546981}, {42., 0.6976713375483357}, {47., 0.7695457305383542}}], {{ LineBox[{{22., 1.2253597444056685`}, {22., 0.3020423997125869}}], LineBox[{ Offset[{1.5, 0}, {22., 1.2253597444056685`}], Offset[{-1.5, 0}, {22., 1.2253597444056685`}]}], LineBox[{ Offset[{1.5, 0}, {22., 0.3020423997125869}], Offset[{-1.5, 0}, {22., 0.3020423997125869}]}]}, { LineBox[{{26., 1.1444596038202446`}, {26., 0.36772396643962707`}}], LineBox[{ Offset[{1.5, 0}, {26., 1.1444596038202446`}], Offset[{-1.5, 0}, {26., 1.1444596038202446`}]}], LineBox[{ Offset[{1.5, 0}, {26., 0.36772396643962707`}], Offset[{-1.5, 0}, {26., 0.36772396643962707`}]}]}, { LineBox[{{28., 1.0626802549615288`}, {28., 0.34435172842716955`}}], LineBox[{ Offset[{1.5, 0}, {28., 1.0626802549615288`}], Offset[{-1.5, 0}, {28., 1.0626802549615288`}]}], LineBox[{ Offset[{1.5, 0}, {28., 0.34435172842716955`}], Offset[{-1.5, 0}, {28., 0.34435172842716955`}]}]}, { LineBox[{{29., 1.1678092851044646`}, {29., 0.47200451723207604`}}], LineBox[{ Offset[{1.5, 0}, {29., 1.1678092851044646`}], Offset[{-1.5, 0}, {29., 1.1678092851044646`}]}], LineBox[{ Offset[{1.5, 0}, {29., 0.47200451723207604`}], Offset[{-1.5, 0}, {29., 0.47200451723207604`}]}]}, { LineBox[{{30., 1.1305347800310712`}, {30., 0.45915186793877044`}}], LineBox[{ Offset[{1.5, 0}, {30., 1.1305347800310712`}], Offset[{-1.5, 0}, {30., 1.1305347800310712`}]}], LineBox[{ Offset[{1.5, 0}, {30., 0.45915186793877044`}], Offset[{-1.5, 0}, {30., 0.45915186793877044`}]}]}, { LineBox[{{40., 0.8915122404049848}, {40., 0.39331494410441137`}}], LineBox[{ Offset[{1.5, 0}, {40., 0.8915122404049848}], Offset[{-1.5, 0}, {40., 0.8915122404049848}]}], LineBox[{ Offset[{1.5, 0}, {40., 0.39331494410441137`}], Offset[{-1.5, 0}, {40., 0.39331494410441137`}]}]}, { LineBox[{{42., 0.9350410425049545}, {42., 0.46030163259171686`}}], LineBox[{ Offset[{1.5, 0}, {42., 0.9350410425049545}], Offset[{-1.5, 0}, {42., 0.9350410425049545}]}], LineBox[{ Offset[{1.5, 0}, {42., 0.46030163259171686`}], Offset[{-1.5, 0}, {42., 0.46030163259171686`}]}]}, { LineBox[{{47., 0.9816120345267103}, {47., 0.557479426549998}}], LineBox[{ Offset[{1.5, 0}, {47., 0.9816120345267103}], Offset[{-1.5, 0}, {47., 0.9816120345267103}]}], LineBox[{ Offset[{1.5, 0}, {47., 0.557479426549998}], Offset[{-1.5, 0}, {47., 0.557479426549998}]}]}}}, { RGBColor[0.6, 0.24, 0.4428931686004542], LineBox[{{26., 1.7653929994607047`}, {28., 1.6399366241947284`}, {29., 1.7147421834125431`}, {30., 1.673515822072929}, {40., 1.843973569812377}, {42., 1.9403019012286435`}, {47., 1.9816214000914343`}}], {{ LineBox[{{26., 2.106725492426047}, {26., 1.4240605064953624`}}], LineBox[{ Offset[{1.5, 0}, {26., 2.106725492426047}], Offset[{-1.5, 0}, {26., 2.106725492426047}]}], LineBox[{ Offset[{1.5, 0}, {26., 1.4240605064953624`}], Offset[{-1.5, 0}, {26., 1.4240605064953624`}]}]}, { LineBox[{{28., 1.9537469024781244`}, {28., 1.3261263459113324`}}], LineBox[{ Offset[{1.5, 0}, {28., 1.9537469024781244`}], Offset[{-1.5, 0}, {28., 1.9537469024781244`}]}], LineBox[{ Offset[{1.5, 0}, {28., 1.3261263459113324`}], Offset[{-1.5, 0}, {28., 1.3261263459113324`}]}]}, { LineBox[{{29., 2.0179117147079593`}, {29., 1.411572652117127}}], LineBox[{ Offset[{1.5, 0}, {29., 2.0179117147079593`}], Offset[{-1.5, 0}, {29., 2.0179117147079593`}]}], LineBox[{ Offset[{1.5, 0}, {29., 1.411572652117127}], Offset[{-1.5, 0}, {29., 1.411572652117127}]}]}, { LineBox[{{30., 1.9655414321546514`}, {30., 1.3814902119912065`}}], LineBox[{ Offset[{1.5, 0}, {30., 1.9655414321546514`}], Offset[{-1.5, 0}, {30., 1.9655414321546514`}]}], LineBox[{ Offset[{1.5, 0}, {30., 1.3814902119912065`}], Offset[{-1.5, 0}, {30., 1.3814902119912065`}]}]}, { LineBox[{{40., 2.060769174529352}, {40., 1.627177965095402}}], LineBox[{ Offset[{1.5, 0}, {40., 2.060769174529352}], Offset[{-1.5, 0}, {40., 2.060769174529352}]}], LineBox[{ Offset[{1.5, 0}, {40., 1.627177965095402}], Offset[{-1.5, 0}, {40., 1.627177965095402}]}]}, { LineBox[{{42., 2.146795190397568}, {42., 1.733808612059719}}], LineBox[{ Offset[{1.5, 0}, {42., 2.146795190397568}], Offset[{-1.5, 0}, {42., 2.146795190397568}]}], LineBox[{ Offset[{1.5, 0}, {42., 1.733808612059719}], Offset[{-1.5, 0}, {42., 1.733808612059719}]}]}, { LineBox[{{47., 2.2572441262133234`}, {47., 1.7059986739695452`}}], LineBox[{ Offset[{1.5, 0}, {47., 2.2572441262133234`}], Offset[{-1.5, 0}, {47., 2.2572441262133234`}]}], LineBox[{ Offset[{1.5, 0}, {47., 1.7059986739695452`}], Offset[{-1.5, 0}, {47., 1.7059986739695452`}]}]}}}}, {}}, {}}, AspectRatio -> NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> True, AxesLabel -> { FormBox["\"Z\"", TraditionalForm], FormBox["\"\[Sigma]\"", TraditionalForm]}, AxesOrigin -> {20., 0}, GridLines -> Automatic, GridLinesStyle -> Directive[ RGBColor[0.9, 0.9, 0.9]], ImageSize -> Full, Method -> {}, PlotRange -> {{20, 50}, {0, 2.5}}, PlotRangeClipping -> True, PlotRangePadding -> {Automatic, Automatic}], TemplateBox[{ "\"\\!\\(\\*SubscriptBox[\\(\[Sigma]\\), \\(\[Alpha]\\)]\\)\"", "\"\\!\\(\\*SubscriptBox[\\(\[Sigma]\\), \\(\[Beta]\\)]\\)\""}, "PointLegend", DisplayFunction -> (StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, 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 -> "Times"}, Background -> Automatic, StripOnInput -> False]& ), Editable -> True, InterpretationFunction :> (RowBox[{"PointLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.24720000000000014`", ",", "0.24`", ",", "0.6`"}], "]"}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.6`", ",", "0.24`", ",", "0.4428931686004542`"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}], ",", RowBox[{"LegendMarkers", "\[Rule]", "False"}]}], "}"}]}], "]"}]& )]}, "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.580451274615082*^9, 3.5804512995940533`*^9}, { 3.5804513543637457`*^9, 3.580451370823258*^9}, {3.580451482927791*^9, 3.580451597899193*^9}, {3.5804516296143723`*^9, 3.5804516400450287`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Style", "[", RowBox[{ RowBox[{"Grid", "[", RowBox[{ RowBox[{"Join", "[", RowBox[{ RowBox[{"{", RowBox[{"{", RowBox[{ "\"\\"", ",", "\"\<\!\(\*SubscriptBox[\(\[Sigma]\), \(\[Alpha]\)]\)\>\"", ",", "\"\<\!\(\*SubscriptBox[\(\[Sigma]\), \(\[Beta]\)]\)\>\""}], "}"}], "}"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"With", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Z", "=", RowBox[{"First", "[", RowBox[{ SubscriptBox["En", "\[Alpha]"], "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[Sigma]a", "=", RowBox[{"First", "[", RowBox[{"sa", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[CapitalDelta]\[Sigma]a", "=", RowBox[{"Last", "[", RowBox[{"sa", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ",", RowBox[{"\[Sigma]b", "=", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "2"}], ",", RowBox[{"First", "[", RowBox[{"sb", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], "]"}]}], "]"}]}], ",", RowBox[{"\[CapitalDelta]\[Sigma]b", "=", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "2"}], ",", RowBox[{"Last", "[", RowBox[{"sb", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], "]"}]}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"Z", ",", RowBox[{"PlusMinus", "[", RowBox[{ RowBox[{"NumberForm", "[", RowBox[{"\[Sigma]a", ",", RowBox[{"{", RowBox[{"2", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"NumberForm", "[", RowBox[{"\[CapitalDelta]\[Sigma]a", ",", RowBox[{"{", RowBox[{"2", ",", "2"}], "}"}]}], "]"}]}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "2"}], ",", RowBox[{"PlusMinus", "[", RowBox[{ RowBox[{"NumberForm", "[", RowBox[{"\[Sigma]b", ",", RowBox[{"{", RowBox[{"3", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"NumberForm", "[", RowBox[{"\[CapitalDelta]\[Sigma]b", ",", RowBox[{"{", RowBox[{"2", ",", "2"}], "}"}]}], "]"}]}], "]"}], ",", "\"\<\>\""}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"Length", "[", "sa", "]"}]}], "}"}]}], "]"}]}], "]"}], ",", RowBox[{"Alignment", "\[Rule]", RowBox[{"{", RowBox[{"Center", ",", "Automatic"}], "}"}]}], ",", RowBox[{"Dividers", "\[Rule]", RowBox[{"{", RowBox[{"None", ",", RowBox[{"{", RowBox[{"2", "\[Rule]", RowBox[{"GrayLevel", "[", "0.7", "]"}]}], "}"}]}], "}"}]}], ",", RowBox[{"Spacings", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"2", "\[Rule]", "2"}], ",", RowBox[{"3", "\[Rule]", "2"}]}], "}"}], ",", RowBox[{"{", RowBox[{"2", "\[Rule]", "0.75"}], "}"}]}], "}"}]}]}], "]"}], ",", "\"\\"", ",", RowBox[{"FontSize", "\[Rule]", "11"}]}], "]"}]], "Input", CellChangeTimes->{{3.580451764347727*^9, 3.580451997788866*^9}, { 3.580452112406464*^9, 3.58045219931326*^9}, {3.580452258781191*^9, 3.580452346443873*^9}, {3.580452385503578*^9, 3.5804525170541983`*^9}}], Cell[BoxData[ StyleBox[ TagBox[GridBox[{ {"\<\"Z\"\>", "\<\"\\!\\(\\*SubscriptBox[\\(\[Sigma]\\), \\(\[Alpha]\\)]\ \\)\"\>", "\<\"\\!\\(\\*SubscriptBox[\\(\[Sigma]\\), \\(\[Beta]\\)]\\)\"\>"}, {"22", RowBox[{ TagBox[ InterpretationBox["\<\"0.76\"\>", 0.7637010720591277, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.46\"\>", 0.46165867234654084`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], "\<\"\"\>"}, {"26", RowBox[{ TagBox[ InterpretationBox["\<\"0.76\"\>", 0.7560917851299358, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.39\"\>", 0.38836781869030873`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.77\"\>", 1.7653929994607047`, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.34\"\>", 0.34133249296534224`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"28", RowBox[{ TagBox[ InterpretationBox["\<\"0.70\"\>", 0.7035159916943492, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.36\"\>", 0.35916426326717965`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.64\"\>", 1.6399366241947284`, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.31\"\>", 0.31381027828339614`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"29", RowBox[{ TagBox[ InterpretationBox["\<\"0.82\"\>", 0.8199069011682703, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.35\"\>", 0.3479023839361942, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.71\"\>", 1.7147421834125431`, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.30\"\>", 0.30316953129541624`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"30", RowBox[{ TagBox[ InterpretationBox["\<\"0.79\"\>", 0.7948433239849209, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.34\"\>", 0.3356914560461504, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.67\"\>", 1.673515822072929, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.29\"\>", 0.29202561008172245`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"40", RowBox[{ TagBox[ InterpretationBox["\<\"0.64\"\>", 0.6424135922546981, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.25\"\>", 0.2490986481502867, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.84\"\>", 1.843973569812377, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.22\"\>", 0.21679560471697512`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"42", RowBox[{ TagBox[ InterpretationBox["\<\"0.70\"\>", 0.6976713375483357, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.24\"\>", 0.2373697049566188, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.94\"\>", 1.9403019012286435`, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.21\"\>", 0.2064932891689246, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]}, {"47", RowBox[{ TagBox[ InterpretationBox["\<\"0.77\"\>", 0.7695457305383542, AutoDelete->True], NumberForm[#, {2, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.21\"\>", 0.21206630398835613`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}], RowBox[{ TagBox[ InterpretationBox["\<\"1.98\"\>", 1.9816214000914343`, AutoDelete->True], NumberForm[#, {3, 2}]& ], "\[PlusMinus]", TagBox[ InterpretationBox["\<\"0.28\"\>", 0.27562272612188915`, AutoDelete->True], NumberForm[#, {2, 2}]& ]}]} }, AutoDelete->False, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Automatic}}}, GridBoxDividers->{ "Columns" -> {{None}}, "RowsIndexed" -> {2 -> GrayLevel[0.7]}}, GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings->{ "ColumnsIndexed" -> {2 -> 2, 3 -> 2}, "RowsIndexed" -> {2 -> 0.75}}], "Grid"], "DialogStyle", StripOnInput->False, FontSize->11]], "Output", CellChangeTimes->{{3.5804519242790527`*^9, 3.580451933182029*^9}, { 3.5804519746978703`*^9, 3.580451998269339*^9}, {3.5804521167841663`*^9, 3.580452199752586*^9}, 3.5804522982024803`*^9, 3.580452347002678*^9, { 3.5804524207370577`*^9, 3.5804524512010193`*^9}, {3.5804524854544563`*^9, 3.580452517308381*^9}}] }, Open ]] }, WindowSize->{740, 701}, WindowMargins->{{337, Automatic}, {14, Automatic}}, FrontEndVersion->"9.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 20, \ 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 147, 2, 28, "Input"], Cell[CellGroupData[{ Cell[729, 26, 2263, 70, 149, "Input"], Cell[2995, 98, 7440, 188, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10472, 291, 2054, 62, 133, "Input"], Cell[12529, 355, 6547, 165, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19113, 525, 197, 4, 28, "Input"], Cell[19313, 531, 405, 10, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19755, 546, 996, 26, 70, "Input"], Cell[20754, 574, 882, 21, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21673, 600, 898, 24, 71, "Input"], Cell[22574, 626, 852, 19, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[23463, 650, 1700, 46, 107, "Input"], Cell[25166, 698, 1031, 26, 63, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26234, 729, 1553, 44, 107, "Input"], Cell[27790, 775, 850, 22, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28677, 802, 1191, 33, 47, "Input"], Cell[29871, 837, 1648, 45, 80, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31556, 887, 1087, 31, 56, "Input"], Cell[32646, 920, 1465, 40, 80, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34148, 965, 1274, 31, 77, "Input"], Cell[35425, 998, 11261, 234, 410, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[46723, 1237, 4098, 102, 189, "Input"], Cell[50824, 1341, 6075, 193, 165, "Output"] }, Open ]] } ] *) (* End of internal cache information *)