(***************************  filename: nbdata.m  ****************************)

<<Miscellaneous`PhysicalConstants`
<<Miscellaneous`Units`

$MaxPrecision = Infinity
$MaxExtraPrecision = Infinity

USimplify[x_] := Simplify[x, Assumptions -> {Meter > 0, Centimeter > 0,
                                             Second > 0, Year > 0}]

$TextStyle = {"Courier", FontSize -> 18}

ic = StyleForm["c", FontSlant-> "Italic"]
ik = StyleForm["k", FontSlant-> "Italic"]
it = StyleForm["t", FontSlant-> "Italic"]
iA = StyleForm["A", FontSlant-> "Italic"]
iB = StyleForm["B", FontSlant-> "Italic"]
iH = StyleForm["H", FontSlant-> "Italic"]
iQ = StyleForm["Q", FontSlant-> "Italic"]
iR = StyleForm["R", FontSlant-> "Italic"]
iS = StyleForm["S", FontSlant-> "Italic"]
iX = StyleForm["X", FontSlant-> "Italic"]
iY = StyleForm["Y", FontSlant-> "Italic"]
iZ = StyleForm["Z", FontSlant-> "Italic"]

Off[General::spell1]


(***  In order to utilize Mathematica's high precision arithmetic,  ***)
(***  enter all data as rational numbers (no decimal points).       ***)

(* The basic inputs are H0, Hmax, and Q0. H0 and Hmax are entered here.      *)
(* Q0 is prompted for as the evaluation of CosmicEvolution.nb proceeds.      *)

    H0 = 72 (Kilo Meter/Second)/(Mega Parsec) (* H0 := H(t0); t0 is the      *)
                                              (* present epoch in seconds.   *)
  Hmax = 5 10^60 H0                           (*                             *)

Switch[Atype, Apos, Q0 = Input["Enter a target between 0 and 1 for the
acceleration parameter Q0 at the present epoch, e.g., \"1/100\" or \"1/2\" or
\"99/100\" (enter as a rational number -- no decimal point; without
quotes):  "],
              Aneg, Q0 = Input["Enter a target less than 1 for the
acceleration parameter Q0 at the present epoch, e.g., \"0\" or \"-1\" (enter
as a rational number -- no decimal point; without quotes):  "] ]


      k = -1/Centimeter^2

      c = Convert[SpeedOfLight, Centimeter/Second]

  H0sec = Convert[H0, 1/Second]

Hmaxsec = Convert[Hmax, 1/Second]

   H0cm = H0sec/c

 Hmaxcm = Hmaxsec/c


Print[StringForm[" "]]

Print[StringForm["`1` = ``,     `3` = ``,     \!\(c\^2\) `5` = ``",
                  ik, N[k],     ic, N[c],     ik, N[c^2 k]         ]]

Print[StringForm[" "]]

Print[StringForm["\!\(H\_0\) = `` = ``,     \!\(H\_0\)/`3` = ``",
                  H0, N[H0sec]        ,     ic, N[H0cm]          ]]

Print[StringForm[" "]]

Print[StringForm["\!\(H\_max\) = `` = ``,     \!\(H\_max\)/`3` = ``",
                  N[Hmax], N[Hmaxsec]   ,     ic, N[Hmaxcm]          ]]

Print[StringForm[" "]]

Switch[Atype, Azero,
       Print[StringForm["\!\(H\_max\)/\!\(H\_0\) = ``", N[Hmax/H0]]], _, 
       Print[StringForm["\!\(Q\_0\) = ``,     \!\(H\_max\)/\!\(H\_0\) = ``",
                         N[Q0]          ,     N[Hmax/H0]                    ]]]

Print[StringForm[" "]]

totaltimeused = 0

newtimeused = TimeUsed[] - totaltimeused

totaltimeused = totaltimeused + newtimeused

Print[StringForm["Memory in use = ``",
                  N[MemoryInUse[]/2^20,7] megabytes]]

(**************************  End of file nbdata.m  ***************************)
