code: def _prop0(self, rho, T):
          
          rhoc = self._constants.get("rhoref", self.rhoc)
          Tc = self._constants.get("Tref", self.Tc)
          delta = rho/rhoc
          tau = Tc/T
          ideal = self._phi0(tau, delta)
          fio = ideal["fio"]
          fiot = ideal["fiot"]
          fiott = ideal["fiott"]
  
          propiedades = _fase()
          propiedades.h = self.R*T*(1+tau*fiot)
          propiedades.s = self.R*(tau*fiot-fio)
          propiedades.cv = -self.R*tau**2*fiott
          propiedades.cp = self.R*(-tau**2*fiott+1)
          propiedades.alfap = 1/T
          propiedades.betap = rho
          return propiedades
docs: Ideal gas properties