library(mosaic)
drug_remaining <- function(dose,duration,time_constant) {
  dose* exp(- duration / time_constant)
}
t = 0:20
  plotFun(drug_remaining(dose = 100, time_constant = 4, duration = t) ~ t, t.lim=range (0,20))