summaryrefslogtreecommitdiff
path: root/623/CH27
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /623/CH27
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '623/CH27')
-rwxr-xr-x623/CH27/EX5.5.1/U5_C5_1.sce22
-rwxr-xr-x623/CH27/EX5.5.10/U5_C5_10.sce11
-rwxr-xr-x623/CH27/EX5.5.11/U5_C5_11.sce10
-rwxr-xr-x623/CH27/EX5.5.12/U5_C5_12.sce14
-rwxr-xr-x623/CH27/EX5.5.13/U5_C5_13.sce15
-rwxr-xr-x623/CH27/EX5.5.14/U5_C5_14.sce14
-rwxr-xr-x623/CH27/EX5.5.15/U5_C5_15.sce19
-rwxr-xr-x623/CH27/EX5.5.16/U5_C5_16.sce16
-rwxr-xr-x623/CH27/EX5.5.17/U5_C5_17.sce12
-rwxr-xr-x623/CH27/EX5.5.18/U5_C5_18.sce14
-rwxr-xr-x623/CH27/EX5.5.19/U5_C5_19.sce12
-rwxr-xr-x623/CH27/EX5.5.2/U5_C5_2.sce21
-rwxr-xr-x623/CH27/EX5.5.20/U5_C5_20.sce11
-rwxr-xr-x623/CH27/EX5.5.21/U5_C5_21.sce18
-rwxr-xr-x623/CH27/EX5.5.3/U5_C5_3.sce28
-rwxr-xr-x623/CH27/EX5.5.4/U5_C5_4.sce10
-rwxr-xr-x623/CH27/EX5.5.5/U5_C5_5.sce26
-rwxr-xr-x623/CH27/EX5.5.6/U5_C5_6.sce9
-rwxr-xr-x623/CH27/EX5.5.7/U5_C5_7.sce12
-rwxr-xr-x623/CH27/EX5.5.8/U5_C5_8.sce29
-rwxr-xr-x623/CH27/EX5.5.9/U5_C5_9.sce10
21 files changed, 333 insertions, 0 deletions
diff --git a/623/CH27/EX5.5.1/U5_C5_1.sce b/623/CH27/EX5.5.1/U5_C5_1.sce
new file mode 100755
index 000000000..68c846bbb
--- /dev/null
+++ b/623/CH27/EX5.5.1/U5_C5_1.sce
@@ -0,0 +1,22 @@
+//function for calculating the wave number
+function[wave]=F(j)
+ wave=B*j*(j+1);
+endfunction
+
+//variable initialization
+r=1.21*10^-10; //internuclear distance (meter)
+m=2.7*10^-26; //mass of oxygen atom (kg)
+h=6.626*10^-34; //Plank's constant (joule second)
+c=3*10^8; //speed of light (meter/second)
+
+//(a) moment of inertia
+mu=m/2; //reduced mass (kg)
+I=mu*r^2; //moment of inertia (kg m^2)
+
+//(b) rotational constant
+B=h/(8*%pi^2*I*c); //rotational constant (m-1)
+
+//(c) wave number
+waveno=F(1)-F(0); //wave no. of the line corresponding to the transition J=0 to J=1 (m-1)
+
+printf("\n(a) I = %.3e kg m^2\n(b) B = %.1f m-1\n(c) wave number = %.0f m-1",I,B,waveno);
diff --git a/623/CH27/EX5.5.10/U5_C5_10.sce b/623/CH27/EX5.5.10/U5_C5_10.sce
new file mode 100755
index 000000000..373af0166
--- /dev/null
+++ b/623/CH27/EX5.5.10/U5_C5_10.sce
@@ -0,0 +1,11 @@
+//variable initialization
+h=6.6*10^-34; //Plank's constant (joule second)
+mu=1.62*10^-27; //reduced mass of HCl molecule (kg)
+c=3*10^8; //speed of light (m/s)
+v=2.886*10^5; //wave no. of absorption line in infrared spectrum (m-1)
+
+//calculation of amplitude of vibration
+k=4*(%pi*c*v)^2*mu; //force constant of HCl molecule (N/m)
+amp=sqrt((h*c*v)/k)*10^10; //amplitude of vibration in the ground state (Å)
+
+printf("\namplitude of vibration = %.2f Å",amp);
diff --git a/623/CH27/EX5.5.11/U5_C5_11.sce b/623/CH27/EX5.5.11/U5_C5_11.sce
new file mode 100755
index 000000000..c71463c42
--- /dev/null
+++ b/623/CH27/EX5.5.11/U5_C5_11.sce
@@ -0,0 +1,10 @@
+//variable initialization
+v1=214330 //fundamental band for CO molecule (m-1)
+v2=425970 //first overtone for CO molecule (m-1)
+A=[1 -2;2 -6]; //coefficient matrix
+b=[v1;v2]; //right hand side matrix
+
+//calculation
+x=inv(A)*b; //values of omega and x*omega (m-1)
+
+printf("\nωe = %.0f m-1\nxe*ωe = %.0f m-1",x(1),x(2));
diff --git a/623/CH27/EX5.5.12/U5_C5_12.sce b/623/CH27/EX5.5.12/U5_C5_12.sce
new file mode 100755
index 000000000..0241fe411
--- /dev/null
+++ b/623/CH27/EX5.5.12/U5_C5_12.sce
@@ -0,0 +1,14 @@
+//variable initialization
+v1=288600 //intense absorption (m-1)
+v2=566800 //intense absorption (m-1)
+v3=834700 //intense absorption (m-1)
+A=[1 -2;2 -6]; //coefficient matrix
+b=[v1;v2]; //right hand side matrix
+mu=1.61*10^-27; //reduced mass (kg)
+c=3*10^8; //speed of light (m/s)
+
+//calculation
+x=inv(A)*b; //values of omega and x*omega (m-1)
+k=4*(%pi*c*x(1))^2*mu; //force constant (N/m)
+
+printf("\nωe = %.0f m-1\nxe*ωe = %.0f m-1\nforce constant = %.1f N/m",x(1),x(2),k);
diff --git a/623/CH27/EX5.5.13/U5_C5_13.sce b/623/CH27/EX5.5.13/U5_C5_13.sce
new file mode 100755
index 000000000..b5b30c979
--- /dev/null
+++ b/623/CH27/EX5.5.13/U5_C5_13.sce
@@ -0,0 +1,15 @@
+//variable initialization
+v1=8.657*10^13; //frequency of rotation absorption spectrum (Hz)
+v2=8.483*10^13; //frequency of rotation absorption spectrum (Hz)
+h=6.6*10^-34; //Plank's constant (joule second)
+mu=1.544*10^-27; //Recuced mass of CH molecule (kg)
+
+//(i) equilibrium separation
+I=h/(2*%pi^2*(v1-v2)); //Moment of inertia (kg m^2)
+r=sqrt(I/mu); //equilibrium internuclear distance (m)
+
+//(ii) force constant of molecule
+v0=(v1+v2)/2; //Central frequency (Hz)
+k=4*mu*(%pi*v0)^2; //Force constant of CH molecule (N/m)
+
+printf("\n (i) equilibrium separation = %.2e meter\n (ii) force constant = %.0f N/m",r,k);
diff --git a/623/CH27/EX5.5.14/U5_C5_14.sce b/623/CH27/EX5.5.14/U5_C5_14.sce
new file mode 100755
index 000000000..ac613b118
--- /dev/null
+++ b/623/CH27/EX5.5.14/U5_C5_14.sce
@@ -0,0 +1,14 @@
+//variable initialization
+k=448 //force constant of CH molecule (N/m)
+mu=4.002*10^-27; //reduced mass of CH molecule (kg)
+r=0.112*10^-9; //internuclear distance (m)
+h=6.6*10^-34; //Plank's constant (joule second)
+
+//Calculation of peak frequencies
+v0=(1/(2*%pi))*sqrt(k/mu); //central frequency (s-1)
+I=mu*r^2; //moment of inertia of molecule (kg m^2)
+x=h/(4*%pi^2*I); //additional frequency (s-1)
+v1=v0+x; //peak frequency (Hz)
+v2=v0-x; //peak frequency (Hz)
+
+printf("\n Peak frequencies = %.3e Hz, %.3e Hz",v1,v2);
diff --git a/623/CH27/EX5.5.15/U5_C5_15.sce b/623/CH27/EX5.5.15/U5_C5_15.sce
new file mode 100755
index 000000000..1303183a7
--- /dev/null
+++ b/623/CH27/EX5.5.15/U5_C5_15.sce
@@ -0,0 +1,19 @@
+//variable initialization
+v1=2174.07 //peak wave number (cm-1)
+v2=2166.35 //peak wave number (cm-1)
+h=6.6*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+mu=1.145*10^-26; //Reduced mass of CO molecule (kg)
+
+//(a) central frequency
+B=(v1-v2)/4; //Rotational constant (cm-1)
+v0=(v1+v2)/2; //Central frequency (cm-1)
+
+//(b) internuclear distance
+I=h/(8*%pi^2*B*100*c); //moment of inertia of molecule (kg m^2)
+r=sqrt(I/mu)*10^10; //equilibrium internuclear distance (Å)
+
+//(c) force constant
+k=4*mu*(%pi*c*v0*100)^2; //force constant (N/m)
+
+printf("\n(a) central frequency = %.2f cm-1\n(b) internuclear distance = %.2f Å\n(c) force constant = %.0f N/m",v0,r,k);
diff --git a/623/CH27/EX5.5.16/U5_C5_16.sce b/623/CH27/EX5.5.16/U5_C5_16.sce
new file mode 100755
index 000000000..3b46c5cf0
--- /dev/null
+++ b/623/CH27/EX5.5.16/U5_C5_16.sce
@@ -0,0 +1,16 @@
+//variable initialization
+mu=3.142*10^-27; //reduced mass of the molecule (kg)
+r=1.288*10^-10; //internuclear distance (m)
+h=6.6*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+v0=201100 //central frequency (m-1)
+
+//Calculation
+I=mu*r^2; //moment of inertia of molecule (kg m^2)
+B=h/(8*%pi^2*I*c); //Rotational constant (m-1)
+vR0=v0+(2*B); //wave no. of 1st line of R-branch (m-1)
+vR1=v0+(4*B); //wave no. of 2nd line of R-branch (m-1)
+vP1=v0-(2*B); //wave no. of 1st line of P-branch (m-1)
+vP2=v0-(4*B); //wave no. of 2nd line of P-branch (m-1)
+
+printf("\n V_R(0) = %.0f m-1\n V_R(1) = %.0f m-1\n V_P(1) = %.0f m-1\n V_P(2) = %.0f m-1",vR0,vR1,vP1,vP2);
diff --git a/623/CH27/EX5.5.17/U5_C5_17.sce b/623/CH27/EX5.5.17/U5_C5_17.sce
new file mode 100755
index 000000000..24ddadff4
--- /dev/null
+++ b/623/CH27/EX5.5.17/U5_C5_17.sce
@@ -0,0 +1,12 @@
+//variable initialization
+mu=1.62*10^-27; //reduced mass of HCl molecule (kg)
+r=1.293*10^-10; //internuclear distance (m)
+h=6.6*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+
+//Calculation of separation between lines R(0) and P(1) of the fundamental band of HCl 35
+I=mu*r^2; //moment of inertia of molecule (kg m^2)
+B=h/(8*%pi^2*I*c); //Rotational constant (m-1)
+sep=4*B; //separation b/w lines R(0) and P(1) (m-1)
+
+printf("\nΔν = %.0f m-1",sep);
diff --git a/623/CH27/EX5.5.18/U5_C5_18.sce b/623/CH27/EX5.5.18/U5_C5_18.sce
new file mode 100755
index 000000000..d9db6f9ad
--- /dev/null
+++ b/623/CH27/EX5.5.18/U5_C5_18.sce
@@ -0,0 +1,14 @@
+//variable initialization
+a=214.6*100; //(m-1)
+b=0.6*100; //(m-1)
+h=6.6*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+no=1/(%e); //number of molecules in state with respect to ground state
+k=1.38*10^-23; //Boltzmann constant (J K-1)
+
+//Calculation
+deltaE=h*c*(a-2*b); //difference in the energies of state 0 and state 1 (J)
+T1=deltaE/k; //temperature at which number of molecules in state 1 is 1/e times of state 0 (K)
+T2=deltaE/(k*log(10)); //temperature at which number of molecules in state 1 is 10% of state 0 (K)
+
+printf("n(i) T = %.0f K\n(ii) T = %.0f K",T1,T2);
diff --git a/623/CH27/EX5.5.19/U5_C5_19.sce b/623/CH27/EX5.5.19/U5_C5_19.sce
new file mode 100755
index 000000000..4aa796298
--- /dev/null
+++ b/623/CH27/EX5.5.19/U5_C5_19.sce
@@ -0,0 +1,12 @@
+//variable initialization
+vexc=4358*10^-10; //wavelength of exciting line (m)
+vsto=4458*10^-10; //wavelength of Stokes line (m)
+
+//calculation of wavelength of Anti-stokes line
+vbar_exc=1/vexc; //wave number of exciting line (m-1)
+vbar_sto=1/vsto; //wave number of Stokes line (m-1)
+delta_vbar=vbar_exc-vbar_sto; //Raman shift (m-1)
+vbar_antistoke=vbar_exc+delta_vbar; //Wave number of Anti-Stokes line (m-1)
+lembda_antistoke=(1/vbar_antistoke)*10^10; //Wavelength of Anti-Stokes line (Å)
+
+printf("\nwavelength of Anti-stokes line = %.1f Å",lembda_antistoke);
diff --git a/623/CH27/EX5.5.2/U5_C5_2.sce b/623/CH27/EX5.5.2/U5_C5_2.sce
new file mode 100755
index 000000000..5e712bd26
--- /dev/null
+++ b/623/CH27/EX5.5.2/U5_C5_2.sce
@@ -0,0 +1,21 @@
+//function for calculating the energy level
+function[energy]=F(j)
+ energy=a*j*(j+1);
+endfunction
+
+//variable initialization
+m=1.6738*10^-27; //mass of hydrogen atom (kg)
+r=0.74*10^-10; //intermolecular distance of hydrogen molecule (meter)
+h=1.054*10^-34; //Planck's constant (joule second)
+e=1.6*10^-19; //Charge of electron (coulombs)
+
+//calculation of rotational energy levels
+mu=m/2; //reduced mass of hydrogen atom (kg)
+I=mu*r^2; //moment of inertia of molecule (kg meter^2)
+a=h^2/(2*I*e); //constant (eV)
+E0=F(0); //energy of level 0 (eV)
+E1=F(1); //energy of level 1 (eV)
+E2=F(2); //energy of level 2 (eV)
+E3=F(3); //energy of level 3 (eV)
+
+printf("\nE0 = %.0f\nE1 = %.2e eV\nE2 = %.2e eV\nE3 = %.2e eV",E0,E1,E2,E3);
diff --git a/623/CH27/EX5.5.20/U5_C5_20.sce b/623/CH27/EX5.5.20/U5_C5_20.sce
new file mode 100755
index 000000000..9bf6d7c62
--- /dev/null
+++ b/623/CH27/EX5.5.20/U5_C5_20.sce
@@ -0,0 +1,11 @@
+//variable initialization
+h=6.62*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+x=62.4*100; //(m-1)
+y=41.6*100; //(m-1)
+
+//calculation of the moment of inertia of HCl molecule
+B=y/4; //Rotational constant of HCl (m-1)
+I=h/(8*%pi^2*B*c); //Moment of inertia (kg m^2)
+
+printf("\n I = %.1e kg m^2",I);
diff --git a/623/CH27/EX5.5.21/U5_C5_21.sce b/623/CH27/EX5.5.21/U5_C5_21.sce
new file mode 100755
index 000000000..05d4aac6d
--- /dev/null
+++ b/623/CH27/EX5.5.21/U5_C5_21.sce
@@ -0,0 +1,18 @@
+//function for calculating the vibrational energy of O2 molecule
+function[energy]=F(v)
+ energy=(((v+.5)*a)-(((v+.5)^2)*b))*h*c;
+endfunction
+
+//variable initialization
+h=6.62*10^-34; //Plank's constant (joule second)
+c=3*10^8; //Speed of light (m/s)
+a=1580.36*100; //value of ωe (m-1)
+b=12.07*100; //value of ωexe (m-1)
+
+//Calculation of zero point energy
+E0=F(0); //Zero point energy of the molecule (J)
+
+//Calculation of vibrational Raman shift
+shift=(F(1)-F(0))/(h*c); //Expected vibrational Raman shift (m-1)
+
+printf("\nZero point energy = %.3e J\nExpected vibrational Raman shift = %.0f m-1",E0,shift);
diff --git a/623/CH27/EX5.5.3/U5_C5_3.sce b/623/CH27/EX5.5.3/U5_C5_3.sce
new file mode 100755
index 000000000..b3f1583f5
--- /dev/null
+++ b/623/CH27/EX5.5.3/U5_C5_3.sce
@@ -0,0 +1,28 @@
+//variable initialization
+u=1.68*10^-27; //mass of hydrogen atom (kg)
+m1=16; //mass of oxygen atom in terms u
+m2=1; //mass of hydrogen atom in terms of u
+I=1.48*10^-47; //moment of inertia of OH-radical (kg m^2)
+h_bar=1.054*10^-34; //Planck's constant (joule second)
+j=5; //energy level of OH-radical
+c=3*10^8; //speed of light (meter/second)
+h=6.626*10^-34; //Plank's constant (joule second)
+
+//(a) internuclear distance
+mu=((m1*m2)/(m1+m2))*u; //reduced mass of the molecule (kg)
+r=(sqrt(I/mu))*10^10; //internuclear distance of molecule (Å)
+
+//(b) angular momentum
+P=h_bar*sqrt(j*(j+1)); //angular momentum of molecule (joule second)
+
+//(c) angular velocity
+omega=P/I; //angular velocity of molecule (radian/second)
+
+//(d) wave number
+B=h/(8*%pi^2*I*c); //rotational constant (m-1)
+no=2*B*(j+1); //wave no. of line corresponding to transition j=5 to j=6 (m-1)
+
+//(e) energy absorbed
+E=c*h*no; //energy absorbed in the transition j=6 to j=5 (joule)
+
+printf("\n(a) r = %.3f Å\n(b) J = %.2e joule second\n(c) ω = %.2e radian/second\n(d) wave number = %.2e m-1\n(e) E = %.1e joule",r,P,omega,no,E);
diff --git a/623/CH27/EX5.5.4/U5_C5_4.sce b/623/CH27/EX5.5.4/U5_C5_4.sce
new file mode 100755
index 000000000..233af4fe0
--- /dev/null
+++ b/623/CH27/EX5.5.4/U5_C5_4.sce
@@ -0,0 +1,10 @@
+//variable initialization
+h=6.63*10^-34; //Plank's constant (joule second)
+v=1.153*10^11; //Frequency of absorption line (Hz)
+mu=11.38*10^-27; //Recuced mass of the molecule (kg)
+
+//Calculation of the internuclear distance
+I=h/(4*%pi^2*v); //moment of inertia of CO molecule (kg m^2)
+r=sqrt(I/mu)*10^10; //Internuclear distance (Å)
+
+printf("\n Internuclear distance = %.2f Å",r);
diff --git a/623/CH27/EX5.5.5/U5_C5_5.sce b/623/CH27/EX5.5.5/U5_C5_5.sce
new file mode 100755
index 000000000..063f671c4
--- /dev/null
+++ b/623/CH27/EX5.5.5/U5_C5_5.sce
@@ -0,0 +1,26 @@
+//variable initialization
+mu=1.62*10^-27; //Reduced mass of HCL (kg)
+c=3*10^8; //Velocity of light (m/s)
+h=6.62*10^-34; //Plank's constant (joule second)
+v1_P=2906.3 //Wave no. of P branch (cm-1)
+v2_P=2927.5 //Wave no. of P branch (cm-1)
+v3_P=2948.7 //Wave no. of P branch (cm-1)
+v4_P=2969.9 //Wave no. of P branch (cm-1)
+v1_R=3012.2 //Wave no. of R branch (cm-1)
+v2_R=3033.4 //Wave no. of R branch (cm-1)
+v3_R=3054.6 //Wave no. of R branch (cm-1)
+v4_R=3078.8 //Wave no. of R branch (cm-1)
+
+//(i) Equilibrium internuclear distance
+delta_v=v2_P-v1_P; //Separation of successive line of P and R branch (cm-1)
+B=delta_v/2; //rotational constant (cm-1)
+I=h/(8*%pi^2*B*10^2*c); //Moment of inertia (kg m^2)
+r=sqrt(I/mu)*10^10; //Equilibrium internuclear distance (Å)
+
+//(ii) Force constant
+v0=(v4_P+v1_R)/2; //Equlibrium frequency (cm-1)
+k=4*%pi^2*mu*c^2*v0^2*10^4; //Force constant of HCl (N/m)
+
+printf("\n(i) Equilibrium internuclear distance = %.2f Å\n(ii) Force constant = %.0f N/m",r,k);
+
+//Note: the answer of (ii) part is wrong in the book
diff --git a/623/CH27/EX5.5.6/U5_C5_6.sce b/623/CH27/EX5.5.6/U5_C5_6.sce
new file mode 100755
index 000000000..5ed6e1bd6
--- /dev/null
+++ b/623/CH27/EX5.5.6/U5_C5_6.sce
@@ -0,0 +1,9 @@
+//variable initialization
+mu=8.37*10^-28; //Reducec mass of hydrogen molecule (kg)
+h=6.58*10^-16; //Plank's constant (eV s)
+E0=0.273 //Ground state vibrational energy of hydrogen molecule (eV)
+
+//calculation of force constant of the molecule
+k=mu*((2*E0)/h)^2; //force constant of hydrogen molecule (N/m)
+
+printf("\n Force constant = %.0f N/m",k);
diff --git a/623/CH27/EX5.5.7/U5_C5_7.sce b/623/CH27/EX5.5.7/U5_C5_7.sce
new file mode 100755
index 000000000..273cfbc1c
--- /dev/null
+++ b/623/CH27/EX5.5.7/U5_C5_7.sce
@@ -0,0 +1,12 @@
+//variable initialization
+m1=1; //molar mass of H atom (amu)
+m2=35; //molar mass of Cl atom (amu)
+u=1.68*10^-27; //atomic mass unit (kg)
+v=2885.9*100; //wave no. of line (m-1)
+c=3*10^8; //Velocity of light (m/s)
+
+//calculation of force constant
+mu=((m1*m2)/(m1+m2))*u; //reduced mass of HCl molecule (kg)
+k=4*(%pi*c*v)^2*mu; //force constant of HCl molecule (N/m)
+
+printf("\n force constant = %.2f N/m",k);
diff --git a/623/CH27/EX5.5.8/U5_C5_8.sce b/623/CH27/EX5.5.8/U5_C5_8.sce
new file mode 100755
index 000000000..4036a404f
--- /dev/null
+++ b/623/CH27/EX5.5.8/U5_C5_8.sce
@@ -0,0 +1,29 @@
+//function for calculating the various vibrational energy levels of CO molecule
+function[energy]=F(V)
+ energy=((V+.5)*h*v)/e;
+endfunction
+
+//function for converting eV to cm-1
+function[energy]=G(V)
+ energy=((V+.5)*h*v*8065)/e;
+endfunction
+
+//variable initialization
+m1=12; //molar mass of C atom (amu)
+m2=16; //molar mass of O atom (amu)
+u=1.68*10^-27; //atomic mass unit (kg)
+k=1870 //force constant of CO molecule (N/m)
+h=6.6*10^-34; //Plank's constant (joule second)
+e=1.602*10^-19; //charge of electron (Coulomb)
+
+//calculation of energy levels
+mu=((m1*m2)/(m1+m2))*u; //reduced mass of CO molecule (kg)
+v=(1/(2*%pi))*sqrt(k/mu); //frequency of vibration of CO molecule (Hz)
+e1=F(0); //energy of 1st level of CO molecule (eV)
+E1=G(0); //energy of 1st level of CO molecule (cm-1)
+e2=F(1); //energy of 2nd level of CO molecule (eV)
+E2=G(1); //energy of 2nd level of CO molecule (cm-1)
+e3=F(2); //energy of 3rd level of CO molecule (eV)
+E3=G(2); //energy of 3rd level of CO molecule (cm-1)
+
+printf("\nE = %.3f eV, %.3f eV, %.3f eV..........................\n = %.1f cm-1, %.1f cm-1, %.1f cm-1.................",e1,e2,e3,E1,E2,E3);
diff --git a/623/CH27/EX5.5.9/U5_C5_9.sce b/623/CH27/EX5.5.9/U5_C5_9.sce
new file mode 100755
index 000000000..ccd1e3f5e
--- /dev/null
+++ b/623/CH27/EX5.5.9/U5_C5_9.sce
@@ -0,0 +1,10 @@
+//variable initialization
+mu=1.61*10^-27; //reduced mass of HCl molecule (kg)
+c=3*10^8; //speed of light (m/s)
+lembda=3.465*10^-6; //wavelength of infrared radiation (m)
+
+//calculation of force constant
+v=c/lembda; //frequency of radiation (s-1)
+k=4*(%pi*v)^2*mu; //force constant of HCl molecule (N/m)
+
+printf("\nforce constant = %.1f N/m",k);