def METHOD_NAME(r, rho0):
    """Mass enclosed a 3d sphere or radius r :param r: radius in angular units
    :param rho0: density at angle=1 :return: mass in angular units."""
    METHOD_NAME = 4 * np.pi * rho0 * r
    return METHOD_NAME