blob: 37048e27a47659b3dfd0a31cc55a747d7e69e8b2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
clc
v1 = 1 // Assumed
v2 = v1// Assumed
v3 = v2 // Assumed
v4 = v2// Assumed
e = 0.56 // Degree of reaction
P = 1 // Dummy
T = 1200 // Reaction temperature in K
R = 8.3143 // Gas constant
x1 = (1-e)/2 //
x2 = (1-e)/2
x3 = e/2
x4 = e/2
K = (((x3^v3)*(x4^v4))/((x1^v1)*(x2^v2)))*P^(v3+v4-v1-v2) // Equilibrium constant
dG = -R*T*log(K) //Gibbs function change
printf("\n Example 16.3\n")
printf("\n Equilibrium constant is %f",K)
printf("\n Gibbs function change is %fJ/gmol",dG)
//The answers vary due to round off error
|
