lam <- c(0.7, 0.7, 0.7, 0.7, 0.7, 0.7)
a <- c(0, 0, 0, 0, 0, 0, 0, 0, 0)

