docs: def bm3_k(p, v0, k0, k0p):
      """
      calculate bulk modulus, wrapper for cal_k_bm3
      cannot handle uncertainties
  
      :param p: pressure
      :param v0: volume at reference conditions
      :param k0: bulk modulus at reference conditions
      :param k0p: pressure derivative of bulk modulus at different conditions
      :return: bulk modulus at high pressure
      """
      return cal_k_bm3(p, [v0, k0, k0p])