函数[z] = prsv (T, P, Tc, Pc, w, k1, x, kij)
nrosust =长度(Tc);
R = 8.314772;
为我= 1:nrosust
k0 (i) = 0.378893 + 1.4897153 * 0.17131848 w - * w。^ 2 + 0.0196554 * w。^ 3;
k (i) = k0 (i) + k1 * (1 + (T / Tc (i)) ^ 0.5) * (0.7 - (T / Tc (i)));
α(i) = (1 + k(我)* (1 - ((T / Tc (i)) ^ 0.5))) ^ 2;
(我)= (0.457235 * (R ^ 2 * Tc (i) ^ 2 / Pc (i))) *α(我);
b (i) = 0.07779607 * (R * Tc(我)/ Pc (i));
结束
sum_a = 0;
为我= 1:nrosust
为j = 1: nrosust
aij (i, j) =((我)* (j)) ^ 0.5 * (1-kij (i, j));
sum_a = sum_a + x(我)* (j) * aij (i, j);
结束
结束
b = (x。* b)总和;
= sum_a * P / R * T ^ 2;
B = B * P / (R (T);
波尔= [1 1 B + B a - 3 * ^ 2 * B - a * B + B B ^ ^ 2 + 3];
Z =根(Pol);
z = max (z);
