(********************  filename: ProcEqS.m (S := R/Rmin)  ********************)

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

T1[t_] := USimplify[c^2 A/3 Second^2] S[t]

T2[t_] := (N[USimplify[2 c^2 B Second^2],1000]/N[(3 Rmin^6),1000])/S[t]^5

state = First[NDSolve`ProcessEquations[
              {Sdot'[t] == T1[t] + T2[t],
                  S'[t] == Sdot[t],
                   S[0] == 1, Sdot[0] == 0},
              {S, Sdot}, t,
              AccuracyGoal -> 50, PrecisionGoal -> 50,
              WorkingPrecision -> 100, MaxSteps -> Infinity]]


Print[StringForm[" "]]

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

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