code: def _Region4(P, x):
      
      T = _TSat_P(P)
      if T > 623.15:
          rhol = 1./_Backward3_sat_v_P(P, T, 0)
          P1 = _Region3(rhol, T)
          rhov = 1./_Backward3_sat_v_P(P, T, 1)
          P2 = _Region3(rhov, T)
      else:
          P1 = _Region1(T, P)
          P2 = _Region2(T, P)
  
      propiedades = {}
      propiedades["T"] = T
      propiedades["P"] = P
      propiedades["v"] = P1["v"]+x*(P2["v"]-P1["v"])
      propiedades["h"] = P1["h"]+x*(P2["h"]-P1["h"])
      propiedades["s"] = P1["s"]+x*(P2["s"]-P1["s"])
      propiedades["cp"] = None
      propiedades["cv"] = None
      propiedades["w"] = None
      propiedades["alfav"] = None
      propiedades["kt"] = None
      propiedades["region"] = 4
      propiedades["x"] = x
      return propiedades
docs: Basic equation for region 4
  
      Parameters
      ----------
      P : float
          Pressure, [MPa]
      x : float
          Vapor quality, [-]
  
      Returns
      -------
      prop : dict
          Dict with calculated properties. The available properties are:
  
              * T: Saturated temperature, [K]
              * P: Saturated pressure, [MPa]
              * x: Vapor quality, [-]
              * v: Specific volume, [m³/kg]
              * h: Specific enthalpy, [kJ/kg]
              * s: Specific entropy, [kJ/kgK]