# NOT RUN {
peak <- tgrid(0,6,0.2)
sparse <- tgrid(0,24,4)
day1 <- c(peak,sparse)
design <- c(day1, day1+72, day1+240)
# }
# NOT RUN {
mod <- mrgsolve::house()
out <- mod %>% ev(amt=1000, ii=24, addl=10) %>% mrgsim(tgrid=design)
plot(out,CP~., type='b')
# }
Run the code above in your browser using DataLab