function PVT
% 已知P和T，計算V
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/26 $

clear all
clc

P = 9.33;   % atm
T = 300.2;  % K
n = 2;      % mol
a = 4.17;
b = 0.0371;
R = 0.08206;
V0 = n*R*T/P
[V,fval] = fzero(@PVTeq,V0,[],P,T,n,a,b,R)

% ------------------------------------------------------------------
function f = PVTeq(V,P,T,n,a,b,R)
f = (P + a*n^2/V^2) * (V-n*b) - n*R*T;