def rho_YY(rho_SL_k, rho_SL_k_comp, rho_SL_l, rho_SL_l_comp, rho_SS_kl,
    rho_LL_kl, rho_YL_kl, rho_YL_lk):
    """
    rho_SL_k, rho_SL_k_comp:  G2F parameters for asset k
    rho_SL_l, rho_SL_l_comp:  G2F parameters for asset l
    rho_SS_kl, rho_LL_kl, rho_YL_kl, rho_YL_lk:  cross terms k-l
    """
    den = rho_SL_k_comp * rho_SL_l_comp
    num = (rho_SS_kl - rho_SL_k * rho_SL_l * rho_LL_kl - rho_SL_k *
        rho_SL_l_comp * rho_YL_lk - rho_SL_l * rho_SL_k_comp * rho_YL_kl)
    return num / den
