R- deSolve中微分方程的求解

Question

我有一个等式，如下所示；

dN/dt = N(t)G(t)

G(t)由以下等式给出：dG/dt = a * G我如何在R中使用deSolve包中的ode函数来解决这个问题？

Answer 1

正如dario已经提到的，这个问题缺乏一些细节。尽管如此，让我们尝试一下答案。如果我们假设a < 0，该模型看起来像Gompertz growth的ode公式：

dN/dt = N * G
dG/dt = a * G

然后，这可以通过以下方式解决：

library(deSolve)

model <- function(t, y, p) {
  with(as.list(c(y, p)), {
    dN <-  N * G
    dG <-  a * G
    list(c(dN, dG))
  })
}

y      <- c(N = 1, G = 1)
parms  <- c(a = -0.1)
times  <- seq(0, 100)
out <- ode(y, times, model, parms)
plot(out)