Actually I figured it out. The issue was I got in the process so confused I did not realise that time t_star is provided by the function as solution.x. I instead decided to make my own vector like this
t_star = linspace(0,5,length(sol.x)).
Don't ask me how this came to my mind more naturaly than just taking sol.x.... But that could work in case the sol.x would have consistent sampling.
t_star = linspace(0,5,length(sol.x)).
Don't ask me how this came to my mind more naturaly than just taking sol.x.... But that could work in case the sol.x would have consistent sampling.