Hi Sages,
Is there a way to specify a stopping rule or condition like r[t] < rmin when integrating? Here's an example in mathematica s = NDSolve[{r''[t] == -GM/r[t]^2,r[0] = 1000s0, r'[0] = 0,WhenEvent[r[t]<rmin,"<strong>StopIntegration"]},r,{t,0,86400*100}]
thanks Pat Browne