lam <- c(1, 1, 1, 1, 1, 1)

