Complete in Matlab, MATLAB, MATLAB

 

% Given data

T = 350; % Temperature in Kelvin

P = 9.4573; % Pressure in bar

Tc = 425.1; % Critical temperature in Kelvin

Pc = 37.96; % Critical pressure in bar

R = 0.08314; % Universal gas constant in bar L / (mol K)

omega_small = 0.200; % Acentric factor

% Calculate reduced temperature and reduced pressure

Tr = T / Tc;

Pr = P / Pc;

% Calculate a

a = (1 + (0.37464 + 1.54226 * omega_small - 0.26992 * omega_small^2) * (1 - sqrt(Tr)))^2;

%Calculate q

psi=0.45724;

omega_big =0.07780;

a_T=psi*a*(R^2*Tc^2/Pc);

b=omega_big*R*Tc/Pc;

q=(a_T)/(b*R*T);

%Calculate beta

beta=(b*P)/(R*T);

%Parameters

epsilon=-0.414;

alpha_small=2.214;

E=epsilon*beta;

O=alpha_small*beta;

q_B=q*beta;

%Finding Z by finding the roots

Z=1+beta-(q_B*(Z-beta)/(Z+E)*(Z+O));

% Select the largest root for vapor phase and the smallest root for liquid phase

V_vapor=Z1*R*T/P;

V_liquid=Z3*R*T/P;