t <- c(0:35); d <- 0.85; k <- 1.90; E <- (0.142*t+1.095)*(t+17.8)*k*d; plot(t, E);