blob: 9116102744f8b48e232b331cda6399ec894697a9 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
//Example 7.2
clear;
clc;
//Given
R=0.08205;//gas constant in dm^3 atm K^-1 mol^-1
b=0.0391;//Van der Waals constant in dm^3 mol^-1
T=1273;//Temperature in K
P=1000;//pressure in atm
//To calculate the fugacity coefficient
k=(b*P)/(R*T);//k=log(f/P)
f=P*exp(k);//fugacity coefficient in atm
mprintf('Fugacity coefficient = %f atm',f);
//end
|
