function v3 = v_redlich(vol);
global R T P Tc Pc

a = 0.42748*R^2*Tc^2.5/Pc;
b = 0.08664*R*Tc/Pc;

v3 = R*T - (P + a/(T^0.5*vol*(vol+b)))*(vol - b);
%-----------