def METHOD_NAME(size=1):
    """Generate a random sampling of points on the surface of the unit sphere.
    Parameters
    ----------
    size : int
        The number of points to generate.
    Returns
    -------
    rep : `~astropy.coordinates.UnitSphericalRepresentation`
        The random points.
    """
    rng = np.random  # can maybe switch to this being an input later - see #11628
    lon = rng.uniform(0, _TWOPI, size) * u.rad
    lat = np.arcsin(rng.uniform(-1, 1, size=size)) * u.rad
    return UnitSphericalRepresentation(lon, lat)