updated the code

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