常微分方程（ODEs）求解：如何避免出现负值？

我现在正在做一个空间动态的蚊虫动力学模型：给每个空间网格单元都单独部署一套常微分方程（ODEs）系统，简单来说就是每个景观单元对应一个独立的ODE模型。模型里的三个核心变量——易感蚊虫数量（Sv）、暴露蚊虫数量（Se）和感染蚊虫数量（St）——会在每个时间步完成数值更新，而且这套ODE模型还和离散的动物运动主体模型做了耦合，实现空间层面的动态交互。

下面是我目前在用的ODE模型运行的示例代码片段：

library(deSolve) 
mod1 <- function(out_tab, time_step, var){ 
  Sv <- out_tab[time_step...

内容的提问来源于stack exchange，提问作者Nell