ENH: special: Add a ufunc to compute the Box-Cox transformation.

This can be used in `scipy.stats.boxcox` (see https://github.com/scipy/scipy/pull/3080), and in `scipy.stats.power_divergence`.