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 ~ time, type = 'b')
}
Run the code above in your browser using DataLab