library(ggplot2)
des1CMTbolus <- function(t,y,p){
  with(as.list(c(p,y)), {
    K      <- CL/V
    DADT1  <- - K*A1
    CP     <- A1/V
    list(c(DADT1),c(CP=CP))
  })
}
pars   <- c(CL=1,V=1)
init   <- c(A1=0)
events <- dose_func(cmt=1,value=10,times=0)
times  <- seq(0,24,length.out=200)
out    <- as.data.frame(deSolve::lsoda(init,times,des1CMTbolus, pars,events=list(data=events)))
out    <- reshape2::melt(out,id='time')
ggplot(out,aes(time,value)) + geom_line() + facet_wrap(~variable,scales='free')
