非线性微分方程数值解系统

Question

我正在尝试求解非线性微分方程组。我的系统有15个必须估计的参数，我正在对其进行优化。有些参数会使我的系统僵化，而有些则不会。因为ODE求解器会因错误而终止，所以这可能会使我的优化结果成为垃圾。解决此类问题的最佳方法是什么？

如果我用隐式方法解决它，也许我的系统不会僵硬，可以吗？如果是，那么您建议的最佳数值方法是什么。

非常感谢。