# NOT RUN {
## IASA model
# Parameters and initial conditions.
pars_solve_iasa = c(
b1 = 21871, b2 = 4374,
df1 = 0.104, dm1 = 0.098, df2 = 0.125, dm2 = 0.118,
sf1 = 0.069, sf2 = 0.05, sm1 = 0.028, sm2 = 0.05,
k1 = 98050, k2 = 8055, h1 = 1, h2 = 0.5,
a = 0.054, alpha = 0.1, v = 0.2, z = 0.1)
# Set ranges 10 % greater and lesser than the
# point estimates.
rg_solve_iasa <- SetRanges(pars_solve_iasa)
# }
Run the code above in your browser using DataLab