logllk_div_time_two function

Compute loglikelihood of divergence times for a(t) = c/(1-t)^2