def moffat_r(r, alpha, beta):
    """
    Moffat profile

    :param r: radial coordinate
    :param alpha:
    :param beta:
    :return:
    """
    return 2. * (beta -1) / alpha ** 2 * (1 + (r/alpha) ** 2) ** (-beta)