“…The dashed curve with diamonds is, instead, the result of a natural continuation with secant prediction, recall Figure 1 (bottom-right), always implementd in Matlab. Following [32], the correction step is made by using the Broyden's update, integrals and external IVPs are solved simultaneously via the embedded Runge-Kutta pair DOPRI54 (see, e.g., [22]: basically, it is also behind ode45) and the maturation age is directly obtained as an output of the latter, which is indeed capable of event detection, i.e., it automatically detects when (10) is satisfied during the integration of (8). As for the previous approach, integrals, IVPs and maturation age are computed from scratch at every continuation step, independently of the same quantities computed at the previous step.…”