pop_A <- c(15000, 20000, 10000)
deaths_A <- c(360, 400, 280)
pop_B <- c(40000, 52000, 8000)
deaths_B <- c(1000, 1040, 240)
CDR_A <- sum(deaths_A) / sum(pop_A) * 1000
CDR_B <- sum(deaths_B) / sum(pop_B) * 1000
cat("Crude Death Rate for Town A:", CDR_A, "\n")
cat("Crude Death Rate for Town B:", CDR_B, "\n")