From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 2258/CH2/EX2.1/2_1.sce | 9 +++++++++ 2258/CH2/EX2.10/2_10.sce | 15 +++++++++++++++ 2258/CH2/EX2.11/2_11.sce | 14 ++++++++++++++ 2258/CH2/EX2.12/2_12.sce | 10 ++++++++++ 2258/CH2/EX2.13/2_13.sce | 15 +++++++++++++++ 2258/CH2/EX2.14/2_14.sce | 17 +++++++++++++++++ 2258/CH2/EX2.15/2_15.sce | 17 +++++++++++++++++ 2258/CH2/EX2.16/2_16.sce | 12 ++++++++++++ 2258/CH2/EX2.17/2_17.sce | 11 +++++++++++ 2258/CH2/EX2.2/2_2.sce | 9 +++++++++ 2258/CH2/EX2.3/2_3.sce | 17 +++++++++++++++++ 2258/CH2/EX2.4/2_4.sce | 18 ++++++++++++++++++ 2258/CH2/EX2.5/2_5.sce | 13 +++++++++++++ 2258/CH2/EX2.6/2_6.sce | 15 +++++++++++++++ 2258/CH2/EX2.7/2_7.sce | 10 ++++++++++ 2258/CH2/EX2.8/2_8.sce | 17 +++++++++++++++++ 2258/CH2/EX2.9/2_9.sce | 10 ++++++++++ 17 files changed, 229 insertions(+) create mode 100755 2258/CH2/EX2.1/2_1.sce create mode 100755 2258/CH2/EX2.10/2_10.sce create mode 100755 2258/CH2/EX2.11/2_11.sce create mode 100755 2258/CH2/EX2.12/2_12.sce create mode 100755 2258/CH2/EX2.13/2_13.sce create mode 100755 2258/CH2/EX2.14/2_14.sce create mode 100755 2258/CH2/EX2.15/2_15.sce create mode 100755 2258/CH2/EX2.16/2_16.sce create mode 100755 2258/CH2/EX2.17/2_17.sce create mode 100755 2258/CH2/EX2.2/2_2.sce create mode 100755 2258/CH2/EX2.3/2_3.sce create mode 100755 2258/CH2/EX2.4/2_4.sce create mode 100755 2258/CH2/EX2.5/2_5.sce create mode 100755 2258/CH2/EX2.6/2_6.sce create mode 100755 2258/CH2/EX2.7/2_7.sce create mode 100755 2258/CH2/EX2.8/2_8.sce create mode 100755 2258/CH2/EX2.9/2_9.sce (limited to '2258/CH2') 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 -- cgit