code: def vap(x, a, b, c):
      
      return np.exp(a+b/x+c*np.log(x))
docs: Vapor pressure model
      
      Parameters
      ----------
      x: int
      a: float
      b: float
      c: float
  
      Returns
      -------
      float
          np.exp(a+b/x+c*np.log(x))