initial commit / add all books

17 files changed, 229 insertions, 0 deletions

diff --git a/2258/CH2/EX2.1/2_1.sce b/2258/CH2/EX2.1/2_1.sce
new file mode 100755
index 000000000..182ca4144
--- /dev/null
+++ b/2258/CH2/EX2.1/2_1.sce
@@ -0,0 +1,9 @@
+clc();

+clear;

+// To calculate the Fermi function

+// given that E-Ef = kT

+// fermi function FE = 1/(1+exp((E-Ef)/kT)

+// therefore FE = 1/(1+exp(kT/kT));

+// FE = 1/(1+exp(1))

+FE=1/(1+exp(1));

+printf("fermi function is %f",FE);

diff --git a/2258/CH2/EX2.10/2_10.sce b/2258/CH2/EX2.10/2_10.sce
new file mode 100755
index 000000000..554327a23
--- /dev/null
+++ b/2258/CH2/EX2.10/2_10.sce
@@ -0,0 +1,15 @@
+clc();

+clear;

+// To calculate the number of states per unit volume

+m=9.1*10^-31;   //mass in kg

+h=6.626*10^-34;

+A=(8*m)^(3/2);

+B=%pi/(2*h^3);

+EfeV=3.10;    //fermi energy in eV

+Ef=EfeV*1.6*10^-19;    //fermi energy in J

+EFeV=EfeV+0.02;    //energy after interval in eV

+EF=EFeV*1.6*10^-19;    //energy after interval in J

+function Q=f(E),Q=A*B*sqrt(E),endfunction

+I=intg(Ef,EF,f)

+printf("number of energy states per unit volume is");

+disp(I);

diff --git a/2258/CH2/EX2.11/2_11.sce b/2258/CH2/EX2.11/2_11.sce
new file mode 100755
index 000000000..62df8945d
--- /dev/null
+++ b/2258/CH2/EX2.11/2_11.sce
@@ -0,0 +1,14 @@
+clc();

+clear;

+// To calculate the mean free path of electron

+T=300;   //temperature in K

+n=8.5*10^28;   //density per m^3

+rho=1.69*10^-8;   //resistivity in ohm/m^3

+me=9.11*10^-31;    //mass of electron in kg

+e=1.6*10^-19;    //charge in coulomb

+KB=1.38*10^-23;    //boltzmann constant in J/k

+lamda=sqrt(3*KB*me*T)/(n*(e^2)*rho);

+printf("mean free path of electron in m is");

+disp(lamda);

+

+//answer given in the book is wrong

diff --git a/2258/CH2/EX2.12/2_12.sce b/2258/CH2/EX2.12/2_12.sce
new file mode 100755
index 000000000..eee3a1345
--- /dev/null
+++ b/2258/CH2/EX2.12/2_12.sce
@@ -0,0 +1,10 @@
+clc();

+clear;

+// To calculate the relaxation time of conduction electrons

+rho=1.43*10^-8;    //resistivity in ohm-m

+n=6.5*10^28;   //electron/m^3

+m=9.11*10^-34;   //mass in kg

+e=1.6*10^-19;   //charge in coulomb

+tow=m/(n*(e^2)*rho);

+printf("relaxation time of conduction electrons in sec is");

+disp(tow);

diff --git a/2258/CH2/EX2.13/2_13.sce b/2258/CH2/EX2.13/2_13.sce
new file mode 100755
index 000000000..85222dc9c
--- /dev/null
+++ b/2258/CH2/EX2.13/2_13.sce
@@ -0,0 +1,15 @@
+clc();

+clear;

+// To calculate the mobility and average time of collision of electrons

+d=8.92*10^3;   //density in kg/m^3

+rho=1.73*10^-8;   //resistivity in ohm-m

+m=9.1*10^-31;    //mass in kg

+M=63.5;    //atomic weight

+e=1.6*10^-19;    //charge in coulomb

+A=6.02*10^26;    //avagadro number

+n=(d*A)/M;

+mew=1/(rho*n*e);

+tow=m/(n*(e^2)*rho);

+printf("mobility of electrons in Copper is %f m/Vs",mew);

+printf("average time of collision of electrons in copper in sec is");

+disp(tow);

diff --git a/2258/CH2/EX2.14/2_14.sce b/2258/CH2/EX2.14/2_14.sce
new file mode 100755
index 000000000..ebef1a637
--- /dev/null
+++ b/2258/CH2/EX2.14/2_14.sce
@@ -0,0 +1,17 @@
+clc();

+clear;

+// To calculate the order of magnitude of velocity of molecules

+MH=1.008*2*1.67*10^-27;    //mass in kg

+T=30;   //temperature in C

+T=T+273;    //temperature in K

+KB=1.38*10^-23;    //boltzmann constant in J/k

+KE=(3/2)*KB*T;   //kinetic energy in J

+KEeV=KE*6.24*10^18;   //kinetic energy in eV

+cbar=sqrt((3*KB*T)/MH);

+printf("average kinetic energy in J is");

+disp(KE);

+printf("average kinetic energy in eV is");

+disp(KEeV);

+printf("velocity of molecules is %f m/s",cbar);

+

+//answers for average kinetic energy in eV and velocity of electrons given in the book are wrong

diff --git a/2258/CH2/EX2.15/2_15.sce b/2258/CH2/EX2.15/2_15.sce
new file mode 100755
index 000000000..321554ccc
--- /dev/null
+++ b/2258/CH2/EX2.15/2_15.sce
@@ -0,0 +1,17 @@
+clc();

+clear;

+// To calculate the velocity of an electron and proton

+Ee=10;   //electron kinetic energy in eV

+EeeV=Ee*1.6*10^-19;   //electron kinetic energy in J

+Ep=10;   //proton kinetic energy in eV

+EpeV=Ep*1.6*10^-19;   //proton kinetic energy in J

+me=9.1*10^-31;   //mass of electron in kg

+mp=1.67*10^-27;   //mass of proton in kg

+cebar=sqrt((2*EeeV)/me);

+cpbar=sqrt((2*EpeV)/mp);

+printf("velocity of electron in m/s is");

+disp(cebar);

+printf("velocity of proton in m/s is");

+disp(cpbar);

+

+//answers given in the book are wrong

diff --git a/2258/CH2/EX2.16/2_16.sce b/2258/CH2/EX2.16/2_16.sce
new file mode 100755
index 000000000..da59c6275
--- /dev/null
+++ b/2258/CH2/EX2.16/2_16.sce
@@ -0,0 +1,12 @@
+clc();

+clear;

+// To calculate the drift velocity of free electrons

+A=10;   //area of cross section in mm^2

+A=A*10^-6;   //area of cross section in m^2

+i=100;   //current in amp

+n=8.5*10^28;    //number of electrons per mm^3

+e=1.6*10^-19;    //electron charge in coulumb

+vd=1/(n*A*e);

+printf("drift velocity is %f m/s",vd);

+

+//answer given in the book is wrong

diff --git a/2258/CH2/EX2.17/2_17.sce b/2258/CH2/EX2.17/2_17.sce
new file mode 100755
index 000000000..129ec9057
--- /dev/null
+++ b/2258/CH2/EX2.17/2_17.sce
@@ -0,0 +1,11 @@
+clc();

+clear;

+// To calculate the thermal conductivity of copper

+tow=3*10^-14;    //relaxation time in sec

+n=8*10^28;    //density of electrons per m^3

+KB=1.38*10^-23;   //boltzmann constant in J/k

+T=0;   //temperature in C

+T=T+273;   //temperature in K

+m=9.1*10^-31;   //mass of electron in kg

+sigma_T=((3*n*tow*(KB^2)*T)/(2*m));

+printf("thermal conductivity of copper is %f ohm-1",sigma_T);

diff --git a/2258/CH2/EX2.2/2_2.sce b/2258/CH2/EX2.2/2_2.sce
new file mode 100755
index 000000000..182ca4144
--- /dev/null
+++ b/2258/CH2/EX2.2/2_2.sce
@@ -0,0 +1,9 @@
+clc();

+clear;

+// To calculate the Fermi function

+// given that E-Ef = kT

+// fermi function FE = 1/(1+exp((E-Ef)/kT)

+// therefore FE = 1/(1+exp(kT/kT));

+// FE = 1/(1+exp(1))

+FE=1/(1+exp(1));

+printf("fermi function is %f",FE);

diff --git a/2258/CH2/EX2.3/2_3.sce b/2258/CH2/EX2.3/2_3.sce
new file mode 100755
index 000000000..51852afec
--- /dev/null
+++ b/2258/CH2/EX2.3/2_3.sce
@@ -0,0 +1,17 @@
+clc();

+clear;

+// To calculate the temperature

+FE=10/100;    //fermi function is 10%

+Ef=5.5;   //fermi energy of silver in eV

+k=1.38*10^-23;

+E=Ef+(Ef/100);

+//FE=1/(1+exp((E-Ef)/(k*T)))

+//therefore 1/FE = 1+exp((E-Ef)/(k*T))

+//therefore (1/FE)-1 = exp((E-Ef)/(k*T))

+//therefore log((1/FE)-1) = (E-Ef)/(k*T)

+//therefore T = (E-Ef)/(k*log((1/FE)-1))

+//let X=E-Ef;    

+X=E-Ef;     //energy in eV

+X=X*1.6*10^-19;    //energy in J

+T = (X/(k*log((1/FE)-1)));

+printf("temperature is %f K",T);

diff --git a/2258/CH2/EX2.4/2_4.sce b/2258/CH2/EX2.4/2_4.sce
new file mode 100755
index 000000000..9fecc5793
--- /dev/null
+++ b/2258/CH2/EX2.4/2_4.sce
@@ -0,0 +1,18 @@
+clc();

+clear;

+// To calculate the temperature

+//let X=E-Ef

+X=0.5;   //E-Ef=0.5 in eV

+X=X*1.6*10^-19;   //X in J

+FE=1/100;    //fermi function is 1% 

+k=1.38*10^-23;

+//FE=1/(1+exp(X/(k*T)))

+//therefore 1/FE = 1+exp(X/(k*T))

+//therefore (1/FE)-1 = exp(X/(k*T))

+//therefore log((1/FE)-1) = X/(k*T)

+//but log(x) = 2.303*log10(x)

+//therefore T = X/(k*log((1/FE)-1))

+//but log(x)=2.303*log10(x)

+//therefore T = X/(k*2.303*log10((1/FE)-1))

+T = X/(k*2.303*log10((1/FE)-1));

+printf("temperature is %f K",T);

diff --git a/2258/CH2/EX2.5/2_5.sce b/2258/CH2/EX2.5/2_5.sce
new file mode 100755
index 000000000..a579261a5
--- /dev/null
+++ b/2258/CH2/EX2.5/2_5.sce
@@ -0,0 +1,13 @@
+clc();

+clear;

+// To calculate the density and mobility of electrons in silver

+rho_s=10.5*10^3;   //density in kg/m^3

+NA=6.02*10^26;    //avagadro number per kmol

+MA=107.9;    

+n=(rho_s*NA)/MA;

+sigma=6.8*10^7;

+e=1.6*10^-19;   //charge in coulomb

+mew=sigma/(n*e);

+printf("density of electrons is");

+disp(n);

+printf("mobility of electrons in silver is %f m^2/Vs",mew);

diff --git a/2258/CH2/EX2.6/2_6.sce b/2258/CH2/EX2.6/2_6.sce
new file mode 100755
index 000000000..b807d2cee
--- /dev/null
+++ b/2258/CH2/EX2.6/2_6.sce
@@ -0,0 +1,15 @@
+clc();

+clear;

+// To calculate the mobility and average time of collision of electrons

+d=8.92*10^3;   //density in kg/m^3

+rho=1.73*10^-8;   //resistivity in ohm-m

+m=9.1*10^-31;    //mass in kg

+w=63.5;    //atomic weight

+e=1.6*10^-19;    //charge in coulomb

+A=6.02*10^26;    //avagadro number

+n=(d*A)/w;

+mew=1/(rho*n*e);

+tow=m/(n*(e^2)*rho);

+printf("mobility of electrons in Copper is %f m/Vs",mew);

+printf("average time of collision of electrons in copper in sec is");

+disp(tow);

diff --git a/2258/CH2/EX2.7/2_7.sce b/2258/CH2/EX2.7/2_7.sce
new file mode 100755
index 000000000..fd6ce9d9b
--- /dev/null
+++ b/2258/CH2/EX2.7/2_7.sce
@@ -0,0 +1,10 @@
+clc();

+clear;

+// To calculate the relaxation time of conduction electrons

+rho=1.54*10^-8;    //resistivity in ohm-m

+n=5.8*10^28;   //electron/m^3

+m=9.108*10^-31;   //mass in kg

+e=1.602*10^-19;   //charge in coulomb

+tow=m/(n*(e^2)*rho);

+printf("relaxation time of conduction electrons in sec is");

+disp(tow);

diff --git a/2258/CH2/EX2.8/2_8.sce b/2258/CH2/EX2.8/2_8.sce
new file mode 100755
index 000000000..51852afec
--- /dev/null
+++ b/2258/CH2/EX2.8/2_8.sce
@@ -0,0 +1,17 @@
+clc();

+clear;

+// To calculate the temperature

+FE=10/100;    //fermi function is 10%

+Ef=5.5;   //fermi energy of silver in eV

+k=1.38*10^-23;

+E=Ef+(Ef/100);

+//FE=1/(1+exp((E-Ef)/(k*T)))

+//therefore 1/FE = 1+exp((E-Ef)/(k*T))

+//therefore (1/FE)-1 = exp((E-Ef)/(k*T))

+//therefore log((1/FE)-1) = (E-Ef)/(k*T)

+//therefore T = (E-Ef)/(k*log((1/FE)-1))

+//let X=E-Ef;    

+X=E-Ef;     //energy in eV

+X=X*1.6*10^-19;    //energy in J

+T = (X/(k*log((1/FE)-1)));

+printf("temperature is %f K",T);

diff --git a/2258/CH2/EX2.9/2_9.sce b/2258/CH2/EX2.9/2_9.sce
new file mode 100755
index 000000000..d87a1dcc6
--- /dev/null
+++ b/2258/CH2/EX2.9/2_9.sce
@@ -0,0 +1,10 @@
+clc();

+clear;

+// To calculate the Fermi distribution function

+// given that E-Ef = kT

+// fermi function FpE = 1/(1+exp((E-Ef)/kT)

+// therefore FpE = 1/(1+exp(kT/kT));

+// FpE = 1/(1+exp(1))

+FpE=1/(1+exp(1));

+printf("fermi function is %f",FpE);

+//the presence of electron at that energy level is not certain