我有一个由微分方程系统给出的模型

与 和。该系统的工作原理如下：将 ode 与初始值直到达到某个值（例如），获取发生这种情况的时间并重新开始初始值 。所有这一切，直到遇到一些停止条件。$y=(y_1,y_2,y_3)$$f:\mathbb{R}^3 \to \mathbb{R}^3$$(y_1^0,y_2^0,y_3^0)$$y_1$$y_1=0$$t_f$$(y_1^0 ,y_2^0, y_3(t_f))$

我的问题是：用什么软件可以模拟这样的系统？最好是开源的，但欢迎任何建议。