# NOT RUN {
## See the function gemDualLinearProgramming.
A <- matrix(c(
0, 0, 0, 1,
8, 6, 1, 0,
4, 2, 1.5, 0,
2, 1.5, 0.5, 0
), 4, 4, TRUE)
B <- matrix(c(
60, 30, 20, 0,
0, 0, 0, 0,
0, 0, 0, 0,
0, 0, 0, 0
), 4, 4, TRUE)
S0Exg <- {
S0Exg <- matrix(NA, 4, 4)
S0Exg[2:4, 4] <- c(48, 20, 8)
S0Exg
}
ge <- sdm2(
A = A, B = B, S0Exg = S0Exg,
maxIteration = 1,
numberOfPeriods = 1000,
ts = TRUE
)
matplot(ge$ts.q, type = "l")
ge2 <- sdm2(
A = A, B = B, S0Exg = S0Exg,
maxIteration = 1,
numberOfPeriods = 1000,
ts = TRUE,
policy = makePolicyMeanValue(150)
)
matplot(ge2$ts.q, type = "l")
# }
Run the code above in your browser using DataLab