n <- c(56, 63, 121, 123) r <- c(6, 9, 18, 7) H <- 4