peak <- tgrid(0,6,0.2)
sparse <- tgrid(0,24,4)
day1 <- c(peak,sparse)
design <- c(day1, day1+72, day1+240)
if (FALSE) {
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