anoxicmin = Min*(1-O2/(O2+ks))* SO4/(SO4+ks2)
dOM =Flux - oxicmin - anoxicmin
dO2 =-oxicmin -2*rox...*HS*(O2/(O2+ks)) + D*(BO2-O2)
dSO4 = -0.5*anoxicmin +rox*HS*(O2/(O2+ks)) +D*(BSO4-SO4)
dHS =0.5*...anoxicmin -rox*HS*(O2/(O2+ks)) + D*(BHS-HS)
list(c(dOM, dO2, dSO4, dHS), SumS = SO4+HS)
})
}
#参数的设置...pars<- c(D = 1, Flux = 100, r = 0.1, rox = 1,
ks = 1, ks2 = 1, BO2 = 100, BSO4 = 10000, BHS