MAINT: an easier way to make ranges from slices
MAINT: split distributions.py into three files
MAINT: remove globals() as suggested by @evr-br
MAINT: make ``python setup.py bdist_wheel`` work.
TST: [stats, discrete] tighten cdf-ppf tests
TST: [stats, discrete] tighten cdf-pmf tests
TST: [stats, discrete] remove non-running tests
TST: [stats, discrete] remove sample tests

They are too noisy: only work for decimal=0.
BUG: integrate: odeint handled the banded jacobian incorrectly.

When `col_deriv` is False, the C wrapper that calls the user's python code that
computes the jacobian matrix uses a macro, MATRIXC2F(jac, data, n, m), to copy
the data into the array that is used by the Fortran solver.  In that macro.
`jac` is the destination, `data` is the source, and `n` and `m` are the number
of rows and columns of the matrix.  The wrapper had the last two arguments
reversed.  In the case of a full jacobian, this did not matter, because the
matrix is square.  In the banded case, the number of rows and columns are not
the same.  This means that the only way to correctly use a banded jacobian was
to put the bands in the columns of the matrix, and use `col_deriv=True` so
MATRIXC2F was not used.  That would result in the data being in the correct
Fortran order, since the Fortran code expects the bands to be in the rows.

This commit fixes the use of MATRIXC2F.  Also, the docstring of odeint is
updated to reflect the correct API:  the bands should be put in the rows
of the matrix, not the columns (unless `col_deriv=True` is used).
BUG: special: comb returned incorrect result when given large np.int64 arguments.

When given `exact=True`, cast the arguments of `comb` to python ints to ensure
that the result is a python integer.  Otherwise a call such as
    comb(np.int64(70), np.int64(30), exact=True)
returns an incorrect result, because it does the calculations with np.int64s.
The true result is too large for a 64 bit integer.
Merge pull request #3146 from WarrenWeckesser/bug-misc-comb

BUG: special: comb returned incorrect result when given large np.int64 arguments
ENH: special: Add a ufunc to compute the Box-Cox transformation.
Merge pull request #3150 from WarrenWeckesser/enh-special-boxcox

ENH: special: Add a ufunc to compute the Box-Cox transformation.
move factorial* to scipy.special
update docs
Merge pull request #3145 from WarrenWeckesser/bug-odeint-banded-jac

BUG: integrate: odeint handled the banded jacobian incorrectly.
Merge pull request #3135 from rgommers/bdist-wheel

MAINT: make ``python setup.py bdist_wheel`` work.
Merge pull request #3123 from perimosocordiae/patch-2

MAINT: sparse: an easier way to make ranges from slices
Merge pull request #3151 from richardtsai/move_factorial_funcs_to_special

MAINT: Move scipy.misc.factorial* to scipy.special
DOC: complete stats.boxcox docstring.  Add add boxcox functions to refguide.

Also do a PEP8 cleanup of `boxcox`.