congruent_hbds_model function

Generate a congruent homogenous-birth-death-sampling model.