I am trying to implement some numerical algorithm for an optimal control problem (still).
It involves double cycle. In the inner one two DE systems are solved sequentially and the result is used to update values of some objective function.
This part uses remeber option thanks to pagan and is ok.
However, in the outer cycle I use the values of this objective function to obtain solution for yet another DE and update ICs for the inner cycle.