calc_burden function

Calculate case burden rate