def cal_k_bm3(p, k):
    """
    calculate bulk modulus

    :param p: pressure
    :param k: [v0, k0, k0p]
    :return: bulk modulus at high pressure
    """
    v = cal_v_bm3(p, k)
    return cal_k_bm3_from_v(v, k)