H <- 10
N <- 5
D <- 7
C <- 9

r <- N / D
d <- H * (1 + r - 2 * r ^ C) / (1 - r)
cat(d, "\n", sep = "")