function v2 = v_vander(vol);
global R T P Tc Pc

a = 27*(R^2)*(Tc^2)/(64*Pc);

b = R*Tc/(8*Pc);

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