Как включить изменение параметров во времени в дифференциальных уравнениях

Question

Я создаю динамическую модель биопленки и хочу включить итерации определенных параметров, которые меняются со временем.

Я создал свою модель, чтобы решить ее с помощью пакета deSolve, теперь я не знаю, гдеЯ должен ввести изменение параметров для решения проблемы.

Пример роста водорослей

Biofilm <- function(t, state, parameters) {
  with(as.list(c(state, parameters)), {
    dX <-  (umax * B) * (1/(1+ kinv*B))* (exp(Beta*(T-T0)))*(I/I+kI) -(Cdet * u*(B - B0)) 
           - (kflood* u*(B - B0))-(Cdet * u*(B - B0)) - (kflood* u*(B - B0)) - LGr
     list(c(dX))
  })
}

parameters <- c(B0=0, I = 1, u = 0.37, umax = 1, kinv = 0.9, kflood = 1.6, Cdet = 0.5, Beta = -0.08, T0 = 20, kI = 1, T = 7, LGr=1)
state      <- c(B = 1)
times      <- seq(0, 55, by = 0.041)

out <- ode(y=state, times, Biofilm, parameters)

plot(out)

Я получил логистическую кривую, как и ожидалось, но в моем уравнении мне нужно предоставить I и T, которые меняются со временем, чтобы получить синусоидальный ответ.