library(deSolve)

model <- function(time, Tr, k) {
  dTr <- -(k[1]*k[2]*k[3]*(Tr - k[5]))/(k[4]*k[2]*k[3])
  list(dTr)
}

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

k = 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 = k)
print(solution)