summaryrefslogtreecommitdiff
path: root/3401
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /3401
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '3401')
-rwxr-xr-x3401/CH1/EX1.1/Ex1_1.sce6
-rwxr-xr-x3401/CH1/EX1.3/Ex1_3.sce6
-rwxr-xr-x3401/CH10/EX10.1/Ex10_1.sce9
-rwxr-xr-x3401/CH10/EX10.10/Ex10_10.sce12
-rwxr-xr-x3401/CH10/EX10.12/Ex10_12.sce33
-rwxr-xr-x3401/CH10/EX10.2/Ex10_2.sce16
-rwxr-xr-x3401/CH10/EX10.3/Ex10_3.sce11
-rwxr-xr-x3401/CH10/EX10.4/Ex10_4.sce32
-rwxr-xr-x3401/CH10/EX10.5/Ex10_5.sce29
-rwxr-xr-x3401/CH10/EX10.7/Ex10_7.sce16
-rwxr-xr-x3401/CH10/EX10.8/Ex10_8.sce17
-rwxr-xr-x3401/CH10/EX10.9/Ex10_9.sce11
-rwxr-xr-x3401/CH11/EX11.1/Ex11_1.sce13
-rwxr-xr-x3401/CH11/EX11.10/Ex11_10.sce18
-rwxr-xr-x3401/CH11/EX11.11/Ex11_11.sce8
-rwxr-xr-x3401/CH11/EX11.2/Ex11_2.sce14
-rwxr-xr-x3401/CH11/EX11.3/Ex11_3.sce15
-rwxr-xr-x3401/CH11/EX11.4/Ex11_4.sce21
-rwxr-xr-x3401/CH11/EX11.5/Ex11_5.sce21
-rwxr-xr-x3401/CH11/EX11.6/Ex11_6.sce22
-rwxr-xr-x3401/CH11/EX11.7/Ex11_7.sce36
-rwxr-xr-x3401/CH11/EX11.8/Ex11_8.sce15
-rwxr-xr-x3401/CH11/EX11.9/Ex11_9.sce19
-rwxr-xr-x3401/CH2/EX2.1/Ex2_1.sce10
-rwxr-xr-x3401/CH2/EX2.2/Ex2_2.sce10
-rwxr-xr-x3401/CH2/EX2.3/Ex2_3.sce24
-rwxr-xr-x3401/CH2/EX2.4/Ex2_4.sce12
-rwxr-xr-x3401/CH2/EX2.5/Ex2_5.sce12
-rwxr-xr-x3401/CH3/EX3.1/Ex3_1.sce14
-rwxr-xr-x3401/CH3/EX3.2/Ex3_2.sce18
-rwxr-xr-x3401/CH3/EX3.3/Ex3_3.sce7
-rwxr-xr-x3401/CH3/EX3.4/Ex3_4.sce8
-rwxr-xr-x3401/CH3/EX3.5/Ex3_5.sce8
-rwxr-xr-x3401/CH3/EX3.6/Ex3_6.sce9
-rwxr-xr-x3401/CH3/EX3.7/Ex3_7.sce16
-rwxr-xr-x3401/CH4/EX4.1/Ex4_1.sce16
-rwxr-xr-x3401/CH4/EX4.10/Ex4_10.sce16
-rwxr-xr-x3401/CH4/EX4.11/Ex4_11.sce19
-rwxr-xr-x3401/CH4/EX4.12/Ex4_12.sce24
-rwxr-xr-x3401/CH4/EX4.13/Ex4_13.sce14
-rwxr-xr-x3401/CH4/EX4.14/Ex4_14.sce16
-rwxr-xr-x3401/CH4/EX4.2/Ex4_2.sce16
-rwxr-xr-x3401/CH4/EX4.3/Ex4_3.sce15
-rwxr-xr-x3401/CH4/EX4.4/Ex4_4.sce12
-rwxr-xr-x3401/CH4/EX4.5/Ex4_5.sce15
-rwxr-xr-x3401/CH4/EX4.6/Ex4_6.sce11
-rwxr-xr-x3401/CH4/EX4.7/Ex4_7.sce12
-rwxr-xr-x3401/CH4/EX4.9/Ex4_9.sce14
-rwxr-xr-x3401/CH5/EX5.1/Ex5_1.sce21
-rwxr-xr-x3401/CH5/EX5.2/Ex5_2.sce27
-rwxr-xr-x3401/CH5/EX5.3/Ex5_3.sce27
-rwxr-xr-x3401/CH5/EX5.4/Ex5_4.sce11
-rwxr-xr-x3401/CH5/EX5.5/Ex5_5.sce13
-rwxr-xr-x3401/CH5/EX5.6/Ex5_6.sce10
-rwxr-xr-x3401/CH5/EX5.7/Ex5_7.sce15
-rwxr-xr-x3401/CH6/EX6.10/Ex6_10.sce13
-rwxr-xr-x3401/CH6/EX6.5/Ex6_5.sce17
-rwxr-xr-x3401/CH6/EX6.6/Ex6_6.sce22
-rwxr-xr-x3401/CH6/EX6.8/Ex6_8.sce24
-rwxr-xr-x3401/CH7/EX7.1/Ex7_1.sce15
-rwxr-xr-x3401/CH7/EX7.2/Ex7_2.sce14
-rwxr-xr-x3401/CH7/EX7.3/Ex7_3.sce12
-rwxr-xr-x3401/CH7/EX7.4/Ex7_4.sce16
-rwxr-xr-x3401/CH7/EX7.5/Ex7_5.sce14
-rwxr-xr-x3401/CH7/EX7.6/Ex7_6.sce14
-rwxr-xr-x3401/CH8/EX8.1/Ex8_1.sce15
-rwxr-xr-x3401/CH8/EX8.2/Ex8_2.sce16
-rwxr-xr-x3401/CH8/EX8.3/Ex8_3.sce23
-rwxr-xr-x3401/CH8/EX8.4/Ex8_4.sce17
-rwxr-xr-x3401/CH8/EX8.5/Ex8_5.sce14
-rwxr-xr-x3401/CH8/EX8.6/Ex8_6.sce15
-rwxr-xr-x3401/CH8/EX8.7/ex8_7.sce18
-rwxr-xr-x3401/CH9/EX9.1/Ex9_1.sce27
-rwxr-xr-x3401/CH9/EX9.2/Ex9_2.sce22
-rwxr-xr-x3401/CH9/EX9.3/Ex9_3.sce15
-rwxr-xr-x3401/CH9/EX9.4/Ex9_4.sce13
-rwxr-xr-x3401/CH9/EX9.5/Ex9_5.sce27
-rwxr-xr-x3401/CH9/EX9.6/Ex9_6.sce19
-rwxr-xr-x3401/CH9/EX9.7/Ex9_7.sce14
-rwxr-xr-x3401/CH9/EX9.8/Ex9_8.sce29
80 files changed, 1317 insertions, 0 deletions
diff --git a/3401/CH1/EX1.1/Ex1_1.sce b/3401/CH1/EX1.1/Ex1_1.sce
new file mode 100755
index 000000000..56827aa29
--- /dev/null
+++ b/3401/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,6 @@
+clc
+
+a=5*10^-8// a=5A = 5*10^-8cm
+n=2
+d=n/a^3
+disp(d,"the value of d in atoms per cm^3 is")
diff --git a/3401/CH1/EX1.3/Ex1_3.sce b/3401/CH1/EX1.3/Ex1_3.sce
new file mode 100755
index 000000000..14f637b58
--- /dev/null
+++ b/3401/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,6 @@
+clc
+
+a1=5*10^-8// a=5A = 5*10^-8cm
+n=2// number of atoms is 2
+d=n/(a1*a1*2^0.5)
+disp(d,"the value of d in atoms per cm^2 is")
diff --git a/3401/CH10/EX10.1/Ex10_1.sce b/3401/CH10/EX10.1/Ex10_1.sce
new file mode 100755
index 000000000..864eef9a7
--- /dev/null
+++ b/3401/CH10/EX10.1/Ex10_1.sce
@@ -0,0 +1,9 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+y=0.9967
+//NB/NC=z
+z=(1/y)-1
+disp(z,"NB/NE value is=")
diff --git a/3401/CH10/EX10.10/Ex10_10.sce b/3401/CH10/EX10.10/Ex10_10.sce
new file mode 100755
index 000000000..1de0d5066
--- /dev/null
+++ b/3401/CH10/EX10.10/Ex10_10.sce
@@ -0,0 +1,12 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+alphaF=0.99
+alphaR=0.20
+Ic=1//mA
+IB=0.050//mA
+
+Vcesat=k*T*log(((Ic*(1-alphaR)+IB)*alphaF)/((alphaF*IB-(1-alphaF)*Ic)*alphaR))
+disp(Vcesat,"VCEsat in V is=")
diff --git a/3401/CH10/EX10.12/Ex10_12.sce b/3401/CH10/EX10.12/Ex10_12.sce
new file mode 100755
index 000000000..8629a6fa0
--- /dev/null
+++ b/3401/CH10/EX10.12/Ex10_12.sce
@@ -0,0 +1,33 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+IE=1*10^-3//A
+Cje=1*10^-12//F
+xB=0.5*10^-4//cm
+Dn=25// cm2/s
+xdc=2.4*10^-4//cm
+vs=10^7
+rc=20//ohm
+Cu=0.1*10^-12//F
+Cs=0.1*10^-12//F
+beta=100
+
+re1=(k*T)/IE
+disp(re1,"re1 in ohm is=")
+taue=re1*Cje
+disp(taue,"taue in sec is=")
+taub=xB^2/(2*Dn)
+disp(taub,"taub in sec is=")
+taub1=xdc/vs
+disp(taub1,"taub in sec is=")
+tauc=rc*(Cu+Cs)
+disp(tauc,"tauc in sec is=")
+tauec=taue+taub+taub1+tauc
+disp(tauec,"tauec in sec is=")
+ftau=1/(2*%pi*tauec)
+disp(ftau,"ftau in GHz is=")
+fbeta=ftau/beta
+disp(fbeta,"fbeta in Hz is=")
+
diff --git a/3401/CH10/EX10.2/Ex10_2.sce b/3401/CH10/EX10.2/Ex10_2.sce
new file mode 100755
index 000000000..07597db1c
--- /dev/null
+++ b/3401/CH10/EX10.2/Ex10_2.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+DB=10//cm^2/s
+tauB0=10^-7//sec
+alphatau=0.9967
+//(xB/LB)=z
+z=acosh(1/alphatau)
+disp(z,"xB/LB is=")
+
+LB=sqrt(DB*tauB0)
+disp(LB,"LB in cm is=")
+xB=z*10^-4
+disp(xB,"xB in meter is=")
diff --git a/3401/CH10/EX10.3/Ex10_3.sce b/3401/CH10/EX10.3/Ex10_3.sce
new file mode 100755
index 000000000..33d4bccf0
--- /dev/null
+++ b/3401/CH10/EX10.3/Ex10_3.sce
@@ -0,0 +1,11 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Js0=10^-11//A/cm^2
+Jr0=10^-8//A/cm^2
+delta=0.9967
+VBE=2*k*T*log(delta*10^3/(1-delta))
+disp(VBE,"VBE in V is=")
+
diff --git a/3401/CH10/EX10.4/Ex10_4.sce b/3401/CH10/EX10.4/Ex10_4.sce
new file mode 100755
index 000000000..0b1e78ed8
--- /dev/null
+++ b/3401/CH10/EX10.4/Ex10_4.sce
@@ -0,0 +1,32 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+DE=10//cm^2/s
+DB=25//cm^2/s
+xB=0.70*10^-6//m
+xE=0.50*10^-6//m
+NE=1*10^18//cm^-3
+NB=1*10^16//cm^-3
+VBE=0.65//V
+tauB0=5*10^-7//s
+tauE0=1*10^-7//s
+Jr0=5*10^-8//A/cm62
+pE0=2.25*10^2//cm^-3
+nB0=2.25*10^4//cm^-3
+LE=10^-3//cm
+LB=3.54*10^-3//cm
+
+gamma1=1/(1+(pE0*DE*LB*tanh(0.0198))/(nB0*DB*LE*tanh(0.050)))
+disp(gamma1,"gamma1 is=")
+alphatau=1/cosh(xB/LB)
+disp(alphatau,"alphatau is=")
+Js0=(e*DB*nB0)/(LB*tanh(xB/LB))
+disp(Js0,"Js0 is A/cm^2")
+delta=1/(1+(Jr0/Js0)*exp(-VBE/(2*0.0259)))
+disp(delta,"delta is=")
+a=gamma1*alphatau*delta
+disp(a,"a is=")
+beta1=a/(1-a)
+disp(beta1,"beta1 is=")// ans varies cause of long no.of digits
diff --git a/3401/CH10/EX10.5/Ex10_5.sce b/3401/CH10/EX10.5/Ex10_5.sce
new file mode 100755
index 000000000..83b7f4dfe
--- /dev/null
+++ b/3401/CH10/EX10.5/Ex10_5.sce
@@ -0,0 +1,29 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+NB=5*10^16//cm^-3
+NC=2*10^15//cm^-3
+epsilons=11.7
+eps=8.85*10^-14 //V
+ni=1.5*10**10 //cm^-3
+x=0.70*10^-6
+a=9.96*10^-12///solving the equation
+Vbi=((k*T))*log((NB*NC)/ni^2)
+disp(Vbi,"Vbi in V is=")
+
+//for
+VCB=2 //V
+xdB1=sqrt(a*(Vbi+VCB))
+disp(xdB1,"xdB in meter is=")
+xB1=x-xdB1
+disp(xB1,"xB in meter is=")// textbook ans are wrong
+
+//for
+VCB=10 //V
+xdB2=sqrt(a*(Vbi+VCB))
+disp(xdB2,"xdB in meter is=")
+xB2=x-xdB2
+disp(xB2,"xB in meter is=")// textbook ans are wrong
+
diff --git a/3401/CH10/EX10.7/Ex10_7.sce b/3401/CH10/EX10.7/Ex10_7.sce
new file mode 100755
index 000000000..6c3bd0f68
--- /dev/null
+++ b/3401/CH10/EX10.7/Ex10_7.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300// K
+k=1.3806*10^-23// JK^-1
+e=1.6*10^-19// eV
+ni=1.5*10^10// cm^-3
+NE=10^18// cm^-3
+pE01=(ni*ni)/NE// neglecting bandgap
+disp(pE01,"pE01 the value of pE01 in cm^-3 is")
+NE=10^19// cm^-3
+pE03=(ni*ni)/NE// neglecting bandgap
+disp(pE03,"pE03 the value of pE03 in cm^-3 is")
+pE02=((ni*ni)/NE)*exp(0.030/(k*(T/e)))// including bandgap
+disp(pE02,"pE02 the value of pE02 in cm^-3 is")
+pE04=((ni*ni)/NE)*exp(0.1/(k*(T/e)))// including bandgap
+disp(pE04,"pE04 the value of pE04 in cm^-3 is")
diff --git a/3401/CH10/EX10.8/Ex10_8.sce b/3401/CH10/EX10.8/Ex10_8.sce
new file mode 100755
index 000000000..f18313bb7
--- /dev/null
+++ b/3401/CH10/EX10.8/Ex10_8.sce
@@ -0,0 +1,17 @@
+clc
+
+T=300// K
+k=1.3806*10^-23// JK^-1
+e=1.6*10^-19// eV
+epsilon=8.85*10^-14 //V
+eps=11.7
+NB=10^16//cm^-3
+Vpt=25//V
+WB=0.5*10^-4
+
+
+//Vpt=(e*WB^2*NB*(NC+NB))/(2*epsilon*eps*NC)
+a=(Vpt*2*epsilon*eps)/(e*WB^2*NB)
+disp(a)
+NC=NB/(a-1)
+disp(NC,"NC in cm^-3 is=")
diff --git a/3401/CH10/EX10.9/Ex10_9.sce b/3401/CH10/EX10.9/Ex10_9.sce
new file mode 100755
index 000000000..e8dd0ef64
--- /dev/null
+++ b/3401/CH10/EX10.9/Ex10_9.sce
@@ -0,0 +1,11 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+NB=10^17//cm^-3
+beta=100
+BVCEO=15//V
+
+BVCBO=(beta)^(1/3)*BVCEO
+disp(BVCBO,"in V is=")
diff --git a/3401/CH11/EX11.1/Ex11_1.sce b/3401/CH11/EX11.1/Ex11_1.sce
new file mode 100755
index 000000000..437815f72
--- /dev/null
+++ b/3401/CH11/EX11.1/Ex11_1.sce
@@ -0,0 +1,13 @@
+clc
+
+T=300//K
+Na=10^16// cm^-3
+ni=1.5*10^10 // cm^-3
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+e=1.6*10^-19// eV
+Vt=0.0259 // V
+phifp=Vt*log(Na/ni)
+disp(phifp,"the value of phisp in V is")
+xdT=(4*(epsilons*epsilon0)*phifp/(e*Na))^0.5
+disp(xdT,"the value of xdT in meter is")
diff --git a/3401/CH11/EX11.10/Ex11_10.sce b/3401/CH11/EX11.10/Ex11_10.sce
new file mode 100755
index 000000000..b219bf235
--- /dev/null
+++ b/3401/CH11/EX11.10/Ex11_10.sce
@@ -0,0 +1,18 @@
+clc
+
+T=300// K
+Na=3*10^16// cm^-3
+tox=500*10^-8// tox=500A =500*10^-8 cm
+k=1.3806*10^-23// JK^-1
+VSB=1// V
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+epsilonox=3.9
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+phifp=(((k*T)/e)*log(Na/ni))// phifp=V1*log(Na/ni)=((k*T)/e)*log(Na/ni)
+disp(phifp,"the value of phifp in V is")
+Cox=(epsilon0*epsilonox)/tox
+disp(Cox,"the value of Cox in F/cm^2 is")
+deltaVT=(sqrt(2*e*(epsilon0*epsilons)*Na))/Cox*((sqrt(2*phifp+VSB))-(sqrt(2*phifp)))
+disp(deltaVT,"the value of deltaVT in V is")
diff --git a/3401/CH11/EX11.11/Ex11_11.sce b/3401/CH11/EX11.11/Ex11_11.sce
new file mode 100755
index 000000000..3a6d56f17
--- /dev/null
+++ b/3401/CH11/EX11.11/Ex11_11.sce
@@ -0,0 +1,8 @@
+clc
+
+micron=400// cm^2/V-s
+L=4*10^-6// m
+VT=1// V
+VGS=3// V
+ftau=(micron*(VGS-VT))/(2*%pi*L*L)
+disp(ftau,"the value of fr is")
diff --git a/3401/CH11/EX11.2/Ex11_2.sce b/3401/CH11/EX11.2/Ex11_2.sce
new file mode 100755
index 000000000..d40646132
--- /dev/null
+++ b/3401/CH11/EX11.2/Ex11_2.sce
@@ -0,0 +1,14 @@
+clc
+
+phim=3.20// v
+x=3.25// v
+Eg=1.11//eV
+Na=10^14// cm^-3
+k=1.3806*10^-23// JK^-1
+T=300// K
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+phifp=(((k*T)/e)*log(Na/ni))
+disp(phifp,"the value of phifp in V is")
+phims=phim-(x+(Eg/2)+phifp)
+disp(phims,"work function difference in V is ")
diff --git a/3401/CH11/EX11.3/Ex11_3.sce b/3401/CH11/EX11.3/Ex11_3.sce
new file mode 100755
index 000000000..3d64fe573
--- /dev/null
+++ b/3401/CH11/EX11.3/Ex11_3.sce
@@ -0,0 +1,15 @@
+clc
+
+Nd=10^16// cm^-3
+tox=500*10^-8// a=500A =500*10^-8 cm
+Qss=10^11// electronic charge per cm
+phims=-1.1// V
+e=1.6*10^-19// eV
+epsilonox=3.9
+epsilon0=8.85*10^-14// F/m
+C=((epsilonox*epsilon0)/tox)
+disp(C,"the value of C in F/cm^2 is")
+Qss=10^11*e
+disp(Qss,"the value of Qss in C/cm^2 is")
+VFB=phims-(Qss/C)
+disp(VFB,"the value of VFB in V is" )
diff --git a/3401/CH11/EX11.4/Ex11_4.sce b/3401/CH11/EX11.4/Ex11_4.sce
new file mode 100755
index 000000000..15c9f9796
--- /dev/null
+++ b/3401/CH11/EX11.4/Ex11_4.sce
@@ -0,0 +1,21 @@
+clc
+
+T=300// K
+k=1.3806*10^-23// JK^-1
+Na=3*10^16// cm^-3
+Qss=10^11// cm^-2
+VTN=0.65// V
+phims=-1.13// V
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+epsilonox=3.9
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+phifp=(((k*T)/e)*log(Na/ni))
+disp(phifp,"the value of phifp in V is")
+xdT=(4*(epsilons*epsilon0)*phifp/(e*Na))^0.5
+disp(xdT,"the value of xdT in meter is")
+QSD=e*Na*xdT// [QSD(max)]=e*Na*xdT
+disp(QSD,"the value of QSD in C.cm^2 is")
+tox=((VTN-phims-2*phifp)*(epsilonox*epsilon0))/(QSD-(Qss*e))
+disp(tox,"in meter is")
diff --git a/3401/CH11/EX11.5/Ex11_5.sce b/3401/CH11/EX11.5/Ex11_5.sce
new file mode 100755
index 000000000..cc8a0ed98
--- /dev/null
+++ b/3401/CH11/EX11.5/Ex11_5.sce
@@ -0,0 +1,21 @@
+clc
+
+T=300// K
+k=1.3806*10^-23// JK^-1
+Na=10^14// cm^-3
+Qss=10^10// cm^-2
+tox=500*10^-8// a=500A =500*10^-8 cm
+phims=-0.83// V
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+epsilonox=3.9
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+phifp=(((k*T)/e)*log(Na/ni))
+disp(phifp,"the value of phifp in V is")
+xdT=(4*(epsilons*epsilon0)*phifp/(e*Na))^0.5
+disp(xdT,"the value of xdT in microm is")
+QSD=e*Na*xdT// [QSD(max)]=e*Na*xdT
+disp(QSD,"the value of QSD in C/cm^2")
+VTN=(QSD-(Qss*e))*(tox/(epsilonox*epsilon0))+phims+2*phifp// VTN=(QSD(max)-Qss)*(tox/epsilonox)+phims+2*phifp
+disp(VTN,"the value of VTN in V is")
diff --git a/3401/CH11/EX11.6/Ex11_6.sce b/3401/CH11/EX11.6/Ex11_6.sce
new file mode 100755
index 000000000..b37e64181
--- /dev/null
+++ b/3401/CH11/EX11.6/Ex11_6.sce
@@ -0,0 +1,22 @@
+clc
+
+T=300// K
+k=1.3806*10^-23// JK^-1
+Qss=10^10// cm^-2
+tox=650*10^-8// tox=650A =650*10^-8 cm
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+epsilonox=3.9
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+Nd=2.5*10^14// cm^-3
+phifn=(((k*T)/e)*log(Nd/ni))//phifn=V1*log(Nd/ni)=((k*T)/e)*log(Nd/ni)
+disp(phifn,"the value of phifn in V is")
+xdT=(4*(epsilons*epsilon0)*phifn/(e*Nd))^0.5
+disp(xdT,"the value of xdT in meter is")
+QSD=e*Nd*xdT// [QSD(max)]=e*Na*xdT
+disp(QSD,"the value of QSD in C/cm^2")
+phims=-0.35// V
+VTP=(-QSD-(Qss*e))*(tox/(epsilonox*epsilon0))+phims-(2*phifn)// VTP=(-QSD(max)-Qss)*(tox/epsilonox)+phims+2*phifn
+disp(VTP,"the value of VTP in V is")
+
diff --git a/3401/CH11/EX11.7/Ex11_7.sce b/3401/CH11/EX11.7/Ex11_7.sce
new file mode 100755
index 000000000..e93729faa
--- /dev/null
+++ b/3401/CH11/EX11.7/Ex11_7.sce
@@ -0,0 +1,36 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+espox=3.9
+esp0=8.85*10^-14
+esp=11.7
+tox=550*10^-8//550 Armstrong
+Na=10^16//cm^-3
+ni=1.5*10^10 //cm^-3
+
+Cox=(espox*esp0)/tox
+disp(Cox,"oxide capacitance in F/cm^2 is= ")
+
+phi=(k*T)*log(Na/ni)
+disp(phi,"mini capacitance in V is= ")
+
+xdt=sqrt((4*esp*esp0*phi)/(e*Na))
+disp(xdt,"in cm^-4 is= ")
+
+Cmin=(espox*esp0)/(tox+(espox/esp)*xdt)
+disp(Cmin,"in F/cm^2 is =")
+
+a=Cmin/Cox
+disp(a,"ratio of Cmin to Cox is= ")
+
+Cfb=(espox*esp0)/(tox+(espox/esp)*sqrt((k*T*esp*esp0)/(e*Na)))
+disp(Cfb,"Cfb in F/cm^2 is=")
+
+b=Cfb/Cox
+disp(b,"ratio of Cfb to Cox is= ")
+
+
+
+
diff --git a/3401/CH11/EX11.8/Ex11_8.sce b/3401/CH11/EX11.8/Ex11_8.sce
new file mode 100755
index 000000000..12985cf6f
--- /dev/null
+++ b/3401/CH11/EX11.8/Ex11_8.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Cox=6.9*10^-8//F/cm2
+esp0=8.85*10^-14
+Vtau=0.65//V
+VGS=5//V
+L=1.25*10^-4//cm
+u=650//cm^2/Vs
+IDsat=4*10^-3//A
+
+W=(IDsat*2*L)/(u*Cox*(VGS-Vtau)^2)
+disp(W,"W in meter is=")
diff --git a/3401/CH11/EX11.9/Ex11_9.sce b/3401/CH11/EX11.9/Ex11_9.sce
new file mode 100755
index 000000000..283c07e34
--- /dev/null
+++ b/3401/CH11/EX11.9/Ex11_9.sce
@@ -0,0 +1,19 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Cox=6.9*10^-8//F/cm2
+esp0=8.85*10^-14
+Vtau=0.65//V
+VGS1=1.5//V
+VGS2=2.5//V
+VDS=0.10//V
+L=2*10^-6//cm
+u=650//cm^2/Vs
+ID1=35*10^-6//A
+ID2=75*10^-6//A
+W=15*10^-6//m
+
+un=(ID2-ID1)*L/(W*Cox*(VGS2-VGS1)*VDS)
+disp(un,"un in cm^2/Vs is=")
diff --git a/3401/CH2/EX2.1/Ex2_1.sce b/3401/CH2/EX2.1/Ex2_1.sce
new file mode 100755
index 000000000..91f0cbc73
--- /dev/null
+++ b/3401/CH2/EX2.1/Ex2_1.sce
@@ -0,0 +1,10 @@
+clc
+
+lambda=0.708*10^-8// cm
+h=6.625*10^-34// J*s Plank's constant
+c=3*10^10// cm/s
+e=1.6*10^-19// eV
+E=h*c/lambda// E=hv=hc/lambda
+disp(E,"the value of E in J is")
+E=E/e
+disp(E,"the value of E in eV is")
diff --git a/3401/CH2/EX2.2/Ex2_2.sce b/3401/CH2/EX2.2/Ex2_2.sce
new file mode 100755
index 000000000..f84fd017f
--- /dev/null
+++ b/3401/CH2/EX2.2/Ex2_2.sce
@@ -0,0 +1,10 @@
+clc
+
+m=9.11*10^-31// kg*m/s
+v=10^5//m/s
+h=6.625*10^-34//js
+
+p=m*v
+disp(p,"momentum is")
+lambda=h/p
+disp(lambda,"de broglie wavelength in meter is")
diff --git a/3401/CH2/EX2.3/Ex2_3.sce b/3401/CH2/EX2.3/Ex2_3.sce
new file mode 100755
index 000000000..72edd006c
--- /dev/null
+++ b/3401/CH2/EX2.3/Ex2_3.sce
@@ -0,0 +1,24 @@
+clc
+
+a=5*10^-8// a=5A = 5*10^-8cm
+h=1.054*10^-34// J*s Plank's constant
+m=9.11*10^-31// kg*m/s
+e=1.6*10^-19// eV
+
+n=1
+En=(h^2*n^2*%pi^2)/(2*m*a^2)
+disp(En,"the value of En in J")
+En=(En/e)
+disp(En,"the value of En in eV")
+
+n=2
+E2=(h^2*n^2*%pi^2)/(2*m*a^2)
+disp(En,"the value of E2 in J")
+E2=(E2/e)
+disp(E2,"the value of E2 in eV")
+
+n=3
+E3=(h^2*n^2*%pi^2)/(2*m*a^2)
+disp(E3,"the value of E3 in J")
+E3=(E3/e)
+disp(E3,"the value of E3 in eV")
diff --git a/3401/CH2/EX2.4/Ex2_4.sce b/3401/CH2/EX2.4/Ex2_4.sce
new file mode 100755
index 000000000..cafb61e77
--- /dev/null
+++ b/3401/CH2/EX2.4/Ex2_4.sce
@@ -0,0 +1,12 @@
+clc
+
+v=10^5// m/s
+m=9.11*10^-31// kg*m/s
+e=1.6*10^-19// eV
+h=1.054*10^-34
+E=0.5*m*v*v
+disp(E,"the value of E in J is eV")
+E1=E/e // value of E in eV
+disp(E1,"eV=")
+d=sqrt((h*h)/(2*m*E))
+disp(d,"the value of d in m is ")
diff --git a/3401/CH2/EX2.5/Ex2_5.sce b/3401/CH2/EX2.5/Ex2_5.sce
new file mode 100755
index 000000000..1d4b6749a
--- /dev/null
+++ b/3401/CH2/EX2.5/Ex2_5.sce
@@ -0,0 +1,12 @@
+clc
+
+E=2// eV
+V0=20// eV
+a=3*10^-10// a=3A = 3*10^-10 m
+m=9.11*10^-31// kg*m/s
+e=1.6*10^-19// eV
+h=1.054*10^-34// J*s
+K2=((2*m*(V0-E)*e)/(h*h))^0.5
+disp(K2,"the value of K2 in m^-1 is")
+T=16*(E/V0)*(1-E/V0)*exp(-2*K2*a)
+disp(T,"the value of T is")
diff --git a/3401/CH3/EX3.1/Ex3_1.sce b/3401/CH3/EX3.1/Ex3_1.sce
new file mode 100755
index 000000000..789ac8c6c
--- /dev/null
+++ b/3401/CH3/EX3.1/Ex3_1.sce
@@ -0,0 +1,14 @@
+clc
+
+v1=10^5 //m per sec
+deltav=0.01 //m per sec
+m=9.11*10^-31// kg*m/s
+e=1.6*10**-9 //Coulombs
+
+//deltaE=0.5*m*(v2**2-v1**2)
+//deltav=v2-v1.............deltav<<v1
+
+deltaE=m*v1*deltav
+disp(deltaE,"in J is")
+deltaE=deltaE/e
+disp(deltaE,"in eV is")
diff --git a/3401/CH3/EX3.2/Ex3_2.sce b/3401/CH3/EX3.2/Ex3_2.sce
new file mode 100755
index 000000000..0e5b916e5
--- /dev/null
+++ b/3401/CH3/EX3.2/Ex3_2.sce
@@ -0,0 +1,18 @@
+clc
+
+P=10
+a=5*10^-10 // a=5 Armstrong
+h=1.054*10^-34
+m=9.11*10^-31 //kg
+
+//alpha*a=%pi
+//sqrt((2*m*E2)/h^2)*a=%pi
+E2=(%pi^2*h^2)/(2*m*a^2)
+disp(E2,'E2= %f J\n')
+E2=E2*6.2415*10^18
+disp(E2,'E2= %f eV')
+
+E1=1.053 //eV
+AE=E2-E1
+disp(AE,'AE= %f eV')
+
diff --git a/3401/CH3/EX3.3/Ex3_3.sce b/3401/CH3/EX3.3/Ex3_3.sce
new file mode 100755
index 000000000..dcb5e60fe
--- /dev/null
+++ b/3401/CH3/EX3.3/Ex3_3.sce
@@ -0,0 +1,7 @@
+clc
+
+m=9.11*10^-31 //kg
+E=1.6*10^-19 //C
+h=6.625*10^-34 //J sec
+N=(4*%pi*(2*m)^(3/2)*2*E^(3/2))/(h^3*3)
+disp(N,'E2= %f per meter^3\n')
diff --git a/3401/CH3/EX3.4/Ex3_4.sce b/3401/CH3/EX3.4/Ex3_4.sce
new file mode 100755
index 000000000..7029cd472
--- /dev/null
+++ b/3401/CH3/EX3.4/Ex3_4.sce
@@ -0,0 +1,8 @@
+clc
+
+y=10
+Ni=10
+gi=10
+//(gi-Ni)!=1
+z=(factorial(gi)/factorial(Ni)*factorial(gi-Ni))
+printf('z=%f\n',z)
diff --git a/3401/CH3/EX3.5/Ex3_5.sce b/3401/CH3/EX3.5/Ex3_5.sce
new file mode 100755
index 000000000..5eb25e7ff
--- /dev/null
+++ b/3401/CH3/EX3.5/Ex3_5.sce
@@ -0,0 +1,8 @@
+clc
+gi=10
+Ni=9
+x=gi-Ni
+
+//factorial(gi-Ni)=1
+z=(factorial(gi)/factorial(Ni)*factorial(gi-Ni))
+printf('z=%f\n',z)
diff --git a/3401/CH3/EX3.6/Ex3_6.sce b/3401/CH3/EX3.6/Ex3_6.sce
new file mode 100755
index 000000000..9bf6b31fc
--- /dev/null
+++ b/3401/CH3/EX3.6/Ex3_6.sce
@@ -0,0 +1,9 @@
+clc
+
+T=300 //K
+//a=E-Ef
+//a=3*k*T
+//fF(E)=1/(1+exp(E-Ef/kT))
+z=exp(3)
+y=1/(1+z) //y=fF(E)
+printf('fF(E)=%f\n',y)
diff --git a/3401/CH3/EX3.7/Ex3_7.sce b/3401/CH3/EX3.7/Ex3_7.sce
new file mode 100755
index 000000000..dc2193807
--- /dev/null
+++ b/3401/CH3/EX3.7/Ex3_7.sce
@@ -0,0 +1,16 @@
+clc
+
+Ef=6.25 //eV
+E=5.95 //eV
+p=0.01
+k= 8.617 *10^-5 //eV K^-1
+
+//p=1-fF(E)
+//fF(E)=1/(1+exp(E-Ef/kT))
+//p=1-1/(1+exp(E-Ef/kT)) equation 1
+x=E-Ef
+
+y=(1/(1-p))-1 // solving the above equtaion 1
+
+T=x/(k*log(y))
+printf('T=%f K\n',T)
diff --git a/3401/CH4/EX4.1/Ex4_1.sce b/3401/CH4/EX4.1/Ex4_1.sce
new file mode 100755
index 000000000..ed6e6fe84
--- /dev/null
+++ b/3401/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300 //K
+Nc=2.8*10^19 //cm^-3
+k=8.617 *10^-5 //eV K^-1
+//E=Ec
+//a=Ec-Ef
+a=-0.25
+//fF(E)=1/(1+exp(Ec-Ef/kT)) =exp(-(Ec-Ef)/kT)
+//b=k*T
+//disp(b)
+c=k*T
+fFE=exp(a/c)
+disp(fFE,"fF(E) value is=")
+n0=Nc*exp(-a/(k*T))
+disp(n0,"n0 value in cm^-3 is")
diff --git a/3401/CH4/EX4.10/Ex4_10.sce b/3401/CH4/EX4.10/Ex4_10.sce
new file mode 100755
index 000000000..293d94b13
--- /dev/null
+++ b/3401/CH4/EX4.10/Ex4_10.sce
@@ -0,0 +1,16 @@
+clc
+
+Nd=5*10**13 //cm^-3
+Na=0
+ni=2.4*10**13 //cm^-3
+T=300 //K
+k=8.617*10^-5 //eV/K
+
+n0=((Nd-Na)/2)+sqrt((((Nd-Na)/2)^2)+ni^2)
+disp(n0,"n0 in cm^-3 is")
+
+p0=(ni^2)/n0
+disp(p0,"p0 in cm^-3")
+
+
+
diff --git a/3401/CH4/EX4.11/Ex4_11.sce b/3401/CH4/EX4.11/Ex4_11.sce
new file mode 100755
index 000000000..c8b2cc85c
--- /dev/null
+++ b/3401/CH4/EX4.11/Ex4_11.sce
@@ -0,0 +1,19 @@
+clc
+
+T=300 //K
+Nd=3*10**15 //cm**-3
+Na=10**16 //cm**-3
+k=8.617*10^-5 //eV/K
+ni=1.5*10**16 //cm^-3
+
+a=((Na-Nd)/2)
+disp(a)
+b=sqrt(((Na-Nd)/2)^2+(ni^2))
+disp(b)
+
+p0=a+b
+disp(p0,"p0 in cm^-3 is=") //textbook ans is wrpng
+
+n0=(ni^2)/p0
+disp(n0,"n0 in cm^-3 is=") //textbook ans is wrpng
+
diff --git a/3401/CH4/EX4.12/Ex4_12.sce b/3401/CH4/EX4.12/Ex4_12.sce
new file mode 100755
index 000000000..cec464bc8
--- /dev/null
+++ b/3401/CH4/EX4.12/Ex4_12.sce
@@ -0,0 +1,24 @@
+clc
+
+T1=550 //K
+T2=300 //K
+Nc=2.8*10**19 //cm**-3
+Ne=1.04*10**19 //cm**-3
+k=8.617*10^-5 //eV/K
+Eg=1.12
+
+//ni^2=Nc*Ne*exp(-Eg/(k*T))
+
+ni=sqrt(Nc*Ne*(T1/T2)^3*exp((-Eg/(k*T2))*(T2/T1)))
+disp(ni ,"ni in cm^-3 is = ")
+
+//n0=1.05*Nd
+//n0=(Nd/2)+sqrt((Nd/2)^2+ni^2) equation 1
+
+a=(((1.05*2)-1)/2)^2 //simplifying and sloving the about equaton 1
+b=ni^2
+c=((4*a)-1)/4
+Nd=sqrt(b/c)
+disp(Nd,"Nd value in cm^-3 is=")
+
+
diff --git a/3401/CH4/EX4.13/Ex4_13.sce b/3401/CH4/EX4.13/Ex4_13.sce
new file mode 100755
index 000000000..acce3a3de
--- /dev/null
+++ b/3401/CH4/EX4.13/Ex4_13.sce
@@ -0,0 +1,14 @@
+clc
+
+T=300 //K
+Na=10**16 //cm^-3
+Nc=2.8*10**19 //cm**-3
+k=8.617*10^-5 //eV/K
+
+//a=Ec-Ef
+a=0.20 //eV
+//Ec-Ef=k*T*log(Ne/(Nd-Na))
+//Nd-Na=Nc*exp(-(Ec-Ef)/k*T)
+
+Nd=Na+Nc*exp(-a/(k*T))
+disp(Nd,"Nd in cm^-3 is= ")
diff --git a/3401/CH4/EX4.14/Ex4_14.sce b/3401/CH4/EX4.14/Ex4_14.sce
new file mode 100755
index 000000000..e32eda8aa
--- /dev/null
+++ b/3401/CH4/EX4.14/Ex4_14.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300 //K
+k=8.617*10^-5 //eV/K
+ni=1.5*10**10 //cm^-3
+//Ef-Ea=a
+a=3*k*T
+//Ea-Ev=b
+b=0.045 //eV
+Eg=1.12 //eV
+
+//Efi-Ef=(Eg/2)-(Ea-Ev)-(Ef-Ea)=kT*log(Na/ni)
+c=(Eg/2)-(b)-(a)
+Na=ni*exp(c/(k*T))
+disp(Na,"Na in cm^-3 is= ")
+
diff --git a/3401/CH4/EX4.2/Ex4_2.sce b/3401/CH4/EX4.2/Ex4_2.sce
new file mode 100755
index 000000000..90c83ebdd
--- /dev/null
+++ b/3401/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,16 @@
+clc
+
+T1=400 //K
+T2=300 //K
+N=1.04*10^19 //cm^-3
+k=8.617*10^-5 //eV/K
+
+Nv=N*(T1/T2)^(3/2)
+disp(Nv,'Nv=%f cm^-3\n\n')
+
+a=k*T2*(T1/T2)
+disp(a,'kT=%f eV\n\n')
+
+//p0=Nv*exp(-(Ef-Ev)/kT)
+p0=Nv*exp(-0.27/a)
+disp(p0,'p0=%f cm^-3 \n')
diff --git a/3401/CH4/EX4.3/Ex4_3.sce b/3401/CH4/EX4.3/Ex4_3.sce
new file mode 100755
index 000000000..03db78419
--- /dev/null
+++ b/3401/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,15 @@
+clc
+
+T1=300 //K
+T2=450 //K
+Nc=4.7*10^17 //cm^-3
+N=7*10^18 //cm^-3
+k=8.617*10^-5 //eV/K
+a=k*T1*(T2/T1)
+disp(a,'kT=%f eV\n\n')
+
+ni=sqrt(Nc*N*exp(-1.42/(k*T1)))
+disp(ni,'ni=%f cm^-3\n\n')
+
+ni=sqrt(Nc*N*(T2/T1)^3*exp(-1.42/a))
+disp(ni,'ni=%f cm^-3\n\n')
diff --git a/3401/CH4/EX4.4/Ex4_4.sce b/3401/CH4/EX4.4/Ex4_4.sce
new file mode 100755
index 000000000..cd7b0cff8
--- /dev/null
+++ b/3401/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,12 @@
+clc
+
+T=300 //K
+mn=1.08 //m0
+mp=0.56 //m0
+k=8.617*10^-5 //eV/K
+
+//Efi-Emidgap=(3/4)k*T*log(mp/mn)
+//a=Efi-Emidgap
+
+a=(0.75)*k*T*log(mp/mn)
+printf('Efi-Emidgap=%f eV',a) //textbook ans is wrong
diff --git a/3401/CH4/EX4.5/Ex4_5.sce b/3401/CH4/EX4.5/Ex4_5.sce
new file mode 100755
index 000000000..54a4893bb
--- /dev/null
+++ b/3401/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300 //K
+Nn=2.8*10^19 //cm^-3
+Np=1.04*10^19 //cm^-3
+//a=Ef-Ev
+an=0.25 //eV
+ap=0.87 //eV
+k=8.617*10^-5 //eV/K
+
+n0=Nn*exp(-an/(k*T))
+disp(n0,"n0 in cm^-3 is=")
+
+p0=Np*exp(-ap/(k*T))
+disp(p0,"p0 in cm^-3 is=")
diff --git a/3401/CH4/EX4.6/Ex4_6.sce b/3401/CH4/EX4.6/Ex4_6.sce
new file mode 100755
index 000000000..d95296064
--- /dev/null
+++ b/3401/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,11 @@
+clc
+
+nf=2
+T=300 //K
+Nc=2.8*10**19 //cm^-3
+
+//n0=(2/sqrt(%pi))*Nc*F12(nf)
+//x=F12(nf)=2.3
+x=2.3
+n0=(2/sqrt(%pi))*Nc*x
+disp(n0,"n0 in cm^-3 is= ")
diff --git a/3401/CH4/EX4.7/Ex4_7.sce b/3401/CH4/EX4.7/Ex4_7.sce
new file mode 100755
index 000000000..6a7a69c6e
--- /dev/null
+++ b/3401/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,12 @@
+clc
+
+Nd=10**16 //cm^-3
+Ne=2.8*10**19 //cm^-3
+T=300 //K
+
+//(nd/(n0+nd))=z=1/(1+(Ne/2*Nd)*exp(-(Ec-Ed)/kT))
+//y=Ec-Ed
+y=0.045
+k=8.617*10^-5 //eV/K
+z=1/(1+(Ne/(2*Nd))*exp(-y/(k*T)))
+disp(z,"the donor states is=")
diff --git a/3401/CH4/EX4.9/Ex4_9.sce b/3401/CH4/EX4.9/Ex4_9.sce
new file mode 100755
index 000000000..1f48a5ba2
--- /dev/null
+++ b/3401/CH4/EX4.9/Ex4_9.sce
@@ -0,0 +1,14 @@
+clc
+
+Nd=10**16 //cm^-3
+Na=0
+ni=1.5*10**10 //cm^-3
+T=300 //K
+k=8.617*10^-5 //eV/K
+
+n0=((Nd-Na)/2)+sqrt((((Nd-Na)/2)^2)+ni^2)
+disp(n0,"n0 in cm^-3 is")
+
+p0=(ni^2)/n0
+disp(p0,"p0 in cm^-3")
+
diff --git a/3401/CH5/EX5.1/Ex5_1.sce b/3401/CH5/EX5.1/Ex5_1.sce
new file mode 100755
index 000000000..aac57710d
--- /dev/null
+++ b/3401/CH5/EX5.1/Ex5_1.sce
@@ -0,0 +1,21 @@
+clc
+
+Nd=10**16 //cm^-3
+Na=0
+ni=1.8*10**6 //cm^-3
+T=300 //K
+k=8.617*10^-5 //eV/K
+E=10 //V/cm^2
+e=1.6*10**-19
+un=8500
+
+n0=((Nd-Na)/2)+sqrt((((Nd-Na)/2)^2)+ni^2)
+disp(n0,"n0 in cm^-3 is")
+
+p0=(ni^2)/n0
+disp(p0,"p0 in cm^-3")
+
+//Jdrf=e*(un*n0+up*p0)*E= e*un*Nd*E
+Jdrf=e*un*Nd*E
+disp(Jdrf,"Jdrf in A/cm^2 is= ")
+
diff --git a/3401/CH5/EX5.2/Ex5_2.sce b/3401/CH5/EX5.2/Ex5_2.sce
new file mode 100755
index 000000000..e96e8a2cb
--- /dev/null
+++ b/3401/CH5/EX5.2/Ex5_2.sce
@@ -0,0 +1,27 @@
+clc
+
+Na=10**17 //cm**-3
+T=300 //K
+k=8.617*10^-5 //eV/K
+E=10 //V/cm^2
+e=1.6*10**-19 //C
+
+
+//sigma=e*un*n0=e*un*(Nd-Na)
+ //if
+Nd=2*10**17 //cm**-3
+sigma=8.16//ohm/cm
+un=sigma/(e*(Nd-Na))
+disp(un,"un in cm^2/Vs is=")
+
+//if
+Nd=5*10**17 //cm^-3
+sigma=20.8//ohm/cm
+un=sigma/(e*(Nd-Na))
+disp(un,"un in cm^2/Vs is=")
+
+//if
+Nd=3.5*10**17 //cm^-3
+sigma=16//ohm/cm
+un=sigma/(e*(Nd-Na))
+disp(un,"un in cm^2/Vs is=")
diff --git a/3401/CH5/EX5.3/Ex5_3.sce b/3401/CH5/EX5.3/Ex5_3.sce
new file mode 100755
index 000000000..9208174fe
--- /dev/null
+++ b/3401/CH5/EX5.3/Ex5_3.sce
@@ -0,0 +1,27 @@
+clc
+
+V=5 //v
+R=10*10^3 //ohm
+J=50 //A/cm^2
+E=100
+Na=1.25*10**16 //cm^-3
+Nd=5*10**15 //cm^-3
+e=1.6*10**-19 //C
+up=410 //cm**2/Vs
+
+I=V/R
+disp(I,"I current in Ampere is=")
+
+A=I/J
+disp(A,"A cross sectional area in cm^2 is=")
+
+L=V/E
+disp(L,"L length of resistor in cm is= ")
+
+sigma=L/(R*A)
+disp(sigma,"sigma conductivity in per ohm cm is=")
+
+//sigma=e*up*p0=e*up*(Na-Nd)
+sigma=e*up*(Na-Nd)
+disp(sigma,"sigma conductivity in per ohm cm is=")
+
diff --git a/3401/CH5/EX5.4/Ex5_4.sce b/3401/CH5/EX5.4/Ex5_4.sce
new file mode 100755
index 000000000..ccb7bacde
--- /dev/null
+++ b/3401/CH5/EX5.4/Ex5_4.sce
@@ -0,0 +1,11 @@
+clc
+
+T=300 // ..K
+Dn=225 //cm^2/s
+e=1.6*10^-19// C
+deltax=0.10 //cm
+deltan=(1*10**18-7*10**17)//cm^-3
+
+//Jnxdif=e*Dx*derivative (n,x)=e*Dn*(deltan/deltax)
+Jnxdif=e*Dn*(deltan/deltax)
+disp(Jnxdif,"diffusion current density in A/cm^2 is= ")
diff --git a/3401/CH5/EX5.5/Ex5_5.sce b/3401/CH5/EX5.5/Ex5_5.sce
new file mode 100755
index 000000000..7a6017f5d
--- /dev/null
+++ b/3401/CH5/EX5.5/Ex5_5.sce
@@ -0,0 +1,13 @@
+clc
+
+T=300 //K
+k=8.617*10^-5 //eV/K
+//derivative(Ndx,x)=a
+a=-10**19 //cm**-4
+Ndx=(10^16-10^19) //cm^3
+l=1
+
+//Ex=-(k*T/l)*(1/Ndx)*derivative(Ndx,x)
+Ex=-(k*T/l)*(1/Ndx)*a*10^3
+disp(Ex,"induced electric field in V/cm is= ")
+
diff --git a/3401/CH5/EX5.6/Ex5_6.sce b/3401/CH5/EX5.6/Ex5_6.sce
new file mode 100755
index 000000000..9ef7daf5e
--- /dev/null
+++ b/3401/CH5/EX5.6/Ex5_6.sce
@@ -0,0 +1,10 @@
+clc
+
+T=300 //K
+u=1000 //cm^2/s
+k=8.617*10^-5 //eV/K
+e=1.6*10**-19 //C
+
+//D=((k*T)/e)*u
+D=k*T*u
+disp(D,"diffusion coefficient in cm^2/s is =")
diff --git a/3401/CH5/EX5.7/Ex5_7.sce b/3401/CH5/EX5.7/Ex5_7.sce
new file mode 100755
index 000000000..11be3e258
--- /dev/null
+++ b/3401/CH5/EX5.7/Ex5_7.sce
@@ -0,0 +1,15 @@
+clc
+
+I=10^-3 //A
+Bz=5*10**-2 //500gauss
+e=1.6810*10**-19 //C
+Vh=-6.25*10^-3 //V
+Vx=12.5 //V
+W=10**-4//m
+d=10**-5 //m
+
+u=-(I*Bz)/(e*Vh*d)
+disp(u,"electron concemtration in m^-3 is= ")
+
+un=(I*I)/(e*Bz*Vx*W*d)
+disp(un,"electron mobility in /Vs is=")// textbook ans is wrong
diff --git a/3401/CH6/EX6.10/Ex6_10.sce b/3401/CH6/EX6.10/Ex6_10.sce
new file mode 100755
index 000000000..b13f654d8
--- /dev/null
+++ b/3401/CH6/EX6.10/Ex6_10.sce
@@ -0,0 +1,13 @@
+clc
+
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Dp=10//cm^2/sec
+Lp=31.6*10**-4 //m
+g1taup0=10^14 //cm^-3
+deltap0=10^13 //cm6-3
+
+//deltap0=g1taup0*[g/((Dp/Lp)+s)]
+s=(Dp/Lp)*((g1taup0/deltap0)-1)
+disp(s,"surface recombination velocity in cm per sec is= ")
+
diff --git a/3401/CH6/EX6.5/Ex6_5.sce b/3401/CH6/EX6.5/Ex6_5.sce
new file mode 100755
index 000000000..efb787d4e
--- /dev/null
+++ b/3401/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,17 @@
+clc
+
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+un=1200
+Nd=10^16 //cm^-3
+esp0=8.85*10^-14
+espr=11.7
+
+sigma=e*un*Nd
+disp(sigma,"conductivity in per ohm cm is= ")
+
+esp=espr*esp0
+disp(esp,"permittivity of silicon in F/cm")
+
+taud=esp/sigma
+disp(taud,"dielectric relaxtion time constant in sec is= ")
diff --git a/3401/CH6/EX6.6/Ex6_6.sce b/3401/CH6/EX6.6/Ex6_6.sce
new file mode 100755
index 000000000..1e0de9083
--- /dev/null
+++ b/3401/CH6/EX6.6/Ex6_6.sce
@@ -0,0 +1,22 @@
+clc
+
+T=300//K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+n0=10^15 //cm^-3
+p0=10^5 //cm^-3
+ni=10^10 //cm^-3
+deltan=10**13 //cm**-3
+deltap=10**13 //cm**-3
+
+//Ef-Efi=a
+a=(k*T)*log(n0/ni)
+disp(a,"fermi level for thermal equlibrium in eV is= ")
+
+//Efn-Efi=b
+b=(k*T)*log((n0+deltan)/ni)
+disp(b,"quasi fermi level for electrons in eV is= ")
+
+//Efi-Efp=c
+c=(k*T)*log((p0+deltap)/ni)
+disp(c,"quasi fermi level for holes in eV is= ")
diff --git a/3401/CH6/EX6.8/Ex6_8.sce b/3401/CH6/EX6.8/Ex6_8.sce
new file mode 100755
index 000000000..885a33b29
--- /dev/null
+++ b/3401/CH6/EX6.8/Ex6_8.sce
@@ -0,0 +1,24 @@
+clc
+
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+x=0
+taup0=10**-6//ses
+taup01=10**-7 //sec
+deltapb=10**14 //cm**-3
+Dn=10 //cm^2/sec
+Dp=10 //cm^2/sec
+B=-9*10^13
+
+deltap=deltapb*(taup01/taup0)
+disp(deltap,"deltap in cm^-3 is=")
+
+g=deltap/taup0
+disp(g,"g generation in cm^-3s^-1 is= ")
+
+//deltapx=10^14*(1-0.9*exp(-x/Lp))
+Lp=sqrt(Dp*taup0)
+disp(Lp,"Lp in meter is=")
+
+deltapx=10^14*(1-0.9*exp(-x/Lp))
+disp(deltapx,"distance from the surface is =")
diff --git a/3401/CH7/EX7.1/Ex7_1.sce b/3401/CH7/EX7.1/Ex7_1.sce
new file mode 100755
index 000000000..b4a1ebd6a
--- /dev/null
+++ b/3401/CH7/EX7.1/Ex7_1.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300//K
+Na=1*10^18// Na=L*10^18 cm^-3
+Nd=1*10^15// cm^-3
+ni=1.5*10^10// cm^-3
+e=1.6*10^-19// eV
+k=1.3806*10^-23// JK^-1
+Vbi=(((k*T)/e)*log(Na*Nd/ni^2))
+disp(Vbi,"the value of Vbi in V is")
+
+//changing the Na value
+Na=10^16// Na=l*10^16 cm^1
+Vbi=(((k*T)/e)*log(Na*Nd/ni^2))
+disp(Vbi,"the value of Vbi in V is")
diff --git a/3401/CH7/EX7.2/Ex7_2.sce b/3401/CH7/EX7.2/Ex7_2.sce
new file mode 100755
index 000000000..f502d447d
--- /dev/null
+++ b/3401/CH7/EX7.2/Ex7_2.sce
@@ -0,0 +1,14 @@
+clc
+
+T=300//K
+Na=10^16// cm^-3
+Nd=10^15// cm^-3
+e=1.6*10^-19// eV
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+Vbi=0.635// V
+W=(((2*(epsilons*epsilon0)*Vbi)/e)*((Na+Nd)/(Na*Nd)))^0.5
+disp(W,"W spacing charge width in meter is")
+xn=0.864*10^-4// xp=0.086 micro m
+Emax=-e*Nd*xn/(epsilons*epsilon0)
+disp(Emax,"the value of Emax in V/cm is")
diff --git a/3401/CH7/EX7.3/Ex7_3.sce b/3401/CH7/EX7.3/Ex7_3.sce
new file mode 100755
index 000000000..381a17a5c
--- /dev/null
+++ b/3401/CH7/EX7.3/Ex7_3.sce
@@ -0,0 +1,12 @@
+clc
+
+T=300//K
+Na=10^16// cm^-3
+Nd=10^15// cm^-3
+e=1.6*10^-19// eV
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+Vbi=0.635// V
+Vr=5// V
+W=(((2*(epsilons*epsilon0)*(Vbi+Vr))/e)*((Na+Nd)/(Na*Nd)))^0.5
+disp(W,"the value of W in m is")
diff --git a/3401/CH7/EX7.4/Ex7_4.sce b/3401/CH7/EX7.4/Ex7_4.sce
new file mode 100755
index 000000000..2b4eb39c5
--- /dev/null
+++ b/3401/CH7/EX7.4/Ex7_4.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300//K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Vr=25 //V
+Emax=3*10^5//V/cm
+Na=10^18//cm^-3
+esp=11.7
+esp0=8.85*10^-14
+
+//Emax=sqrt((2*e*Vr*(Na+Nd))/(esp*(Na+Nd))
+x=((Emax^2)*esp*esp0)/(2*e*Vr) //solving the above equation we get
+
+Nd=(Na*x)/(Na-x)
+disp(Nd,"doping concentration in cm^-3 is =")
diff --git a/3401/CH7/EX7.5/Ex7_5.sce b/3401/CH7/EX7.5/Ex7_5.sce
new file mode 100755
index 000000000..44e7188b7
--- /dev/null
+++ b/3401/CH7/EX7.5/Ex7_5.sce
@@ -0,0 +1,14 @@
+clc
+
+T=300//K
+Na=10^16// cm^-3
+Nd=10^15// cm^-3
+e=1.6*10^-19// eV
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+Vbi=0.635// V
+Vr=5// V
+C=((e*(epsilons*epsilon0)*Na*Nd)/(2*(Vbi+Vr)*(Na+Nd)))^0.5
+A=10^-4// cm^2
+C=C*A
+disp(C,"the value of c with cross section area in F is")
diff --git a/3401/CH7/EX7.6/Ex7_6.sce b/3401/CH7/EX7.6/Ex7_6.sce
new file mode 100755
index 000000000..c9783f95e
--- /dev/null
+++ b/3401/CH7/EX7.6/Ex7_6.sce
@@ -0,0 +1,14 @@
+clc
+
+T=300//K
+ni=1.5*10^10// cm^-3
+Nd=10^15// cm^-3
+Vbi=0.855// V
+e=1.6*10^-19// eV
+k=1.3806*10^-23// JK^-1
+S=1.32*10^15// Slope=1.32*10^-15 (F/cm^2)^-2*(V)-1
+epsilon0=8.85*10^-14// F/m
+epsilons=11.7
+Nd=2/((e*(epsilons*epsilon0))*S)
+Na=((ni*ni)/Nd)*exp((e*Vbi)/(k*T))
+disp(Na,"the value of Na in cm^-3 is")
diff --git a/3401/CH8/EX8.1/Ex8_1.sce b/3401/CH8/EX8.1/Ex8_1.sce
new file mode 100755
index 000000000..2f8674cd6
--- /dev/null
+++ b/3401/CH8/EX8.1/Ex8_1.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300// K
+ni=1.5*10^10 //cm^-3
+k=8.617*10^-5 //eV/K
+Nd=10^16 //cm^-3
+Va=0.60 //V
+
+//pn=pn0*exp(e*Va/k*T)
+pn0=ni^2/Nd
+disp(pn0,"thermal equlibrium minority carrier hole concentration in cm^-3 is= ")
+
+pn=pn0*exp(Va/(k*T))
+disp(pn,"minority carrier hole concentration in cm^-3 is= ")
+
diff --git a/3401/CH8/EX8.2/Ex8_2.sce b/3401/CH8/EX8.2/Ex8_2.sce
new file mode 100755
index 000000000..e99028b56
--- /dev/null
+++ b/3401/CH8/EX8.2/Ex8_2.sce
@@ -0,0 +1,16 @@
+clc
+
+T=300//K
+Na=10**16 //cm^-3
+Nd=10**16 //cm^-3
+ni=1.5*10**10 //cm^-3
+Dn=25 //cm^2/s
+Dp=10 //cm^2/s
+taup0=5*10^-7 //s
+taun0=5*10^-7 //s
+epsilon=11.7
+e=1.6*10**-19 //C
+
+//J=(E*D*np0/Ln)+(E*D*pn0/Lp)
+Js=e*ni^2*[(1/Na)*sqrt(Dn/taun0)+(1/Nd)*sqrt(Dp/taup0)]
+disp(Js,"reverse saturation current density in A/cm^2 is=")
diff --git a/3401/CH8/EX8.3/Ex8_3.sce b/3401/CH8/EX8.3/Ex8_3.sce
new file mode 100755
index 000000000..849184c85
--- /dev/null
+++ b/3401/CH8/EX8.3/Ex8_3.sce
@@ -0,0 +1,23 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Jp=5 //A/cm^2
+Jn=20 //A/cm^2
+Va=0.65 //V
+ni=1.5*10**10 //cm^-3
+Dn=25 //cm^2/s
+Dp=10 //cm^2/s
+taup0=5*10^-7 //s
+taun0=5*10^-7 //s
+
+//Jn=(e*Dn*np0/Ln)*[exp(e*Va/k*T)-1]==e*sqrt(Dn/taun0)*(ni^2/Na)*[exp(e*Va/k*T)-1]
+
+Na=e*sqrt(Dn/taun0)*(ni^2/Jn)*[(exp(Va/(k*T)))-1]
+disp(Na,"Na elctron diffusion current density in cm^-3 is= ")
+
+//Jp=(e*Dp*pn0/Lp)*[exp(e*Va/k*T)-1]==e*sqrt(Dp/taup0)*(ni^2/Nd)
+
+Nd=e*sqrt(Dp/taup0)*(ni^2/Jp)*[(exp(Va/(k*T)))-1]
+disp(Nd,"Nd hole diffusion current density in cm^-3 is= ")
diff --git a/3401/CH8/EX8.4/Ex8_4.sce b/3401/CH8/EX8.4/Ex8_4.sce
new file mode 100755
index 000000000..06b5520d9
--- /dev/null
+++ b/3401/CH8/EX8.4/Ex8_4.sce
@@ -0,0 +1,17 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Va=0.65 //V
+Js=4.15*10^-11// A/cm^2
+Nd=10**16 //cm**-3
+un=1350 //cm^2/Vs
+Jn=3.29 //A/cm^2
+
+J=Js*[exp(Va/(k*T))-1]
+disp(J,"current density in A/cm^2 is= ")
+
+E=Jn/(e*Nd*un)
+disp(E,"electric field in V/cm is= ")
+
diff --git a/3401/CH8/EX8.5/Ex8_5.sce b/3401/CH8/EX8.5/Ex8_5.sce
new file mode 100755
index 000000000..a59de2b1e
--- /dev/null
+++ b/3401/CH8/EX8.5/Ex8_5.sce
@@ -0,0 +1,14 @@
+clc
+
+T1=300 //K
+T2=310 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Va1=0.60 //V
+Va2=0.5827 //V
+E=1.12 //eV
+
+//J=exp(-Eg/(k*T))*exp((e*Va)/(k*T))
+
+a=([e*Va2*k*T1]+[e*Va1*k*T2])/((k*T1)-(k*T2))
+disp(a)
diff --git a/3401/CH8/EX8.6/Ex8_6.sce b/3401/CH8/EX8.6/Ex8_6.sce
new file mode 100755
index 000000000..f13c2a04c
--- /dev/null
+++ b/3401/CH8/EX8.6/Ex8_6.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Vt=0.0259 //v
+lp0=10^-3 //A
+taup0=10^-7 //s
+Idq=1*10**-3//A
+
+Cd=(1*lp0*taup0)/(2*Vt)
+disp(Cd,"diffusion capacitance in F is= ")
+
+vd=(Vt/Idq)
+disp(vd,"diffusion in ohm is= ")
diff --git a/3401/CH8/EX8.7/ex8_7.sce b/3401/CH8/EX8.7/ex8_7.sce
new file mode 100755
index 000000000..c6c994d59
--- /dev/null
+++ b/3401/CH8/EX8.7/ex8_7.sce
@@ -0,0 +1,18 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Na=10**16 //cm^-3
+Nd=10**16 //cm^-3
+ni=1.5*10**10 //cm^-3
+tau0=5*10^-7 //s
+eps=11.7
+//Vbr+Vr=z
+z=5 //V
+
+W=sqrt(((2*eps))*((Na+Nd)/Na*Nd)*z)
+disp(W,"depletion width in cm is= ")// textbook ans is wrong
+
+Jgen=(e*ni*W)/(2*tau0)
+disp(Jgen,"generation current density in A/cm^2 is= ")// textbook ans is wrong
diff --git a/3401/CH9/EX9.1/Ex9_1.sce b/3401/CH9/EX9.1/Ex9_1.sce
new file mode 100755
index 000000000..4808cf386
--- /dev/null
+++ b/3401/CH9/EX9.1/Ex9_1.sce
@@ -0,0 +1,27 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Nd=10**16 //cm^-3
+Nr=2.8*10**19 //cm^-3
+Vbi=0.33
+phim=4.55 //V
+psi=4.01 //V
+esp=11.7
+esp0=8.85*10^-14
+
+hb=phim-psi
+disp(hb,"barrier height in V is= ")
+
+phiu=(k*T)*log(Nr/Nd)
+disp(phiu,"schottky barrier height in V is= ")
+
+Vbt=hb-phiu
+disp(Vbt,"in V is=")
+
+xn=sqrt((2*esp*esp0*Vbi)/(e*Nd))
+disp(xn,"space charge in cm is= ")
+
+Emax=(e*Nd*xn)/(esp*esp0)
+disp(Emax,"maxi electric field in V/cm is= ")
diff --git a/3401/CH9/EX9.2/Ex9_2.sce b/3401/CH9/EX9.2/Ex9_2.sce
new file mode 100755
index 000000000..7858e956c
--- /dev/null
+++ b/3401/CH9/EX9.2/Ex9_2.sce
@@ -0,0 +1,22 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+esp=11.7
+esp0=8.85*10^-14
+Nd=2.7*10**17 //cm^-3
+Na=2.8*10**19 //cm^-3
+Vbi=0.40
+
+//(1/C)^2=2*(Vbi+Vr)/(e*esp*Nd)
+//delta(1/C)^2/deltaVr=z
+z=4.4*10**13
+Nd=2/(e*esp*esp0*z)
+disp(Nd,"semiconductor doping in cm^-3 is= ")
+
+phin=(k*T)*log(Na/Nd)
+disp(phin,"in v is=")
+
+phibn=Vbi+phin
+disp(phibn,"schottky barrier height in V is=")
diff --git a/3401/CH9/EX9.3/Ex9_3.sce b/3401/CH9/EX9.3/Ex9_3.sce
new file mode 100755
index 000000000..cc5df689b
--- /dev/null
+++ b/3401/CH9/EX9.3/Ex9_3.sce
@@ -0,0 +1,15 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+esp=11.7
+esp0=8.85*10^-14
+E=6.8*10**4 //V/cm
+
+deltaphi=sqrt((e*E)/(4*%pi*esp*esp0))
+disp(deltaphi,"schottky barrier lowring in V is =")
+
+xm=sqrt(e/(16*%pi*esp*esp0*E))
+disp(xm,"maxi barrier height in meter is=")
+
diff --git a/3401/CH9/EX9.4/Ex9_4.sce b/3401/CH9/EX9.4/Ex9_4.sce
new file mode 100755
index 000000000..edb448c6d
--- /dev/null
+++ b/3401/CH9/EX9.4/Ex9_4.sce
@@ -0,0 +1,13 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+esp=11.7
+esp0=8.85*10^-14
+phibe=0.67 //V
+Jst=6*10**-5 //A/cm^2
+
+A=exp(phibe/(k*T))*(Jst/T^2)
+disp(A,"Richardson constant in A/K^2-cm^2 is=") //textbook ans is different
+
diff --git a/3401/CH9/EX9.5/Ex9_5.sce b/3401/CH9/EX9.5/Ex9_5.sce
new file mode 100755
index 000000000..85b4ac766
--- /dev/null
+++ b/3401/CH9/EX9.5/Ex9_5.sce
@@ -0,0 +1,27 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+esp=11.7
+esp0=8.85*10^-14
+phibe=0.67 //V
+A=114 //A/K^2-cm^2
+Na=10^18 //cm^-3
+Nd=10^16 //cm^-3
+taup0=10^-7
+taun0=10^-7
+Dp=10 //cm^2/s
+Dn=25 //cm^2/s
+Lp=1.0*10**-3 //cm
+Ln=1.58*10**-3 //cm
+pn0=2.25*10**4 //cm^-3
+np0=2.25*10**2 //cm^-3
+
+
+a=k*T
+Jst=(A*T^2)*exp(-phibe/a)
+disp(Jst,"Richardson constant in A/K^2-cm^2 is=")
+
+Js=((e*Dn*np0)/Ln)+((e*Dp*pn0)/Lp)
+disp(Js,"reverse saturation current density in A/cm^2 is=")
diff --git a/3401/CH9/EX9.6/Ex9_6.sce b/3401/CH9/EX9.6/Ex9_6.sce
new file mode 100755
index 000000000..16fc51ab5
--- /dev/null
+++ b/3401/CH9/EX9.6/Ex9_6.sce
@@ -0,0 +1,19 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+esp=11.7
+esp0=8.85*10^-14
+J=10 //A/K^2-cm^2
+Jst=5.98*10^-5//A/K^2-cm^2
+Js=3.66*10^-11 //A/K^2-cm^2
+
+//J=Jst*[exp(e*Va/k*T)-1]
+//Va=(k*T/e)*log(J/Jst)
+Va=(k*T)*log(J/Jst)
+disp(Va,"forward bias voltage in V is= ")
+
+Va=(k*T)*log(J/Js)
+disp(Va,"forward bias voltage in V is= ")
+
diff --git a/3401/CH9/EX9.7/Ex9_7.sce b/3401/CH9/EX9.7/Ex9_7.sce
new file mode 100755
index 000000000..1ce4f4525
--- /dev/null
+++ b/3401/CH9/EX9.7/Ex9_7.sce
@@ -0,0 +1,14 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+Nd=7*10^18 //cm^3
+esp=11.7
+esp0=8.85*10^-14
+phibn=0.67 //V
+Vbi=0.67
+
+xn=sqrt((2*esp*esp0*Vbi)/(e*Nd))
+disp(xn,"space charge width in cm is= ")
+
diff --git a/3401/CH9/EX9.8/Ex9_8.sce b/3401/CH9/EX9.8/Ex9_8.sce
new file mode 100755
index 000000000..1b81213fd
--- /dev/null
+++ b/3401/CH9/EX9.8/Ex9_8.sce
@@ -0,0 +1,29 @@
+clc
+
+T=300 //K
+k=8.617*10**-5//eV/K
+e=1.6*10**-19 //C
+ni=2.4*10^13 //cm^-3
+Na=7*10^18 //cm^3
+Nd=10^16 //cm^-3
+pp0=6*10^18 //cm^-3
+esp=11.7
+esp0=8.85*10^-14
+xn=4.13
+xp=4.07
+
+//AE1=e*(xn-xp)
+AE=(xn-xp)
+disp(AE,"AE in eV is=")
+
+deltaE=1.43-0.67
+
+AE1=deltaE-AE
+disp(AE1,"AE1 in eV is= ")
+
+pn0=(ni^2)/Nd
+disp(pn0,"pn0 in cm^-3 is=")
+
+Vbi=AE1+(k*T)*log((Nd*pp0)/(pn0*Na))
+disp(Vbi,"Vbi in V is=")
+