ngay_hl <- seq(as.Date("2014-05-13"), by="day", length.out = 30)
ngay_kt <- ngay_hl + 300
dt_bh <- sample(5000000:7000000,30)
d <- data.frame(ngay_hl, ngay_kt, dt_bh)

sum.dt_bh.period <- function(x,y) {
  a <- as.Date("x","%d/%m/%Y")
  b <- as.Date("y","%d/%m/%Y")
  D_temp  <- mutate(d,
  var_M =
  ifelse(ngay_hl %in% a:b & ngay_kt %in% a:b, (ngay_kt-ngay_hl)*dt_bh/as.numeric(ngay_kt-ngay_hl),
  ifelse(ngay_hl %in% a:b & ngay_kt >b, (b-ngay_hl)*dt_bh/as.numeric(ngay_kt-ngay_hl),
  ifelse(ngay_hl < a & ngay_kt %in% a:b, (ngay_kt-a)*dt_bh/as.numeric(ngay_kt-ngay_hl),
  ifelse(ngay_hl < a & ngay_kt > b,(b-a)*dt_bh/as.numeric(ngay_kt-ngay_hl),0)))))
colSums(D_temp[,"var_M"],na.rm = TRUE)
}