Hello Everyone! Why is this giving me negative values? Mp = 2.4353*10^18; m = 1.8*10^13; Rm = (1.8*10^13)/(2.4353*10^18); tfin = 10^7; ysol = ParametricNDSolve[{y''[ t] + (G + 6) Sqrt[(y'[t]^2 + Rm^2 y[t]^2)/6] y'[t] + Rm^2 y[t] == 0, y[0] == -16.8, y'[0] == 0.00033}, y, {t, 0, tfin}, {G}]; HSol[G_] = (Sqrt[((D[y[G][t], t] /. ysol)^2 + Rm^2 (y[G][t] /. ysol)^2)/6]); Table[(16 D[HSol[G], t]/HSol[G]^2) /. t -> tfin, {G, 0, 1, 0.1}]