summaryrefslogtreecommitdiff
path: root/3411/CH10/EX2.8.u2/Ex2_8_u2.sce
blob: bf0d8d95e8e5c53f631ca5ff7a861c1d31cfee4f (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
//Example 2_8_u2
clc();
clear;
//TO Find the energy at which probability of occupancy is 0.9 and density of states and the population density
fe=0.9
k=(1/fe)-1
logk=log(k)
kb=8.62*10^-5
t=1000      //units in K
E=logk*kb*t        //units in eV
ef=7.06      //units in eV
energy=E+ef      //units in eV
printf("The energy of this state is E=%.2f eV",energy)
n=9.1*10^-31
EE=energy           //units in eV
h=4.14*10^-15
ZE=(8*sqrt(2)*%pi*n^1.5*sqrt(EE))/h^3        //units in met^-3 (eV)^-1
ZE=ZE*1.56*10^28                           //units met^-3 (eV)^-1
printf("\nThe density of the states for this energy is Z(E)=")
disp(ZE)
printf("met^-3 (eV)^-1")
ne=ZE*fe   //units in met^-3 (eV)^-1
printf("\nThe population density for this energy is N(E)=")
disp(ne)
printf("met^-3 (eV)^-1")