function cp=licl_birerong(t,x)
T=t+273.15;
tao=T/228-1;
a=[1.43980;-1.24317;-0.12070;0.12825;0.62934;58.5225;-105.6343;47.7948];
b=[88.7891;-120.1958;-16.9264;52.4654;0.10826;0.46988];
cps=b(1)+b(2)*tao^0.02+b(3)*tao^0.04+b(4)*tao^0.06+b(5)*tao^1.8+b(6)*tao^8;
f1=a(6)*tao^0.02+a(7)*tao^0.04+a(8)*tao^0.06;
f2=a(4)+a(5)*x;
cp=cps*(1-f1*f2);