The logic sounds simple. LHS (theta) can be any physical property like critical temperature, critical pressure, boiling point temperature and critical volume; for each of these property the contants a to f are given in table. You need to just substitute molecular weight, M.
M = 44.01 ; % Co2 molecular weight
gama = 1.53 ; % specific gravity of co2
K = [544.4 0.2998 1.0555 -1.3478*10-3 -0.61641 0.0
4.5203*104 -0.8063 1.6015 -1.8078*10-3 -0.3084 0.0
1.206*10-2 0.20378 -1.3036 -2.657*10-3 0.5287 2.6012*10-3
6.77857 0.401673 -1.58262 3.77409*10-3 2.984036 -4.25288*10-3 ] ; % Table 1-4
a = K(:,1) ;
b = K(:,2) ;
c = K(:,3) ;
d = K(:,4) ;
e = K(:,5) ;
f = K(:,6) ;
theta = a.*M.^b.*gama.^c.*exp(d*M+e*gama+f.*M*gama) ;
theta would be a 4x1 array giving you critical temperature, critical pressure, boiling point temperature and critical volume respectively.
