initial commit / add all books

17 files changed, 259 insertions, 0 deletions

diff --git a/3557/CH3/EX3.11/Ex3_11.sce b/3557/CH3/EX3.11/Ex3_11.sce
new file mode 100644
index 000000000..4a1d91191
--- /dev/null
+++ b/3557/CH3/EX3.11/Ex3_11.sce
@@ -0,0 +1,14 @@
+//Example 3.11//

+u=1;

+u1=1;

+v=1;

+v1=1;

+w=0;

+w1=1;

+a=(u*u1)+(v*v1)+(w*w1)

+//mprintf("a = %i",a)

+b=(sqrt((u^2)+(v^2)+(w^2)))*(sqrt((u1^2)+(v1^2)+(w1^2)))

+//mprintf("b = %i",b)

+c=acosd(a/b)

+mprintf("c = %f degree ",c)

+

diff --git a/3557/CH3/EX3.12/Ex3_12.sce b/3557/CH3/EX3.12/Ex3_12.sce
new file mode 100644
index 000000000..c4ec37ffa
--- /dev/null
+++ b/3557/CH3/EX3.12/Ex3_12.sce
@@ -0,0 +1,3 @@
+//Example 3.12//

+//As the  problem  is in the statement in the book

+mprintf("As the  problem  is in the statement in the book it cannot be solved using scilab")

diff --git a/3557/CH3/EX3.13/Ex3_13.sce b/3557/CH3/EX3.13/Ex3_13.sce
new file mode 100644
index 000000000..3e33591dc
--- /dev/null
+++ b/3557/CH3/EX3.13/Ex3_13.sce
@@ -0,0 +1,3 @@
+//Example 3.13//

+//As the  problem  is in the statement in the book

+mprintf("As the  problem  cannot be solved using scilab This problem is same as sample problem 3.10 ")

diff --git a/3557/CH3/EX3.14/Ex3_14.sce b/3557/CH3/EX3.14/Ex3_14.sce
new file mode 100644
index 000000000..3ed84b091
--- /dev/null
+++ b/3557/CH3/EX3.14/Ex3_14.sce
@@ -0,0 +1,15 @@
+//Example 3.14 (a)//
+a=2;//given
+rw=0.137;//nm // atomic radius of Tungsten
+r=a*rw
+mprintf("r = %f nm",r)
+r1=1/(r) //Taking inverse of r
+mprintf("\nr1 = %f atoms/nm",r1)
+//Example 3.14 (b)
+b=0.143;// atomic radius of Aluminium
+a1=(4*b)/(sqrt(2))  //Face centered cubic
+mprintf("\n a1 = %f nm",a1)
+r2=sqrt(3)*a1; //body diagonal length
+mprintf("\n r2 = %f nm",r2)
+r3=1/(r2); //linear density     //Taking inverse of r2 i.e body diagonal length
+mprintf("\n r3 = %f atoms/nm",r3)
diff --git a/3557/CH3/EX3.15/Ex3_15.sce b/3557/CH3/EX3.15/Ex3_15.sce
new file mode 100644
index 000000000..767b9065f
--- /dev/null
+++ b/3557/CH3/EX3.15/Ex3_15.sce
@@ -0,0 +1,32 @@
+//Example 3.15 (a)//

+rW=0.137;//nm //atomic radius of tungsten (From appendix 2)

+a=(4*rW)/(sqrt(3))//Body centered cubic

+mprintf("a = %f nm",a)

+l=sqrt(2)*a; // face diagonal length

+mprintf("\n l = %f nm",l)

+

+//The area of the  (111) plane within yhe unit cell

+c=sqrt(3);//given

+d=2;//given

+h=(c/d)*l

+//mprintf("h = %f ",h)

+A=(1/2)*l*h

+mprintf("\nA = %f nm^2",A)

+c1=3;//atoms

+d1=1/6;//atoms

+ad=(c1*d1)/A

+mprintf("\nad = %f atoms/nm^2",ad)

+

+//(b)

+// Following the calculations of sample problem 3.14b we find that the length of the body diagonal is

+b=0.143;// atomic radius of Aluminium

+a1=(4*b)/(sqrt(2))  //Face centered cubic

+//mprintf("\n a1 = %f nm",a1)

+l1=sqrt(2)*a1;

+mprintf("\nl1 = %f nm",l1)

+//the area of the (111) plane within the unit cell is

+A1=(1/2)*l1*(c/d)*l1

+mprintf("\nA1 = %f nm^2",A1)

+e1=(1/2);//atoms

+ad2=((c1*d1)+(c1*e1))/A1

+mprintf("\nad2 %f atoms/nm^2",ad2)

diff --git a/3557/CH3/EX3.16/Ex3_16.sce b/3557/CH3/EX3.16/Ex3_16.sce
new file mode 100644
index 000000000..cc9a48863
--- /dev/null
+++ b/3557/CH3/EX3.16/Ex3_16.sce
@@ -0,0 +1,15 @@
+//Example 3.16//

+// Following the calculations of sample problem 3.3 we find that the length of the body diagonal is

+rMg=0.078;//nm // Ionic radius of Magnesium (From Appendix 2)

+rO=0.132;//nm // Ionic radius of Oxygen (From Appendix 2)

+a=2*rMg+2*rO

+//mprintf("a = %f nm",a)

+l=sqrt(3)*a

+mprintf("l = %f nm",l)

+c=1;// Mg^2+

+i=c/l//nm

+mprintf("\n i = %f Mg^2+/nm",i)

+//similarly

+i2=c/l

+mprintf("\n i2 = %f O2-/nm",i2)

+mprintf("\n(1.37Mg2+ + 1.37O2-)/nm")

diff --git a/3557/CH3/EX3.17/Ex3_17.sce b/3557/CH3/EX3.17/Ex3_17.sce
new file mode 100644
index 000000000..01c39c5d6
--- /dev/null
+++ b/3557/CH3/EX3.17/Ex3_17.sce
@@ -0,0 +1,19 @@
+//Example 3.17//

+

+// Following the calculations of sample problem 3.3 we find that the length of the body diagonal is

+rMg=0.078;//nm // Ionic radius of Magnesium (From Appendix 2)

+rO=0.132;//nm // Ionic radius of Oxygen (From Appendix 2)

+a=2*rMg+2*rO

+//mprintf("a = %f nm",a)

+l=sqrt(2)*a

+mprintf("l = %f nm",l)

+d=(sqrt(3)*l)/2 //height

+//mprintf("\nd = %f",d)

+A=(1/2)*l*d //planar area

+mprintf("\n A = %f nm^-2",A)

+c=2;//ions

+id=c/A; //ionic density for Mg2+

+mprintf("\n id = %f nm^-2  (ionic density for Mg2+)",id)

+id1=c/A;//ionic density for O2-

+mprintf("\n id1 = %f nm^-2  (ionic density for O2-)",id1)

+mprintf("\n 13.1(Mg^2+ or O^2-)/nm^2")

diff --git a/3557/CH3/EX3.18/Ex3_18.sce b/3557/CH3/EX3.18/Ex3_18.sce
new file mode 100644
index 000000000..84c99f241
--- /dev/null
+++ b/3557/CH3/EX3.18/Ex3_18.sce
@@ -0,0 +1,9 @@
+//Example 3.18//

+rsi=0.117;//nm //atomic radius of silicon (From appendix 2)

+a=8;//given // (a is obtain by cross multiplication)

+l= a*rsi//nm //body diagonal length

+mprintf("l = %f nm",l)

+//the linear density

+b=2;//atoms //From the figure 3.23 there are two atoms per lattice point therefore we choose value 2 atoms in linear density

+ld=b/l

+mprintf("\n ld = %f atoms/nm",ld)

diff --git a/3557/CH3/EX3.19/Ex3_19.sce b/3557/CH3/EX3.19/Ex3_19.sce
new file mode 100644
index 000000000..077d69580
--- /dev/null
+++ b/3557/CH3/EX3.19/Ex3_19.sce
@@ -0,0 +1,12 @@
+//Example 3.19//

+

+e=0.117;//nm //atomic radius of silicon (From Appendix 2)

+a=(8/sqrt(3))*e

+mprintf("a= %f nm",a)

+s=sqrt(2)*a

+mprintf("\n s= %f nm",s)

+i=2;//atoms //From the figure 3.23 there are two atoms per lattice point therefore we choose value 2 atoms in planar density

+A=(1/2)*s*(sqrt(3)/2)*s //area of traingle

+mprintf("\n A = %f nm^2",A)

+p=i/A;//planar density

+mprintf("\n p = %f atoms/nm^2",p)

diff --git a/3557/CH3/EX3.2/Ex3_2.sce b/3557/CH3/EX3.2/Ex3_2.sce
new file mode 100644
index 000000000..9ce3b10ee
--- /dev/null
+++ b/3557/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,11 @@
+//Example 3.2//

+rCu=0.128;//nm //atomic radius copper (From appendix 2)

+a=(4/sqrt(2))*rCu

+mprintf("a = %f nm",a)

+//The density of the unit cells is

+a1=4;// atoms

+b1=63.55;//gram //atomic mass of copper

+c1=0.6023*10^24;//atoms// Avogardo's number

+d=10^7;//nm/cm

+p=(a1/a^3)*(b1/c1)*d^3

+mprintf("\n p = %f g/cm^3",p)

diff --git a/3557/CH3/EX3.20/Ex3_20.sce b/3557/CH3/EX3.20/Ex3_20.sce
new file mode 100644
index 000000000..ead2cc6b1
--- /dev/null
+++ b/3557/CH3/EX3.20/Ex3_20.sce
@@ -0,0 +1,29 @@
+//Example 3.20//

+c=1;//centimeter // opposite side of a triangle

+e=3;//centimeter // adjacent side of a triangle

+a=atand(c/e)// (As tan = oppposite side/adjacent side)

+mprintf("a = %f degree ",a)

+a1=180;//degree

+b1=2;//given

+theta= (a1-a)/b1

+mprintf("\ntheta = %f degree",theta)

+//Braggs law

+rMg=0.078;//nm // Ionic radius of Magnesium (From Appendix 2)

+rO=0.132;//nm // Ionic radius of Oxygen (From Appendix 2)

+a2=2*rMg+2*rO

+//mprintf("a2 = %f nm",a2)

+h=1; //spacing between adjacent plane

+k=1;//spacing between adjacent plane

+l=1;//spacing between adjacent plane

+d=(a2)/sqrt(h^2+k^2+l^2)

+mprintf("\n d = %f nm",d)

+//substituting to obtain lamda for n=1

+//for n=1

+l1=b1*d*sind(theta)

+mprintf("\n l1= %f nm",l1)

+//for n=2

+l2=(b1*d*sind(theta))/b1

+mprintf("\n l2 = %f nm",l2)

+//for n=3

+l3=(b1*d*sind(theta))/e;

+mprintf("\n l3 = %f nm",l3)

diff --git a/3557/CH3/EX3.21/Ex3_21.sce b/3557/CH3/EX3.21/Ex3_21.sce
new file mode 100644
index 000000000..433d29ee1
--- /dev/null
+++ b/3557/CH3/EX3.21/Ex3_21.sce
@@ -0,0 +1,31 @@
+//Example 3.21//

+a= 0.404;//nm //lattice parameter

+a1=1;//given

+b1=1;//given

+c1=1;//given

+b=sqrt(a1+b1+c1)

+d111=a/b

+mprintf("d111 = %f = nm",d111)

+a2=2;//given

+b2=0;//given

+c2=0;//given

+d200=a/sqrt(a2^2+b2+c2);

+mprintf("\n d200 = %f nm",d200)

+a3=2;//given

+b3=2;//given

+c3=0;//given

+d220=a/sqrt(a3^2+b3^2+c3);

+mprintf("\n d220 = %f nm",d220)

+l=0.1542;//nm// from the figure 3.39

+thetha111=asind(l/(a2*d111))

+mprintf("\nthetha111 = %f degree",thetha111)

+t111=a2*thetha111

+mprintf("\nt111 = %f degree",t111)

+thetha200=asind(l/(a2*d200))

+mprintf("\nthetha200 = %f degree",thetha200)

+t200=a2*thetha200

+mprintf("\nt200 = %f degree",t200)

+thetha220=asind(l/(a2*d220))

+mprintf("\nthetha220 = %f degree",thetha220)

+t220=a2*thetha220

+mprintf("\nt220 = %f degree",t220)

diff --git a/3557/CH3/EX3.3/Ex3_3.sce b/3557/CH3/EX3.3/Ex3_3.sce
new file mode 100644
index 000000000..61e2f9b38
--- /dev/null
+++ b/3557/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,13 @@
+//Example 3.3//

+rMg=0.078;//nm // Ionic radius of Magnesium (From Appendix 2)

+rO=0.132;//nm // Ionic radius of Oxygen (From Appendix 2)

+a=2*rMg+2*rO

+mprintf("a = %f nm",a)

+Vu=(a)^3;//nm

+mprintf("\nVu = %f nm^3",Vu)

+b=4;//by formula

+c=4/3;//By formula

+volume=((b*c)*%pi*(rMg)^3)+((b*c)*%pi*(rO)^3)

+mprintf("\nvolume = %f nm^3",volume)

+IPF=volume/Vu;

+mprintf("\nIPF = %f ",IPF)

diff --git a/3557/CH3/EX3.4/Ex3_4.sce b/3557/CH3/EX3.4/Ex3_4.sce
new file mode 100644
index 000000000..2b6a439df
--- /dev/null
+++ b/3557/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,9 @@
+//Example 3.4//

+a=24.31;//gram //atomic mass of magnesium

+b=16.00;//gram // atomic mass of oxygen

+c=0.6023*10^24;//Avogardo's number

+v=0.0741;//nm^3 //unit cell volume

+d=10^7;//nm/cm

+e=4;//Number of electrons

+p=((((e*a)+(e*b))/(c))/(v))*d^3

+mprintf("p = %f g/cm^3",p)

diff --git a/3557/CH3/EX3.5/Ex3_5.sce b/3557/CH3/EX3.5/Ex3_5.sce
new file mode 100644
index 000000000..99679293d
--- /dev/null
+++ b/3557/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,21 @@
+//Example 3.5//

+a=0.741;//nm //unit cell dimensions

+b=0.494;//nm //unit cell dimensions

+c=0.255;//nm //unit cell dimensions

+v=a*b*c

+mprintf("v = %f nm^3",v)

+a1=12.01;//gram //atomic mass of carbon

+b1=1.008;//gram // atomic mass of Hydogen

+c1=0.6023*10^24;//atoms // Avogardo's number

+d1=2;//Number of electrons

+e1=4;//Number of electrons

+m=((d1*a1)+(e1*b1))/c1

+mprintf("\nm = (%e n)g",m)

+//Therefore, the unit cell density is,

+d=10^7;//nm/cm

+p=(m/v)*d^3

+mprintf("\n p = %f g/cm^3 (As answer in the textbook is calculated  wrong)",p)

+//solving for n gives

+n=2

+//Aa a result, there are

+mprintf("\n4(=2n)C atoms + 8(=4n)H atoms per unit cell.")

diff --git a/3557/CH3/EX3.6/Ex3_6.sce b/3557/CH3/EX3.6/Ex3_6.sce
new file mode 100644
index 000000000..7110ff1e6
--- /dev/null
+++ b/3557/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,12 @@
+//Example 3.6//

+a=2;//body diagonal

+b=4;//body diagonal

+c=a*b ;//(using cross multiplication)

+//mprintf("c= %i ",c)

+d=sqrt(3);

+Vu=(c/d)^3

+mprintf("Vu = %f rSi^3",Vu)

+Va=c*(4/3)*%pi

+mprintf("\n Va = %f rSi^3",Va)

+APF=Va/Vu;

+mprintf("\nAPF = %f ",APF)

diff --git a/3557/CH3/EX3.7/Ex3_7.sce b/3557/CH3/EX3.7/Ex3_7.sce
new file mode 100644
index 000000000..e56452257
--- /dev/null
+++ b/3557/CH3/EX3.7/Ex3_7.sce
@@ -0,0 +1,11 @@
+//Example 3.7//

+a=98.5;// Unit cell volume

+b=0.117;//nm //nanometer //atomic radius of Silicon

+V=a*b^3

+mprintf("V = %f nm^3",V)

+a1=8;//atoms

+c=28.09;//gram //atomic mass of silicon

+d=0.6023*10^24;//atoms //Avogardo's number

+e=10^7;//nm/cm

+P=(a1/V)*(c/d)*(e^3)

+mprintf("\nP = %f g/cm^3",P)