@staticmethod
def beta_r(r, r_ani):
    """
        anisotropy as a function of radius

        :param r: 3d radius
        :param r_ani: anisotropy radius
        :return: beta
        """
    return r ** 2 / (r_ani ** 2 + r ** 2)
