sample_c_one function

Sample divergence function parameter c for a(t) = c / (1-t) through Gibbs sampler