library(deSolve)

model <- function(time, Tr, parms) {
#  pa <- parms$pa
#  Ca <- parms$Ca
#  V <- parms$V
#  Text <- parms$Text

   dTr <- parms[1]*parms[2]*Tr
#  dTr <- -(qv * pa * Ca * (Tr - Text)) / (V * pa * Ca)
  list(dTr,parms[2])
}

data <- read.csv("input.csv")
head(data)

#parms = list(qv = 54, pa = 1.204, Ca = 1006, V = 3000,Text=data$temp)


parms <- c(growth_rate=0.03, 
       input_values=data$temp)

solution <- ode(y = 100, times = data$time, func = model, parms = parms)
print(solution)