How to store value of variable at specific time steps while MATLAB simulation runs?

Question

I am simulating a dynamic system using the ode15s(..) MATLAB solver:

tspan = 0:0.01:8;
[t,theta] = ode15s(@ode_simulation, tspan, theta)
<