In [2]:

    

import CoolProp
import CoolProp.CoolProp as CP
from math import sqrt
import numpy as np
import matplotlib.pyplot as plt
import scipy.optimize
%matplotlib inline


    

In [9]:

    

Tc = 590.7
rhoc = 351
b_0 =3.66766530
b_1 =-0.210687796
b_2 =-0.781330239
b_3 =0.130979005
b_4 =2.09502491
b_5 =6.28891793
k_1 =-3.94616949
k_2 =5.48487930

T = np.array([290,450])
tau = Tc/T
c_0 = b_0+1
c_1 = b_1/(-0.5)
c_2 = b_2/(-1.0/6.0)
c_3 = b_3/(-1.0/12.0)
c_4 = b_4
c_5 = b_5
R = 8.314472/60.05196
cp = R*(c_0+c_1/tau+c_2/tau**2+c_3/tau**3+c_5*(c_4*tau)**2*np.exp(-c_4*tau)/(1-np.exp(-c_4*tau))**2)
print cp


    

    




[ 1.0342473   1.45158917]

In [ ]: