f_gam function

Calculate serial interval mixture density assuming underlying gamma distribution