summaryrefslogtreecommitdiff
path: root/3411/CH5
diff options
context:
space:
mode:
Diffstat (limited to '3411/CH5')
-rw-r--r--3411/CH5/EX5.10.u1/Ex5_10_u1.sce14
-rw-r--r--3411/CH5/EX5.10.u1/Ex5_10_u1.txt2
-rw-r--r--3411/CH5/EX5.11.u1/Ex5_11_u1.sce74
-rw-r--r--3411/CH5/EX5.11.u1/Ex5_11_u1.txt13
-rw-r--r--3411/CH5/EX5.12.u1/Ex5_12_u1.sce51
-rw-r--r--3411/CH5/EX5.12.u1/Ex5_12_u1.txt10
-rw-r--r--3411/CH5/EX5.13.u1/Ex5_13_u1.sce96
-rw-r--r--3411/CH5/EX5.13.u1/Ex5_13_u1.txt18
-rw-r--r--3411/CH5/EX5.14.u1/Ex5_14_u1.sce10
-rw-r--r--3411/CH5/EX5.15.u1/Ex5_15_u1.sce8
-rw-r--r--3411/CH5/EX5.16.u1/Ex5_16_u1.sce18
-rw-r--r--3411/CH5/EX5.16.u1/Ex5_16_u1.txt4
-rw-r--r--3411/CH5/EX5.9.u1/Ex5_9_u1.sce14
-rw-r--r--3411/CH5/EX5.9.u1/Ex5_9_u1.txt4
14 files changed, 255 insertions, 81 deletions
diff --git a/3411/CH5/EX5.10.u1/Ex5_10_u1.sce b/3411/CH5/EX5.10.u1/Ex5_10_u1.sce
index 8689098f1..578eda983 100644
--- a/3411/CH5/EX5.10.u1/Ex5_10_u1.sce
+++ b/3411/CH5/EX5.10.u1/Ex5_10_u1.sce
@@ -2,13 +2,13 @@
clc();
clear;
//To calculate the lattice constant
-h=6.63*10^-34 //Plancks Constant
-m=1.804*10^-27
-KB=1.38*10^-23
-T=300
-lamda=h/sqrt(3*m*KB*T) //units in mts
+h=6.63*10^-34 //units in m^2 kg s^-1
+m=1.804*10^-27 //units in Kgs
+KB=1.38*10^-23 //units in m^2 kg s^-2 K^-1
+T=300 //units in K
+lamda=h/sqrt(3*m*KB*T) //units in meters
n=2
-a=(sqrt(3)*lamda)/2 //units in mts
+a=(sqrt(3)*lamda)/2 //units in meters
printf("Lattice constant a=");
disp(a);
-printf("mts")
+printf("meters")
diff --git a/3411/CH5/EX5.10.u1/Ex5_10_u1.txt b/3411/CH5/EX5.10.u1/Ex5_10_u1.txt
index dc254f563..5d5133ca9 100644
--- a/3411/CH5/EX5.10.u1/Ex5_10_u1.txt
+++ b/3411/CH5/EX5.10.u1/Ex5_10_u1.txt
@@ -1,3 +1,3 @@
Lattice constant a=
1.213D-10
-mts \ No newline at end of file
+meters \ No newline at end of file
diff --git a/3411/CH5/EX5.11.u1/Ex5_11_u1.sce b/3411/CH5/EX5.11.u1/Ex5_11_u1.sce
index e98d92ef0..264e905ef 100644
--- a/3411/CH5/EX5.11.u1/Ex5_11_u1.sce
+++ b/3411/CH5/EX5.11.u1/Ex5_11_u1.sce
@@ -2,6 +2,78 @@
clc();
clear;
//To determine the unitcell and its dimensions
+//Experimental data
+//We have relation sin^(theta)=(lamda/2*a)^2 and (h^2+k^2+l^2)=j*(lamda/2*a)^2
+theta21=12.1 //units in degrees
+theta22=17.1 //units in degrees
+theta23=21 //units in degrees
+theta24=24.3 //units in degrees
+theta25=27.2 //units in degrees
+theta26=29.9 //units in degrees
+theta28=34.7 //units in degrees
+theta29=36.9 //units in degrees
+theta210=38.9 //units in degrees
+theta211=40.9 //units in degrees
+theta212=42.8 //units in degrees
+theta1=theta21/2 //units in degrees
+theta2=theta22/2 //units in degrees
+theta3=theta23/2 //units in degrees
+theta4=theta24/2 //units in degrees
+theta5=theta25/2 //units in degrees
+theta6=theta26/2 //units in degrees
+theta8=theta28/2 //units in degrees
+theta9=theta29/2 //units in degrees
+theta10=theta210/2 //units in degrees
+theta11=theta211/2 //units in degrees
+theta12=theta212/2 //units in degrees
+//sin^2(theta) values
+sin1=(sin(theta1*%pi/180))^2
+sin2=(sin(theta2*%pi/180))^2
+sin3=(sin(theta3*%pi/180))^2
+sin4=(sin(theta4*%pi/180))^2
+sin5=(sin(theta5*%pi/180))^2
+sin6=(sin(theta6*%pi/180))^2
+sin8=(sin(theta8*%pi/180))^2
+sin9=(sin(theta9*%pi/180))^2
+sin10=(sin(theta10*%pi/180))^2
+sin11=(sin(theta11*%pi/180))^2
+sin12=(sin(theta12*%pi/180))^2
+//sin^2(theta)/0.0111 value
+temp1=sin1/sin1
+temp2=sin2/sin1
+temp3=sin3/sin1
+temp4=sin4/sin1
+temp5=sin5/sin1
+temp6=sin6/sin1
+temp8=sin8/sin1
+temp9=sin9/sin1
+temp10=sin10/sin1
+temp11=sin11/sin1
+temp12=sin12/sin1
+//(h,k,l) values are determined such that the sum h^2+k^2+l^2=temp value in that manner hence we have to select the (h,k,l) values
+//(h,k,l) values
+hkl1=100 //As h^2+k^2+l^2=1
+hkl2=110 //As h^2+k^2+l^2=2
+hkl3=111 //As h^2+k^2+l^2=3
+hkl4=200 //As h^2+k^2+l^2=4
+hkl5=210 //As h^2+k^2+l^2=5
+hkl6=211 //As h^2+k^2+l^2=6
+hkl8=220 //As h^2+k^2+l^2=8
+hkl9=300 //As h^2+k^2+l^2=9
+hkl10=310 //As h^2+k^2+l^2=10
+hkl11=311 //As h^2+k^2+l^2=11
+hkl12=222 //As h^2+k^2+l^2=12
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta21,hkl1,temp1 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta22,hkl2,temp2 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta23,hkl3,temp3 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta24,hkl4,temp4 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is%.1f\n",theta25,hkl5,temp5 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta26,hkl6,temp6 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta28,hkl8,temp8 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta29,hkl9,temp9 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta210,hkl10,temp10 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta211,hkl11,temp11 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %.1f\n",theta212,hkl12,temp12 )
lamda=71 //units in pm
-a=lamda/(2*sqrt(0.0111)) //units in pm
+a=lamda/(2*sqrt(sin1)) //units in pm
printf("The unitcell and its dimensions are %dpm",a)
diff --git a/3411/CH5/EX5.11.u1/Ex5_11_u1.txt b/3411/CH5/EX5.11.u1/Ex5_11_u1.txt
index d5e90c7e6..52c0e8f14 100644
--- a/3411/CH5/EX5.11.u1/Ex5_11_u1.txt
+++ b/3411/CH5/EX5.11.u1/Ex5_11_u1.txt
@@ -1 +1,12 @@
- The unitcell and its dimensions are 336pm \ No newline at end of file
+ unit cell Dimensions when 2*theta=12.1 is (100) where sin^2(theta)/0.0111 is 1
+unit cell Dimensions when 2*theta=17.1 is (110) where sin^2(theta)/0.0111 is 2.0
+unit cell Dimensions when 2*theta=21.0 is (111) where sin^2(theta)/0.0111 is 3.0
+unit cell Dimensions when 2*theta=24.3 is (200) where sin^2(theta)/0.0111 is 4.0
+unit cell Dimensions when 2*theta=27.2 is (210) where sin^2(theta)/0.0111 is5.0
+unit cell Dimensions when 2*theta=29.9 is (211) where sin^2(theta)/0.0111 is 6.0
+unit cell Dimensions when 2*theta=34.7 is (220) where sin^2(theta)/0.0111 is 8.0
+unit cell Dimensions when 2*theta=36.9 is (300) where sin^2(theta)/0.0111 is 9.0
+unit cell Dimensions when 2*theta=38.9 is (310) where sin^2(theta)/0.0111 is 10.0
+unit cell Dimensions when 2*theta=40.9 is (311) where sin^2(theta)/0.0111 is 11.0
+unit cell Dimensions when 2*theta=42.8 is (222) where sin^2(theta)/0.0111 is 12.0
+The unitcell and its dimensions are 336pm \ No newline at end of file
diff --git a/3411/CH5/EX5.12.u1/Ex5_12_u1.sce b/3411/CH5/EX5.12.u1/Ex5_12_u1.sce
index 35bb1c6fd..a07b9e963 100644
--- a/3411/CH5/EX5.12.u1/Ex5_12_u1.sce
+++ b/3411/CH5/EX5.12.u1/Ex5_12_u1.sce
@@ -2,13 +2,56 @@
clc();
clear;
//To determine the cubic structure of element and lattice constant and to identify element
-ratio=(sin(20*%pi/180)/sin(29*%pi/180))^2
+//Diffraction data
+theta21=40 //units in degrees
+theta22=58 //units in degrees
+theta23=73 //units in degrees
+theta24=86.8 //units in degrees
+theta25=100.4 //units in degrees
+theta26=114.7 //units in degrees
+theta1=theta21/2 //units in degrees
+theta2=theta22/2 //units in degrees
+theta3=theta23/2 //units in degrees
+theta4=theta24/2 //units in degrees
+theta5=theta25/2 //units in degrees
+theta6=theta26/2 //units in degrees
+//sin^2(theta) values
+sin1=(sin(theta1*%pi/180))^2
+sin2=(sin(theta2*%pi/180))^2
+sin3=(sin(theta3*%pi/180))^2
+sin4=(sin(theta4*%pi/180))^2
+sin5=(sin(theta5*%pi/180))^2
+sin6=(sin(theta6*%pi/180))^2
+//sin^2(theta)/0.111 value
+temp1=sin1/sin1
+temp2=sin2/sin1
+temp3=sin3/sin1
+temp4=sin4/sin1
+temp5=sin5/sin1
+temp6=sin6/sin1
+//(h,k,l) values are determined such that the sum h^2+k^2+l^2=temp value in that manner hence we have to select the (h,k,l) values
+//(h,k,l) values
+hkl1=100 //As h^2+k^2+l^2=1
+hkl2=110 //As h^2+k^2+l^2=2
+hkl3=111 //As h^2+k^2+l^2=3
+hkl4=200 //As h^2+k^2+l^2=4
+hkl5=210 //As h^2+k^2+l^2=5
+hkl6=211 //As h^2+k^2+l^2=6
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta21,hkl1,temp1 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta22,hkl2,temp2 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta23,hkl3,temp3 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta24,hkl4,temp4 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is%d\n",theta25,hkl5,temp5 )
+printf("unit cell Dimensions when 2*theta=%.1f is (%d) where sin^2(theta)/0.0111 is %d\n",theta26,hkl6,temp6 )
+
+ratio=sin1/sin2
printf("The ratio of sin(theta)^2 values for first and second angles is %.2f\n Hence the crystal structure is bcc\n",ratio)
lamda=0.154 //units in nm
-h=1
-k=1
+//As we have used ratio of angles of 2*theta=40 degrees and 58 degrees above we use h=1,k=1,l=0 and a^2=(lamda/2)*sqrt(sqrt(h^2+k^2+l^2)/sin^2(theta))
+h=2
+k=0
l=0
theta=20 //units in degrees
a=(lamda/2)*(sqrt(sqrt(h^2+k^2+l^2)/sin(theta*(%pi/180))^2)) //units in nm
printf("Lattice constant a=%.3fnm \n And the element is tungsten Since Tungsten has lattice constant of %.3fnm and crystallizes in bcc structure",a,a)
-//Given in textbook lattice constant as a=0.318nm but the correct answer is a=0.268nm
+//Given in textbook to find lattice constant h=1,k=1,l=1 but the correct answer is h=2,k=0,l=0
diff --git a/3411/CH5/EX5.12.u1/Ex5_12_u1.txt b/3411/CH5/EX5.12.u1/Ex5_12_u1.txt
index 59d8efafb..1e2945a29 100644
--- a/3411/CH5/EX5.12.u1/Ex5_12_u1.txt
+++ b/3411/CH5/EX5.12.u1/Ex5_12_u1.txt
@@ -1,4 +1,10 @@
+unit cell Dimensions when 2*theta=40.0 is (100) where sin^2(theta)/0.0111 is 1
+unit cell Dimensions when 2*theta=58.0 is (110) where sin^2(theta)/0.0111 is 2
+unit cell Dimensions when 2*theta=73.0 is (111) where sin^2(theta)/0.0111 is 3
+unit cell Dimensions when 2*theta=86.8 is (200) where sin^2(theta)/0.0111 is 4
+unit cell Dimensions when 2*theta=100.4 is (210) where sin^2(theta)/0.0111 is5
+unit cell Dimensions when 2*theta=114.7 is (211) where sin^2(theta)/0.0111 is 6
The ratio of sin(theta)^2 values for first and second angles is 0.50
Hence the crystal structure is bcc
-Lattice constant a=0.268nm
- And the element is tungsten Since Tungsten has lattice constant of 0.268nm and crystallizes in bcc structure \ No newline at end of file
+Lattice constant a=0.318nm
+ And the element is tungsten Since Tungsten has lattice constant of 0.318nm and crystallizes in bcc structure \ No newline at end of file
diff --git a/3411/CH5/EX5.13.u1/Ex5_13_u1.sce b/3411/CH5/EX5.13.u1/Ex5_13_u1.sce
index a3dc6d537..c94889846 100644
--- a/3411/CH5/EX5.13.u1/Ex5_13_u1.sce
+++ b/3411/CH5/EX5.13.u1/Ex5_13_u1.sce
@@ -2,34 +2,74 @@
clc();
clear;
//To determine the crystal structure and indices of plane and lattice parameter of the material
-theta1=20.7
-peak1=sin((theta1/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 1 is %d\n",peak1)
-theta2=28.72
-peak2=sin((theta2/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 2 is %.1f\n",peak2)
-theta3=35.36
-peak3=sin((theta3/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 3 is %.1f\n",peak3)
-theta4=41.07
-peak4=sin((theta4/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 4 is %.1f\n",peak4)
-theta5=46.19
-peak5=sin((theta5/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 5 is %.1f\n",peak5)
-theta6=50.9
-peak6=sin((theta6/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 6 is %.1f\n",peak6)
-theta7=55.28
-peak7=sin((theta7/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 7 is %.1f\n",peak7)
-theta8=59.42
-peak8=sin((theta8/2)*%pi/180)^2/0.0308
-printf("Indices of the plane corresponding to peak 8 is %.1f\n",peak8)
+theta21=20.7 //units in degrees
+theta22=28.72 //units in degrees
+theta23=35.36 //units in degrees
+theta24=41.07 //units in degrees
+theta25=46.19 //units in degrees
+theta26=50.90 //units in degrees
+theta28=55.28 //units in degrees
+theta29=59.4 //units in degrees
+
+theta1=theta21/2 //units in degrees
+theta2=theta22/2 //units in degrees
+theta3=theta23/2 //units in degrees
+theta4=theta24/2 //units in degrees
+theta5=theta25/2 //units in degrees
+theta6=theta26/2 //units in degrees
+theta8=theta28/2 //units in degrees
+theta9=theta29/2 //units in degrees
+//sin^2(theta) values
+sin1=(sin(theta1*%pi/180))^2
+sin2=(sin(theta2*%pi/180))^2
+sin3=(sin(theta3*%pi/180))^2
+sin4=(sin(theta4*%pi/180))^2
+sin5=(sin(theta5*%pi/180))^2
+sin6=(sin(theta6*%pi/180))^2
+sin8=(sin(theta8*%pi/180))^2
+sin9=(sin(theta9*%pi/180))^2
+//sin^2(theta)/0.0308 values
+temp1=sin1/sin1
+temp2=sin2/sin1
+temp3=sin3/sin1
+temp4=sin4/sin1
+temp5=sin5/sin1
+temp6=sin6/sin1
+temp8=sin8/sin1
+temp9=sin9/sin1
+
+h2k2l21=temp1*2
+
+h2k2l22=temp2*2
+h2k2l23=temp3*2
+h2k2l24=temp4*2
+h2k2l25=temp5*2
+h2k2l26=temp6*2
+h2k2l28=temp8*2
+h2k2l29=temp9*2
+//(h,k,l) values are determined such that the sum h^2+k^2+l^2=temp value in that manner hence we have to select the (h,k,l) values
+//(h,k,l) values
+hkl1=110 //As h^2+k^2+l^2=2
+hkl2=200 //As h^2+k^2+l^2=4
+hkl3=211 //As h^2+k^2+l^2=6
+hkl4=220 //As h^2+k^2+l^2=8
+hkl5=310 //As h^2+k^2+l^2=10
+hkl6=232 //As h^2+k^2+l^2=12
+hkl8=321 //As h^2+k^2+l^2=14
+hkl9=400 //As h^2+k^2+l^2=16
+
+printf("unit cell Dimensions for peak 1 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is %.2f\n",theta21,hkl1,ceil(h2k2l21) )
+printf("unit cell Dimensions for peak 2 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is %.2f\n",theta22,hkl2,ceil(h2k2l22) )
+printf("unit cell Dimensions for peak 3 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is%.2f\n",theta23,hkl3,ceil(h2k2l23))
+printf("unit cell Dimensions for peak 4 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is %.2f\n",theta24,hkl4,ceil(h2k2l24))
+printf("unit cell Dimensions for peak 5 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is%.2f\n",theta25,hkl5,ceil(h2k2l25))
+printf("unit cell Dimensions for peak 6 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is%.2f\n",theta26,hkl6,ceil(h2k2l26))
+printf("unit cell Dimensions for peak 7 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is %.2f\n",theta28,hkl8,ceil(h2k2l28))
+printf("unit cell Dimensions for peak 8 when 2*theta=%.1f is (%d) where sin^2(theta)/0.0308 is %.2f\n",theta29,hkl9,ceil(h2k2l29))
+
printf("The material corresonds to bcc structure\n")
+//Consider peak no 8 where theta=29.71
lamda=0.07107 //units in nm
-theta=29.71 //units in degrees
-d400=lamda/(2*sin(theta*(%pi/180))) //units in nm
-hkl=16
-a=d400*sqrt(hkl) //units in nm
+d400=lamda/(2*sin(theta9*(%pi/180))) //units in nm
+a=d400*sqrt(ceil(h2k2l29)) //units in nm
printf("Lattice parameter of the material a=%.4fnm",a)
diff --git a/3411/CH5/EX5.13.u1/Ex5_13_u1.txt b/3411/CH5/EX5.13.u1/Ex5_13_u1.txt
index ca39815c5..01e496b7c 100644
--- a/3411/CH5/EX5.13.u1/Ex5_13_u1.txt
+++ b/3411/CH5/EX5.13.u1/Ex5_13_u1.txt
@@ -1,10 +1,10 @@
-Indices of the plane corresponding to peak 1 is 1
-Indices of the plane corresponding to peak 2 is 2.0
-Indices of the plane corresponding to peak 3 is 3.0
-Indices of the plane corresponding to peak 4 is 4.0
-Indices of the plane corresponding to peak 5 is 5.0
-Indices of the plane corresponding to peak 6 is 6.0
-Indices of the plane corresponding to peak 7 is 7.0
-Indices of the plane corresponding to peak 8 is 8.0
+unit cell Dimensions for peak 1 when 2*theta=20.7 is (110) where sin^2(theta)/0.0308 is 2
+unit cell Dimensions for peak 2 when 2*theta=28.7 is (200) where sin^2(theta)/0.0308 is 4.0
+unit cell Dimensions for peak 3 when 2*theta=35.4 is (211) where sin^2(theta)/0.0308 is 6.0
+unit cell Dimensions for peak 4 when 2*theta=41.1 is (220) where sin^2(theta)/0.0308 is 8.0
+unit cell Dimensions for peak 5 when 2*theta=46.2 is (310) where sin^2(theta)/0.0308 is10.0
+unit cell Dimensions for peak 6 when 2*theta=50.9 is (232) where sin^2(theta)/0.0308 is 12.0
+unit cell Dimensions for peak 7 when 2*theta=55.3 is (321) where sin^2(theta)/0.0308 is 14.0
+unit cell Dimensions for peak 8 when 2*theta=59.4 is (400) where sin^2(theta)/0.0308 is 15.94
The material corresonds to bcc structure
-Lattice parameter of the material a=0.2868nm \ No newline at end of file
+Lattice parameter of the material a=0.2863nm \ No newline at end of file
diff --git a/3411/CH5/EX5.14.u1/Ex5_14_u1.sce b/3411/CH5/EX5.14.u1/Ex5_14_u1.sce
index 847088d47..5c247d9c7 100644
--- a/3411/CH5/EX5.14.u1/Ex5_14_u1.sce
+++ b/3411/CH5/EX5.14.u1/Ex5_14_u1.sce
@@ -9,9 +9,9 @@ l=1
d=a/sqrt(h^2+k^2+l^2) //units in nm
theta=28.5 //units in degrees
lamda=2*d*sin(theta*(%pi/180)) //units in nm
-h=6.63*10^-34 //Plancks Constant
-m=1.67*10^-27
-KB=1.38*10^-23
-lamda=lamda*10^-9 //units in mts
-T=h^2/(3*m*KB*lamda^2)
+h=6.63*10^-34 //units in m^2 kg s^-1
+m=1.67*10^-27 //units in Kgs
+KB=1.38*10^-23 //units in m^2 kg s^-2 K^-1
+lamda=lamda*10^-9 //units in meters
+T=h^2/(3*m*KB*lamda^2) //units in K
printf("The effective temprature of neutrons is T=%dK",T)
diff --git a/3411/CH5/EX5.15.u1/Ex5_15_u1.sce b/3411/CH5/EX5.15.u1/Ex5_15_u1.sce
index 4030c57de..c3326bb73 100644
--- a/3411/CH5/EX5.15.u1/Ex5_15_u1.sce
+++ b/3411/CH5/EX5.15.u1/Ex5_15_u1.sce
@@ -2,10 +2,10 @@
clc();
clear;
//To calculate the Braggs angle
-h=6.63*10^-34 //Plancks Constant
-m=9.1*10^-31
-e=1.6*10^-19
-v=80
+h=6.63*10^-34 //units in m^2 kg s^-1
+m=9.1*10^-31 //units in Kgs
+e=1.6*10^-19 //units in coulombs
+v=80 //units in volts
lamda=h/sqrt(2*m*e*v) //units in mts
lamda=lamda*10^9 //units in nm
a=0.35 //units in nm
diff --git a/3411/CH5/EX5.16.u1/Ex5_16_u1.sce b/3411/CH5/EX5.16.u1/Ex5_16_u1.sce
index 98483f50c..d0535a5fd 100644
--- a/3411/CH5/EX5.16.u1/Ex5_16_u1.sce
+++ b/3411/CH5/EX5.16.u1/Ex5_16_u1.sce
@@ -2,8 +2,8 @@
clc();
clear;
//To calculate the difference between the samples
-d=0.2552
-a=d*sqrt(2)
+d=0.2552 //units in nm
+a=d*sqrt(2) //units in nm
lamda=0.152 //units in nm
theta=21 //units in degrees
//For sample A
@@ -11,14 +11,16 @@ d111=lamda/(2*sin(theta*%pi/180)) //units in nm
h=1
k=1
l=1
-a=d111*sqrt(h^2+k^2+l^2) //units in nm
-printf("For sample A a=%.4f nm",a)
+a1=d111*sqrt(h^2+k^2+l^2) //units in nm
+printf("For sample A a=%.4f nm",a1)
//For sample B
-theta=21.38
+theta=21.38 //units in degrees
d111=lamda/(2*sin(theta*%pi/180)) //units in nm
h=1
k=1
l=1
-a=d111*sqrt(h^2+k^2+l^2) //units in nm
-printf("\nFor sample B a=%.4f nm",a)
-printf("\n Sample B is pure high purity copper as lattice parameter of A is 1.75percent greater than that of pure copper")
+a2=d111*sqrt(h^2+k^2+l^2) //units in nm
+change=((a1-a2)/a2)*100
+printf("\nFor sample B a=%.4f nm",a2)
+printf("\n Sample B is pure high purity copper as lattice parameter of A is %.2f percent greater than that of pure copper",change)
+//Given in text book change in lattice parameter is 1.75% greater but it is 1.73%
diff --git a/3411/CH5/EX5.16.u1/Ex5_16_u1.txt b/3411/CH5/EX5.16.u1/Ex5_16_u1.txt
index 47695e7f2..d99dff05c 100644
--- a/3411/CH5/EX5.16.u1/Ex5_16_u1.txt
+++ b/3411/CH5/EX5.16.u1/Ex5_16_u1.txt
@@ -1,3 +1,3 @@
- For sample A a=0.3673 nm
+For sample A a=0.3673 nm
For sample B a=0.3611 nm
- Sample B is pure high purity copper as lattice parameter of A is 1.75percent greater than that of pure copper \ No newline at end of file
+ Sample B is pure high purity copper as lattice parameter of A is 1.73 percent greater than that of pure copper \ No newline at end of file
diff --git a/3411/CH5/EX5.9.u1/Ex5_9_u1.sce b/3411/CH5/EX5.9.u1/Ex5_9_u1.sce
index 7c8c82e26..b200c2b79 100644
--- a/3411/CH5/EX5.9.u1/Ex5_9_u1.sce
+++ b/3411/CH5/EX5.9.u1/Ex5_9_u1.sce
@@ -2,14 +2,14 @@
clc();
clear;
//To determine the interplanar spacing
-h=6.63*10^-34 //Plancks Constant
-m=9.1*10^-31
-e=1.6*10^-19
-v=844
-lamda=h/sqrt(2*m*e*v) //units in mts
+h=6.63*10^-34 //units in m^2 kg s^-1
+m=9.1*10^-31 //units in Kgs
+e=1.6*10^-19 //units in coulombs
+v=844 //units in Volts
+lamda=h/sqrt(2*m*e*v) //units in meters
n=1
theta=58 //units in degrees
-d=(n*lamda)/(2*sin(theta*(%pi/180))) //units in mts
+d=(n*lamda)/(2*sin(theta*(%pi/180))) //units in meters
printf("The interplanar spacing d=")
disp(d)
-printf("mts")
+printf("meters")
diff --git a/3411/CH5/EX5.9.u1/Ex5_9_u1.txt b/3411/CH5/EX5.9.u1/Ex5_9_u1.txt
index 05ff349b8..5a0c1075e 100644
--- a/3411/CH5/EX5.9.u1/Ex5_9_u1.txt
+++ b/3411/CH5/EX5.9.u1/Ex5_9_u1.txt
@@ -1,3 +1,3 @@
- The interplanar spacing d=
+The interplanar spacing d=
2.493D-11
-mts \ No newline at end of file
+meters \ No newline at end of file
generated by cgit (git 2.34.1) at 2024-12-20 11:49:08 +0000