POLYMATH Report          NLE 
 Explicit Equations 2016-May-31 

Calculated values of explicit variables
    Variable Value
1 T 250
2 LCP_CO2 8.82E+04
3 LDN_CO2 23.7577
4 LTC_CO2 0.136225
5 LVS_CO2 0.000155
6 VP_CO2 1.79E+06

Problem source text
# S. 16 - Complex consecutive calculations
# Physical property calculator for liquid CO2
# Verified Solution: LCP = 8.82e+4, LDN = 23.7577
# LTC =0.136225 , LVS = 0.000155, VP = 1.79E+6
# Source of the correlations: DIPPR database
LCP_CO2 = -8304300 + 104370 * T - 433.33 * T ^ 2 + 0.60052 * T ^ 3 # J/kmol*K
LDN_CO2 = 2.768 / 0.26212 ^ (1 + (1 - T / 304.21) ^ 0.2908) # kmol/m^3
LTC_CO2 = 0.4406 - 0.0012175 * T # W/m*K
LVS_CO2 = exp(18.775 - 402.92 / T - 4.6854 * ln(T) - 6.9171E-26 * T ^ 10) # Pa*s
VP_CO2 = exp(140.54 - 4735 / T - 21.268 * ln(T) + 0.040909 * T) # Pa
# Valid temperature range: 217 K - 304 K
T= 250 # Temperature (K)

Matlab formatted problem
Create m file called PolyExplicit.m and paste the following text into it.
% S. 16 - Complex consecutive calculations
% Physical property calculator for liquid CO2
% Verified Solution: LCP = 8.82e+4, LDN = 23.7577
% LTC =0.136225 , LVS = 0.000155, VP = 1.79E+6
% Source of the correlations: DIPPR database
%Temperature (K)
T = 250;
%J/kmol*K
LCP_CO2 = -8304300 + 104370 * T - (433.33 * T ^ 2) + 0.60052 * T ^ 3;
%kmol/m^3
LDN_CO2 = 2.768 / (0.26212 ^ (1 + (1 - (T / 304.21)) ^ 0.2908));
%W/m*K
LTC_CO2 = 0.4406 - (0.0012175 * T);
%Pa*s
LVS_CO2 = exp(18.775 - (402.92 / T) - (4.6854 * log(T)) - (6.9171E-26 * T ^ 10));
%Pa
VP_CO2 = exp(140.54 - (4735 / T) - (21.268 * log(T)) + 0.040909 * T);

General Settings
number of explicit equations: 6
reporting digits 8