# NOT RUN {
A <- matrix(c(
0.05, 0.05, 0.1,
0.1, 0, 0.1,
0, 0.15, 0.05
), 3, 3, TRUE)
S <- diag(3)
# a market-clearing price vector
p <- c(0.6, 0.9, 1)
result <- F_Z(A, p, S)
# Each sales rate is equal to 1
result$q
# the purchase matrix
A %*% diag(result$z)
# a non-market-clearing price vector
p <- c(1, 1, 1)
result <- F_Z(A, p, S)
# Some sales rates don't equal 1
result$q
# the purchase matrix
A %*% diag(result$z)
# }
Run the code above in your browser using DataLab