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

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

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[    116967,       2608]
NotebookOptionsPosition[    107081,       2446]
NotebookOutlinePosition[    107853,       2471]
CellTagsIndexPosition[    107810,       2468]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[" FRW Models ", "Title",
 FontSize->24,ExpressionUUID->"6ab56dfe-2a2d-4adc-9284-94516474d8bf"],

Cell[TextData[{
 "This notebook numerically integrates the modified FRW equation \n\n",
 Cell[BoxData[
  FormBox[
   TemplateBox[<|"boxes" -> FormBox[
       RowBox[{
         FractionBox["1", "2"], 
         SuperscriptBox[
          RowBox[{"(", 
            FractionBox[
             RowBox[{
               StyleBox["d", "TI"], 
               OverscriptBox[
                StyleBox["a", "TI"], "~"]}], 
             RowBox[{
               StyleBox["d", "TI"], 
               OverscriptBox[
                StyleBox["t", "TI"], "~"]}]], ")"}], "2"], "+", 
         SubscriptBox[
          StyleBox["U", "TI"], 
          RowBox[{
            StyleBox["e", "TI"], 
            StyleBox["f", "TI"], 
            StyleBox["f", "TI"]}]], "\[LongEqual]", 
         FractionBox[
          SubscriptBox["\[CapitalOmega]", 
           StyleBox["k", "TI"]], "2"]}], TraditionalForm], "errors" -> {}, 
     "input" -> 
     "\\frac{1}{2} \\left( \\frac{d\\tilde{a}}{d\\tilde{t}} \\right)^2 \
+{\\mathcal U}_{eff} = \\frac{\\Omega_k}{2}", "state" -> "Boxes"|>,
    "TeXAssistantTemplate"], TraditionalForm]],ExpressionUUID->
  "822401ad-f808-4776-8779-34559545d832"],
 "\n\nwith\n\n",
 Cell[BoxData[
  FormBox[
   TemplateBox[<|"boxes" -> FormBox[
       RowBox[{
         SubscriptBox[
          StyleBox["U", "TI"], 
          RowBox[{
            StyleBox["e", "TI"], 
            StyleBox["f", "TI"], 
            StyleBox["f", "TI"]}]], "\[LongEqual]", "-", 
         FractionBox["1", "2"], 
         RowBox[{"(", 
           RowBox[{
             SubscriptBox["\[CapitalOmega]", "\[CapitalLambda]"], 
             SuperscriptBox[
              OverscriptBox[
               StyleBox["a", "TI"], "~"], "2"], "+", 
             FractionBox[
              SubscriptBox["\[CapitalOmega]", 
               StyleBox["m", "TI"]], 
              OverscriptBox[
               StyleBox["a", "TI"], "~"]], "+", 
             FractionBox[
              SubscriptBox["\[CapitalOmega]", 
               StyleBox["r", "TI"]], 
              SuperscriptBox[
               OverscriptBox[
                StyleBox["a", "TI"], "~"], "2"]]}], ")"}]}], TraditionalForm],
      "errors" -> {}, "input" -> 
     "U_{eff} = - \\frac{1}{2} \\left( \\Omega_\\Lambda \
\\tilde{a}^2+\\frac{\\Omega_m}{\\tilde{a}} + \\frac{\\Omega_r}{\\tilde{a}^2} \
\\right)", "state" -> "Boxes"|>,
    "TeXAssistantTemplate"], TraditionalForm]],ExpressionUUID->
  "9e56f951-e277-430e-abc6-abfa563bdd36"],
 "\n\nto compute a specific homogeneous isotropic cosmological model for the \
inputs (1) the Hubble constant ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["H", 
    RowBox[{"0", " "}]], TraditionalForm]],ExpressionUUID->
  "3639216a-3184-44a4-a337-c255b3882439"],
 ", and the three \[CapitalOmega]'s for radiation, matter, and cosmological \
constant. Note bene: The program assumes these parameters are such that the \
universe started with a big bang, which we know may not always be the case.  \
But it appears to be for our universe."
}], "Text",
 CellChangeTimes->{{3.734870189567423*^9, 3.73487022526422*^9}, {
  3.924003241961752*^9, 3.924003274114379*^9}, {3.924003389420176*^9, 
  3.924003396562475*^9}, {3.924003494391067*^9, 3.924003601293056*^9}, {
  3.986976005808714*^9, 
  3.9869760117279577`*^9}},ExpressionUUID->"b8391488-a377-44e7-b60f-\
e71e9d98e4cb"],

Cell[CellGroupData[{

Cell["Clearing the variables used:", "Subsection",ExpressionUUID->"e6aedb34-b49a-400c-af4e-68999a2f63d2"],

Cell["\<\
Clear all the variables that will be used in the calculation:\
\>", "Text",
 CellChangeTimes->{
  3.7348702385853167`*^9},ExpressionUUID->"e4806ddc-21e9-42f4-9dd9-\
6885efc3084b"],

Cell[BoxData[
 RowBox[{"Clear", "[", 
  RowBox[{
  "\[CapitalOmega]r", ",", "\[CapitalOmega]m", ",", 
   "\[CapitalOmega]\[CapitalLambda]", ",", "\[CapitalOmega]k", ",", "h", ",", 
   "Th", ",", " ", "x", ",", "y", ",", "a0", ",", "t0", ",", 
   "\[CapitalOmega]crit", ",", "yend", ",", "ymax"}], "]"}]], "Input",
 CellChangeTimes->{{3.544186480361396*^9, 3.544186501559106*^9}, {
  3.734870064257477*^9, 3.7348700965798798`*^9}, {3.9240037407156*^9, 
  3.924003741155175*^9}},
 CellLabel->"In[43]:=",ExpressionUUID->"23463587-9071-455e-9b3f-05e542a62f83"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Parameters of a FRW model:", "Subsection",ExpressionUUID->"deb930ed-d043-4f30-a67d-2d9a6b775af9"],

Cell[TextData[{
 "Four parameters specify a FRW cosmological model. First, there are the \
three",
 StyleBox[" \[CapitalOmega] '",
  FontSlant->"Italic"],
 "s :"
}], "Text",
 CellChangeTimes->{{3.544186507848159*^9, 3.544186514689322*^9}, 
   3.924003658411167*^9},ExpressionUUID->"ad8bab4c-b46f-41b2-a51a-\
fea95c48d9df"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalOmega]r", "=", 
  RowBox[{".1", "*", 
   RowBox[{"10", "^", 
    RowBox[{"(", 
     RowBox[{"-", "5"}], ")"}]}]}]}]], "Input",
 CellChangeTimes->{{3.543849613990481*^9, 3.543849615711028*^9}, 
   3.544186528259637*^9, {3.92408233265228*^9, 3.924082334957397*^9}, {
   3.924189709082595*^9, 3.924189709809947*^9}, {3.986975922404738*^9, 
   3.986975925899373*^9}},
 CellLabel->"In[44]:=",ExpressionUUID->"30cc047b-3a96-4273-a0eb-8c275e05154d"],

Cell[BoxData["1.0000000000000002`*^-6"], "Output",
 CellChangeTimes->{
  3.543849617066078*^9, 3.544185876242299*^9, 3.544186682357336*^9, 
   3.607337257405223*^9, 3.924004120789707*^9, 3.924004279312621*^9, 
   3.924077939874192*^9, {3.924078007073792*^9, 3.9240780329237843`*^9}, 
   3.924078087828835*^9, 3.9240781342817383`*^9, 3.924078165930067*^9, 
   3.9240782304957247`*^9, 3.92408215263313*^9, {3.924082336179925*^9, 
   3.9240823608142767`*^9}, {3.924189727281798*^9, 3.924189851712995*^9}, 
   3.986975961481249*^9, 3.986979671975518*^9, 3.9869797030168467`*^9},
 CellLabel->"Out[44]=",ExpressionUUID->"89ac8e4d-921f-4b33-aba2-bf6339dd19ef"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalOmega]m", "=", ".3"}]], "Input",
 CellChangeTimes->{{3.543849620690419*^9, 3.543849622321785*^9}, 
   3.544186528826647*^9, {3.5441866652498627`*^9, 3.544186669765218*^9}, {
   3.544186756170513*^9, 3.544186758873932*^9}, {3.924077923800068*^9, 
   3.924077923903597*^9}, {3.9240780260425053`*^9, 3.924078026779439*^9}, {
   3.924078154653441*^9, 3.924078156173443*^9}, {3.9240823395875273`*^9, 
   3.92408233996104*^9}, {3.924189715282794*^9, 3.924189847854536*^9}, {
   3.986975930354418*^9, 3.986975939477872*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"f252d9d9-6a33-4bc7-9c28-de6b3e4e97c8"],

Cell[BoxData["0.3`"], "Output",
 CellChangeTimes->{
  3.543849628090397*^9, 3.5441858775317907`*^9, 3.544186683355575*^9, 
   3.5441867598627977`*^9, 3.60733725919232*^9, 3.924004120802409*^9, 
   3.924004279319478*^9, 3.924077939879917*^9, {3.924078007091679*^9, 
   3.924078032929825*^9}, 3.924078087834594*^9, {3.9240781342875433`*^9, 
   3.92407816593591*^9}, 3.924078230501398*^9, 3.924082152639534*^9, 
   3.924082360830448*^9, {3.924189727287821*^9, 3.924189851719092*^9}, 
   3.9869759615044937`*^9, 3.9869796719806356`*^9, 3.98697970302071*^9},
 CellLabel->"Out[45]=",ExpressionUUID->"df38846c-2aab-452b-8345-a4504613c6bd"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[CapitalOmega]\[CapitalLambda]", "=", ".7"}]], "Input",
 CellChangeTimes->{
  3.544186529274458*^9, {3.5441866734054537`*^9, 3.5441866742612543`*^9}, 
   3.73487010520903*^9, {3.924077917726976*^9, 3.9240779189820423`*^9}, {
   3.924078076374337*^9, 3.924078077862118*^9}, {3.924082344566*^9, 
   3.924082349990781*^9}, {3.924189758273526*^9, 3.924189758881691*^9}, {
   3.986976028714679*^9, 3.9869760294020348`*^9}, {3.986979695294531*^9, 
   3.986979696258795*^9}, {3.986979726489603*^9, 
   3.9869797461843767`*^9}},ExpressionUUID->"d6a4e5a2-0971-402d-a52c-\
3d50751f6ea5"],

Cell[BoxData["0.2`"], "Output",
 CellChangeTimes->{
  3.543849629401236*^9, 3.5441858786010838`*^9, 3.5441866844101877`*^9, 
   3.6073372620372543`*^9, 3.9240041208176126`*^9, 3.9240042793374777`*^9, 
   3.924077939884907*^9, {3.924078007097678*^9, 3.924078032946207*^9}, 
   3.924078087840117*^9, 3.924078134295467*^9, 3.924078165954981*^9, 
   3.924078230507018*^9, 3.924082152654806*^9, 3.92408236083605*^9, {
   3.9241897273056726`*^9, 3.9241898517301893`*^9}, 3.986975961508101*^9, 
   3.986979671992787*^9, 3.986979703023638*^9},
 CellLabel->"Out[46]=",ExpressionUUID->"4f82c729-2b0b-4960-8cb1-a47b4f54b8f6"]
}, Open  ]],

Cell[TextData[{
 "The last parameter is the Hubble constant ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["H", "0"], TraditionalForm]],ExpressionUUID->
  "3c15b498-041c-4834-966f-cb701b70f3c1"],
 ", or equivalently the Hubble time ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["T", "0"], TraditionalForm]],ExpressionUUID->
  "80fce731-0403-4e78-924e-2c867d72fd1f"],
 "=1",
 StyleBox["/",
  FontSlant->"Italic"],
 Cell[BoxData[
  FormBox[
   SubscriptBox["H", "0"], TraditionalForm]],ExpressionUUID->
  "cccb818f-66fd-47f7-afe1-52f66a0b978f"],
 " which we denote here by ",
 StyleBox["Th. ",
  FontWeight->"Bold"],
 "  ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["T", "0"], TraditionalForm]],ExpressionUUID->
  "cd8c3419-009f-4377-bf4b-144d1677be1c"],
 " has the dimensions of time.  A billion years (Gyr)  is a convenient unit \
of time for cosmology and ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["T", "0"], TraditionalForm]],ExpressionUUID->
  "0ba0597f-ed50-4e22-8f1e-bac2e072f811"],
 StyleBox["=9.788 ",
  FontSlant->"Italic"],
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SuperscriptBox["h", 
      RowBox[{"-", "1"}]], " ", "Gyr"}], ",", " ", 
    RowBox[{
     RowBox[{"where", "  ", "h"}], "=", 
     SubscriptBox["H", "0"]}]}], TraditionalForm]],ExpressionUUID->
  "b2dc386f-08be-4cea-a277-6d168c67207b"],
 "/[100 (km/s)/Mpc]. Thus, for ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["H", "0"], TraditionalForm]],ExpressionUUID->
  "ce447692-55c2-4329-982a-ae0f17cf70ad"],
 "=67.3",
 StyleBox["(km/sec)/Mpc",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox[",  ",
  FontSlant->"Italic"]
}], "Text",
 CellChangeTimes->{{3.6072723847630587`*^9, 
  3.607272385964484*^9}},ExpressionUUID->"3ae14679-7710-4cb7-a0ba-\
12b2a2f3a7c3"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"h", "=", ".7"}]], "Input",
 CellChangeTimes->{{3.607272389925571*^9, 3.607272391677503*^9}, {
  3.924077997761505*^9, 3.92407802228087*^9}, {3.986975946768263*^9, 
  3.986975948074192*^9}},
 CellLabel->"In[47]:=",ExpressionUUID->"8ee56d3f-3d2f-4baa-af9a-f9f58d70cc79"],

Cell[BoxData["0.7`"], "Output",
 CellChangeTimes->{
  3.543849634174951*^9, 3.544185890239428*^9, 3.5441866889985228`*^9, 
   3.607337267016509*^9, 3.924004120823676*^9, 3.924004279347785*^9, 
   3.924077939889864*^9, {3.924078007102768*^9, 3.924078032952948*^9}, 
   3.9240780878568287`*^9, 3.9240781343108807`*^9, 3.924078165972435*^9, 
   3.9240782305241537`*^9, 3.9240821526602907`*^9, 3.924082360849195*^9, {
   3.92418972731434*^9, 3.92418985174833*^9}, 3.9869759615110607`*^9, 
   3.9869796720020447`*^9, 3.986979703029059*^9},
 CellLabel->"Out[47]=",ExpressionUUID->"dd55df7d-7e22-4b8c-9f4e-701a32803d6a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Th", "=", 
  RowBox[{"9.788", "/", "h"}]}]], "Input",
 CellLabel->"In[48]:=",ExpressionUUID->"c00bea6e-86a4-48c3-891f-9bfbaadb391c"],

Cell[BoxData["13.982857142857144`"], "Output",
 CellChangeTimes->{
  3.543849642240851*^9, 3.544185892953197*^9, 3.5441866911040783`*^9, 
   3.607337269794415*^9, 3.924004120841691*^9, 3.924004279365653*^9, 
   3.924077939895309*^9, {3.92407800711061*^9, 3.924078032971128*^9}, 
   3.9240780878635798`*^9, 3.924078134316928*^9, 3.924078165978362*^9, 
   3.924078230531205*^9, 3.924082152665372*^9, 3.92408236085483*^9, {
   3.9241897273287477`*^9, 3.924189851754624*^9}, 3.986975961514104*^9, 
   3.986979672017254*^9, 3.986979703040361*^9},
 CellLabel->"Out[48]=",ExpressionUUID->"3a70c617-6fd6-4078-9e7e-7b54a16e1be8"]
}, Open  ]],

Cell[TextData[{
 "We define ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "cfc0d154-f04e-410b-ab82-f363963aaf94"],
 " for curvature. It is related to the other ",
 StyleBox["\[CapitalOmega]",
  FontSlant->"Italic"],
 "'s by:"
}], "Text",
 CellChangeTimes->{
  3.5441865299978027`*^9, {3.924003701353318*^9, 
   3.924003724316926*^9}},ExpressionUUID->"24e96611-0854-4827-b7ea-\
4c4bf1f88fe4"],

Cell[BoxData[
 RowBox[{"\[CapitalOmega]k", ":=", 
  RowBox[{
  "1", "-", "\[CapitalOmega]r", "-", "\[CapitalOmega]m", "-", 
   "\[CapitalOmega]\[CapitalLambda]"}]}]], "Input",
 CellChangeTimes->{{3.544186530850996*^9, 3.54418653246659*^9}, 
   3.734870116642231*^9, {3.924003730160471*^9, 3.924003731009303*^9}},
 CellLabel->"In[49]:=",ExpressionUUID->"d13c5cb9-d8b4-45ab-aced-923155930d27"],

Cell["For the parameters above ", "Text",ExpressionUUID->"eb1643e8-f0cf-42d8-a1b8-6fb73e47c416"],

Cell[CellGroupData[{

Cell[BoxData["\[CapitalOmega]k"], "Input",
 CellChangeTimes->{
  3.544186533001226*^9, {3.924003733207151*^9, 3.924003733644986*^9}},
 CellLabel->"In[50]:=",ExpressionUUID->"2896bc8b-dfd3-4241-9fe7-f9285d66cba8"],

Cell[BoxData["0.49999899999999997`"], "Output",
 CellChangeTimes->{
  3.543849661006405*^9, 3.544186706406258*^9, 3.5441867669958267`*^9, 
   3.6073372766462793`*^9, 3.924004120860054*^9, 3.924004279375374*^9, 
   3.924077939915461*^9, {3.9240780071285963`*^9, 3.924078032990347*^9}, 
   3.92407808788342*^9, 3.924078134328587*^9, 3.9240781659984417`*^9, 
   3.924078230539895*^9, 3.924082152687841*^9, 3.924082360863113*^9, {
   3.924189727338809*^9, 3.9241898517779713`*^9}, 3.986975961537969*^9, 
   3.986979672022389*^9, 3.9869797030560102`*^9},
 CellLabel->"Out[50]=",ExpressionUUID->"0713fa54-e78a-4d80-8bca-1943608e6ee7"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrating the FRW equation in dimensionless variables:", "Subsection",ExpressionUUID->"9200126f-1431-4dca-b993-724b62b14593"],

Cell[TextData[{
 "Its convenient to rewrite the FRW equation in terms of  dimensionless \
variables as above.  Here it is convenient to write ",
 StyleBox["x",
  FontSlant->"Italic"],
 " for ",
 Cell[BoxData[
  FormBox[
   OverscriptBox["t", "~"], TraditionalForm]],ExpressionUUID->
  "422466bb-5ed6-4435-81c3-61c46ca237ec"],
 ", and",
 StyleBox[" y ",
  FontSlant->"Italic"],
 "for ",
 Cell[BoxData[
  FormBox[
   OverscriptBox["a", "~"], TraditionalForm]],ExpressionUUID->
  "86ccecdb-c6a9-4b14-a44a-e584e9753bcd"],
 ". The FRW equation is then the same as that of a non-relativistic particle \
moving in a potential U",
 StyleBox["(y) ",
  FontSlant->"Italic"],
 "where ",
 StyleBox["y ",
  FontSlant->"Italic"],
 "is the particle's position  and ",
 StyleBox["x ",
  FontSlant->"Italic"],
 "is the time. "
}], "Text",
 CellChangeTimes->{{3.924003762966796*^9, 
  3.924003803539091*^9}},ExpressionUUID->"f84c4137-99f8-49c0-b8a6-\
9e737ee38167"],

Cell[BoxData[
 RowBox[{
  RowBox[{"U", "[", "y_", "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"1", "/", "2"}], ")"}], 
   RowBox[{"(", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "\[CapitalOmega]r"}], "/", 
      RowBox[{"y", "^", "2"}]}], " ", "-", 
     RowBox[{"\[CapitalOmega]m", "/", "y"}], "-", 
     RowBox[{"\[CapitalOmega]\[CapitalLambda]", " ", 
      RowBox[{"y", "^", "2"}]}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.544186533633973*^9, 3.54418653503204*^9}, 
   3.734870125540658*^9},
 CellLabel->"In[51]:=",ExpressionUUID->"d1b3b8b5-568f-4ade-9bb3-a03a15e62fb4"],

Cell[TextData[{
 "That is, the FRW equation becomes:\n\n       ",
 Cell[BoxData[
  FormBox[
   SuperscriptBox[
    StyleBox[
     RowBox[{"(", 
      RowBox[{"dy", "/", "dx"}], ")"}],
     FontSlant->"Italic"], "2"], TraditionalForm]],ExpressionUUID->
  "c0b7b0c2-7ff2-469c-946d-b5305251183d"],
 "= ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "c31dbeae-300e-4467-a953-a95bdf4ff826"],
 "(",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "r"], TraditionalForm]],ExpressionUUID->
  "5c36519b-04dd-4d02-a3ad-8c760111db02"],
 ", ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    RowBox[{
     SubscriptBox["\[CapitalOmega]", "m"], ","}], " "}], TraditionalForm]],
  ExpressionUUID->"8d643f48-defb-4d27-b999-592c1a62046c"],
 " ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "\[CapitalLambda]"], TraditionalForm]],
  ExpressionUUID->"01487ce3-8f8d-490e-9d46-5d7e6cb86fce"],
 ")-2 U",
 StyleBox["(y; ",
  FontSlant->"Italic"],
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", 
    RowBox[{"r", " "}]], TraditionalForm]],ExpressionUUID->
  "f3ff7bfd-e9b6-464e-845d-40a391d91eff"],
 ", ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "m"], TraditionalForm]],ExpressionUUID->
  "beb87be3-2f82-40b4-a0a6-7005e0533df2"],
 ", ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "\[CapitalLambda]"], TraditionalForm]],
  ExpressionUUID->"2d73aaac-aa5a-4c68-b23a-329b4a9e2ef9"],
 ")  . "
}], "Text",
 CellChangeTimes->{{3.924003823638357*^9, 3.9240038332179527`*^9}, {
  3.9240041873162117`*^9, 
  3.924004193656588*^9}},ExpressionUUID->"63ef8e17-9bfd-48f6-b1e0-\
b89c5c223bba"],

Cell[TextData[{
 "We plot this ",
 StyleBox["U(y) ",
  FontSlant->"Italic"],
 "out to a value ",
 StyleBox["y=yend, ",
  FontSlant->"Italic"],
 "also showing a horizontal line at the value of ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "46217160-df75-47c3-a61b-5577058261e6"],
 "/2. (You might have to change ",
 StyleBox["yend ",
  FontWeight->"Bold"],
 StyleBox["and edit ",
  FontVariations->{"CompatibilityType"->0}],
 " ",
 StyleBox["PlotRange ",
  FontWeight->"Bold"],
 StyleBox["in the ",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["Show  ",
  FontWeight->"Bold",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["statement below to get an informative plot.) ",
  FontVariations->{"CompatibilityType"->0}]
}], "Text",
 CellChangeTimes->{{3.924003844283738*^9, 
  3.924003853337681*^9}},ExpressionUUID->"8b5360dc-577f-4c17-9a22-\
5ecaf4f3952d"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"yend", "=", "4"}]], "Input",
 CellChangeTimes->{{3.924078128357442*^9, 3.924078128901462*^9}, {
  3.92407822263478*^9, 3.9240782231774178`*^9}},
 CellLabel->"In[52]:=",ExpressionUUID->"995d8353-ca9f-4707-b6d5-4fdfc91f9b03"],

Cell[BoxData["4"], "Output",
 CellChangeTimes->{
  3.5438496833080263`*^9, 3.544185944364217*^9, 3.5441867764716797`*^9, 
   3.607337284371716*^9, 3.924004120877771*^9, 3.9240042793941708`*^9, 
   3.924077939927125*^9, {3.924078007136654*^9, 3.924078033011364*^9}, 
   3.924078087895585*^9, 3.924078134345646*^9, 3.924078166006835*^9, 
   3.9240782305605803`*^9, 3.924082152705525*^9, 3.9240823608843517`*^9, {
   3.9241897273639727`*^9, 3.924189851798718*^9}, 3.986975961586488*^9, 
   3.986979672037847*^9, 3.9869797030704803`*^9},
 CellLabel->"Out[52]=",ExpressionUUID->"e23a12d4-07b5-48c3-a217-5ca93a973c2d"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"plotu", ":=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"U", "[", "y", "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"y", ",", "0", ",", "yend"}], "}"}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\"\<y\>\"", ",", "\"\<U\>\""}], "}"}]}], ",", 
    RowBox[{"DisplayFunction", "\[Rule]", "Identity"}]}], "]"}]}]], "Input",
 CellLabel->"In[53]:=",ExpressionUUID->"1b57ef22-2ef1-4acb-ad18-7778d8d853a5"],

Cell[BoxData[
 RowBox[{"plotomc", ":=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]k", "/", "2"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"y", ",", "0.", ",", "yend"}], "}"}], ",", 
    RowBox[{"DisplayFunction", "->", "Identity"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.544186545352792*^9, {3.924003863771332*^9, 3.924003864515663*^9}},
 CellLabel->"In[54]:=",ExpressionUUID->"dd64f4dc-f8de-411d-9282-04f016c07a05"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plotu", ",", "plotomc", ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "yend"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", 
        RowBox[{"-", "3"}]}], "}"}]}], "}"}]}], ",", 
   RowBox[{"DisplayFunction", "->", "$DisplayFunction"}]}], "]"}]], "Input",
 CellLabel->"In[55]:=",ExpressionUUID->"8dc0139a-4316-4e32-a66e-fca498f5ca84"],

Cell[BoxData[
 GraphicsBox[{
   InterpretationBox[{
     TagBox[{{{}, {}, 
        TagBox[
         {RGBColor[0.24, 0.6, 0.8], AbsoluteThickness[2], Opacity[1.], 
          LineBox[CompressedData["
1:eJwV13k41N0XAHBLi4ShJIRS9KYkUpRwJluylCVJorRQJEpJpUjpTRGVpCwh
lHpDsm9nEpEwM/Z1MPYxw1gzlH7f31/f5/M89/ne7Zxz71U87WVzToCPj8+B
n4/v/9+dMvOjBfAJ0j8sCPstWUqxft7ufmhzOigkNpzpcllCUTXflCZbmQ7v
ClVeuf8SpAgJXhwadMsAnun2EzpPBSknCxY/45FM2P24eDCFJkCxWLOJrpxC
+BezkHpNgKJz+QD38UwmOH6jHrFXEKCsUYlQOxb1GebzlmlV+vFTBIKzLUsH
P8PY2ygjx638FG5P60Ul7Sw4oXJLOpnBR/kZrfiR25IFSd+HHv44xEe5K/Tl
nxCZbIjZul5xc+Qibj3w9GGWdzbIv72Q2COxiPXBXiMdldkwsTlXsvf5H1QS
VP2o5psD7hZOSX1vf+P3P8lqjfW5MKTpyxodm0cv3Xvhf1TywD/BQzg/ZB6l
b7lwN9/NA/otfYbPP/N4YU4+68aOfMgwL/+yz4OHIlNRu9eHFoAwH9d6UnwO
czSuRZn2FYD5kVEFScovdPa2/XVZpxDOCRvft7r8CzM4pILy4UKQYEmfrGyb
RZvhh7ruxsXgpWB/9j1lBhc2u8U9jy2GmoTv/U/vzGDyOePF4qliaFgf7Sms
P4MzTH4KKakEBoPfJVz+Oo3RXTcNcxZLYdd8b1RB6xTul3NIZtghLFc9eTMv
cQpZx7WXCn1CkH5wadeAxxQqsG+c+bSdAqpwpj1NcAq9Lmgp83VSQPNkiWMv
TGL5w5SmyOVfYXVhkNFZiUmUfi/5QEXzK+iMRjEa+iaQMjg5aPPoK/BXBSY4
PZ7AVWcz3qfsKQOZV9VSF/q56HpPwUHnbBkor3i3J66Yi0VJYSuo4WUg8l+R
XdgLLp7t9XCfGywDPclanbCDXMxx3qJq/uIbbM97RCb5jOMxh8QM7kQ5yFZ9
NnviwcHxjiPCSbIV8Kf4xV1bTQ4+cBY6Z2tYAfp3M7snFtj45YyXbG5kBWz9
/e9YajgbxS7pBd/S/g5uEC1ylTKKFUFtx5fdrgSRba8fihmx8MTSsJy8lEro
jhmiuq9l4eS/ZPELdZXA8G9NcBkdwQ1h7yp+rq+CbsqL7FUvRtD/5TX1p2VV
cL9NoFljfBh3/iexTE6oGt4EO1ZLFw/hD7WKU7Xq1XAp7pZAy8shPPXZr+iO
QzU0rHRQXn51CMNyu717P1QD2/EfiZztQzhM+dSRavkTrstaXpJIHcT4poOf
NZ7XQEi57ldIG8Ddx/6s7CuqgWfSq0fJjwewpj3TNbK/BhiXVt6N9hxAXvda
uV+7a2HZZV97Y80BPMIaeFDcWgsqaZ5K3mX9uPLv3RMHFKgg72E0NT3Wh9EX
7J1VTajQq54Ssr+5D5UaVU9JXKKCxqBqZl1JH+qntZzpKKFCaaXFSGpYHxak
ReXkqdHAR2rd593qfcixZReb19HA0CihkxTIxF0rWvT0mTQQGhkO/uPOxJul
X0vVZ2lwr+Loa4WjTFyu8pKyRoEOV9TM9edUmai4uL+825MOCaki+bOdvWj3
/uVPH9F6kK9csI407sVYpyALV8V6cB00+eih2YvMVZ61x3bXw1MHq7SLir3o
dduAqudUDx+dWP4jf3owxIZTv/xTPQxw9E5R8nuw5LdBe4xFAwTauPT47e5B
Zeux4W+PG2FcQeX6gQgGygYc0HFLbIRn8nkBh64wUPxTwmPhvEZg/t4dHH2E
gfNCtmrWzEawlrA7NinDQBolz4expwm6kqfD3r3rwpvqdxfnBppgb2nfh2t1
nejl3H44dqEJvoQnXrbN7sSzoZqJINEMzw7MtJx/3YmHhwcMH+g2w5zfmPAe
t05USjALWf28GSooXwztBTuxjiQpqUZuAUkvGwVZ4w78pud5jm7XAqTrly+r
qnVgvsf33KseLVBGCwryXNuBbyv9HIqiWuC17ir2DVY7+gV0xZtyWuBi+ov4
9GftuHE8VeXM61aIOu65ZNNoG16v2wvR020gu6/S93RBK4qp5Qz8EGqHOcs/
w2nvWzE5TD10Qa4d4M695aToVqRb/NPmbNwOv90fG2pcb8VtPyV9lKPa4Z6V
c0GdVit2V46nZml3gI2+tcyeohY88DVFrPZmJzBdijfrtjQjY4NizuKTTjhh
MJNS+aMZrwbGOqq/7QSu++rxC8XNmAjP3z//2QmnDk2oTCY240JJoMFxuS7I
nDJ1F7nUjBmFjr5DJV2QkJvTYyDcjGuzJRj8/N1w5/zyvGTrJjwRunCkX6ob
zCVjFenGTZh0duDnd9VuUKzRTlyl04SqawoKHh/rBnan577GjU1o4HvqxZrM
bpDZxO9nOdOIWrH2qZE7emBc3cJ+Nq4RXZxTnd4zekDsk5XTlj8NeN9s3atX
Iz1AkxmJz5hqwFStiMZH0z3gWFceeYDVgGyxm+aewr0wf2btlrSWBvRFy707
tXrh/AO1Hzu+NGDohhnJ4rBemKneR9np2YC5TMMaqi4TNFUN/BzZ9dhsenci
7AATPon59Q731+NseqmUhQ0TTpeM2N/rqkftW3tdqtyYUP9FfV13XT3mr94x
S3nGhP7t52t2ZNVjgZHshqxhJghUbv+ldKsei1O5PpGRffDsAl24Z009doqo
vbJ50wd+M+3TLFI9/r7iUSr+oQ9U+75n8QnXox4MCoVhHyRtE/C0WKRjaUtH
XDCrD5r/uO4+M0RHFKqs9CX3Q/bbcuu7RXQsc49b58DuJ+Lpe/F6Dzp+GfHw
cJ7tBwPmVFChKx2Tz+sUneEbAKlMBT2n03QMdm05fklyABLatdwKHehoelri
9T3dASjfNfao3pSOVIdg6fTQAdieThW9soWOlFbb81+iBmBadk3pv0p0zLTf
mJ+fMACC11zpyRvo+MwO7b9lD0AxY9ffaWk6HrXmRbV2DoDkxPow9go6dpp6
rhHcPgjiEVQrDQ4Nayv3nRPSHgSP45ZHC0ZoWGoinCO6fxDOfKnSMx2k4Ruj
d0ek7QYhqNI98mY3DU+Te59vvz0It59rhqvX03BY227VsdpBUG9NPdGRT8OZ
zXqkT55DEBikO7XwiIb6TIPKR35DsHGLnYnMQxo+iDMNOH9vCPRTXmXoBdNQ
SvLI+KZXQ2AWkpISE0BDLX4P6uvyITh3nnHw41Ua+nVEhz9cNwzxTlncVado
SImKP+C6eRiUBdUU451oKGST/NdQYxgWQ4QmdzjS8FVVhtdfk2GYPKbu4HaU
hoU53w/7XhmGhw5H7yhZ0vB3+DTp7I9hsDCPtK3fR0Mj8/nK/Y3DUK8oeK58
Lw1Dl/EFru8ehrUBYZ+KtWkof3slt316GJ7Ea/xETWI+7htp1htGwKwlP4yk
SsMAI6sIuD4CO3v556nyNKz6a2cqHzQC4/kC/22Uo6F4kSPfQugI1Fyr0Lol
S8MEDTfv3KQREDhgR9FbS8xH4Y7V9roRGJZbKSwkQUMB3kfxdcosmA9T/PHf
UhquqzWYgh0sCJyrqvVeQsNdiW1NZ/eyQH6HceJeQRqeO7g8Jt2SBdVfvZLb
+Ij+X51WNvBlQerBj3Ypc1Ts9eQtdwtkQenX2+43Zqk4vz+C9fgRC8o4wl1H
p6moyirJaIpnwZ1Yo+LtXCqG68joXKhkgfakTajZMBXfi2XKPaGzIM6nPO/8
IBUpTJO/WR0sKPRfFIzop+Lko6vlv8dZkOyi2svroeLRDuqhCOlRaFW507u3
nYpeGa4aORtHgWUyZZzQSsWH9/6sblcdBSkb0lNSCxULtm1r27R/FIQOzJUK
NVJRzv/BmbwLo6Aws5Bzto6Ku63kTTp9RuHJ35+PJWqpeEgpewv/nVF4Rq8S
r/xJxYCaXo7Zs1Ho57cUtfhBxegEP5pX7CgkGXpdX1dFxcyrpC+RqaNwa3zq
2tR3KjLl9PwYhaMwS69YU1JOxQVuw3HBCmK83xbS0r9RUbLCXW8LdRTUfc7l
vy+joolntMCVvlHg+yO6NZdCxbSSqdClImzgc7NN311CxbKnj7y2SrFBzuK+
i3sxFdvPKdoc3sAGt32WJ94VUVFE7PDaV7vYEFGyVFynkIrKzAFeiT4brmQW
FocXUFE/17+TacqGastn0Zx8Kno7f0hUdWLDnv2RFV/zqBiyc/99azc2ZF39
IbyHcNKyVlffy2zIGdNyy82lYmH7pYMxt9jAcPjVqEu4IX2pKiWYDbPfth/6
mUNFdlCs2EA4G6TDl9a6EF5irzmx4jUbhk9nmfMRlt9W3aCWzIZRee/y1Gwq
av09lWubzoYllUEadoQtLClBxvlsmPcLYYsQPhOz/rB2GRv+uV+yv+YLFW+O
3FmnUsOG6HN1254TjtBmDMk2E+sxfiXFhXBqsF62SA8bknQupGkTLm6IDVwc
YUN+nO/ONYTrFX9bcKfYIBXntm8+i4ojXo4yzD9scOgh4SBhvtLCgYblHChV
NstvJywlIptVIcGB5Vpzis2EVY/fuJO3jgO+rQK/WgkbvG81S1PmQMt3h519
hI/Naq+N2cGB3MSpumnCXkYv+0L3ciCqv7ZSlOj/wbPZjDuGHDjUzl6rRji2
x87f25IDQU8dC+0IZ6nlmJ6258AG+rYP9wlX+UuuOeLCAWFRX2YhYUa1T6+x
Bweqjzs4zxGekW74pH2NAy6uC9J6xHqJuO28qRLAAVW3h6IhhDfmPDVZF8IB
hqSmdifhQ9ZW3YuxHBCKqpGKJvbj7JuMj9xUDvw3R/7GR+znLY6YHzOTA0oS
jKfehN+H1Ip/L+eACV9f0DkiHkpbVLvy6jgwXSgQyyLcqByaltZKzJdRWO5L
xBN/mZlBGJsDLzvf7Uwm4k1a/INYwCwHNN5WXDIm4lHNeUWHN98Y5L60fM8h
fHy+0ueI5Bj4tWXw2xLx7G36D9lEYQwetSVJixPx/m/UA5E9W8Ygw0FjXSPh
7J3GKet0x+D0g2d0b6RidWDyZVGTMej0+c//IJEvPXWC+n8Pj4GDFZdf5SsV
RT3KmplnxsDtd1nIbyK/XN+C0IfHY3Co8MLaGSI//SfiG2NejIFZVaDcEiJ/
n8NiQtgb4n8v2G0yRH5TOor3Xv4yBpvSwNOByH8ZSR2PPR1j0Kjvfm8LjYo1
93fVfd86DkoHk0YFOqioQ/op8HT3OATNnmPGdBLr98pF25FM+KrwdV0GFe+l
P0kYtxsH26J7rpG9RPvW4SvSd8fB1SGsJf7/9W9b3Fr3lnGI4ZKYMkQ9lcrV
tNjFHAeWkcTtrfNUvE+uDlxkj8PRv+s8DH4T9ebDwtWD/FywD9u71oiox49x
C6NLiguTH/U9coj6bTASlLnMgAuOjGlvIzEa3lCrq/Y4wIWLK0WsH5BomOkj
M0Cz4EJIsHH+T3EaKixmyMTYc8Hcb2S1y2oa8lZ1Be3w5MI/i6sqc6SJ9rra
dseiuQAd952XKtFwfcQoL22MC2E2i3079Ijz3CUgLWuKC6FSS8e+6tPQWXO1
Q9EcFxQkA5qOkmkY26yTXyMwATUCG7RDDGkoLf/Id1xqAmz2FsqLm9Fw1Yct
U7tgAux6mg1S7Gm4tPwcG8Mn4Ntficg+4nxOfsGLqYqcAJ0M0491vsR56hZm
Tn81AXzuz68V+9EwSDjnA/PtBARp3lZI8KfhotWSC0vzJ0A08NWbf+/TcK4r
acCsZwISnNRsmyNpOPqrm9GkPgnXn4i4uefRMNctvMxMaxJ0c+8KthfQMLBF
PxX3TYKNRX6CRTFxf8iL80wzmQQW63eLzldifX0df/ufmASrw0aS+37S8PV0
s4zSw0lI4uOL/UHcX0wnao9c6Z4Evv/meWtWEvcxVlG12JMpED0Yt1XFnY6X
Xh7e5P58CgLuFF0le9JR26jvVkX0FHRuXTh+3JuO1XHC2/3fToHW3mVzL3zp
OGHlED6SPwXZfdEOqvfoqJ83a1veNwUaPP/hxlg6tt3T6Lqxdxqko5pbxOrp
KCb/bry/fxoCblVvmDWsx9ALOa4JSrNwZ3zNwXqNBpwlf+w/H/ALJDY6WG20
b8SrO8uzpb7NAYd3VGfXwybsTi9/TFech7p3xQOaU8R7wL0xRuXmAoR+ZgwI
///9sX/5omvRb9hvvNV8ZqENDyQpun1euwg9dq9XHSztwOvyJzf2nvgLrCA5
Es+/C1XZditbLPnIvNjYQq5ZN64mcYXK5/nI7LhlLoafenAJxmRsiuUnm5n/
cEjM6EWPXUVtK8wFyP5CBeGURiZ+PHg6oXRWgCwp9/NTDH8/zozKig9GCZIV
chKbvysPYI7tAcFIoyVkc568o4zXIAaYBr/JHlxCTjQbcN+YPoSKXWatb4OX
kqsjLtp+nhhGHe37kTI7l5GNXj1MqDvBwhMKVB5/wzLyx9AT+rTCUZypDmN6
311OpisWdrsacDB6ifNVGyUhsrGVw/2LRWO4mpGTJFAuRNa5dmZ30zouLr19
98jeyyvI3kpHhubGuciX5pzBXCVMLrq27FhA0QReeZ0RtqpMmGxI/0ISj5/E
km0PhSQurSSznB+4RQRNoYTZFdUgMRFyyIbkCQWnaezcTmu4USxClq0xkl44
PIOZ8aDU7yxKtq/bytqzZxadMrQk5BdFyfU+t1eJKv7CwODEu50fxMjLePvb
3TfN4adh+RgxCxI58PAS7SEZHj6MCz7FOEQib3zsmgIKPDxrw1FOtyaR1SIC
b0Vv5KFccXGmpT2JnKLb5my+jYePwh2/h54mkePWCgvk6vLwvNarCeEbJLKP
a31B/Ekebrwvabo8lUQuCznEKknh4eIef9GW9yTyo5snDm36wMM2Tl996kcS
WaxDRDYknYcR9l+cTD6TyKOHbYbt8nj4d6uNT3AxiSze7HRquoqHnfTweMEG
Epl5e7DEiM3D/Ae/zjQ0kcjO6bUxmVweRu47qfK2lUR+vWDRKzfDQ/MUtWwD
Bol876TE29k/PCz0q/1xd4RE3hX4hj+DNI9R23eHW7NJZMWq3XIKkvN4mRl7
RHGcGI/Myfgw6XncYnGxmzJNIpOtzrd5Ks7jEv7G5IhfJLLhtMkzhvI89uTs
cz81TyKXyuVmHd46j8Xub3eo/yGR5VRydL6qzePL9Stn/v4lkfluGG7V1JzH
/wFtmZQq
           "]]},
         Annotation[#, "Charting`Private`Tag#1"]& ]}, {}},
      {"WolframDynamicHighlight", <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>}], 
     StyleBox[
      DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {}, 
        Slot["HighlightElements"], 
        Slot["LayoutOptions"], 
        Slot["Meta"], 
        Charting`HighlightActionFunction["DynamicHighlight", {{{{}, {}, 
            Annotation[{
              Directive[
               Opacity[1.], 
               RGBColor[0.24, 0.6, 0.8], 
               AbsoluteThickness[2]], 
              Line[CompressedData["
1:eJwV13k41N0XAHBLi4ShJIRS9KYkUpRwJluylCVJorRQJEpJpUjpTRGVpCwh
lHpDsm9nEpEwM/Z1MPYxw1gzlH7f31/f5/M89/ne7Zxz71U87WVzToCPj8+B
n4/v/9+dMvOjBfAJ0j8sCPstWUqxft7ufmhzOigkNpzpcllCUTXflCZbmQ7v
ClVeuf8SpAgJXhwadMsAnun2EzpPBSknCxY/45FM2P24eDCFJkCxWLOJrpxC
+BezkHpNgKJz+QD38UwmOH6jHrFXEKCsUYlQOxb1GebzlmlV+vFTBIKzLUsH
P8PY2ygjx638FG5P60Ul7Sw4oXJLOpnBR/kZrfiR25IFSd+HHv44xEe5K/Tl
nxCZbIjZul5xc+Qibj3w9GGWdzbIv72Q2COxiPXBXiMdldkwsTlXsvf5H1QS
VP2o5psD7hZOSX1vf+P3P8lqjfW5MKTpyxodm0cv3Xvhf1TywD/BQzg/ZB6l
b7lwN9/NA/otfYbPP/N4YU4+68aOfMgwL/+yz4OHIlNRu9eHFoAwH9d6UnwO
czSuRZn2FYD5kVEFScovdPa2/XVZpxDOCRvft7r8CzM4pILy4UKQYEmfrGyb
RZvhh7ruxsXgpWB/9j1lBhc2u8U9jy2GmoTv/U/vzGDyOePF4qliaFgf7Sms
P4MzTH4KKakEBoPfJVz+Oo3RXTcNcxZLYdd8b1RB6xTul3NIZtghLFc9eTMv
cQpZx7WXCn1CkH5wadeAxxQqsG+c+bSdAqpwpj1NcAq9Lmgp83VSQPNkiWMv
TGL5w5SmyOVfYXVhkNFZiUmUfi/5QEXzK+iMRjEa+iaQMjg5aPPoK/BXBSY4
PZ7AVWcz3qfsKQOZV9VSF/q56HpPwUHnbBkor3i3J66Yi0VJYSuo4WUg8l+R
XdgLLp7t9XCfGywDPclanbCDXMxx3qJq/uIbbM97RCb5jOMxh8QM7kQ5yFZ9
NnviwcHxjiPCSbIV8Kf4xV1bTQ4+cBY6Z2tYAfp3M7snFtj45YyXbG5kBWz9
/e9YajgbxS7pBd/S/g5uEC1ylTKKFUFtx5fdrgSRba8fihmx8MTSsJy8lEro
jhmiuq9l4eS/ZPELdZXA8G9NcBkdwQ1h7yp+rq+CbsqL7FUvRtD/5TX1p2VV
cL9NoFljfBh3/iexTE6oGt4EO1ZLFw/hD7WKU7Xq1XAp7pZAy8shPPXZr+iO
QzU0rHRQXn51CMNyu717P1QD2/EfiZztQzhM+dSRavkTrstaXpJIHcT4poOf
NZ7XQEi57ldIG8Ddx/6s7CuqgWfSq0fJjwewpj3TNbK/BhiXVt6N9hxAXvda
uV+7a2HZZV97Y80BPMIaeFDcWgsqaZ5K3mX9uPLv3RMHFKgg72E0NT3Wh9EX
7J1VTajQq54Ssr+5D5UaVU9JXKKCxqBqZl1JH+qntZzpKKFCaaXFSGpYHxak
ReXkqdHAR2rd593qfcixZReb19HA0CihkxTIxF0rWvT0mTQQGhkO/uPOxJul
X0vVZ2lwr+Loa4WjTFyu8pKyRoEOV9TM9edUmai4uL+825MOCaki+bOdvWj3
/uVPH9F6kK9csI407sVYpyALV8V6cB00+eih2YvMVZ61x3bXw1MHq7SLir3o
dduAqudUDx+dWP4jf3owxIZTv/xTPQxw9E5R8nuw5LdBe4xFAwTauPT47e5B
Zeux4W+PG2FcQeX6gQgGygYc0HFLbIRn8nkBh64wUPxTwmPhvEZg/t4dHH2E
gfNCtmrWzEawlrA7NinDQBolz4expwm6kqfD3r3rwpvqdxfnBppgb2nfh2t1
nejl3H44dqEJvoQnXrbN7sSzoZqJINEMzw7MtJx/3YmHhwcMH+g2w5zfmPAe
t05USjALWf28GSooXwztBTuxjiQpqUZuAUkvGwVZ4w78pud5jm7XAqTrly+r
qnVgvsf33KseLVBGCwryXNuBbyv9HIqiWuC17ir2DVY7+gV0xZtyWuBi+ov4
9GftuHE8VeXM61aIOu65ZNNoG16v2wvR020gu6/S93RBK4qp5Qz8EGqHOcs/
w2nvWzE5TD10Qa4d4M695aToVqRb/NPmbNwOv90fG2pcb8VtPyV9lKPa4Z6V
c0GdVit2V46nZml3gI2+tcyeohY88DVFrPZmJzBdijfrtjQjY4NizuKTTjhh
MJNS+aMZrwbGOqq/7QSu++rxC8XNmAjP3z//2QmnDk2oTCY240JJoMFxuS7I
nDJ1F7nUjBmFjr5DJV2QkJvTYyDcjGuzJRj8/N1w5/zyvGTrJjwRunCkX6ob
zCVjFenGTZh0duDnd9VuUKzRTlyl04SqawoKHh/rBnan577GjU1o4HvqxZrM
bpDZxO9nOdOIWrH2qZE7emBc3cJ+Nq4RXZxTnd4zekDsk5XTlj8NeN9s3atX
Iz1AkxmJz5hqwFStiMZH0z3gWFceeYDVgGyxm+aewr0wf2btlrSWBvRFy707
tXrh/AO1Hzu+NGDohhnJ4rBemKneR9np2YC5TMMaqi4TNFUN/BzZ9dhsenci
7AATPon59Q731+NseqmUhQ0TTpeM2N/rqkftW3tdqtyYUP9FfV13XT3mr94x
S3nGhP7t52t2ZNVjgZHshqxhJghUbv+ldKsei1O5PpGRffDsAl24Z009doqo
vbJ50wd+M+3TLFI9/r7iUSr+oQ9U+75n8QnXox4MCoVhHyRtE/C0WKRjaUtH
XDCrD5r/uO4+M0RHFKqs9CX3Q/bbcuu7RXQsc49b58DuJ+Lpe/F6Dzp+GfHw
cJ7tBwPmVFChKx2Tz+sUneEbAKlMBT2n03QMdm05fklyABLatdwKHehoelri
9T3dASjfNfao3pSOVIdg6fTQAdieThW9soWOlFbb81+iBmBadk3pv0p0zLTf
mJ+fMACC11zpyRvo+MwO7b9lD0AxY9ffaWk6HrXmRbV2DoDkxPow9go6dpp6
rhHcPgjiEVQrDQ4Nayv3nRPSHgSP45ZHC0ZoWGoinCO6fxDOfKnSMx2k4Ruj
d0ek7QYhqNI98mY3DU+Te59vvz0It59rhqvX03BY227VsdpBUG9NPdGRT8OZ
zXqkT55DEBikO7XwiIb6TIPKR35DsHGLnYnMQxo+iDMNOH9vCPRTXmXoBdNQ
SvLI+KZXQ2AWkpISE0BDLX4P6uvyITh3nnHw41Ua+nVEhz9cNwzxTlncVado
SImKP+C6eRiUBdUU451oKGST/NdQYxgWQ4QmdzjS8FVVhtdfk2GYPKbu4HaU
hoU53w/7XhmGhw5H7yhZ0vB3+DTp7I9hsDCPtK3fR0Mj8/nK/Y3DUK8oeK58
Lw1Dl/EFru8ehrUBYZ+KtWkof3slt316GJ7Ea/xETWI+7htp1htGwKwlP4yk
SsMAI6sIuD4CO3v556nyNKz6a2cqHzQC4/kC/22Uo6F4kSPfQugI1Fyr0Lol
S8MEDTfv3KQREDhgR9FbS8xH4Y7V9roRGJZbKSwkQUMB3kfxdcosmA9T/PHf
UhquqzWYgh0sCJyrqvVeQsNdiW1NZ/eyQH6HceJeQRqeO7g8Jt2SBdVfvZLb
+Ij+X51WNvBlQerBj3Ypc1Ts9eQtdwtkQenX2+43Zqk4vz+C9fgRC8o4wl1H
p6moyirJaIpnwZ1Yo+LtXCqG68joXKhkgfakTajZMBXfi2XKPaGzIM6nPO/8
IBUpTJO/WR0sKPRfFIzop+Lko6vlv8dZkOyi2svroeLRDuqhCOlRaFW507u3
nYpeGa4aORtHgWUyZZzQSsWH9/6sblcdBSkb0lNSCxULtm1r27R/FIQOzJUK
NVJRzv/BmbwLo6Aws5Bzto6Ku63kTTp9RuHJ35+PJWqpeEgpewv/nVF4Rq8S
r/xJxYCaXo7Zs1Ho57cUtfhBxegEP5pX7CgkGXpdX1dFxcyrpC+RqaNwa3zq
2tR3KjLl9PwYhaMwS69YU1JOxQVuw3HBCmK83xbS0r9RUbLCXW8LdRTUfc7l
vy+joolntMCVvlHg+yO6NZdCxbSSqdClImzgc7NN311CxbKnj7y2SrFBzuK+
i3sxFdvPKdoc3sAGt32WJ94VUVFE7PDaV7vYEFGyVFynkIrKzAFeiT4brmQW
FocXUFE/17+TacqGastn0Zx8Kno7f0hUdWLDnv2RFV/zqBiyc/99azc2ZF39
IbyHcNKyVlffy2zIGdNyy82lYmH7pYMxt9jAcPjVqEu4IX2pKiWYDbPfth/6
mUNFdlCs2EA4G6TDl9a6EF5irzmx4jUbhk9nmfMRlt9W3aCWzIZRee/y1Gwq
av09lWubzoYllUEadoQtLClBxvlsmPcLYYsQPhOz/rB2GRv+uV+yv+YLFW+O
3FmnUsOG6HN1254TjtBmDMk2E+sxfiXFhXBqsF62SA8bknQupGkTLm6IDVwc
YUN+nO/ONYTrFX9bcKfYIBXntm8+i4ojXo4yzD9scOgh4SBhvtLCgYblHChV
NstvJywlIptVIcGB5Vpzis2EVY/fuJO3jgO+rQK/WgkbvG81S1PmQMt3h519
hI/Naq+N2cGB3MSpumnCXkYv+0L3ciCqv7ZSlOj/wbPZjDuGHDjUzl6rRji2
x87f25IDQU8dC+0IZ6nlmJ6258AG+rYP9wlX+UuuOeLCAWFRX2YhYUa1T6+x
Bweqjzs4zxGekW74pH2NAy6uC9J6xHqJuO28qRLAAVW3h6IhhDfmPDVZF8IB
hqSmdifhQ9ZW3YuxHBCKqpGKJvbj7JuMj9xUDvw3R/7GR+znLY6YHzOTA0oS
jKfehN+H1Ip/L+eACV9f0DkiHkpbVLvy6jgwXSgQyyLcqByaltZKzJdRWO5L
xBN/mZlBGJsDLzvf7Uwm4k1a/INYwCwHNN5WXDIm4lHNeUWHN98Y5L60fM8h
fHy+0ueI5Bj4tWXw2xLx7G36D9lEYQwetSVJixPx/m/UA5E9W8Ygw0FjXSPh
7J3GKet0x+D0g2d0b6RidWDyZVGTMej0+c//IJEvPXWC+n8Pj4GDFZdf5SsV
RT3KmplnxsDtd1nIbyK/XN+C0IfHY3Co8MLaGSI//SfiG2NejIFZVaDcEiJ/
n8NiQtgb4n8v2G0yRH5TOor3Xv4yBpvSwNOByH8ZSR2PPR1j0Kjvfm8LjYo1
93fVfd86DkoHk0YFOqioQ/op8HT3OATNnmPGdBLr98pF25FM+KrwdV0GFe+l
P0kYtxsH26J7rpG9RPvW4SvSd8fB1SGsJf7/9W9b3Fr3lnGI4ZKYMkQ9lcrV
tNjFHAeWkcTtrfNUvE+uDlxkj8PRv+s8DH4T9ebDwtWD/FywD9u71oiox49x
C6NLiguTH/U9coj6bTASlLnMgAuOjGlvIzEa3lCrq/Y4wIWLK0WsH5BomOkj
M0Cz4EJIsHH+T3EaKixmyMTYc8Hcb2S1y2oa8lZ1Be3w5MI/i6sqc6SJ9rra
dseiuQAd952XKtFwfcQoL22MC2E2i3079Ijz3CUgLWuKC6FSS8e+6tPQWXO1
Q9EcFxQkA5qOkmkY26yTXyMwATUCG7RDDGkoLf/Id1xqAmz2FsqLm9Fw1Yct
U7tgAux6mg1S7Gm4tPwcG8Mn4Ntficg+4nxOfsGLqYqcAJ0M0491vsR56hZm
Tn81AXzuz68V+9EwSDjnA/PtBARp3lZI8KfhotWSC0vzJ0A08NWbf+/TcK4r
acCsZwISnNRsmyNpOPqrm9GkPgnXn4i4uefRMNctvMxMaxJ0c+8KthfQMLBF
PxX3TYKNRX6CRTFxf8iL80wzmQQW63eLzldifX0df/ufmASrw0aS+37S8PV0
s4zSw0lI4uOL/UHcX0wnao9c6Z4Evv/meWtWEvcxVlG12JMpED0Yt1XFnY6X
Xh7e5P58CgLuFF0le9JR26jvVkX0FHRuXTh+3JuO1XHC2/3fToHW3mVzL3zp
OGHlED6SPwXZfdEOqvfoqJ83a1veNwUaPP/hxlg6tt3T6Lqxdxqko5pbxOrp
KCb/bry/fxoCblVvmDWsx9ALOa4JSrNwZ3zNwXqNBpwlf+w/H/ALJDY6WG20
b8SrO8uzpb7NAYd3VGfXwybsTi9/TFech7p3xQOaU8R7wL0xRuXmAoR+ZgwI
///9sX/5omvRb9hvvNV8ZqENDyQpun1euwg9dq9XHSztwOvyJzf2nvgLrCA5
Es+/C1XZditbLPnIvNjYQq5ZN64mcYXK5/nI7LhlLoafenAJxmRsiuUnm5n/
cEjM6EWPXUVtK8wFyP5CBeGURiZ+PHg6oXRWgCwp9/NTDH8/zozKig9GCZIV
chKbvysPYI7tAcFIoyVkc568o4zXIAaYBr/JHlxCTjQbcN+YPoSKXWatb4OX
kqsjLtp+nhhGHe37kTI7l5GNXj1MqDvBwhMKVB5/wzLyx9AT+rTCUZypDmN6
311OpisWdrsacDB6ifNVGyUhsrGVw/2LRWO4mpGTJFAuRNa5dmZ30zouLr19
98jeyyvI3kpHhubGuciX5pzBXCVMLrq27FhA0QReeZ0RtqpMmGxI/0ISj5/E
km0PhSQurSSznB+4RQRNoYTZFdUgMRFyyIbkCQWnaezcTmu4USxClq0xkl44
PIOZ8aDU7yxKtq/bytqzZxadMrQk5BdFyfU+t1eJKv7CwODEu50fxMjLePvb
3TfN4adh+RgxCxI58PAS7SEZHj6MCz7FOEQib3zsmgIKPDxrw1FOtyaR1SIC
b0Vv5KFccXGmpT2JnKLb5my+jYePwh2/h54mkePWCgvk6vLwvNarCeEbJLKP
a31B/Ekebrwvabo8lUQuCznEKknh4eIef9GW9yTyo5snDm36wMM2Tl996kcS
WaxDRDYknYcR9l+cTD6TyKOHbYbt8nj4d6uNT3AxiSze7HRquoqHnfTweMEG
Epl5e7DEiM3D/Ae/zjQ0kcjO6bUxmVweRu47qfK2lUR+vWDRKzfDQ/MUtWwD
Bol876TE29k/PCz0q/1xd4RE3hX4hj+DNI9R23eHW7NJZMWq3XIKkvN4mRl7
RHGcGI/Myfgw6XncYnGxmzJNIpOtzrd5Ks7jEv7G5IhfJLLhtMkzhvI89uTs
cz81TyKXyuVmHd46j8Xub3eo/yGR5VRydL6qzePL9Stn/v4lkfluGG7V1JzH
/wFtmZQq
               "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
         "HighlightElements" -> <|
           "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
          "LayoutOptions" -> <|
           "PanelPlotLayout" -> <||>, 
            "PlotRange" -> {{0, 
              4}, {-2.627182088666209, -0.2476451555125296}}, 
            "Frame" -> {{False, False}, {False, False}}, 
            "AxesOrigin" -> {0, -0.11544865922621433`}, 
            "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
            "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), 
            "DefaultStyle" -> {
              Directive[
               Opacity[1.], 
               RGBColor[0.24, 0.6, 0.8], 
               AbsoluteThickness[2]]}, 
            "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
                (Identity[#]& )[
                 Part[#, 1]], 
                (Identity[#]& )[
                 Part[#, 2]]}& ), 
              "ScalingFunctions" -> {{Identity, Identity}, {
                Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
            False|>, 
          "Meta" -> <|
           "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
            Plot, "GroupHighlight" -> False|>|>]]& )[<|
        "HighlightElements" -> <|
          "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
         "LayoutOptions" -> <|
          "PanelPlotLayout" -> <||>, 
           "PlotRange" -> {{0, 4}, {-2.627182088666209, -0.2476451555125296}},
            "Frame" -> {{False, False}, {False, False}}, 
           "AxesOrigin" -> {0, -0.11544865922621433`}, 
           "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
           "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), 
           "DefaultStyle" -> {
             Directive[
              Opacity[1.], 
              RGBColor[0.24, 0.6, 0.8], 
              AbsoluteThickness[2]]}, 
           "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
               (Identity[#]& )[
                Part[#, 1]], 
               (Identity[#]& )[
                Part[#, 2]]}& ), 
             "ScalingFunctions" -> {{Identity, Identity}, {
               Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
           False|>, 
         "Meta" -> <|
          "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
           Plot, "GroupHighlight" -> False|>|>],
       ImageSizeCache->{{4.503599627370496*^15, -4.503599627370496*^15}, {
        4.503599627370496*^15, -4.503599627370496*^15}}],
      Selectable->False]},
    Annotation[{{{{}, {}, 
        Annotation[{
          Directive[
           Opacity[1.], 
           RGBColor[0.24, 0.6, 0.8], 
           AbsoluteThickness[2]], 
          Line[CompressedData["
1:eJwV13k41N0XAHBLi4ShJIRS9KYkUpRwJluylCVJorRQJEpJpUjpTRGVpCwh
lHpDsm9nEpEwM/Z1MPYxw1gzlH7f31/f5/M89/ne7Zxz71U87WVzToCPj8+B
n4/v/9+dMvOjBfAJ0j8sCPstWUqxft7ufmhzOigkNpzpcllCUTXflCZbmQ7v
ClVeuf8SpAgJXhwadMsAnun2EzpPBSknCxY/45FM2P24eDCFJkCxWLOJrpxC
+BezkHpNgKJz+QD38UwmOH6jHrFXEKCsUYlQOxb1GebzlmlV+vFTBIKzLUsH
P8PY2ygjx638FG5P60Ul7Sw4oXJLOpnBR/kZrfiR25IFSd+HHv44xEe5K/Tl
nxCZbIjZul5xc+Qibj3w9GGWdzbIv72Q2COxiPXBXiMdldkwsTlXsvf5H1QS
VP2o5psD7hZOSX1vf+P3P8lqjfW5MKTpyxodm0cv3Xvhf1TywD/BQzg/ZB6l
b7lwN9/NA/otfYbPP/N4YU4+68aOfMgwL/+yz4OHIlNRu9eHFoAwH9d6UnwO
czSuRZn2FYD5kVEFScovdPa2/XVZpxDOCRvft7r8CzM4pILy4UKQYEmfrGyb
RZvhh7ruxsXgpWB/9j1lBhc2u8U9jy2GmoTv/U/vzGDyOePF4qliaFgf7Sms
P4MzTH4KKakEBoPfJVz+Oo3RXTcNcxZLYdd8b1RB6xTul3NIZtghLFc9eTMv
cQpZx7WXCn1CkH5wadeAxxQqsG+c+bSdAqpwpj1NcAq9Lmgp83VSQPNkiWMv
TGL5w5SmyOVfYXVhkNFZiUmUfi/5QEXzK+iMRjEa+iaQMjg5aPPoK/BXBSY4
PZ7AVWcz3qfsKQOZV9VSF/q56HpPwUHnbBkor3i3J66Yi0VJYSuo4WUg8l+R
XdgLLp7t9XCfGywDPclanbCDXMxx3qJq/uIbbM97RCb5jOMxh8QM7kQ5yFZ9
NnviwcHxjiPCSbIV8Kf4xV1bTQ4+cBY6Z2tYAfp3M7snFtj45YyXbG5kBWz9
/e9YajgbxS7pBd/S/g5uEC1ylTKKFUFtx5fdrgSRba8fihmx8MTSsJy8lEro
jhmiuq9l4eS/ZPELdZXA8G9NcBkdwQ1h7yp+rq+CbsqL7FUvRtD/5TX1p2VV
cL9NoFljfBh3/iexTE6oGt4EO1ZLFw/hD7WKU7Xq1XAp7pZAy8shPPXZr+iO
QzU0rHRQXn51CMNyu717P1QD2/EfiZztQzhM+dSRavkTrstaXpJIHcT4poOf
NZ7XQEi57ldIG8Ddx/6s7CuqgWfSq0fJjwewpj3TNbK/BhiXVt6N9hxAXvda
uV+7a2HZZV97Y80BPMIaeFDcWgsqaZ5K3mX9uPLv3RMHFKgg72E0NT3Wh9EX
7J1VTajQq54Ssr+5D5UaVU9JXKKCxqBqZl1JH+qntZzpKKFCaaXFSGpYHxak
ReXkqdHAR2rd593qfcixZReb19HA0CihkxTIxF0rWvT0mTQQGhkO/uPOxJul
X0vVZ2lwr+Loa4WjTFyu8pKyRoEOV9TM9edUmai4uL+825MOCaki+bOdvWj3
/uVPH9F6kK9csI407sVYpyALV8V6cB00+eih2YvMVZ61x3bXw1MHq7SLir3o
dduAqudUDx+dWP4jf3owxIZTv/xTPQxw9E5R8nuw5LdBe4xFAwTauPT47e5B
Zeux4W+PG2FcQeX6gQgGygYc0HFLbIRn8nkBh64wUPxTwmPhvEZg/t4dHH2E
gfNCtmrWzEawlrA7NinDQBolz4expwm6kqfD3r3rwpvqdxfnBppgb2nfh2t1
nejl3H44dqEJvoQnXrbN7sSzoZqJINEMzw7MtJx/3YmHhwcMH+g2w5zfmPAe
t05USjALWf28GSooXwztBTuxjiQpqUZuAUkvGwVZ4w78pud5jm7XAqTrly+r
qnVgvsf33KseLVBGCwryXNuBbyv9HIqiWuC17ir2DVY7+gV0xZtyWuBi+ov4
9GftuHE8VeXM61aIOu65ZNNoG16v2wvR020gu6/S93RBK4qp5Qz8EGqHOcs/
w2nvWzE5TD10Qa4d4M695aToVqRb/NPmbNwOv90fG2pcb8VtPyV9lKPa4Z6V
c0GdVit2V46nZml3gI2+tcyeohY88DVFrPZmJzBdijfrtjQjY4NizuKTTjhh
MJNS+aMZrwbGOqq/7QSu++rxC8XNmAjP3z//2QmnDk2oTCY240JJoMFxuS7I
nDJ1F7nUjBmFjr5DJV2QkJvTYyDcjGuzJRj8/N1w5/zyvGTrJjwRunCkX6ob
zCVjFenGTZh0duDnd9VuUKzRTlyl04SqawoKHh/rBnan577GjU1o4HvqxZrM
bpDZxO9nOdOIWrH2qZE7emBc3cJ+Nq4RXZxTnd4zekDsk5XTlj8NeN9s3atX
Iz1AkxmJz5hqwFStiMZH0z3gWFceeYDVgGyxm+aewr0wf2btlrSWBvRFy707
tXrh/AO1Hzu+NGDohhnJ4rBemKneR9np2YC5TMMaqi4TNFUN/BzZ9dhsenci
7AATPon59Q731+NseqmUhQ0TTpeM2N/rqkftW3tdqtyYUP9FfV13XT3mr94x
S3nGhP7t52t2ZNVjgZHshqxhJghUbv+ldKsei1O5PpGRffDsAl24Z009doqo
vbJ50wd+M+3TLFI9/r7iUSr+oQ9U+75n8QnXox4MCoVhHyRtE/C0WKRjaUtH
XDCrD5r/uO4+M0RHFKqs9CX3Q/bbcuu7RXQsc49b58DuJ+Lpe/F6Dzp+GfHw
cJ7tBwPmVFChKx2Tz+sUneEbAKlMBT2n03QMdm05fklyABLatdwKHehoelri
9T3dASjfNfao3pSOVIdg6fTQAdieThW9soWOlFbb81+iBmBadk3pv0p0zLTf
mJ+fMACC11zpyRvo+MwO7b9lD0AxY9ffaWk6HrXmRbV2DoDkxPow9go6dpp6
rhHcPgjiEVQrDQ4Nayv3nRPSHgSP45ZHC0ZoWGoinCO6fxDOfKnSMx2k4Ruj
d0ek7QYhqNI98mY3DU+Te59vvz0It59rhqvX03BY227VsdpBUG9NPdGRT8OZ
zXqkT55DEBikO7XwiIb6TIPKR35DsHGLnYnMQxo+iDMNOH9vCPRTXmXoBdNQ
SvLI+KZXQ2AWkpISE0BDLX4P6uvyITh3nnHw41Ua+nVEhz9cNwzxTlncVado
SImKP+C6eRiUBdUU451oKGST/NdQYxgWQ4QmdzjS8FVVhtdfk2GYPKbu4HaU
hoU53w/7XhmGhw5H7yhZ0vB3+DTp7I9hsDCPtK3fR0Mj8/nK/Y3DUK8oeK58
Lw1Dl/EFru8ehrUBYZ+KtWkof3slt316GJ7Ea/xETWI+7htp1htGwKwlP4yk
SsMAI6sIuD4CO3v556nyNKz6a2cqHzQC4/kC/22Uo6F4kSPfQugI1Fyr0Lol
S8MEDTfv3KQREDhgR9FbS8xH4Y7V9roRGJZbKSwkQUMB3kfxdcosmA9T/PHf
UhquqzWYgh0sCJyrqvVeQsNdiW1NZ/eyQH6HceJeQRqeO7g8Jt2SBdVfvZLb
+Ij+X51WNvBlQerBj3Ypc1Ts9eQtdwtkQenX2+43Zqk4vz+C9fgRC8o4wl1H
p6moyirJaIpnwZ1Yo+LtXCqG68joXKhkgfakTajZMBXfi2XKPaGzIM6nPO/8
IBUpTJO/WR0sKPRfFIzop+Lko6vlv8dZkOyi2svroeLRDuqhCOlRaFW507u3
nYpeGa4aORtHgWUyZZzQSsWH9/6sblcdBSkb0lNSCxULtm1r27R/FIQOzJUK
NVJRzv/BmbwLo6Aws5Bzto6Ku63kTTp9RuHJ35+PJWqpeEgpewv/nVF4Rq8S
r/xJxYCaXo7Zs1Ho57cUtfhBxegEP5pX7CgkGXpdX1dFxcyrpC+RqaNwa3zq
2tR3KjLl9PwYhaMwS69YU1JOxQVuw3HBCmK83xbS0r9RUbLCXW8LdRTUfc7l
vy+joolntMCVvlHg+yO6NZdCxbSSqdClImzgc7NN311CxbKnj7y2SrFBzuK+
i3sxFdvPKdoc3sAGt32WJ94VUVFE7PDaV7vYEFGyVFynkIrKzAFeiT4brmQW
FocXUFE/17+TacqGastn0Zx8Kno7f0hUdWLDnv2RFV/zqBiyc/99azc2ZF39
IbyHcNKyVlffy2zIGdNyy82lYmH7pYMxt9jAcPjVqEu4IX2pKiWYDbPfth/6
mUNFdlCs2EA4G6TDl9a6EF5irzmx4jUbhk9nmfMRlt9W3aCWzIZRee/y1Gwq
av09lWubzoYllUEadoQtLClBxvlsmPcLYYsQPhOz/rB2GRv+uV+yv+YLFW+O
3FmnUsOG6HN1254TjtBmDMk2E+sxfiXFhXBqsF62SA8bknQupGkTLm6IDVwc
YUN+nO/ONYTrFX9bcKfYIBXntm8+i4ojXo4yzD9scOgh4SBhvtLCgYblHChV
NstvJywlIptVIcGB5Vpzis2EVY/fuJO3jgO+rQK/WgkbvG81S1PmQMt3h519
hI/Naq+N2cGB3MSpumnCXkYv+0L3ciCqv7ZSlOj/wbPZjDuGHDjUzl6rRji2
x87f25IDQU8dC+0IZ6nlmJ6258AG+rYP9wlX+UuuOeLCAWFRX2YhYUa1T6+x
Bweqjzs4zxGekW74pH2NAy6uC9J6xHqJuO28qRLAAVW3h6IhhDfmPDVZF8IB
hqSmdifhQ9ZW3YuxHBCKqpGKJvbj7JuMj9xUDvw3R/7GR+znLY6YHzOTA0oS
jKfehN+H1Ip/L+eACV9f0DkiHkpbVLvy6jgwXSgQyyLcqByaltZKzJdRWO5L
xBN/mZlBGJsDLzvf7Uwm4k1a/INYwCwHNN5WXDIm4lHNeUWHN98Y5L60fM8h
fHy+0ueI5Bj4tWXw2xLx7G36D9lEYQwetSVJixPx/m/UA5E9W8Ygw0FjXSPh
7J3GKet0x+D0g2d0b6RidWDyZVGTMej0+c//IJEvPXWC+n8Pj4GDFZdf5SsV
RT3KmplnxsDtd1nIbyK/XN+C0IfHY3Co8MLaGSI//SfiG2NejIFZVaDcEiJ/
n8NiQtgb4n8v2G0yRH5TOor3Xv4yBpvSwNOByH8ZSR2PPR1j0Kjvfm8LjYo1
93fVfd86DkoHk0YFOqioQ/op8HT3OATNnmPGdBLr98pF25FM+KrwdV0GFe+l
P0kYtxsH26J7rpG9RPvW4SvSd8fB1SGsJf7/9W9b3Fr3lnGI4ZKYMkQ9lcrV
tNjFHAeWkcTtrfNUvE+uDlxkj8PRv+s8DH4T9ebDwtWD/FywD9u71oiox49x
C6NLiguTH/U9coj6bTASlLnMgAuOjGlvIzEa3lCrq/Y4wIWLK0WsH5BomOkj
M0Cz4EJIsHH+T3EaKixmyMTYc8Hcb2S1y2oa8lZ1Be3w5MI/i6sqc6SJ9rra
dseiuQAd952XKtFwfcQoL22MC2E2i3079Ijz3CUgLWuKC6FSS8e+6tPQWXO1
Q9EcFxQkA5qOkmkY26yTXyMwATUCG7RDDGkoLf/Id1xqAmz2FsqLm9Fw1Yct
U7tgAux6mg1S7Gm4tPwcG8Mn4Ntficg+4nxOfsGLqYqcAJ0M0491vsR56hZm
Tn81AXzuz68V+9EwSDjnA/PtBARp3lZI8KfhotWSC0vzJ0A08NWbf+/TcK4r
acCsZwISnNRsmyNpOPqrm9GkPgnXn4i4uefRMNctvMxMaxJ0c+8KthfQMLBF
PxX3TYKNRX6CRTFxf8iL80wzmQQW63eLzldifX0df/ufmASrw0aS+37S8PV0
s4zSw0lI4uOL/UHcX0wnao9c6Z4Evv/meWtWEvcxVlG12JMpED0Yt1XFnY6X
Xh7e5P58CgLuFF0le9JR26jvVkX0FHRuXTh+3JuO1XHC2/3fToHW3mVzL3zp
OGHlED6SPwXZfdEOqvfoqJ83a1veNwUaPP/hxlg6tt3T6Lqxdxqko5pbxOrp
KCb/bry/fxoCblVvmDWsx9ALOa4JSrNwZ3zNwXqNBpwlf+w/H/ALJDY6WG20
b8SrO8uzpb7NAYd3VGfXwybsTi9/TFech7p3xQOaU8R7wL0xRuXmAoR+ZgwI
///9sX/5omvRb9hvvNV8ZqENDyQpun1euwg9dq9XHSztwOvyJzf2nvgLrCA5
Es+/C1XZditbLPnIvNjYQq5ZN64mcYXK5/nI7LhlLoafenAJxmRsiuUnm5n/
cEjM6EWPXUVtK8wFyP5CBeGURiZ+PHg6oXRWgCwp9/NTDH8/zozKig9GCZIV
chKbvysPYI7tAcFIoyVkc568o4zXIAaYBr/JHlxCTjQbcN+YPoSKXWatb4OX
kqsjLtp+nhhGHe37kTI7l5GNXj1MqDvBwhMKVB5/wzLyx9AT+rTCUZypDmN6
311OpisWdrsacDB6ifNVGyUhsrGVw/2LRWO4mpGTJFAuRNa5dmZ30zouLr19
98jeyyvI3kpHhubGuciX5pzBXCVMLrq27FhA0QReeZ0RtqpMmGxI/0ISj5/E
km0PhSQurSSznB+4RQRNoYTZFdUgMRFyyIbkCQWnaezcTmu4USxClq0xkl44
PIOZ8aDU7yxKtq/bytqzZxadMrQk5BdFyfU+t1eJKv7CwODEu50fxMjLePvb
3TfN4adh+RgxCxI58PAS7SEZHj6MCz7FOEQib3zsmgIKPDxrw1FOtyaR1SIC
b0Vv5KFccXGmpT2JnKLb5my+jYePwh2/h54mkePWCgvk6vLwvNarCeEbJLKP
a31B/Ekebrwvabo8lUQuCznEKknh4eIef9GW9yTyo5snDm36wMM2Tl996kcS
WaxDRDYknYcR9l+cTD6TyKOHbYbt8nj4d6uNT3AxiSze7HRquoqHnfTweMEG
Epl5e7DEiM3D/Ae/zjQ0kcjO6bUxmVweRu47qfK2lUR+vWDRKzfDQ/MUtWwD
Bol876TE29k/PCz0q/1xd4RE3hX4hj+DNI9R23eHW7NJZMWq3XIKkvN4mRl7
RHGcGI/Myfgw6XncYnGxmzJNIpOtzrd5Ks7jEv7G5IhfJLLhtMkzhvI89uTs
cz81TyKXyuVmHd46j8Xub3eo/yGR5VRydL6qzePL9Stn/v4lkfluGG7V1JzH
/wFtmZQq
           "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
     "HighlightElements" -> <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
      "LayoutOptions" -> <|
       "PanelPlotLayout" -> <||>, 
        "PlotRange" -> {{0, 4}, {-2.627182088666209, -0.2476451555125296}}, 
        "Frame" -> {{False, False}, {False, False}}, 
        "AxesOrigin" -> {0, -0.11544865922621433`}, 
        "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
        "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), "DefaultStyle" -> {
          Directive[
           Opacity[1.], 
           RGBColor[0.24, 0.6, 0.8], 
           AbsoluteThickness[2]]}, 
        "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), 
          "ScalingFunctions" -> {{Identity, Identity}, {
            Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> False|>, 
      "Meta" -> <|
       "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
        Plot, "GroupHighlight" -> False|>|>, "DynamicHighlight"]], 
   InterpretationBox[{
     TagBox[{{{}, {}, 
        TagBox[
         {RGBColor[0.24, 0.6, 0.8], AbsoluteThickness[2], Opacity[1.], 
          LineBox[CompressedData["
1:eJxTTMoPSmViYGAwAWIQ7ZX64uKBl6V2l8Sdzv3+f95+gah/W5fkEnsYv5Fj
s3qn5BY4PyJy4foPH4/A+dz/G2Pc5c7D+aqB714c7r4C55efs7Sf8eUmnC++
RfAeI+N9OH/bI+cz520ewfmHsuZKR755Aud/VbPlX5v7HM5n+rlaQFr1FZxv
9j9hW/C6N3D+mRaTc8e03sP58hNe/1z57gOc//r7/XtXDT7B+Ute7T7F1/cZ
zueTXf7+yZMvcH5P5ta0BSrf4PxvDqufZNR/h/NLjI5sETv8A86/v+5I90XF
X3D+76wrszWrfsP52o7s/9J2/4Hz3Rcppm8U/4cIL9l4pYcx/+F8nTeh3Nd9
GRxgfGH+DxxHfiH4LPtnr1eewwjnZ5vsvsnpzQTnr/ZMWrDvG4L/9bWUwLNp
zHD+1mB35ikuLHB+vUfr/C3PEHzFu143FreywvlW5i1TJI3Y4PwYufM/GS8j
+F9P9T4qaGSH82ewxJUEqXAg3H9v6yKmIwg+a21jiGUhJ5zPsDJu/SMhLji/
aNb6XqFDCP5e7Q4OwTxuOF/Qq0iniY8Hzr+je+Fy5R4Ef8M8e5Uncbxwfux6
M0HZfwh+Q+vCxjur+OD8tS9kZ/P58MP5Si0iHuzLEPzp8txf//9H8AGeXsTJ

           "]]},
         Annotation[#, "Charting`Private`Tag#1"]& ]}, {}},
      {"WolframDynamicHighlight", <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>}], 
     StyleBox[
      DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {}, 
        Slot["HighlightElements"], 
        Slot["LayoutOptions"], 
        Slot["Meta"], 
        Charting`HighlightActionFunction["DynamicHighlight", {{{{}, {}, 
            Annotation[{
              Directive[
               Opacity[1.], 
               RGBColor[0.24, 0.6, 0.8], 
               AbsoluteThickness[2]], 
              Line[CompressedData["
1:eJxTTMoPSmViYGAwAWIQ7ZX64uKBl6V2l8Sdzv3+f95+gah/W5fkEnsYv5Fj
s3qn5BY4PyJy4foPH4/A+dz/G2Pc5c7D+aqB714c7r4C55efs7Sf8eUmnC++
RfAeI+N9OH/bI+cz520ewfmHsuZKR755Aud/VbPlX5v7HM5n+rlaQFr1FZxv
9j9hW/C6N3D+mRaTc8e03sP58hNe/1z57gOc//r7/XtXDT7B+Ute7T7F1/cZ
zueTXf7+yZMvcH5P5ta0BSrf4PxvDqufZNR/h/NLjI5sETv8A86/v+5I90XF
X3D+76wrszWrfsP52o7s/9J2/4Hz3Rcppm8U/4cIL9l4pYcx/+F8nTeh3Nd9
GRxgfGH+DxxHfiH4LPtnr1eewwjnZ5vsvsnpzQTnr/ZMWrDvG4L/9bWUwLNp
zHD+1mB35ikuLHB+vUfr/C3PEHzFu143FreywvlW5i1TJI3Y4PwYufM/GS8j
+F9P9T4qaGSH82ewxJUEqXAg3H9v6yKmIwg+a21jiGUhJ5zPsDJu/SMhLji/
aNb6XqFDCP5e7Q4OwTxuOF/Qq0iniY8Hzr+je+Fy5R4Ef8M8e5Uncbxwfux6
M0HZfwh+Q+vCxjur+OD8tS9kZ/P58MP5Si0iHuzLEPzp8txf//9H8AGeXsTJ

               "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
         "HighlightElements" -> <|
           "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
          "LayoutOptions" -> <|
           "PanelPlotLayout" -> <||>, 
            "PlotRange" -> {{0., 4}, {0., 0.49999899999999997`}}, 
            "Frame" -> {{False, False}, {False, False}}, 
            "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 360/GoldenRatio}, 
            "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
            GoldenRatio^(-1), "DefaultStyle" -> {
              Directive[
               Opacity[1.], 
               RGBColor[0.24, 0.6, 0.8], 
               AbsoluteThickness[2]]}, 
            "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
                (Identity[#]& )[
                 Part[#, 1]], 
                (Identity[#]& )[
                 Part[#, 2]]}& ), 
              "ScalingFunctions" -> {{Identity, Identity}, {
                Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
            False|>, 
          "Meta" -> <|
           "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
            Plot, "GroupHighlight" -> False|>|>]]& )[<|
        "HighlightElements" -> <|
          "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
         "LayoutOptions" -> <|
          "PanelPlotLayout" -> <||>, 
           "PlotRange" -> {{0., 4}, {0., 0.49999899999999997`}}, 
           "Frame" -> {{False, False}, {False, False}}, 
           "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 360/GoldenRatio}, 
           "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
           GoldenRatio^(-1), "DefaultStyle" -> {
             Directive[
              Opacity[1.], 
              RGBColor[0.24, 0.6, 0.8], 
              AbsoluteThickness[2]]}, 
           "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
               (Identity[#]& )[
                Part[#, 1]], 
               (Identity[#]& )[
                Part[#, 2]]}& ), 
             "ScalingFunctions" -> {{Identity, Identity}, {
               Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
           False|>, 
         "Meta" -> <|
          "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
           Plot, "GroupHighlight" -> False|>|>],
       ImageSizeCache->{{4.503599627370496*^15, -4.503599627370496*^15}, {
        4.503599627370496*^15, -4.503599627370496*^15}}],
      Selectable->False]},
    Annotation[{{{{}, {}, 
        Annotation[{
          Directive[
           Opacity[1.], 
           RGBColor[0.24, 0.6, 0.8], 
           AbsoluteThickness[2]], 
          Line[CompressedData["
1:eJxTTMoPSmViYGAwAWIQ7ZX64uKBl6V2l8Sdzv3+f95+gah/W5fkEnsYv5Fj
s3qn5BY4PyJy4foPH4/A+dz/G2Pc5c7D+aqB714c7r4C55efs7Sf8eUmnC++
RfAeI+N9OH/bI+cz520ewfmHsuZKR755Aud/VbPlX5v7HM5n+rlaQFr1FZxv
9j9hW/C6N3D+mRaTc8e03sP58hNe/1z57gOc//r7/XtXDT7B+Ute7T7F1/cZ
zueTXf7+yZMvcH5P5ta0BSrf4PxvDqufZNR/h/NLjI5sETv8A86/v+5I90XF
X3D+76wrszWrfsP52o7s/9J2/4Hz3Rcppm8U/4cIL9l4pYcx/+F8nTeh3Nd9
GRxgfGH+DxxHfiH4LPtnr1eewwjnZ5vsvsnpzQTnr/ZMWrDvG4L/9bWUwLNp
zHD+1mB35ikuLHB+vUfr/C3PEHzFu143FreywvlW5i1TJI3Y4PwYufM/GS8j
+F9P9T4qaGSH82ewxJUEqXAg3H9v6yKmIwg+a21jiGUhJ5zPsDJu/SMhLji/
aNb6XqFDCP5e7Q4OwTxuOF/Qq0iniY8Hzr+je+Fy5R4Ef8M8e5Uncbxwfux6
M0HZfwh+Q+vCxjur+OD8tS9kZ/P58MP5Si0iHuzLEPzp8txf//9H8AGeXsTJ

           "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
     "HighlightElements" -> <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
      "LayoutOptions" -> <|
       "PanelPlotLayout" -> <||>, 
        "PlotRange" -> {{0., 4}, {0., 0.49999899999999997`}}, 
        "Frame" -> {{False, False}, {False, False}}, "AxesOrigin" -> {0, 0}, 
        "ImageSize" -> {360, 360/GoldenRatio}, "Axes" -> {True, True}, 
        "LabelStyle" -> {}, "AspectRatio" -> GoldenRatio^(-1), "DefaultStyle" -> {
          Directive[
           Opacity[1.], 
           RGBColor[0.24, 0.6, 0.8], 
           AbsoluteThickness[2]]}, 
        "HighlightLabelingFunctions" -> <|"CoordinatesToolOptions" -> ({
            (Identity[#]& )[
             Part[#, 1]], 
            (Identity[#]& )[
             Part[#, 2]]}& ), 
          "ScalingFunctions" -> {{Identity, Identity}, {
            Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> False|>, 
      "Meta" -> <|
       "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
        Plot, "GroupHighlight" -> False|>|>, "DynamicHighlight"]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox[
     TagBox["\"y\"", HoldForm], TraditionalForm], 
    FormBox[
     TagBox["\"U\"", HoldForm], TraditionalForm]},
  AxesOrigin->{0, -0.11544865922621433`},
  DefaultBaseStyle->{"PlotGraphics", "Graphics"},
  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]],
  ImagePadding->All,
  Method->{
   "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotInteractivity:>True,
  PlotRange->{{0, 4}, {0, -3}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.543849693769854*^9, 3.544185953978735*^9, 3.54418678040519*^9, 
   3.6073372898207607`*^9, 3.9240041209449177`*^9, 3.924004279444892*^9, 
   3.924077940031032*^9, {3.924078007193041*^9, 3.924078033089471*^9}, 
   3.924078087943061*^9, 3.924078134384706*^9, 3.924078166060299*^9, 
   3.924078230611949*^9, 3.924082152759513*^9, 3.9240823609335127`*^9, {
   3.9241897274192963`*^9, 3.9241898518527107`*^9}, 3.986975961620541*^9, 
   3.986979672363201*^9, 3.9869797031052713`*^9},
 CellLabel->"Out[55]=",ExpressionUUID->"9ce7da3d-f19b-4696-bd9e-48d537d89aa9"]
}, Open  ]],

Cell[TextData[{
 "To solve for a big bang FRW model (one that starts at ",
 StyleBox["a=0 ",
  FontSlant->"Italic"],
 "or ",
 StyleBox["y=0)  ",
  FontSlant->"Italic"],
 "is equivalent to carrying out the integral:\n\n       ",
 StyleBox["x(y)=",
  FontSlant->"Italic"],
 Cell[BoxData[
  FormBox[
   RowBox[{
    SubsuperscriptBox["\[Integral]", "0", "y"], " ", 
    RowBox[{
     SuperscriptBox[
      RowBox[{"[", 
       RowBox[{
        RowBox[{
         SubscriptBox["\[CapitalOmega]", "k"], "(", " ", 
         RowBox[{
          SubscriptBox["\[CapitalOmega]", "r"], ",", " ", 
          SubscriptBox["\[CapitalOmega]", "m"], ",", 
          SubscriptBox["\[CapitalOmega]", "v"]}], ")"}], "-", 
        RowBox[{"2", 
         RowBox[{"U", "(", 
          RowBox[{
           RowBox[{"w", ";", " ", 
            SubscriptBox["\[CapitalOmega]", "r"]}], ",", " ", 
           SubscriptBox["\[CapitalOmega]", "m"], ",", " ", 
           SubscriptBox["\[CapitalOmega]", "v"]}], ")"}]}]}], "]"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "/", "2"}], ")"}]], 
     StyleBox[
      RowBox[{
       StyleBox["d",
        FontSlant->"Italic"], "w"}]]}]}], TraditionalForm]],ExpressionUUID->
  "a6cf582d-d738-4362-a9b5-48487af3329a"],
 "                                           (*)\n  "
}], "Text",
 CellChangeTimes->{{3.9240038765685387`*^9, 
  3.9240038768702593`*^9}},ExpressionUUID->"1b00e417-94b2-44d8-aeaa-\
03c916d18945"],

Cell[BoxData[
 RowBox[{
  RowBox[{"x", "[", "y_", "]"}], ":=", 
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"\[CapitalOmega]k", "-", 
       RowBox[{"2", 
        RowBox[{"U", "[", "w", "]"}]}]}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "/", "2"}], ")"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"w", ",", "0", ",", "y"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.544186549753783*^9, {3.924003884398019*^9, 3.924003885056386*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"7ec99045-096b-4996-b719-e2790e99b097"],

Cell["That's the answer. ", "Text",
 CellChangeTimes->{
  3.9240039033526573`*^9},ExpressionUUID->"ba90b5ad-34e6-4a6f-89b5-\
e9771c1e559e"],

Cell[TextData[{
 "For values of ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "3b29b22d-045f-453d-82cf-4b9db6b85d81"],
 " above the maximum of the potential the expansion is unbounded in time;  \
for values that are below the expansion is bounded and the universe \
recollapses.  The value just at the maximum divides these two behaviors. We \
denote that value by ",
 "\[CapitalOmega]",
 StyleBox["crit.",
  FontWeight->"Bold"]
}], "Text",
 CellChangeTimes->{
  3.54418655588161*^9, {3.924003909726821*^9, 
   3.924003909986376*^9}},ExpressionUUID->"8476146b-7e78-4979-9dda-\
1d08c432c551"],

Cell[BoxData[
 RowBox[{"crit", ":=", " ", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], 
     RowBox[{"U", "[", "yy", "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"yy", ",", " ", ".5"}], "}"}]}], "]"}]}]], "Input",
 CellLabel->"In[57]:=",ExpressionUUID->"3740c1c5-696a-4d99-9c68-fb1efb9a576e"],

Cell[BoxData[
 RowBox[{"\[CapitalOmega]crit", ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]\[CapitalLambda]", ">", ".000001"}], ",", 
    RowBox[{"-", 
     RowBox[{"crit", "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "0"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5441865570495*^9, 3.54418655809352*^9}, {
  3.92400415963719*^9, 3.924004170138816*^9}},
 CellLabel->"In[58]:=",ExpressionUUID->"583820f6-dafd-465b-9dd4-97e8c8612159"],

Cell[CellGroupData[{

Cell[BoxData["\[CapitalOmega]crit"], "Input",
 CellChangeTimes->{3.544186559145056*^9},
 CellLabel->"In[59]:=",ExpressionUUID->"01570fca-4386-4bbb-9d3e-58ed4c4aa626"],

Cell[BoxData[
 RowBox[{"-", "0.4952902987449596`"}]], "Output",
 CellChangeTimes->{
  3.543849722480421*^9, 3.544186218600436*^9, 3.5441868050825443`*^9, 
   3.607337302007023*^9, 3.924004121010762*^9, 3.924004279505962*^9, 
   3.9240779401231413`*^9, {3.924078007263979*^9, 3.9240780331557283`*^9}, 
   3.924078088000018*^9, 3.924078134438815*^9, 3.9240781661155376`*^9, 
   3.924078230671218*^9, 3.924082152812028*^9, 3.924082361007127*^9, {
   3.92418972748875*^9, 3.9241898519230556`*^9}, 3.9869759616688967`*^9, 
   3.986979672415781*^9, 3.9869797031318293`*^9},
 CellLabel->"Out[59]=",ExpressionUUID->"3c964217-7dca-444c-a7bf-3e6608ed9230"]
}, Open  ]],

Cell[TextData[{
 "For bounded cases the integration in (*) can extend only to the first \
turning point where the denominator vanishes.  This represents the expansion \
of a universe up to its maximum scale factor ",
 StyleBox["ymax.  ",
  FontSlant->"Italic"],
 "Since the integral is singular at the turning point it is necessary to stop \
it just a bit before by small fraction \[Epsilon] ."
}], "Text",ExpressionUUID->"b8a87e39-2001-4338-801b-07326d4ae205"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"eps", "=", ".0001"}]], "Input",
 CellLabel->"In[60]:=",ExpressionUUID->"dd2d6361-44e0-4c71-a27f-011146a8cb15"],

Cell[BoxData["0.0001`"], "Output",
 CellChangeTimes->{
  3.543849735207077*^9, 3.544186025060561*^9, 3.544186220272471*^9, 
   3.544186818613674*^9, 3.60733730553647*^9, 3.9240041210264273`*^9, 
   3.9240042795112257`*^9, 3.924077940136221*^9, {3.924078007269422*^9, 
   3.924078033173403*^9}, 3.9240780880053053`*^9, 3.924078134455068*^9, 
   3.924078166120699*^9, 3.9240782306872263`*^9, 3.924082152817217*^9, 
   3.924082361012456*^9, {3.924189727505657*^9, 3.924189851942877*^9}, 
   3.986975961672229*^9, 3.986979672419487*^9, 3.9869797031348248`*^9},
 CellLabel->"Out[60]=",ExpressionUUID->"a5f02194-6d1c-4e49-b5a9-af5bcae23f0f"]
}, Open  ]],

Cell[TextData[{
 "The following list of statements looks for the turning point which is the \
smallest positive value of ",
 StyleBox["y ",
  FontSlant->"Italic"],
 "where ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "d21e2517-4dcd-4fec-b549-4d0a2d2ec358"],
 "/2=U(tp). This is relevant only for closed models where ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["\[CapitalOmega]", "k"], TraditionalForm]],ExpressionUUID->
  "94db163c-58bb-4555-b616-abafd38076e9"],
 "/2 < ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SubscriptBox["U", "max"], ".", " "}], TraditionalForm]],ExpressionUUID->
  "45c365c5-9b88-4dab-8d20-c970521a3ce5"],
 " (If the program fails to identify the correct turning point, its possible \
that you might have to insert a  statement ",
 StyleBox["ymax=  ",
  FontWeight->"Bold"],
 StyleBox["to give it the correct value and then run the program again.) ",
  FontVariations->{"CompatibilityType"->0}],
 "\n"
}], "Text",
 CellChangeTimes->{{3.9240039656680326`*^9, 
  3.924003980853455*^9}},ExpressionUUID->"d7ec73d4-64e2-4542-9e89-\
3dbd875a71c6"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"listtprules", "=", 
  RowBox[{"NSolve", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]k", "\[Equal]", 
     RowBox[{"2", 
      RowBox[{"U", "[", "tp", "]"}]}]}], ",", " ", "tp"}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.544186564266329*^9, {3.9240039909430513`*^9, 3.9240039916174507`*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"cf55b747-aa05-4c94-a638-7e3f424bf4bd"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"tp", "\[Rule]", 
     RowBox[{"0.2688936672588881`", "\[VeryThinSpace]", "-", 
      RowBox[{"1.6483035582842214`", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"tp", "\[Rule]", 
     RowBox[{"0.2688936672588881`", "\[VeryThinSpace]", "+", 
      RowBox[{"1.6483035582842214`", " ", "\[ImaginaryI]"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"tp", "\[Rule]", 
     RowBox[{"-", "0.5377840011658397`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"tp", "\[Rule]", 
     RowBox[{"-", "3.3333518520204363`*^-6"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.543849739817165*^9, 3.544186033670064*^9, 3.5441862224027987`*^9, 
   3.544186827200459*^9, 3.607337309763323*^9, 3.9240041210561237`*^9, 
   3.9240042795188923`*^9, 3.924077940150835*^9, {3.924078007286268*^9, 
   3.924078033189554*^9}, 3.9240780880229053`*^9, 3.924078134468679*^9, 
   3.924078166141021*^9, 3.924078230700596*^9, 3.924082152833701*^9, 
   3.92408236102601*^9, {3.9241897275224867`*^9, 3.9241898519535313`*^9}, 
   3.986975961693716*^9, 3.986979672470963*^9, 3.986979703139876*^9},
 CellLabel->"Out[61]=",ExpressionUUID->"4b7ffe62-10ec-49a8-b415-8b5a94847359"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"listtp", "=", 
  RowBox[{"tp", "/.", "listtprules"}]}]], "Input",
 CellLabel->"In[62]:=",ExpressionUUID->"6080aa2a-db61-4e88-abae-2f37d15f5161"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"0.2688936672588881`", "\[VeryThinSpace]", "-", 
    RowBox[{"1.6483035582842214`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.2688936672588881`", "\[VeryThinSpace]", "+", 
    RowBox[{"1.6483035582842214`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"-", "0.5377840011658397`"}], ",", 
   RowBox[{"-", "3.3333518520204363`*^-6"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.543849752519454*^9, 3.5441861152665586`*^9, 3.544186225310961*^9, 
   3.544186829276134*^9, 3.607337317137867*^9, 3.924004121062395*^9, 
   3.924004279534696*^9, 3.924077940160308*^9, {3.924078007298748*^9, 
   3.924078033200556*^9}, 3.924078088030435*^9, 3.9240781344805202`*^9, 
   3.9240781661468477`*^9, 3.924078230706118*^9, 3.9240821528407993`*^9, 
   3.9240823610354*^9, {3.92418972753192*^9, 3.924189851976017*^9}, 
   3.986975961697195*^9, 3.9869796724754868`*^9, 3.986979703154643*^9},
 CellLabel->"Out[62]=",ExpressionUUID->"6acfab07-ed53-4c39-b503-306a5a2567b2"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"listpostp", "=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", 
    RowBox[{"Select", "[", 
     RowBox[{"listtp", ",", 
      RowBox[{
       RowBox[{"#", ">", "0"}], " ", "&"}]}], "]"}], ",", 
    RowBox[{"{", "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.544186565925921*^9, 3.544186566764421*^9}, {
  3.924004001662677*^9, 3.92400400231942*^9}},
 CellLabel->"In[63]:=",ExpressionUUID->"3d1ab90e-aebf-4da3-9510-2a96657115ab"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.543849755464974*^9, 3.5441860441481857`*^9, 3.544186118083769*^9, 
   3.544186227257855*^9, 3.544186831285268*^9, 3.607337318775078*^9, 
   3.924004121074664*^9, 3.9240042795409613`*^9, 3.9240779401749983`*^9, {
   3.924078007313829*^9, 3.92407803321664*^9}, 3.924078088044578*^9, 
   3.924078134494039*^9, 3.9240781661520844`*^9, 3.924078230719664*^9, 
   3.924082152853753*^9, 3.924082361048573*^9, {3.924189727544983*^9, 
   3.924189851982441*^9}, 3.9869759617203703`*^9, 3.986979672484428*^9, 
   3.986979703158341*^9},
 CellLabel->"Out[63]=",ExpressionUUID->"0dc69b4c-9edd-49b0-8842-fc90ba9812e0"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"listpostp", "=", 
  RowBox[{"Sort", "[", "listpostp", "]"}]}]], "Input",
 CellLabel->"In[64]:=",ExpressionUUID->"a4da45c1-75da-4ca9-be97-7fc53ab0d6a2"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.543849759069757*^9, 3.544186053405918*^9, 3.544186118901965*^9, 
   3.5441862282272367`*^9, 3.5441868363762217`*^9, 3.6073373213662977`*^9, 
   3.924004121084227*^9, 3.9240042795620813`*^9, 3.924077940180242*^9, {
   3.924078007319363*^9, 3.924078033221991*^9}, 3.924078088049974*^9, 
   3.924078134499568*^9, 3.9240781661716347`*^9, 3.9240782307315397`*^9, 
   3.924082152859104*^9, 3.924082361053957*^9, {3.92418972755098*^9, 
   3.924189851988007*^9}, 3.986975961724132*^9, 3.986979672487631*^9, 
   3.9869797031680393`*^9},
 CellLabel->"Out[64]=",ExpressionUUID->"a41e1fc8-c4ce-4bb5-bddf-bb9532b50ef1"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", " ", 
   RowBox[{"ymax", ":=", 
    RowBox[{
     RowBox[{"listpostp", "[", 
      RowBox[{"[", "1", "]"}], "]"}], 
     RowBox[{"(", 
      RowBox[{"1", "-", "eps"}], ")"}]}]}], ",", 
   RowBox[{"ymax", ":=", " ", "yend"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.544186568052351*^9, 3.5441865687909327`*^9}, {
  3.92400400961639*^9, 3.92400400993513*^9}},
 CellLabel->"In[65]:=",ExpressionUUID->"0fdc212b-a17f-43b1-ba65-8f89c8ccb18e"],

Cell[CellGroupData[{

Cell[BoxData["ymax"], "Input",
 CellLabel->"In[66]:=",ExpressionUUID->"d8af3ffd-3003-4701-840e-2c0d0fe9a365"],

Cell[BoxData["4"], "Output",
 CellChangeTimes->{
  3.544186120591332*^9, 3.5441868392061977`*^9, 3.924004121102776*^9, 
   3.924004279584198*^9, 3.924077940188068*^9, {3.9240780073400507`*^9, 
   3.924078033243046*^9}, 3.924078088069861*^9, 3.924078134517318*^9, 
   3.924078166195859*^9, 3.924078230748005*^9, 3.924082152867457*^9, 
   3.924082361062097*^9, {3.924189727568307*^9, 3.92418985201167*^9}, 
   3.9869759617443037`*^9, 3.986979672492906*^9, 3.986979703172762*^9},
 CellLabel->"Out[66]=",ExpressionUUID->"191be143-8e9b-4eb0-8cef-2f49962db839"]
}, Open  ]],

Cell[TextData[{
 "In the unbounded cases the integration can extend to the limit of the plot, \
 ",
 StyleBox["yend. ",
  FontSlant->"Italic"]
}], "Text",ExpressionUUID->"d512f051-dd5c-434e-b292-e5cf1120fa38"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plotting the Results", "Subsection",ExpressionUUID->"f8efa3a5-2ecc-4440-bc89-384f35b8cc44"],

Cell[TextData[{
 "A parameter ",
 StyleBox["z ",
  FontSlant->"Italic"],
 "is introduced which runs from ",
 StyleBox["0",
  FontSlant->"Italic"],
 " to ",
 StyleBox["1",
  FontSlant->"Italic"],
 " along the curve to be plotted, either from the big bang to the big crunch \
for bounded models or from the big bang to ",
 StyleBox["yend ",
  FontSlant->"Italic"],
 "for unbounded models.  "
}], "Text",ExpressionUUID->"d6a024c8-d0c6-4d4f-8eb3-c31d29607051"],

Cell[CellGroupData[{

Cell[BoxData["ymax"], "Input",
 CellLabel->"In[67]:=",ExpressionUUID->"0d0d2f30-4955-4528-b355-c666dd0ff455"],

Cell[BoxData["4"], "Output",
 CellChangeTimes->{{3.544186095410017*^9, 3.5441861243087387`*^9}, 
   3.54418623363449*^9, 3.544186846903516*^9, 3.607337327501944*^9, 
   3.924004121108132*^9, 3.924004279590733*^9, 3.924077940194903*^9, {
   3.9240780073542833`*^9, 3.9240780332579327`*^9}, 3.9240780880781317`*^9, 
   3.924078134523041*^9, 3.92407816620156*^9, 3.9240782307531157`*^9, 
   3.924082152874036*^9, 3.924082361069771*^9, {3.924189727574288*^9, 
   3.9241898520324574`*^9}, 3.986975961761402*^9, 3.986979672496673*^9, 
   3.9869797031772423`*^9},
 CellLabel->"Out[67]=",ExpressionUUID->"1b5de4bd-ab8a-4269-9ccb-3c6572e0e793"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"yu", "[", "z_", "]"}], ":=", 
  RowBox[{"yend", " ", "z"}]}]], "Input",
 CellLabel->"In[68]:=",ExpressionUUID->"bb9b0085-e881-4d36-8f51-6934ee6caa7a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"xu", "[", "z_", "]"}], ":=", 
  RowBox[{"x", "[", 
   RowBox[{"yu", "[", "z", "]"}], "]"}]}]], "Input",
 CellLabel->"In[69]:=",ExpressionUUID->"a96938f1-6062-457a-9ef0-6e390803eadc"],

Cell[BoxData[
 RowBox[{
  RowBox[{"yb", "[", "z_", "]"}], ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"z", "<", ".5"}], ",", 
    RowBox[{"ymax", 
     RowBox[{"(", 
      RowBox[{"2", " ", "z"}], ")"}]}], ",", 
    RowBox[{"ymax", " ", 
     RowBox[{"(", 
      RowBox[{"2", 
       RowBox[{"(", 
        RowBox[{"1", "-", "z"}], ")"}]}], ")"}]}]}], "]"}]}]], "Input",
 CellLabel->"In[70]:=",ExpressionUUID->"3a19a1ce-7678-4c27-ba60-b84f42c618f4"],

Cell[BoxData[
 RowBox[{
  RowBox[{"xb", "[", "z_", "]"}], ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"z", "<", ".5"}], ",", 
    RowBox[{"x", "[", 
     RowBox[{"yb", "[", "z", "]"}], "]"}], ",", 
    RowBox[{
     RowBox[{"2", 
      RowBox[{"x", "[", "ymax", "]"}]}], "-", 
     RowBox[{"x", "[", 
      RowBox[{"yb", "[", "z", "]"}], "]"}]}]}], "]"}]}]], "Input",
 CellLabel->"In[71]:=",ExpressionUUID->"91351a5e-c971-4447-8fb2-1e78d3a6b83b"],

Cell[BoxData[
 RowBox[{
  RowBox[{"yp", "[", "z_", "]"}], ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", " ", 
    RowBox[{"yb", "[", "z", "]"}], ",", 
    RowBox[{"yu", "[", "z", "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.544186574654521*^9, 3.544186575319253*^9}, {
  3.924004024596204*^9, 3.9240040255612707`*^9}},
 CellLabel->"In[72]:=",ExpressionUUID->"26313eba-63e8-46f6-bdd6-38b1f862e472"],

Cell[BoxData[
 RowBox[{
  RowBox[{"xp", "[", "z_", "]"}], ":=", 
  RowBox[{"If", "[", 
   RowBox[{
    RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", " ", 
    RowBox[{"xb", "[", "z", "]"}], ",", 
    RowBox[{"xu", "[", "z", "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.54418657634383*^9, 3.544186577086484*^9}, {
  3.924004030167987*^9, 3.924004030493576*^9}},
 CellLabel->"In[73]:=",ExpressionUUID->"9bff297a-87a1-4c90-90e0-ab1fb1f3c99e"],

Cell[BoxData[
 RowBox[{"plotxy", ":=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"xp", "[", "z", "]"}], ",", 
      RowBox[{"yp", "[", "z", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"z", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"AxesLabel", "->", 
     RowBox[{"{", 
      RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
    RowBox[{"DisplayFunction", "->", "Identity"}]}], "]"}]}]], "Input",
 CellLabel->"In[74]:=",ExpressionUUID->"2928ea34-ab48-4cb9-95d7-698db7d8b0b9"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plotxy", ",", 
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Dashing", "[", 
       RowBox[{"{", ".01", "}"}], "]"}], ",", " ", 
      RowBox[{"Line", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "[", "1", "]"}], ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "[", "1", "]"}], ",", "1"}], "}"}]}], "}"}], "]"}]}], 
     "}"}], "]"}], ",", 
   RowBox[{"DisplayFunction", "->", "$DisplayFunction"}]}], "]"}]], "Input",
 CellLabel->"In[75]:=",ExpressionUUID->"c46e22bf-1fb2-45e7-9bb7-4c3e87aa53d5"],

Cell[BoxData[
 GraphicsBox[{
   InterpretationBox[{
     TagBox[{{{}, {}, 
        TagBox[
         {RGBColor[0.24, 0.6, 0.8], AbsoluteThickness[2], Opacity[1.], 
          FaceForm[Opacity[0.3]], LineBox[CompressedData["
1:eJxd1Xs01GkYB/CZMZefcflV25lidyprUlZtZJfNyvOsioRWKKtakbbIsSYh
lA2hspOEmoou65a7rKEakdal0WXUIJnaaWUidnKZRIS2PTt/7XvOc57znOf5
fv59jXaGuP9Eo1AoAR/r3/7/J5efNihy0YGd7zktRw08ITRKaTXkYayZd0EI
w7187m0L+CQqu9TZ5Gf4nYx8ecbjW80+DHoud8mcOGvARry9pn9hFFiF925c
ONdJc/8LjMdJvJhOHpr7BAhyk7i0fdgCUvmqP65nJYGnKLZVW7VVkz8O/R90
A26V+GnyKUCsp/eamQRp8hlgZmVr538tGOYHFaxbtO009PzebS5Zztd4Z8CS
Q3xl9zpU452DPxmf5G2VRWq8i2AVvfRkO/2IxsuFUgtr6+d9R2B15Q6WpzgX
WPyd0qXSBI2fB13q81V2uUkaPx/aS1saZ8KSNX4hbBELvdy+P6XxyyDOVtmy
MPC8xheBo0eEU+jsTNCxeWtawRfB4iXVu2k1mSCKbEuW3xFBzlk37wezLwBj
7KTL8ogqEJqXlpQ/uAQFw8RDmawaVAcCesdqs2Hw5eRjruAG3N1WJb/TWgDR
DxXKypk6SDlIb75lXAkt7h6U4ZFG4LofS9XxrgTfDua23wybIL/G2f9wSiWM
eYmr3Nc0wdmBvs8UE5XA+9EoqCqjCezfPujly0QQu2eoI9q6GXjRN8yrk6vB
6tCvxYyYO1DoEnTWgiOGnLyGLZ8Sd0EtubVe5049xLyzLHVY0AqrdNn3HUQS
ULg6nmr4tR2EDk0hNMojML+X1icc7YLCG5YFURkd8PDq3iQq9Tk02zw6euqm
HOb48S+32r6Abt3Hhvzh57Cgv2XoB5US6L6WfS78F7ApL7CnJLgPdtN3W2ov
fgk+B54oDRcPQNDI9SuLyvsgLeZEskeZCppe8UKfLhuApVp+2c1fDEGJtTen
WqEC11aqsnBwGMqfLvmb7zMEBtZjkg5zNewx9uDImSMgM5WJ9FPeAK90j2Kj
pxpst9d2KZWjUJTjxl139A2IC+65XeaNAWdd6LYfWkeBF3tbGnB4HBpmcb2K
p95C9j75RU7DOzBxTAtOXT0OeoPy8EdGk3D+q/n82tB3YFo0kGga/R7sPiw5
Ij0xAd6x3Fe7a6bgWSl9113xJIR3W2+qmDcDV180HreRvoe87eG63ds/QFUo
9xt25xQ06PhROl0pmGGd15w4PA1FX76fbpik4MqFaWaykRnIVOTnGmdRMfaQ
/+EwKgVpns1SbWcanu89mpmTQEH7wODTdWM0nPxc5hozRUGBjgmj94wWyuYH
3245SEWX/e7v0tfSMS1G7ms2TEVOQGqGqJeOEdIqC4YvDa9PbX6Qk8hAsVf7
PfVjGj7acvK4wUom0l3pnmVrtPDNavkgtY2J8WbipRbXtFAwcLaTH8dCATt0
We8iOq7gBQa68wgc/PqzlrQkOhZN1wtpjQSGd+nli9V0zMoSOK3ap40HcooN
azYzUCgJzH0xh43l51jzxusYqC0Wx8/5g41Jx1rXruMx0dM5fXrWzzr49/e3
elhxTCwIjDGK19dFfbbL1009TNzv+lQSdVMXKdHuCgc7FjrWbTBQ+uih2upZ
R/RpFo5Kv2NyZ/Tw6l4/2xA1C+fllUQ8K9LH4pFly4scCVxLmKbqu5D4bOZ6
nLETgca7Xm1WbCRRn8IVZ20gkFZ/5dOyTST+lHHwcYorgfUHTK64epG4SjFR
u8+DQNs+4zrBThJV3SucV+4g0LJ5gYodRSKLqxN5MYLAzxPmrmflkzhFLd44
N5JA6l9tep0FJK4JE2onRxH417fpbfnFJKqDdy0LP0TgxTezdzhUkLg5zdJk
QzyBhv5kROJNEtML770ePEHgbHt2rlYbiUItJdv/JIHDF1oC2zpINPeJ1+1M
JVA6cWxFzpOPvs+O2rp0AgUVrBp7BYmJe2c6BOcIJIwYsrh+EstKi+ZRMgns
O9Qo3KQi0Sh8YvX+LAKbnxz50WiIREOJ4BvvSwQmpNL660dJbA9+onv/MoH+
qttlqeMkOtxvb7LLJtB+fVyY7+THfEiif0UOgUa5aGM+TaKg+/VL4zwC//tt
ZuE/YRnuxw==
           "]]},
         Annotation[#, "Charting`Private`Tag#1"]& ]}, {}},
      {"WolframDynamicHighlight", <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>}], 
     StyleBox[
      DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {}, 
        Slot["HighlightElements"], 
        Slot["LayoutOptions"], 
        Slot["Meta"], 
        Charting`HighlightActionFunction["DynamicHighlight", {{{{}, {}, 
            Annotation[{
              Directive[
               Opacity[1.], 
               AbsoluteThickness[2], 
               FaceForm[
                Opacity[0.3]], 
               RGBColor[0.24, 0.6, 0.8]], 
              Line[CompressedData["
1:eJxd1Xs01GkYB/CZMZefcflV25lidyprUlZtZJfNyvOsioRWKKtakbbIsSYh
lA2hspOEmoou65a7rKEakdal0WXUIJnaaWUidnKZRIS2PTt/7XvOc57znOf5
fv59jXaGuP9Eo1AoAR/r3/7/J5efNihy0YGd7zktRw08ITRKaTXkYayZd0EI
w7187m0L+CQqu9TZ5Gf4nYx8ecbjW80+DHoud8mcOGvARry9pn9hFFiF925c
ONdJc/8LjMdJvJhOHpr7BAhyk7i0fdgCUvmqP65nJYGnKLZVW7VVkz8O/R90
A26V+GnyKUCsp/eamQRp8hlgZmVr538tGOYHFaxbtO009PzebS5Zztd4Z8CS
Q3xl9zpU452DPxmf5G2VRWq8i2AVvfRkO/2IxsuFUgtr6+d9R2B15Q6WpzgX
WPyd0qXSBI2fB13q81V2uUkaPx/aS1saZ8KSNX4hbBELvdy+P6XxyyDOVtmy
MPC8xheBo0eEU+jsTNCxeWtawRfB4iXVu2k1mSCKbEuW3xFBzlk37wezLwBj
7KTL8ogqEJqXlpQ/uAQFw8RDmawaVAcCesdqs2Hw5eRjruAG3N1WJb/TWgDR
DxXKypk6SDlIb75lXAkt7h6U4ZFG4LofS9XxrgTfDua23wybIL/G2f9wSiWM
eYmr3Nc0wdmBvs8UE5XA+9EoqCqjCezfPujly0QQu2eoI9q6GXjRN8yrk6vB
6tCvxYyYO1DoEnTWgiOGnLyGLZ8Sd0EtubVe5049xLyzLHVY0AqrdNn3HUQS
ULg6nmr4tR2EDk0hNMojML+X1icc7YLCG5YFURkd8PDq3iQq9Tk02zw6euqm
HOb48S+32r6Abt3Hhvzh57Cgv2XoB5US6L6WfS78F7ApL7CnJLgPdtN3W2ov
fgk+B54oDRcPQNDI9SuLyvsgLeZEskeZCppe8UKfLhuApVp+2c1fDEGJtTen
WqEC11aqsnBwGMqfLvmb7zMEBtZjkg5zNewx9uDImSMgM5WJ9FPeAK90j2Kj
pxpst9d2KZWjUJTjxl139A2IC+65XeaNAWdd6LYfWkeBF3tbGnB4HBpmcb2K
p95C9j75RU7DOzBxTAtOXT0OeoPy8EdGk3D+q/n82tB3YFo0kGga/R7sPiw5
Ij0xAd6x3Fe7a6bgWSl9113xJIR3W2+qmDcDV180HreRvoe87eG63ds/QFUo
9xt25xQ06PhROl0pmGGd15w4PA1FX76fbpik4MqFaWaykRnIVOTnGmdRMfaQ
/+EwKgVpns1SbWcanu89mpmTQEH7wODTdWM0nPxc5hozRUGBjgmj94wWyuYH
3245SEWX/e7v0tfSMS1G7ms2TEVOQGqGqJeOEdIqC4YvDa9PbX6Qk8hAsVf7
PfVjGj7acvK4wUom0l3pnmVrtPDNavkgtY2J8WbipRbXtFAwcLaTH8dCATt0
We8iOq7gBQa68wgc/PqzlrQkOhZN1wtpjQSGd+nli9V0zMoSOK3ap40HcooN
azYzUCgJzH0xh43l51jzxusYqC0Wx8/5g41Jx1rXruMx0dM5fXrWzzr49/e3
elhxTCwIjDGK19dFfbbL1009TNzv+lQSdVMXKdHuCgc7FjrWbTBQ+uih2upZ
R/RpFo5Kv2NyZ/Tw6l4/2xA1C+fllUQ8K9LH4pFly4scCVxLmKbqu5D4bOZ6
nLETgca7Xm1WbCRRn8IVZ20gkFZ/5dOyTST+lHHwcYorgfUHTK64epG4SjFR
u8+DQNs+4zrBThJV3SucV+4g0LJ5gYodRSKLqxN5MYLAzxPmrmflkzhFLd44
N5JA6l9tep0FJK4JE2onRxH417fpbfnFJKqDdy0LP0TgxTezdzhUkLg5zdJk
QzyBhv5kROJNEtML770ePEHgbHt2rlYbiUItJdv/JIHDF1oC2zpINPeJ1+1M
JVA6cWxFzpOPvs+O2rp0AgUVrBp7BYmJe2c6BOcIJIwYsrh+EstKi+ZRMgns
O9Qo3KQi0Sh8YvX+LAKbnxz50WiIREOJ4BvvSwQmpNL660dJbA9+onv/MoH+
qttlqeMkOtxvb7LLJtB+fVyY7+THfEiif0UOgUa5aGM+TaKg+/VL4zwC//tt
ZuE/YRnuxw==
               "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
         "HighlightElements" -> <|
           "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
          "LayoutOptions" -> <|
           "PanelPlotLayout" -> <||>, 
            "PlotRange" -> {{0., 3.0806348765583547`}, {0., 4.}}, 
            "Frame" -> {{False, False}, {False, False}}, 
            "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 432}, 
            "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
            Automatic, "DefaultStyle" -> {
              Directive[
               Opacity[1.], 
               AbsoluteThickness[2], 
               FaceForm[
                Opacity[0.3]], 
               RGBColor[0.24, 0.6, 0.8]]}, 
            "HighlightLabelingFunctions" -> <|
             "CoordinatesToolOptions" -> Identity, 
              "ScalingFunctions" -> {{Identity, Identity}, {
                Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
            False|>, 
          "Meta" -> <|
           "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
            ParametricPlot, "GroupHighlight" -> False|>|>]]& )[<|
        "HighlightElements" -> <|
          "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
         "LayoutOptions" -> <|
          "PanelPlotLayout" -> <||>, 
           "PlotRange" -> {{0., 3.0806348765583547`}, {0., 4.}}, 
           "Frame" -> {{False, False}, {False, False}}, 
           "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 432}, 
           "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
           Automatic, "DefaultStyle" -> {
             Directive[
              Opacity[1.], 
              AbsoluteThickness[2], 
              FaceForm[
               Opacity[0.3]], 
              RGBColor[0.24, 0.6, 0.8]]}, 
           "HighlightLabelingFunctions" -> <|
            "CoordinatesToolOptions" -> Identity, 
             "ScalingFunctions" -> {{Identity, Identity}, {
               Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
           False|>, 
         "Meta" -> <|
          "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
           ParametricPlot, "GroupHighlight" -> False|>|>],
       ImageSizeCache->{{4.503599627370496*^15, -4.503599627370496*^15}, {
        4.503599627370496*^15, -4.503599627370496*^15}}],
      Selectable->False]},
    Annotation[{{{{}, {}, 
        Annotation[{
          Directive[
           Opacity[1.], 
           AbsoluteThickness[2], 
           FaceForm[
            Opacity[0.3]], 
           RGBColor[0.24, 0.6, 0.8]], 
          Line[CompressedData["
1:eJxd1Xs01GkYB/CZMZefcflV25lidyprUlZtZJfNyvOsioRWKKtakbbIsSYh
lA2hspOEmoou65a7rKEakdal0WXUIJnaaWUidnKZRIS2PTt/7XvOc57znOf5
fv59jXaGuP9Eo1AoAR/r3/7/J5efNihy0YGd7zktRw08ITRKaTXkYayZd0EI
w7187m0L+CQqu9TZ5Gf4nYx8ecbjW80+DHoud8mcOGvARry9pn9hFFiF925c
ONdJc/8LjMdJvJhOHpr7BAhyk7i0fdgCUvmqP65nJYGnKLZVW7VVkz8O/R90
A26V+GnyKUCsp/eamQRp8hlgZmVr538tGOYHFaxbtO009PzebS5Zztd4Z8CS
Q3xl9zpU452DPxmf5G2VRWq8i2AVvfRkO/2IxsuFUgtr6+d9R2B15Q6WpzgX
WPyd0qXSBI2fB13q81V2uUkaPx/aS1saZ8KSNX4hbBELvdy+P6XxyyDOVtmy
MPC8xheBo0eEU+jsTNCxeWtawRfB4iXVu2k1mSCKbEuW3xFBzlk37wezLwBj
7KTL8ogqEJqXlpQ/uAQFw8RDmawaVAcCesdqs2Hw5eRjruAG3N1WJb/TWgDR
DxXKypk6SDlIb75lXAkt7h6U4ZFG4LofS9XxrgTfDua23wybIL/G2f9wSiWM
eYmr3Nc0wdmBvs8UE5XA+9EoqCqjCezfPujly0QQu2eoI9q6GXjRN8yrk6vB
6tCvxYyYO1DoEnTWgiOGnLyGLZ8Sd0EtubVe5049xLyzLHVY0AqrdNn3HUQS
ULg6nmr4tR2EDk0hNMojML+X1icc7YLCG5YFURkd8PDq3iQq9Tk02zw6euqm
HOb48S+32r6Abt3Hhvzh57Cgv2XoB5US6L6WfS78F7ApL7CnJLgPdtN3W2ov
fgk+B54oDRcPQNDI9SuLyvsgLeZEskeZCppe8UKfLhuApVp+2c1fDEGJtTen
WqEC11aqsnBwGMqfLvmb7zMEBtZjkg5zNewx9uDImSMgM5WJ9FPeAK90j2Kj
pxpst9d2KZWjUJTjxl139A2IC+65XeaNAWdd6LYfWkeBF3tbGnB4HBpmcb2K
p95C9j75RU7DOzBxTAtOXT0OeoPy8EdGk3D+q/n82tB3YFo0kGga/R7sPiw5
Ij0xAd6x3Fe7a6bgWSl9113xJIR3W2+qmDcDV180HreRvoe87eG63ds/QFUo
9xt25xQ06PhROl0pmGGd15w4PA1FX76fbpik4MqFaWaykRnIVOTnGmdRMfaQ
/+EwKgVpns1SbWcanu89mpmTQEH7wODTdWM0nPxc5hozRUGBjgmj94wWyuYH
3245SEWX/e7v0tfSMS1G7ms2TEVOQGqGqJeOEdIqC4YvDa9PbX6Qk8hAsVf7
PfVjGj7acvK4wUom0l3pnmVrtPDNavkgtY2J8WbipRbXtFAwcLaTH8dCATt0
We8iOq7gBQa68wgc/PqzlrQkOhZN1wtpjQSGd+nli9V0zMoSOK3ap40HcooN
azYzUCgJzH0xh43l51jzxusYqC0Wx8/5g41Jx1rXruMx0dM5fXrWzzr49/e3
elhxTCwIjDGK19dFfbbL1009TNzv+lQSdVMXKdHuCgc7FjrWbTBQ+uih2upZ
R/RpFo5Kv2NyZ/Tw6l4/2xA1C+fllUQ8K9LH4pFly4scCVxLmKbqu5D4bOZ6
nLETgca7Xm1WbCRRn8IVZ20gkFZ/5dOyTST+lHHwcYorgfUHTK64epG4SjFR
u8+DQNs+4zrBThJV3SucV+4g0LJ5gYodRSKLqxN5MYLAzxPmrmflkzhFLd44
N5JA6l9tep0FJK4JE2onRxH417fpbfnFJKqDdy0LP0TgxTezdzhUkLg5zdJk
QzyBhv5kROJNEtML770ePEHgbHt2rlYbiUItJdv/JIHDF1oC2zpINPeJ1+1M
JVA6cWxFzpOPvs+O2rp0AgUVrBp7BYmJe2c6BOcIJIwYsrh+EstKi+ZRMgns
O9Qo3KQi0Sh8YvX+LAKbnxz50WiIREOJ4BvvSwQmpNL660dJbA9+onv/MoH+
qttlqeMkOtxvb7LLJtB+fVyY7+THfEiif0UOgUa5aGM+TaKg+/VL4zwC//tt
ZuE/YRnuxw==
           "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
     "HighlightElements" -> <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
      "LayoutOptions" -> <|
       "PanelPlotLayout" -> <||>, 
        "PlotRange" -> {{0., 3.0806348765583547`}, {0., 4.}}, 
        "Frame" -> {{False, False}, {False, False}}, "AxesOrigin" -> {0, 0}, 
        "ImageSize" -> {360, 432}, "Axes" -> {True, True}, "LabelStyle" -> {},
         "AspectRatio" -> Automatic, "DefaultStyle" -> {
          Directive[
           Opacity[1.], 
           AbsoluteThickness[2], 
           FaceForm[
            Opacity[0.3]], 
           RGBColor[0.24, 0.6, 0.8]]}, 
        "HighlightLabelingFunctions" -> <|
         "CoordinatesToolOptions" -> Identity, 
          "ScalingFunctions" -> {{Identity, Identity}, {
            Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> False|>, 
      "Meta" -> <|
       "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
        ParametricPlot, "GroupHighlight" -> False|>|>, "DynamicHighlight"]], 
   {Dashing[{0.01}], 
    LineBox[{{0.8440195317660885, 0}, {0.8440195317660885, 1}}]}},
  Axes->{True, True},
  AxesLabel->{
    FormBox[
     TagBox["\"x\"", HoldForm], TraditionalForm], 
    FormBox[
     TagBox["\"y\"", HoldForm], TraditionalForm]},
  AxesOrigin->{0, 0},
  DefaultBaseStyle->{"PlotGraphics", "Graphics"},
  DisplayFunction->Identity,
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None},
  PlotInteractivity:>True,
  PlotRange->{{0., 3.0806348765583547`}, {0., 4.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.05], 
     Scaled[0.05]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.543849781562581*^9, 3.544186071922852*^9, 3.5441861320291643`*^9, 
   3.544186242062565*^9, 3.54418685958696*^9, 3.6073373437867517`*^9, 
   3.9240041211760187`*^9, 3.924004280563801*^9, 3.9240779406747*^9, {
   3.924078007739085*^9, 3.92407803366186*^9}, 3.924078088357408*^9, 
   3.924078134826106*^9, 3.924078167182392*^9, 3.9240782317368193`*^9, 
   3.924082153856769*^9, 3.924082361405727*^9, {3.924189727997468*^9, 
   3.924189853305759*^9}, 3.986975961868144*^9, 3.986979672891365*^9, 
   3.986979703397184*^9},
 CellLabel->"Out[75]=",ExpressionUUID->"55df9539-1670-4e0d-9938-25dd799129fa"]
}, Open  ]],

Cell["The dashed vertical line is a the time of today. ", "Text",ExpressionUUID->"89487cb2-9e4d-4755-863c-0667bf96fad3"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Rescaling to get ",
 StyleBox["a(t) ",
  FontSlant->"Italic"],
 "as a function of ",
 StyleBox["t:",
  FontSlant->"Italic"]
}], "Subsection",ExpressionUUID->"b05aebc7-3a85-475a-a307-abb4cfeb34c2"],

Cell[TextData[{
 "The Hubble time ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["T", "0"], TraditionalForm]],ExpressionUUID->
  "3e0fb66a-16eb-45da-9519-baa51db1c329"],
 " (",
 StyleBox["T",
  FontWeight->"Bold"],
 StyleBox["h",
  FontWeight->"Bold"],
 StyleBox[") can be used to rescale the ",
  FontVariations->{"CompatibilityType"->0}],
 " dimensionless measure of the scale factor ",
 StyleBox["y ",
  FontSlant->"Italic"],
 StyleBox["and the dimensionless measure of the time ",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["x ",
  FontSlant->"Italic",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox[" that were convenient for computation to give the usual scale \
factor ",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["a(t) ",
  FontSlant->"Italic",
  FontVariations->{"CompatibilityType"->0}],
 StyleBox["as follows:",
  FontVariations->{"CompatibilityType"->0}]
}], "Text",ExpressionUUID->"c31f84d4-6ba7-48c0-8488-8715c7e8c77c"],

Cell[BoxData[
 RowBox[{
  RowBox[{"t", "[", "z_", "]"}], ":=", 
  RowBox[{"Th", " ", 
   RowBox[{"xp", "[", "z", "]"}]}]}]], "Input",
 CellLabel->"In[76]:=",ExpressionUUID->"a8adecf2-baff-4f25-bc88-9030de284db1"],

Cell[BoxData[
 RowBox[{"a0", ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Th", "/", 
     RowBox[{"Abs", "[", "\[CapitalOmega]k", "]"}]}], ")"}], "^", 
   RowBox[{"(", 
    RowBox[{"1", "/", "2"}], ")"}]}]}]], "Input",
 CellChangeTimes->{
  3.544186580488901*^9, {3.924004081111545*^9, 3.924004081659782*^9}},
 CellLabel->"In[77]:=",ExpressionUUID->"fea91ced-534f-4fec-ae6f-9f1d7d43cc1d"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a0", "//", "N"}]], "Input",
 CellLabel->"In[78]:=",ExpressionUUID->"9166b52d-376f-4387-89de-ff60f52333c6"],

Cell[BoxData["5.288267222564942`"], "Output",
 CellChangeTimes->{
  3.543849843359888*^9, 3.544186158011589*^9, 3.544186249720093*^9, 
   3.544186869759287*^9, 3.607337356038678*^9, 3.924004121201255*^9, 
   3.924004280618082*^9, 3.92407794071354*^9, {3.924078007779093*^9, 
   3.924078033703755*^9}, 3.924078088427004*^9, 3.924078134868387*^9, 
   3.92407816724632*^9, 3.9240782318078127`*^9, 3.9240821539228783`*^9, 
   3.924082361452883*^9, {3.9241897280369*^9, 3.924189853374434*^9}, 
   3.986975961886413*^9, 3.986979673004583*^9, 3.9869797034140253`*^9},
 CellLabel->"Out[78]=",ExpressionUUID->"69082cee-2e10-4572-ba9f-a3e669d42be4"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"a", "[", "z_", "]"}], ":=", 
  RowBox[{"a0", " ", 
   RowBox[{"yp", "[", "z", "]"}]}]}]], "Input",
 CellLabel->"In[79]:=",ExpressionUUID->"3e47d984-d093-44c9-a31d-b3a7ef6e01e7"],

Cell[TextData[{
 "We can then plot ",
 StyleBox["a(t)  vs  t .   ",
  FontSlant->"Italic"],
 StyleBox["(Its the same plot as above but on rescaled axes.)",
  FontVariations->{"CompatibilityType"->0}]
}], "Text",ExpressionUUID->"3055344a-aa2f-4051-ba38-b274fa9f8f8e"],

Cell[BoxData[
 RowBox[{"plotat", ":=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"t", "[", "z", "]"}], ",", 
      RowBox[{"a", "[", "z", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"z", ",", "0", ",", "1"}], "}"}], ",", 
    RowBox[{"AxesLabel", "->", 
     RowBox[{"{", 
      RowBox[{"\"\<t (Gyr)\>\"", ",", "\"\<a(t)\>\""}], "}"}]}], ",", 
    RowBox[{"DisplayFunction", "->", "Identity"}]}], "]"}]}]], "Input",
 CellLabel->"In[80]:=",ExpressionUUID->"42fe58a7-d80c-4896-91fb-372e56a7624b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plotat", ",", 
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Dashing", "[", 
       RowBox[{"{", ".01", "}"}], "]"}], ",", " ", 
      RowBox[{"Line", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Th", " ", 
            RowBox[{"x", "[", "1", "]"}]}], ",", "0"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Th", " ", 
            RowBox[{"x", "[", "1", "]"}]}], ",", "a0"}], "}"}]}], "}"}], 
       "]"}]}], "}"}], "]"}], ",", 
   RowBox[{"DisplayFunction", "->", "$DisplayFunction"}]}], "]"}]], "Input",
 CellLabel->"In[81]:=",ExpressionUUID->"bcadcf0a-2bb3-484c-8451-2579314c323c"],

Cell[BoxData[
 GraphicsBox[{
   InterpretationBox[{
     TagBox[{{{}, {}, 
        TagBox[
         {RGBColor[0.24, 0.6, 0.8], AbsoluteThickness[2], Opacity[1.], 
          FaceForm[Opacity[0.3]], LineBox[CompressedData["
1:eJxlznk01XkfB3CailtZ4p4h997hZssly13I+n3j/tJiCaMGYZBIafGQoZIW
ssxRZkyKYoylEN0oIhI1D6EZEYpGzOBas4/omemZ55z73/M953M+53M+n/fr
fNkBx9yCVklJSYX8U//r//ds985l9DmQ6/0yOrXfxZALjqYcraGvJXMyeXj5
p2Ed2WNkfiSeKff+e8I5fcqmYDBKss8kf67s8BaWxRH1iWcZD7J+JLaFB8Jr
xuMl93nEayctKKv3iuS+mJjGjM4kbUonR3qbi5WHS8nGPe5cnUsZkryINLn7
t98YypbkK4hrPZ2+UH1Lkq8hb3vDrPd8KCLn94eMii1rSdKTrygbhzsSr45c
qAwJD0oWSbwn5P663mDceyDxnpIf5keTVuU8kXgtpHDh8dWwoAaSytom4xbc
Ss4xhtrLDBslfhtZN6e1/tzzpxL/BYlKfGnEUm+S+O0k4Neh+062v0j8TpIl
7PKoo/dI/F6y77CnV8+tHmLuMj65tKaPmJ6vml1r9Zr01rdV8KL7CHNeMS4z
+A1Ry02zLfF7S5yztdPQ2keyAhk+mQb9ROow1tL9B4hLhMouqz8GSPU3MSVS
d4eJ2zLnnNqZIaJdPdgvUzND9n3mNt0tmiAf5Z6nZY3OkFjFRL+8kQlyX3EK
biqzJJ/1+NdjzEmyO+wLkW7kLJkz44hkEyeJc9zN6iPcOZJ6ZNUJC78p8sqz
PsO4Yp487Sqfuyk3TSYXy4rs2/4kqq9op7iCWfIF87csdfe/yMXC0NjltAVi
HyfeuXR1DfwCW/wb9T6Sf6tjQFQpB8fr4i8dt0jjZboL51OfEtyczagfY1fj
O922D0U+Kkj89o04hikLWoTfoZwxBh53t3wIL1qP4JVd6XyWBnyP0rrdl+WR
sV3eIytmMzIcXpT77NgIcx/7lLpyLdyPrR2IX1SCsfQrRaUBHbRbdv5slUnH
1+8mtyxY6CHBWnl1VcDn6Opk68uG6mN2OPp2wQZVMH5oP1/YuBXzHbTCfb2b
wONprtStNsYi1eB1LYeBMCGXSG80QalfTqTTdhbyohU2PTTi4mNd4djPiupw
mNJWkDfjIc5/3N/TVwOBdjGxsS58OK29espkOxs6LBvlez4CMI52zbTrbcbW
8Od+Sh6mEOptuc75ezOaWlxMmfvN8DBB4X3QkCZGjve3dDhuwyWzdSe9H2nh
UxN7kPDNsbfkmqdqqjZcA2eabTgWONHwaKwyXAdCyzPyjPWWGOHM3VHYoYtz
nzeoP1u2xC3TOuYuoy2QbxRmyr+zQsga5uhuJT08SzHoTBJZI143dURmRA97
fglspl22weREelpBKwdzzMimCC+Cu5/dPfZXrj58HuY6rsgDx/v3FzEjDJCJ
e+ljJ4Gz3jmXWV9txW3R0uvpLkCsv3j7tYEhfNIs83+3sEVW9rujkcpGYHjQ
aMIMWxQf+C24bdwIhlE9huoztnh70sx30dcYx+2bg/Kd7ZD8bYtv5EtjLCvx
8vUK7DC0W/yNm5UJ5h402rOl7RG03/mQhsgELpUnfi/ytMdB57MyHWwuThfq
1tYU24M+bl4/kMLFm3In5zXSQox22zgN/M1F9X9Upll7hAgJe7FgFcrD4pcp
p91vCBE6rL3rTCcP4lyrnstTQlwpa44so/jI2eBC27CNwgH2Yo2ViI/s5JZW
D3MKd/I0ElTL+TCsKo7MsaCwoO3oulDBh8FSwVmuNYWLnDxxSRUfVxcGj+yz
o5DPc6Uz6vn44Gp994YjhSGqNGz5BR9TvR1jdH8KgaFBGpUTfBwcoed4B1Ao
mbwykTbFx80EttVPgRTmj9ZWhk3z8VLj1lajgxQu/EvZSXv+n/8NPnrvcJhC
3umG6PQVPlTYrRMRkRT+SGV1htMEeF9qpF5zkoL+xp3ZzusFyFydZvEpikL4
9xGHOHICJKTs1UuKobDqWuunQUUB2nYIxVlnKWzOjTZwVRWg+/iI10AchVDN
giUDNQE0nbMrtc5TKC9ob5RlCjCW7q1TepGCXbGu5xN1AXx5wRaz8RSSDdy1
brAFUJhNNxVcotBRFjsdpSkAEYvVohMpqJkU17hrC1BlGCiuTaIQUNEVb6Qr
wH8Bu0TsVw==
           "]]},
         Annotation[#, "Charting`Private`Tag#1"]& ]}, {}},
      {"WolframDynamicHighlight", <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>}], 
     StyleBox[
      DynamicBox[(Charting`HighlightActionBox["DynamicHighlight", {}, 
        Slot["HighlightElements"], 
        Slot["LayoutOptions"], 
        Slot["Meta"], 
        Charting`HighlightActionFunction["DynamicHighlight", {{{{}, {}, 
            Annotation[{
              Directive[
               Opacity[1.], 
               AbsoluteThickness[2], 
               FaceForm[
                Opacity[0.3]], 
               RGBColor[0.24, 0.6, 0.8]], 
              Line[CompressedData["
1:eJxlznk01XkfB3CailtZ4p4h997hZssly13I+n3j/tJiCaMGYZBIafGQoZIW
ssxRZkyKYoylEN0oIhI1D6EZEYpGzOBas4/omemZ55z73/M953M+53M+n/fr
fNkBx9yCVklJSYX8U//r//ds985l9DmQ6/0yOrXfxZALjqYcraGvJXMyeXj5
p2Ed2WNkfiSeKff+e8I5fcqmYDBKss8kf67s8BaWxRH1iWcZD7J+JLaFB8Jr
xuMl93nEayctKKv3iuS+mJjGjM4kbUonR3qbi5WHS8nGPe5cnUsZkryINLn7
t98YypbkK4hrPZ2+UH1Lkq8hb3vDrPd8KCLn94eMii1rSdKTrygbhzsSr45c
qAwJD0oWSbwn5P663mDceyDxnpIf5keTVuU8kXgtpHDh8dWwoAaSytom4xbc
Ss4xhtrLDBslfhtZN6e1/tzzpxL/BYlKfGnEUm+S+O0k4Neh+062v0j8TpIl
7PKoo/dI/F6y77CnV8+tHmLuMj65tKaPmJ6vml1r9Zr01rdV8KL7CHNeMS4z
+A1Ry02zLfF7S5yztdPQ2keyAhk+mQb9ROow1tL9B4hLhMouqz8GSPU3MSVS
d4eJ2zLnnNqZIaJdPdgvUzND9n3mNt0tmiAf5Z6nZY3OkFjFRL+8kQlyX3EK
biqzJJ/1+NdjzEmyO+wLkW7kLJkz44hkEyeJc9zN6iPcOZJ6ZNUJC78p8sqz
PsO4Yp487Sqfuyk3TSYXy4rs2/4kqq9op7iCWfIF87csdfe/yMXC0NjltAVi
HyfeuXR1DfwCW/wb9T6Sf6tjQFQpB8fr4i8dt0jjZboL51OfEtyczagfY1fj
O922D0U+Kkj89o04hikLWoTfoZwxBh53t3wIL1qP4JVd6XyWBnyP0rrdl+WR
sV3eIytmMzIcXpT77NgIcx/7lLpyLdyPrR2IX1SCsfQrRaUBHbRbdv5slUnH
1+8mtyxY6CHBWnl1VcDn6Opk68uG6mN2OPp2wQZVMH5oP1/YuBXzHbTCfb2b
wONprtStNsYi1eB1LYeBMCGXSG80QalfTqTTdhbyohU2PTTi4mNd4djPiupw
mNJWkDfjIc5/3N/TVwOBdjGxsS58OK29espkOxs6LBvlez4CMI52zbTrbcbW
8Od+Sh6mEOptuc75ezOaWlxMmfvN8DBB4X3QkCZGjve3dDhuwyWzdSe9H2nh
UxN7kPDNsbfkmqdqqjZcA2eabTgWONHwaKwyXAdCyzPyjPWWGOHM3VHYoYtz
nzeoP1u2xC3TOuYuoy2QbxRmyr+zQsga5uhuJT08SzHoTBJZI143dURmRA97
fglspl22weREelpBKwdzzMimCC+Cu5/dPfZXrj58HuY6rsgDx/v3FzEjDJCJ
e+ljJ4Gz3jmXWV9txW3R0uvpLkCsv3j7tYEhfNIs83+3sEVW9rujkcpGYHjQ
aMIMWxQf+C24bdwIhlE9huoztnh70sx30dcYx+2bg/Kd7ZD8bYtv5EtjLCvx
8vUK7DC0W/yNm5UJ5h402rOl7RG03/mQhsgELpUnfi/ytMdB57MyHWwuThfq
1tYU24M+bl4/kMLFm3In5zXSQox22zgN/M1F9X9Upll7hAgJe7FgFcrD4pcp
p91vCBE6rL3rTCcP4lyrnstTQlwpa44so/jI2eBC27CNwgH2Yo2ViI/s5JZW
D3MKd/I0ElTL+TCsKo7MsaCwoO3oulDBh8FSwVmuNYWLnDxxSRUfVxcGj+yz
o5DPc6Uz6vn44Gp994YjhSGqNGz5BR9TvR1jdH8KgaFBGpUTfBwcoed4B1Ao
mbwykTbFx80EttVPgRTmj9ZWhk3z8VLj1lajgxQu/EvZSXv+n/8NPnrvcJhC
3umG6PQVPlTYrRMRkRT+SGV1htMEeF9qpF5zkoL+xp3ZzusFyFydZvEpikL4
9xGHOHICJKTs1UuKobDqWuunQUUB2nYIxVlnKWzOjTZwVRWg+/iI10AchVDN
giUDNQE0nbMrtc5TKC9ob5RlCjCW7q1TepGCXbGu5xN1AXx5wRaz8RSSDdy1
brAFUJhNNxVcotBRFjsdpSkAEYvVohMpqJkU17hrC1BlGCiuTaIQUNEVb6Qr
wH8Bu0TsVw==
               "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
         "HighlightElements" -> <|
           "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
          "LayoutOptions" -> <|
           "PanelPlotLayout" -> <||>, 
            "PlotRange" -> {{0., 43.07607738821883}, {0., 
             21.153068890259767`}}, 
            "Frame" -> {{False, False}, {False, False}}, 
            "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 432}, 
            "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
            Automatic, "DefaultStyle" -> {
              Directive[
               Opacity[1.], 
               AbsoluteThickness[2], 
               FaceForm[
                Opacity[0.3]], 
               RGBColor[0.24, 0.6, 0.8]]}, 
            "HighlightLabelingFunctions" -> <|
             "CoordinatesToolOptions" -> Identity, 
              "ScalingFunctions" -> {{Identity, Identity}, {
                Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
            False|>, 
          "Meta" -> <|
           "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
            ParametricPlot, "GroupHighlight" -> False|>|>]]& )[<|
        "HighlightElements" -> <|
          "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
         "LayoutOptions" -> <|
          "PanelPlotLayout" -> <||>, 
           "PlotRange" -> {{0., 43.07607738821883}, {0., 
            21.153068890259767`}}, 
           "Frame" -> {{False, False}, {False, False}}, 
           "AxesOrigin" -> {0, 0}, "ImageSize" -> {360, 432}, 
           "Axes" -> {True, True}, "LabelStyle" -> {}, "AspectRatio" -> 
           Automatic, "DefaultStyle" -> {
             Directive[
              Opacity[1.], 
              AbsoluteThickness[2], 
              FaceForm[
               Opacity[0.3]], 
              RGBColor[0.24, 0.6, 0.8]]}, 
           "HighlightLabelingFunctions" -> <|
            "CoordinatesToolOptions" -> Identity, 
             "ScalingFunctions" -> {{Identity, Identity}, {
               Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> 
           False|>, 
         "Meta" -> <|
          "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
           ParametricPlot, "GroupHighlight" -> False|>|>],
       ImageSizeCache->{{4.503599627370496*^15, -4.503599627370496*^15}, {
        4.503599627370496*^15, -4.503599627370496*^15}}],
      Selectable->False]},
    Annotation[{{{{}, {}, 
        Annotation[{
          Directive[
           Opacity[1.], 
           AbsoluteThickness[2], 
           FaceForm[
            Opacity[0.3]], 
           RGBColor[0.24, 0.6, 0.8]], 
          Line[CompressedData["
1:eJxlznk01XkfB3CailtZ4p4h997hZssly13I+n3j/tJiCaMGYZBIafGQoZIW
ssxRZkyKYoylEN0oIhI1D6EZEYpGzOBas4/omemZ55z73/M953M+53M+n/fr
fNkBx9yCVklJSYX8U//r//ds985l9DmQ6/0yOrXfxZALjqYcraGvJXMyeXj5
p2Ed2WNkfiSeKff+e8I5fcqmYDBKss8kf67s8BaWxRH1iWcZD7J+JLaFB8Jr
xuMl93nEayctKKv3iuS+mJjGjM4kbUonR3qbi5WHS8nGPe5cnUsZkryINLn7
t98YypbkK4hrPZ2+UH1Lkq8hb3vDrPd8KCLn94eMii1rSdKTrygbhzsSr45c
qAwJD0oWSbwn5P663mDceyDxnpIf5keTVuU8kXgtpHDh8dWwoAaSytom4xbc
Ss4xhtrLDBslfhtZN6e1/tzzpxL/BYlKfGnEUm+S+O0k4Neh+062v0j8TpIl
7PKoo/dI/F6y77CnV8+tHmLuMj65tKaPmJ6vml1r9Zr01rdV8KL7CHNeMS4z
+A1Ry02zLfF7S5yztdPQ2keyAhk+mQb9ROow1tL9B4hLhMouqz8GSPU3MSVS
d4eJ2zLnnNqZIaJdPdgvUzND9n3mNt0tmiAf5Z6nZY3OkFjFRL+8kQlyX3EK
biqzJJ/1+NdjzEmyO+wLkW7kLJkz44hkEyeJc9zN6iPcOZJ6ZNUJC78p8sqz
PsO4Yp487Sqfuyk3TSYXy4rs2/4kqq9op7iCWfIF87csdfe/yMXC0NjltAVi
HyfeuXR1DfwCW/wb9T6Sf6tjQFQpB8fr4i8dt0jjZboL51OfEtyczagfY1fj
O922D0U+Kkj89o04hikLWoTfoZwxBh53t3wIL1qP4JVd6XyWBnyP0rrdl+WR
sV3eIytmMzIcXpT77NgIcx/7lLpyLdyPrR2IX1SCsfQrRaUBHbRbdv5slUnH
1+8mtyxY6CHBWnl1VcDn6Opk68uG6mN2OPp2wQZVMH5oP1/YuBXzHbTCfb2b
wONprtStNsYi1eB1LYeBMCGXSG80QalfTqTTdhbyohU2PTTi4mNd4djPiupw
mNJWkDfjIc5/3N/TVwOBdjGxsS58OK29espkOxs6LBvlez4CMI52zbTrbcbW
8Od+Sh6mEOptuc75ezOaWlxMmfvN8DBB4X3QkCZGjve3dDhuwyWzdSe9H2nh
UxN7kPDNsbfkmqdqqjZcA2eabTgWONHwaKwyXAdCyzPyjPWWGOHM3VHYoYtz
nzeoP1u2xC3TOuYuoy2QbxRmyr+zQsga5uhuJT08SzHoTBJZI143dURmRA97
fglspl22weREelpBKwdzzMimCC+Cu5/dPfZXrj58HuY6rsgDx/v3FzEjDJCJ
e+ljJ4Gz3jmXWV9txW3R0uvpLkCsv3j7tYEhfNIs83+3sEVW9rujkcpGYHjQ
aMIMWxQf+C24bdwIhlE9huoztnh70sx30dcYx+2bg/Kd7ZD8bYtv5EtjLCvx
8vUK7DC0W/yNm5UJ5h402rOl7RG03/mQhsgELpUnfi/ytMdB57MyHWwuThfq
1tYU24M+bl4/kMLFm3In5zXSQox22zgN/M1F9X9Upll7hAgJe7FgFcrD4pcp
p91vCBE6rL3rTCcP4lyrnstTQlwpa44so/jI2eBC27CNwgH2Yo2ViI/s5JZW
D3MKd/I0ElTL+TCsKo7MsaCwoO3oulDBh8FSwVmuNYWLnDxxSRUfVxcGj+yz
o5DPc6Uz6vn44Gp994YjhSGqNGz5BR9TvR1jdH8KgaFBGpUTfBwcoed4B1Ao
mbwykTbFx80EttVPgRTmj9ZWhk3z8VLj1lajgxQu/EvZSXv+n/8NPnrvcJhC
3umG6PQVPlTYrRMRkRT+SGV1htMEeF9qpF5zkoL+xp3ZzusFyFydZvEpikL4
9xGHOHICJKTs1UuKobDqWuunQUUB2nYIxVlnKWzOjTZwVRWg+/iI10AchVDN
giUDNQE0nbMrtc5TKC9ob5RlCjCW7q1TepGCXbGu5xN1AXx5wRaz8RSSDdy1
brAFUJhNNxVcotBRFjsdpSkAEYvVohMpqJkU17hrC1BlGCiuTaIQUNEVb6Qr
wH8Bu0TsVw==
           "]]}, "Charting`Private`Tag#1"]}}, {}}, <|
     "HighlightElements" -> <|
       "Label" -> {"XYLabel"}, "Ball" -> {"InterpolatedBall"}|>, 
      "LayoutOptions" -> <|
       "PanelPlotLayout" -> <||>, 
        "PlotRange" -> {{0., 43.07607738821883}, {0., 21.153068890259767`}}, 
        "Frame" -> {{False, False}, {False, False}}, "AxesOrigin" -> {0, 0}, 
        "ImageSize" -> {360, 432}, "Axes" -> {True, True}, "LabelStyle" -> {},
         "AspectRatio" -> Automatic, "DefaultStyle" -> {
          Directive[
           Opacity[1.], 
           AbsoluteThickness[2], 
           FaceForm[
            Opacity[0.3]], 
           RGBColor[0.24, 0.6, 0.8]]}, 
        "HighlightLabelingFunctions" -> <|
         "CoordinatesToolOptions" -> Identity, 
          "ScalingFunctions" -> {{Identity, Identity}, {
            Identity, Identity}}|>, "Primitives" -> {}, "GCFlag" -> False|>, 
      "Meta" -> <|
       "DefaultHighlight" -> {"Dynamic", None}, "Index" -> {}, "Function" -> 
        ParametricPlot, "GroupHighlight" -> False|>|>, "DynamicHighlight"]], 
   {Dashing[{0.01}], 
    LineBox[{{11.801804538466392`, 0}, {11.801804538466392`, 
      5.288267222564942}}]}},
  Axes->{True, True},
  AxesLabel->{
    FormBox[
     TagBox["\"t (Gyr)\"", HoldForm], TraditionalForm], 
    FormBox[
     TagBox["\"a(t)\"", HoldForm], TraditionalForm]},
  AxesOrigin->{0, 0},
  DefaultBaseStyle->{"PlotGraphics", "Graphics"},
  DisplayFunction->Identity,
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  ImagePadding->All,
  ImageSize->{Automatic, 173.5070803668482},
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None},
  PlotInteractivity:>True,
  PlotRange->{{0., 43.07607738821883}, {0., 21.153068890259767`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.05], 
     Scaled[0.05]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.543849851116379*^9, 3.54418616269827*^9, 3.544186258706572*^9, 
   3.544186883365192*^9, 3.60733736390539*^9, 3.924004121235881*^9, 
   3.9240042815606537`*^9, 3.924077941090787*^9, {3.924078008150661*^9, 
   3.9240780340903463`*^9}, 3.924078088562904*^9, 3.92407813501446*^9, 
   3.924078168210853*^9, 3.9240782327466393`*^9, 3.9240821548797398`*^9, 
   3.9240823616317587`*^9, {3.92418972842397*^9, 3.924189854755089*^9}, 
   3.986975961947535*^9, 3.986979673272891*^9, 3.986979703642321*^9},
 CellLabel->"Out[81]=",ExpressionUUID->"8bfc5e6f-4ad4-4d2b-a589-d21d46b559ad"]
}, Open  ]],

Cell["\<\
The vertical dotted line is the present time. The present age (in Gyr)  can \
be evaluated, and, if the evolution is bounded, the total age (in Gyr) and \
the maximum size (in billions of light years.)\
\>", "Text",
 CellChangeTimes->{{3.9240043163830748`*^9, 
  3.924004317054493*^9}},ExpressionUUID->"743a238a-c8b1-489d-a30a-\
af9b2368b862"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"agepresent", "=", 
  RowBox[{"Th", " ", 
   RowBox[{"x", "[", "1", "]"}], " "}]}]], "Input",
 CellLabel->"In[82]:=",ExpressionUUID->"cf88e7e6-2c03-4b1c-ba87-77c91a7aa7f5"],

Cell[BoxData["11.801804538466392`"], "Output",
 CellChangeTimes->{
  3.543849858679595*^9, 3.544186180498046*^9, 3.924004121250045*^9, 
   3.924004281583914*^9, 3.924077941175004*^9, {3.924078008196869*^9, 
   3.924078034145322*^9}, 3.9240780886127653`*^9, 3.924078135023859*^9, 
   3.924078168252315*^9, 3.9240782327876177`*^9, 3.9240821549180927`*^9, 
   3.924082361643025*^9, {3.924189728488822*^9, 3.9241898548137293`*^9}, 
   3.986975961960106*^9, 3.986979673294217*^9, 3.986979703651335*^9},
 CellLabel->"Out[82]=",ExpressionUUID->"6c93cd64-deaa-4ae6-93e3-b6aea127e42c"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", " ", 
   RowBox[{"agebigcrunch", "=", 
    RowBox[{"th", " ", 
     RowBox[{"xp", "[", "1", "]"}]}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.54418658286549*^9, 3.5441865837373753`*^9}, {
  3.924004105190589*^9, 3.924004105965377*^9}},
 CellLabel->"In[83]:=",ExpressionUUID->"67989f1a-c9cd-464a-b020-70e1a324e884"],

Cell[BoxData[
 RowBox[{"If", "[", 
  RowBox[{
   RowBox[{"\[CapitalOmega]k", "<", "\[CapitalOmega]crit"}], ",", " ", 
   RowBox[{"maxradius", "=", 
    RowBox[{"a0", " ", "ymax"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.544186584427169*^9, 3.544186585080593*^9}, {
  3.924004109271515*^9, 3.924004109575737*^9}},
 CellLabel->"In[84]:=",ExpressionUUID->"909068ab-461e-424b-ab7a-631aa06934e8"]
}, Open  ]]
}, Open  ]]
},
ScreenStyleEnvironment->"Working",
WindowToolbars->"EditBar",
WindowSize->{1067, 805},
WindowMargins->{{Automatic, 148}, {Automatic, 0}},
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"Magnification"->1,
"PaperOrientation"->"Portrait",
"PaperSize"->{612, 792},
"PostScriptOutputFile":>FrontEnd`FileName[{$RootDirectory, "home", "hartle", 
   "131book", "math", "FRW"}, "frw1.nb.ps", CharacterEncoding -> "ISO8859-1"]},
Magnification->1.25,
FrontEndVersion->"14.3 for Mac OS X ARM (64-bit) (July 8, 2025)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"1a388225-427b-4868-8602-cc7e31512757"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 99, 1, 87, "Title",ExpressionUUID->"6ab56dfe-2a2d-4adc-9284-94516474d8bf"],
Cell[681, 25, 3335, 84, 374, "Text",ExpressionUUID->"b8391488-a377-44e7-b60f-e71e9d98e4cb"],
Cell[CellGroupData[{
Cell[4041, 113, 105, 0, 66, "Subsection",ExpressionUUID->"e6aedb34-b49a-400c-af4e-68999a2f63d2"],
Cell[4149, 115, 189, 5, 42, "Text",ExpressionUUID->"e4806ddc-21e9-42f4-9dd9-6885efc3084b"],
Cell[4341, 122, 556, 10, 36, "Input",ExpressionUUID->"23463587-9071-455e-9b3f-05e542a62f83"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4934, 137, 103, 0, 66, "Subsection",ExpressionUUID->"deb930ed-d043-4f30-a67d-2d9a6b775af9"],
Cell[5040, 139, 322, 9, 42, "Text",ExpressionUUID->"ad8bab4c-b46f-41b2-a51a-fea95c48d9df"],
Cell[CellGroupData[{
Cell[5387, 152, 475, 10, 36, "Input",ExpressionUUID->"30cc047b-3a96-4273-a0eb-8c275e05154d"],
Cell[5865, 164, 653, 9, 41, "Output",ExpressionUUID->"89ac8e4d-921f-4b33-aba2-bf6339dd19ef"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6555, 178, 629, 9, 36, "Input",ExpressionUUID->"f252d9d9-6a33-4bc7-9c28-de6b3e4e97c8"],
Cell[7187, 189, 632, 9, 41, "Output",ExpressionUUID->"df38846c-2aab-452b-8345-a4504613c6bd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7856, 203, 602, 10, 36, "Input",ExpressionUUID->"d6a4e5a2-0971-402d-a52c-3d50751f6ea5"],
Cell[8461, 215, 614, 9, 41, "Output",ExpressionUUID->"4f82c729-2b0b-4960-8cb1-a47b4f54b8f6"]
}, Open  ]],
Cell[9090, 227, 1752, 57, 96, "Text",ExpressionUUID->"3ae14679-7710-4cb7-a0ba-12b2a2f3a7c3"],
Cell[CellGroupData[{
Cell[10867, 288, 292, 5, 36, "Input",ExpressionUUID->"8ee56d3f-3d2f-4baa-af9a-f9f58d70cc79"],
Cell[11162, 295, 613, 9, 41, "Output",ExpressionUUID->"dd55df7d-7e22-4b8c-9f4e-701a32803d6a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11812, 309, 156, 3, 36, "Input",ExpressionUUID->"c00bea6e-86a4-48c3-891f-9bfbaadb391c"],
Cell[11971, 314, 620, 9, 41, "Output",ExpressionUUID->"3a70c617-6fd6-4078-9e7e-7b54a16e1be8"]
}, Open  ]],
Cell[12606, 326, 452, 14, 42, "Text",ExpressionUUID->"24e96611-0854-4827-b7ea-4c4bf1f88fe4"],
Cell[13061, 342, 391, 7, 36, "Input",ExpressionUUID->"d13c5cb9-d8b4-45ab-aced-923155930d27"],
Cell[13455, 351, 96, 0, 42, "Text",ExpressionUUID->"eb1643e8-f0cf-42d8-a1b8-6fb73e47c416"],
Cell[CellGroupData[{
Cell[13576, 355, 212, 3, 36, "Input",ExpressionUUID->"2896bc8b-dfd3-4241-9fe7-f9285d66cba8"],
Cell[13791, 360, 628, 9, 41, "Output",ExpressionUUID->"0713fa54-e78a-4d80-8bca-1943608e6ee7"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[14468, 375, 133, 0, 66, "Subsection",ExpressionUUID->"9200126f-1431-4dca-b993-724b62b14593"],
Cell[14604, 377, 947, 32, 96, "Text",ExpressionUUID->"f84c4137-99f8-49c0-b8a6-9e737ee38167"],
Cell[15554, 411, 599, 16, 36, "Input",ExpressionUUID->"d1b3b8b5-568f-4ade-9bb3-a03a15e62fb4"],
Cell[16156, 429, 1692, 55, 96, "Text",ExpressionUUID->"63ef8e17-9bfd-48f6-b1e0-b89c5c223bba"],
Cell[17851, 486, 931, 30, 69, "Text",ExpressionUUID->"8b5360dc-577f-4c17-9a22-5ecaf4f3952d"],
Cell[CellGroupData[{
Cell[18807, 520, 247, 4, 36, "Input",ExpressionUUID->"995d8353-ca9f-4707-b6d5-4fdfc91f9b03"],
Cell[19057, 526, 612, 9, 41, "Output",ExpressionUUID->"e23a12d4-07b5-48c3-a217-5ca93a973c2d"]
}, Open  ]],
Cell[19684, 538, 462, 11, 36, "Input",ExpressionUUID->"1b57ef22-2ef1-4acb-ad18-7778d8d853a5"],
Cell[20149, 551, 448, 10, 36, "Input",ExpressionUUID->"dd64f4dc-f8de-411d-9282-04f016c07a05"],
Cell[CellGroupData[{
Cell[20622, 565, 469, 12, 36, "Input",ExpressionUUID->"8dc0139a-4316-4e32-a66e-fca498f5ca84"],
Cell[21094, 579, 35240, 658, 298, "Output",ExpressionUUID->"9ce7da3d-f19b-4696-bd9e-48d537d89aa9"]
}, Open  ]],
Cell[56349, 1240, 1450, 43, 127, "Text",ExpressionUUID->"1b00e417-94b2-44d8-aeaa-03c916d18945"],
Cell[57802, 1285, 602, 17, 36, "Input",ExpressionUUID->"7ec99045-096b-4996-b719-e2790e99b097"],
Cell[58407, 1304, 139, 3, 42, "Text",ExpressionUUID->"ba90b5ad-34e6-4a6f-89b5-e9771c1e559e"],
Cell[58549, 1309, 652, 17, 96, "Text",ExpressionUUID->"8476146b-7e78-4979-9dda-1d08c432c551"],
Cell[59204, 1328, 330, 9, 36, "Input",ExpressionUUID->"3740c1c5-696a-4d99-9c68-fb1efb9a576e"],
Cell[59537, 1339, 468, 10, 36, "Input",ExpressionUUID->"583820f6-dafd-465b-9dd4-97e8c8612159"],
Cell[CellGroupData[{
Cell[60030, 1353, 166, 2, 36, "Input",ExpressionUUID->"01570fca-4386-4bbb-9d3e-58ed4c4aa626"],
Cell[60199, 1357, 646, 10, 41, "Output",ExpressionUUID->"3c964217-7dca-444c-a7bf-3e6608ed9230"]
}, Open  ]],
Cell[60860, 1370, 461, 8, 96, "Text",ExpressionUUID->"b8a87e39-2001-4338-801b-07326d4ae205"],
Cell[CellGroupData[{
Cell[61346, 1382, 134, 2, 36, "Input",ExpressionUUID->"dd2d6361-44e0-4c71-a27f-011146a8cb15"],
Cell[61483, 1386, 635, 9, 41, "Output",ExpressionUUID->"a5f02194-6d1c-4e49-b5a9-af5bcae23f0f"]
}, Open  ]],
Cell[62133, 1398, 1129, 31, 124, "Text",ExpressionUUID->"d7ec73d4-64e2-4542-9e89-3dbd875a71c6"],
Cell[CellGroupData[{
Cell[63287, 1433, 397, 9, 36, "Input",ExpressionUUID->"cf55b747-aa05-4c94-a638-7e3f424bf4bd"],
Cell[63687, 1444, 1244, 25, 44, "Output",ExpressionUUID->"4b7ffe62-10ec-49a8-b415-8b5a94847359"]
}, Open  ]],
Cell[CellGroupData[{
Cell[64968, 1474, 168, 3, 36, "Input",ExpressionUUID->"6080aa2a-db61-4e88-abae-2f37d15f5161"],
Cell[65139, 1479, 1004, 17, 44, "Output",ExpressionUUID->"6acfab07-ed53-4c39-b503-306a5a2567b2"]
}, Open  ]],
Cell[CellGroupData[{
Cell[66180, 1501, 512, 12, 36, "Input",ExpressionUUID->"3d1ab90e-aebf-4da3-9510-2a96657115ab"],
Cell[66695, 1515, 672, 11, 41, "Output",ExpressionUUID->"0dc69b4c-9edd-49b0-8842-fc90ba9812e0"]
}, Open  ]],
Cell[CellGroupData[{
Cell[67404, 1531, 175, 3, 36, "Input",ExpressionUUID->"a4da45c1-75da-4ca9-be97-7fc53ab0d6a2"],
Cell[67582, 1536, 676, 11, 41, "Output",ExpressionUUID->"a41e1fc8-c4ce-4bb5-bddf-bb9532b50ef1"]
}, Open  ]],
Cell[68273, 1550, 549, 13, 36, "Input",ExpressionUUID->"0fdc212b-a17f-43b1-ba65-8f89c8ccb18e"],
Cell[CellGroupData[{
Cell[68847, 1567, 109, 1, 36, "Input",ExpressionUUID->"d8af3ffd-3003-4701-840e-2c0d0fe9a365"],
Cell[68959, 1570, 555, 8, 41, "Output",ExpressionUUID->"191be143-8e9b-4eb0-8cef-2f49962db839"]
}, Open  ]],
Cell[69529, 1581, 209, 5, 42, "Text",ExpressionUUID->"d512f051-dd5c-434e-b292-e5cf1120fa38"]
}, Open  ]],
Cell[CellGroupData[{
Cell[69775, 1591, 97, 0, 66, "Subsection",ExpressionUUID->"f8efa3a5-2ecc-4440-bc89-384f35b8cc44"],
Cell[69875, 1593, 456, 15, 69, "Text",ExpressionUUID->"d6a024c8-d0c6-4d4f-8eb3-c31d29607051"],
Cell[CellGroupData[{
Cell[70356, 1612, 109, 1, 36, "Input",ExpressionUUID->"0d0d2f30-4955-4528-b355-c666dd0ff455"],
Cell[70468, 1615, 635, 9, 41, "Output",ExpressionUUID->"1b5de4bd-ab8a-4269-9ccb-3c6572e0e793"]
}, Open  ]],
Cell[71118, 1627, 185, 4, 36, "Input",ExpressionUUID->"bb9b0085-e881-4d36-8f51-6934ee6caa7a"],
Cell[71306, 1633, 217, 5, 36, "Input",ExpressionUUID->"a96938f1-6062-457a-9ef0-6e390803eadc"],
Cell[71526, 1640, 457, 14, 36, "Input",ExpressionUUID->"3a19a1ce-7678-4c27-ba60-b84f42c618f4"],
Cell[71986, 1656, 455, 13, 36, "Input",ExpressionUUID->"91351a5e-c971-4447-8fb2-1e78d3a6b83b"],
Cell[72444, 1671, 464, 10, 36, "Input",ExpressionUUID->"26313eba-63e8-46f6-bdd6-38b1f862e472"],
Cell[72911, 1683, 461, 10, 36, "Input",ExpressionUUID->"9bff297a-87a1-4c90-90e0-ab1fb1f3c99e"],
Cell[73375, 1695, 547, 14, 62, "Input",ExpressionUUID->"2928ea34-ab48-4cb9-95d7-698db7d8b0b9"],
Cell[CellGroupData[{
Cell[73947, 1713, 682, 19, 62, "Input",ExpressionUUID->"c46e22bf-1fb2-45e7-9bb7-4c3e87aa53d5"],
Cell[74632, 1734, 12839, 257, 561, "Output",ExpressionUUID->"55df9539-1670-4e0d-9938-25dd799129fa"]
}, Open  ]],
Cell[87486, 1994, 120, 0, 42, "Text",ExpressionUUID->"89487cb2-9e4d-4755-863c-0667bf96fad3"]
}, Open  ]],
Cell[CellGroupData[{
Cell[87643, 1999, 214, 7, 66, "Subsection",ExpressionUUID->"b05aebc7-3a85-475a-a307-abb4cfeb34c2"],
Cell[87860, 2008, 964, 29, 69, "Text",ExpressionUUID->"c31f84d4-6ba7-48c0-8488-8715c7e8c77c"],
Cell[88827, 2039, 212, 5, 36, "Input",ExpressionUUID->"a8adecf2-baff-4f25-bc88-9030de284db1"],
Cell[89042, 2046, 391, 10, 36, "Input",ExpressionUUID->"fea91ced-534f-4fec-ae6f-9f1d7d43cc1d"],
Cell[CellGroupData[{
Cell[89458, 2060, 130, 2, 36, "Input",ExpressionUUID->"9166b52d-376f-4387-89de-ff60f52333c6"],
Cell[89591, 2064, 639, 9, 41, "Output",ExpressionUUID->"69082cee-2e10-4572-ba9f-a3e669d42be4"]
}, Open  ]],
Cell[90245, 2076, 212, 5, 36, "Input",ExpressionUUID->"3e47d984-d093-44c9-a31d-b3a7ef6e01e7"],
Cell[90460, 2083, 266, 6, 42, "Text",ExpressionUUID->"3055344a-aa2f-4051-ba38-b274fa9f8f8e"],
Cell[90729, 2091, 554, 14, 62, "Input",ExpressionUUID->"42fe58a7-d80c-4896-91fb-372e56a7624b"],
Cell[CellGroupData[{
Cell[91308, 2109, 753, 21, 62, "Input",ExpressionUUID->"bcadcf0a-2bb3-484c-8451-2579314c323c"],
Cell[92064, 2132, 12981, 260, 238, "Output",ExpressionUUID->"8bfc5e6f-4ad4-4d2b-a589-d21d46b559ad"]
}, Open  ]],
Cell[105060, 2395, 353, 7, 69, "Text",ExpressionUUID->"743a238a-c8b1-489d-a30a-af9b2368b862"],
Cell[CellGroupData[{
Cell[105438, 2406, 195, 4, 36, "Input",ExpressionUUID->"cf88e7e6-2c03-4b1c-ba87-77c91a7aa7f5"],
Cell[105636, 2412, 576, 8, 41, "Output",ExpressionUUID->"6c93cd64-deaa-4ae6-93e3-b6aea127e42c"]
}, Open  ]],
Cell[106227, 2423, 428, 9, 36, "Input",ExpressionUUID->"67989f1a-c9cd-464a-b020-70e1a324e884"],
Cell[106658, 2434, 395, 8, 36, "Input",ExpressionUUID->"909068ab-461e-424b-ab7a-631aa06934e8"]
}, Open  ]]
}, Open  ]]
}
]
*)