def METHOD_NAME():
    reference_cap = np.array(
        [0.643501108793284, 1.570796326794896,
         2.498091544796509, 3.141592653589793])
    reference_regions = np.array([1, 4, 4, 1])
    dim = 2
    N = 10
    s_cap, n_regions = eq.caps(dim, N)
    npt.assert_almost_equal(reference_cap, s_cap)
    npt.assert_almost_equal(reference_regions, n_regions)