# NOT RUN {
times<- seq(2005, 2100, 1/12)
temp_cmin <- 18
temp_cmax <- 26
# Temperature at which performance is at its maximum value.
temp_op <- (temp_cmax+temp_cmin)/3+sqrt(((temp_cmax+temp_cmin)/3)^2-
(temp_cmax*temp_cmin)/3)
ro <- 0.8
# Temperature that occurs in the minimum time of the simulation.
temp_i <- 20
temp <- get_RCP8.5(times)+temp_i
rate <- rate_TPC(temp,ro,temp_cmin,temp_cmax,temp_op)
plot(times,rate, type="l")
# }
Run the code above in your browser using DataLab