summaryrefslogtreecommitdiff
path: root/1019/CH6/EX6.13/Example_6_13.sce
blob: 4f3c9c54b1d35a02c0e5c2f871988db9bb6c5ddc (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
//Example 6.13
clear;
clc;

//Given
T=298;//temperature in K
P=1;//pressure in atm
I=1.9373;//moment of inertia in 10^(-46) kg m^2
R=8.314;//gas constant in J K^-1 mol^-1
h=6.626;//plancks constant in 10^(34) Js
k=1.38;//in (10^(-23)) J K^-1

//To determine the rotational contributions to entropy and free energy
Qr=(8*(%pi^2)*I*k)/(20*(h^2));//rotational transition function
Sr=R*(1+log(Qr*T));//rotational contributions to entropy in J K^-1 mol^-1
Gr=((R*T)-(T*Sr))*0.001;//rotational contributions to free energy in J mol^-1
mprintf('rotational contributions to entropy = %f J K^-1 mol^-1',Sr);
mprintf('\n rotational contributions to free energy = %f kJ mol^-1',Gr);
//end