library(deSolve)

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

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

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

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


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

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