def heatCapacity_Shomate(T, hCP):
    """
    heatCapacity_Shomate(T, hCP)
      NIST vapor, liquid, and solid phases

    heat capacity correlation
      hC (J/mol/K) = A + B*t + C*t^2 + D*t^3 + E/t^2
      t (K) = T/1000.0
    
    Parameters
      T, temperature in Kelvin
      hCP, A=hCP[0], B=hCP[1], C=hCP[2], D=hCP[3], E=hCP[4]
      A, B, C, D, and E are regression coefficients  

    Returns
      heat capacity in J/mol/K at T
    """
    t = T/1000.0
    return hCP[0] + hCP[1]*t + hCP[2]*t**2 + hCP[3]*t**3 + hCP[4]/t**2