facilityR0 function

Calculate basic reproduction number R0