def ball_volume(ndim, radius): """Return the volume of a ball of dimension `ndim` and radius `radius`.""" n = ndim r = radius return np.pi ** (n / 2) / special.gamma(n / 2 + 1) * r ** n