Conditionally execute a function inside a vector field while integrating with an ODE solver

I am integrating some equations that represent an astrodynamics problem. They are defined by a vector field such as:

where x are the states (e.g. position, velocity, mass...). The dynamics that I'm integrating are fed by a control u , that is computed based on the states and time: . This control u is computed inside the vector field, which we could explicitely state as . I am integrating with an ODE solver (e.g. ode45()), and thus the solver computes the control at every time step.

Now, 我万t is to avoid computing the control at every time step, but rather only at certain times, with a given frequency (e.g. let's say I want to compute the control at 0.1 Hz, every 10s). How would you recommend me to approach this?

I hope I was clear enough. Thank you for any help.