sample_c_two function

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