Hello, I have a conceptual question to ask. Say I want to utilize the ode15s function to solve a system of coupled differential equations given a specific set of initial conditions, however, instead of just imputing the initial and final time span and allowing the program to workout the optimal time step, in the back end, I want to discritize the time span into specified time steps and format this into a time for loop.
The reason I ask is because I already tried implementing the ode15s function in the for loop and set the time span to be equivalent to the for loop time step, however this significantly increases the run time. Moreover, I am doing this because there is a separate process I want the program to run at every single time step which is why I want to dicritize in the first place.
I believe the issue arises from the fact that although I am discriticizing the time span in the for loop, the ode solver is calculating all the solutions in between each define time step. I already tried searching online in the forums and resources but the only parameter that is related to the time step that one may specify is the MaxStep size. There is no way to define a minimum step size. What should one do in this scenario? PS - May it have something to do with defining a Jacobian for my system of odes? Or is that completely unrelated? If you are interested in looking at the for loop with the nested ode solver I link it here: Code 
