can you correct this code on matlab and compleat the code to find pressure saturation and fugacity liquid and fugacity for vapor when philiq equal phivap for each tempretaure.$($Tc $=658\mathrm{.}3$; $\%$ Critical temperature in K

Pc $=1\mathrm{.}817$e$6$; $\%$ Critical pressure in Pa

omega $=0\mathrm{.}562$; $\%$ Acentric factor

$\%$ Universal gas constant

R $=8\mathrm{.}314$; $\%$ J$/($mol K$)$

b $=0\mathrm{.}07780*$ R $*$ Tc $/$ Pc;

kappa $=0\mathrm{.}37464+1\mathrm{.}54226*$ omega $-0\mathrm{.}26992*$ omega$^2$;

$\%$ Reduced temperatures

Tr $=[0\mathrm{.}7,0\mathrm{.}8,0\mathrm{.}9,1\mathrm{.}0]$;

$\%$loop

for i$=1$:length$($T$)$

T$=$Tr$($i$)$;

Tr$=$ T $/$ Tc;

alpha $=(1+$ kappa $*(1-$ sqrt$($T $/$ Tc$)))^2$;

a $=0\mathrm{.}45724*$ R$^2*$ Tc$^2/$ Pc $*$ alpha;

$\%$initial guess for psat $($pa$)$

if Tr $==0\mathrm{.}7$

psat$=48418$;

end

if Tr$==0\mathrm{.}8$

psat$=227470$;

end

if Tr$==0\mathrm{.}9$

psat$=711240$;

end

if Tr$==1$

psat $=1704911$;

end

tolerance $=1$e$-4$;

error$=0$;

while error $<$ tolerance

$\%$solve

A$=$ a$*$psat$/($R $*$ T$)^2$;

B $=$ b$*$psat $/($R $*$ T$)$;

coeffs $=[1,-1+$B$,$ A$-(2*$B$)-(3*$B$^2),-($A$*$B$)+($B$^2)+($B$^3)]$;

z$=$roots$($coeffs$)$;

z$=$real $($z$($imag$($z$)==0))$; $\%$keep one real root

z$\_$min $=$ min$($z$)$;

z$\_$max $=$ max$($z$)$;

$\%$calculate fugacity coefficient

philiq $=$ exp$(($z$\_$min$-1)-$log$($z$\_$min$-$B$)-($A$/(2*$sqrt$(2)*$B$))*$log$(($z$\_$min$+(1+$sqrt$(2))*$B$)/($z$\_$min$+(1-$sqrt$(2))*$B$)))$;

phivap $=$ exp$(($z$\_$max$-1)-$log$($z$\_$max$-$B$)-($A$/(2*$sqrt$(2)*$B$))*$log$(($z$\_$max$+(1+$sqrt$(2))*$B$)/($z$\_$max$+(1-$sqrt$(2))*$B$)))$;

$\%$calculate psat new

psatnew $=$ psat$*($philiq$/$phivap$)$;

error $=$ abs$(($philiq$/$phivap$)-1)$;

psat$=$ psatnew ;

end

end

$\%$calculate volumes

Vliq $=$ z$\_$min$*$R$*$T$/$psat;

Vvap $=$ z$\_$max$*$R$*$T$/$psat;$)$