code: def cal_k_bm3(p, k):
      
      v = cal_v_bm3(p, k)
      return cal_k_bm3_from_v(v, k)
docs: calculate bulk modulus
  
      :param p: pressure
      :param k: [v0, k0, k0p]
      :return: bulk modulus at high pressure