#构建模型
model<-function(t,y, pars) {
with(as.list(c(y, pars)),{
Min = r*OM
oxicmin = Min*...anoxicmin -rox*HS*(O2/(O2+ks)) + D*(BHS-HS)
list(c(dOM, dO2, dSO4, dHS), SumS = SO4+HS)
})
}
#参数的设置
pars...= 1, O2 = 1, SO4 = 1, HS = 1)
#计算
print(system.time(
ST2 <- runsteady(y = y, func = model,parms = pars...区别就是为了适应生物学的方程stode和stodes函数有一个关键的参数pos可以保证参数的正性:
print(system.time(
ST <- stode(y = y, func = model, parms =pars