A <- rbind(c(0.7, 0.2), c(0.3, 0.8))
eigen(A)

####################################################

A <- rbind(c(0.7, 0.2), c(1, 0))
b <- c(1, 1)
solve(A, b)

####################################################

A <- rbind(c(0.7, 0.2), c(-0.7, -0.2), c(1, 0))
b <- c(1, -1, 1)
qr.solve(A, b)

####################################################

p <- c(0.7, 0.2, 0.1)
sample(3, 1, prob=p)