Thanks for your work with this package!
I'm hoping to generate some data to test out a collocation approach using different choices of points, which means integrating the ODE and evaluating at irregularly spaced points. It's common in integration methods to allow the user to choose the points to evaluate the solution. E.g. solve_ivp's parameter t_eval.
Desired
It would be great if DynSys.make_trajectory() could handle t_eval.
Alternative
User finds highest dt that divides all their evaluation times, oversamples using make_trajectory(), then determines the indices to thin out to the desired evaluation times.