summaryrefslogtreecommitdiff
path: root/296/CH1
diff options
context:
space:
mode:
Diffstat (limited to '296/CH1')
-rwxr-xr-x296/CH1/EX1.1/eg1_1.sce10
-rwxr-xr-x296/CH1/EX1.2/eg1_2.sce10
-rwxr-xr-x296/CH1/EX1.3/eg1_3.sce8
-rwxr-xr-x296/CH1/EX1.4/eg1_4.sce16
4 files changed, 44 insertions, 0 deletions
diff --git a/296/CH1/EX1.1/eg1_1.sce b/296/CH1/EX1.1/eg1_1.sce
new file mode 100755
index 000000000..947def86d
--- /dev/null
+++ b/296/CH1/EX1.1/eg1_1.sce
@@ -0,0 +1,10 @@
+a = 5; // lattice constant
+b = 0.5*sqrt(a^2 + a^2 +a^2); //separation between nearest atoms
+r = 0.5*b; //radius of each atom
+V = 4/3*%pi*r*r*r; //Volume of each atom
+n = 1+8*0.125; //number of atoms per cube
+pf = V*n/(a*a*a); //packing fraction
+disp("1)")
+disp(a,"lattice constant (in armstrong)=")
+disp(r,"radius of atoms (in armstrong) =")
+disp(pf,"maximum packing fraction =") \ No newline at end of file
diff --git a/296/CH1/EX1.2/eg1_2.sce b/296/CH1/EX1.2/eg1_2.sce
new file mode 100755
index 000000000..a879d15c8
--- /dev/null
+++ b/296/CH1/EX1.2/eg1_2.sce
@@ -0,0 +1,10 @@
+disp("The plane illustrated in Fig. 1-5 has intercepts at 2a, 4b and lc along the three crystal axes. Taking the reciprocals of these intercepts, we get 1/4, 1/2,and 1. These three fractions have the same relationship to each other as the integers 2,1, and 4 (obtained by multiplying each fraction by 4).Thus the plane can be referred to as a (214) plane. The only exception is if the intercept is a fraction of the lattice constant a. In that case, we do not reduce it to the lowest set of integers.")
+x = 0:0.05:2;
+y = 0:0.1:4;
+a=2;
+b=4;
+c=1;
+deff('[z]=fs(x,y)','z = (1-(0.5*x)-(0.25*y))');
+//z = (1-(0.5*x)-(0.25*y));
+fplot3d(x,y,fs);
+
diff --git a/296/CH1/EX1.3/eg1_3.sce b/296/CH1/EX1.3/eg1_3.sce
new file mode 100755
index 000000000..ad956ff42
--- /dev/null
+++ b/296/CH1/EX1.3/eg1_3.sce
@@ -0,0 +1,8 @@
+a = 5.43*10^-8; // lattice constant
+d = (1+4*0.25)/(a*a); //areal density on (100) plane
+n = 2*(6*0.5+8*0.125); //number of atoms per cube
+V = n/(a*a*a); //volume density
+disp("3)")
+disp(a,"lattice constant (in cm)=")
+disp(d,"areal density on (100) plane (in per cm square) =")
+disp(V,"volume density (in per cm cube) =") \ No newline at end of file
diff --git a/296/CH1/EX1.4/eg1_4.sce b/296/CH1/EX1.4/eg1_4.sce
new file mode 100755
index 000000000..19a3d9d40
--- /dev/null
+++ b/296/CH1/EX1.4/eg1_4.sce
@@ -0,0 +1,16 @@
+n = 10^16; // desired density of P atoms
+k = 0.35;
+l = 5000; //initial load of Si in grams
+w =31;//atomic weight of P
+d = 2.33; //density of Si
+i = n/k; //initial concentration of P in melt, assuming C(S)=kC(L)
+V = l/d; //volume of Si
+N = i*V; //number of P atoms
+W = N*w/(6.02*10^23)
+disp("4.a)")
+disp(n,"desired density of P atoms (per cubic centimeter)=")
+disp(i,"initial concentration of P in melt (in per cubic cm )=")
+disp("4.b)")
+disp(V,"Volume of Si (in cubic cm) =")
+disp(N,"number of P atoms =")
+disp(W,"weight of phosphorus to be added(in grams) =") \ No newline at end of file