Added estimation for critical alpha

The critical alpha, when the shape function changes from linear to
non-linear, can be estimated by the time it takes to get to the best
discrepancy divided by the location of that discrepancy in the small
alpha (=0.5) regime.