//clc()
P = 10^6;//Pa
T = 373;//K
Tc = 405.5;//K
Pc = 112.8 * 10 ^ 5//Pa
R = 8.314;
a = 0.4278 * (R^2) * (Tc ^ 2.5)/Pc;
b = 0.0867 * R * Tc / Pc;
//P1 = (R*T/(V - b) - a/((T^0.5)*V*(V + b)))/10^5;
//10^6=((8.314*373)/(V-2.59*10^-5))-8.68/((373^0.5)*V*(V+2.59*10^-5) 
//solving this we get,
V = 3.0 * 10^-3;//m^3/mol
disp("m^3/mol",V,"molar volume of gas = ")