summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/chapter1.ipynb2121
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/chapter2.ipynb1262
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/chapter5.ipynb64
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ1_chapter1.pngbin0 -> 18937 bytes
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ_ch1.pngbin0 -> 21564 bytes
-rw-r--r--Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_percentChangeinDiodeCurrent_chapter2.pngbin0 -> 34465 bytes
-rw-r--r--sample_notebooks/SUMITPRADHAN/chapter6.ipynb1632
7 files changed, 5079 insertions, 0 deletions
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter1.ipynb b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter1.ipynb
new file mode 100644
index 00000000..7de612d4
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter1.ipynb
@@ -0,0 +1,2121 @@
+{
+ "cells": [
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "# Chapter-1 Review of fundamentals of semiconductor"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.5.1 Pg 1-7"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 92,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Forbidden gap for Si : \n",
+ "at 35 degree C = 1.099 eV\n",
+ "at 60 degree C = 1.090 eV\n",
+ "\n",
+ "Forbidden gap for Ge : \n",
+ "at 35 degree C = 0.7163 eV\n",
+ "at 60 degree C = 0.7107 eV\n"
+ ]
+ }
+ ],
+ "source": [
+ "from __future__ import division\n",
+ "#Given : \n",
+ "t1=35##degreeC\n",
+ "t2=60##degreeC\n",
+ "T1=t1+273##K\n",
+ "T2=t2+273##K\n",
+ "print \"Forbidden gap for Si : \"\n",
+ "EG1_Si=1.21-3.6*10**-4*T1##eV\n",
+ "print \"at 35 degree C = %0.3f eV\"%EG1_Si\n",
+ "EG2_Si=1.21-3.6*10**-4*T2##eV\n",
+ "print \"at 60 degree C = %0.3f eV\"%EG2_Si\n",
+ "print \"\\nForbidden gap for Ge : \"\n",
+ "EG1_Ge=0.785-2.23*10**-4*T1##eV\n",
+ "print \"at 35 degree C = %0.4f eV\"%EG1_Ge\n",
+ "EG2_Ge=0.785-2.23*10**-4*T2##eV\n",
+ "print \"at 60 degree C = %0.4f eV\"%EG2_Ge"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.9.1 Pg 1-22"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 93,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "(i) Concentration of electron = 2.88e+21 m**3 \n",
+ "(ii) Drift velocity = 2.167 m/s \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=6*10**-2##m\n",
+ "V=1##Volt\n",
+ "A=10*10**-6##m**2\n",
+ "I=10*10**-3##A\n",
+ "q=1.602*10**-19##Coulomb\n",
+ "mu_n=1300*10**-4##m**2/V-s\n",
+ "E=V/l##V/m\n",
+ "v=mu_n*E##m/s\n",
+ "J=I/A##A/m**2\n",
+ "n=J/(q*mu_n*E)##per m**3\n",
+ "print \"(i) Concentration of electron = %0.2e m**3 \"%n\n",
+ "print \"(ii) Drift velocity = %0.3f m/s \"%v"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.9.2 Pg 1-23"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 94,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron mobility = 0.365 m**2/V-s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=6*10**-2##m\n",
+ "V=12##Volt\n",
+ "v=73##m/s\n",
+ "E=V/l##V/m\n",
+ "mu=v/E##m**2/V-s\n",
+ "print \"Electron mobility = %0.3f m**2/V-s\"%mu"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.10.1 Pg 1-25"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 95,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Concentration of electron = 9.615e+20 per m**3 \n",
+ "Electron velocity = 3.250 m/s \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=4*10**-2##m\n",
+ "A=10*10**-6##m**2\n",
+ "V=1##Volt\n",
+ "I=5*10**-3##A\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu=1300##cm**2/V-s\n",
+ "J=I/A##A/m**2\n",
+ "E=V/l##V/m\n",
+ "n=J/(q*mu*10**-4*E)#\n",
+ "v=mu*10**-4*E##m/s\n",
+ "print \"Concentration of electron = %0.3e per m**3 \"%n\n",
+ "print \"Electron velocity = %0.3f m/s \"%v"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.1 Pg 1-28"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 96,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Conductivity = 5.520e-04 (ohm-m)**-1 \n",
+ "Resistivity = 1811.6 ohm-m\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.5*10**10/10**-6##per m**3\n",
+ "mu_n=1800*10**-4##m**2/V-s\n",
+ "mu_p=500*10**-4##m**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_i=ni*(mu_n+mu_p)*q##(ohm-m)**-1\n",
+ "print \"Conductivity = %0.3e (ohm-m)**-1 \"%sigma_i\n",
+ "rho_i=1/sigma_i##ohm-m\n",
+ "print \"Resistivity = %0.1f ohm-m\"%rho_i"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.2 Pg 1-28"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 97,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Intrinsic carrier concentration = 1.5e+10 per cm**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "from math import sqrt, exp\n",
+ "#Given : \n",
+ "T=300##K\n",
+ "Ao=2.735*10**31##constant for Si\n",
+ "k=86*10**-6##boltzman constant\n",
+ "EGO=1.1##volt(Bandgap energy)\n",
+ "ni=sqrt(Ao*T**3*exp(-EGO/k/T))##per cm**3\n",
+ "print \"Intrinsic carrier concentration = %0.1e per cm**3\"%ni"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.3 Pg 1-29"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 98,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "(i) Current = 3.00 A\n",
+ "(ii) Conductivity = 2.78e+07 (ohm-m)**-1\n",
+ "(iii) velocity of free electrons = 2.08e-02 m/s\n",
+ "(iv) Mobility = 0.193 m**2/V-s \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "A=1*10**-6##m**2\n",
+ "R=3.6*10**-4/10**-2##ohm/m\n",
+ "n=9*10**26##electrons/m**3\n",
+ "J=3*10**6##A/m**2\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "I=J*A##A\n",
+ "print \"(i) Current = %0.2f A\"%I\n",
+ "rho=R*A##ohm-m\n",
+ "sigma=1/rho##(ohm-m)**-1\n",
+ "print \"(ii) Conductivity = %0.2e (ohm-m)**-1\"%sigma\n",
+ "v=J/n/q##m/s\n",
+ "print \"(iii) velocity of free electrons = %0.2e m/s\"%v\n",
+ "mu=sigma/n/q##m**2/V-s\n",
+ "print \"(iv) Mobility = %0.3f m**2/V-s \"%mu"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.4 Pg 1-31"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 99,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Intrinsic concentration at 30 degree C = 1.16e+16 per m**3) \n",
+ "Intrinsic concentration at 100 degree C = 1.221e+18 (per m**3)\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "rho=3*10**5*10**-2##ohm-m\n",
+ "T1=30+273##K\n",
+ "mu_n=0.13##m**2/V-s\n",
+ "mu_p=0.05##m**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "T2=100+273##K\n",
+ "sigma_i=1/rho##(ohm-m)**-1\n",
+ "ni1=sigma_i/q/(mu_n+mu_p)##electrons/m**3\n",
+ "print \"Intrinsic concentration at 30 degree C = %0.2e per m**3) \"%ni1\n",
+ "k=8.62*10**-5##eV/K(Boltzman constant)\n",
+ "EGO=1.21##V(Energy band gap)\n",
+ "Ao=ni1**2/(T1**3*exp(-EGO/k/T1))##constant\n",
+ "ni2=sqrt(Ao*T2**3*exp(-EGO/k/T2))##per cm**3\n",
+ "print \"Intrinsic concentration at 100 degree C = %0.3e (per m**3)\"%ni2"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.5 Pg 1-32"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 100,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Majority Carrier density = 3.720e+20 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=0.1*10**-2##m\n",
+ "R=1.5*10**3##ohm\n",
+ "mu_n=0.14##m**2/V-s\n",
+ "mu_p=0.05##m**2/V-s\n",
+ "A=8*10**-8##m**2\n",
+ "ni=1.5*10**10*10**6## per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "rho_n=R*A/l##ohm-m\n",
+ "sigma_n=1/rho_n##(ohm-m)**-1\n",
+ "ND=sigma_n/mu_n/q##\n",
+ "print \"Majority Carrier density = %0.3e per m**3\"%ND"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.6 Pg 1-32"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 101,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Length of the bar = 0.855 mm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "A=2.5*10**-4##m**2\n",
+ "n=1.5*10**16##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=0.14##m**2/V-s\n",
+ "mu_p=0.05##m**2/V-s\n",
+ "I=1.2*10**-3##A\n",
+ "V=9##Volts\n",
+ "ni=n## per m**3\n",
+ "sigma_i=ni*q*(mu_n+mu_p)##(ohm-m)**-1\n",
+ "rho_i=1/sigma_i##ohm-m\n",
+ "R=V/I##ohm\n",
+ "l=R*A/rho_i##m\n",
+ "print \"Length of the bar = %0.3f mm\"%(l*1000)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.7 Pg 1-34"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 102,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistivity = 2.3148e+05 ohm-cm\n",
+ "Ratio of donor impurity atom to Si atom : 1e-08\n"
+ ]
+ }
+ ],
+ "source": [
+ "from __future__ import division\n",
+ "#Given : \n",
+ "n=5*10**22##per cm**3\n",
+ "mu_n=1300##cm**2/V-s\n",
+ "mu_p=500##cm**2/V-s\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "T=300##K\n",
+ "rho_n=9.5##ohm-cm\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_i=ni*q*(mu_n+mu_p)##(ohm-cm)**-1\n",
+ "rho_i=1/sigma_i##ohm-cm\n",
+ "print \"Resistivity = %0.4e ohm-cm\"%rho_i\n",
+ "sigma_n=1/rho_n##(ohm-cm)**-1\n",
+ "ND=sigma_n/mu_n/q##per m**3\n",
+ "Ratio=ND/n#\n",
+ "print \"Ratio of donor impurity atom to Si atom : %0.e\"%(Ratio)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.8 Pg 1-35"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 103,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistivity of intrinsic Si = 2246.91 ohm-m \n",
+ "Resistivity of doped Si = 9.259e-02 ohm-m\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "n=5*10**22##per cm**3\n",
+ "ni=1.52*10**10*10**6##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=0.135##m**2/V-s\n",
+ "mu_p=0.048##m**2/V-s\n",
+ "impurity=1/10**8##atoms\n",
+ "sigma_i=ni*q*(mu_n+mu_p)##(ohm-cm)**-1\n",
+ "rho_i=1/sigma_i##ohm-cm\n",
+ "print \"Resistivity of intrinsic Si = %0.2f ohm-m \"%rho_i\n",
+ "ND=n*impurity*10**6##per m**3\n",
+ "sigma_n=ND*mu_n*q##(ohm-m)**-1\n",
+ "rho_n=1/sigma_n##ohm-m\n",
+ "print \"Resistivity of doped Si = %0.3e ohm-m\"%rho_n\n",
+ "#Answer in the book is not accurate."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.9 Pg 1-36"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 104,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Ratio of doner atom to Si atom per unit volume : 1e-08\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "rho=9.6*10**-2##ohm-m\n",
+ "mu_n=1300*10**-4##m**2/V-s\n",
+ "sigma_n=1/rho##(ohm-cm)**-1\n",
+ "TotalAtoms=5*10**28##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "ND=sigma_n/mu_n/q##per m**3\n",
+ "ratio=ND/TotalAtoms#\n",
+ "print \"Ratio of doner atom to Si atom per unit volume : %0.e\"%ratio"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.10 Pg 1-37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 105,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistivity of Ge = 45.0 ohm-cm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=2.5*10**13##per cm**3\n",
+ "mu_p=1800##cm**2/V-s\n",
+ "mu_n=3800##cm**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_i=ni*q*(mu_n+mu_p)##(ohm-cm)**-1\n",
+ "rho_i=1/sigma_i##ohm-cm\n",
+ "print \"Resistivity of Ge =\",round(rho_i),\"ohm-cm\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.11.11 Pg 1-37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 106,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron concentration = 1.440e+10 per m**3 \n",
+ "Conductivity of Si = 80 (ohm-m)**-1\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.2*10**16##per m**3\n",
+ "p=10**22##per m**3\n",
+ "mu_p=500*10**-4##cm**2/V-s\n",
+ "mu_n=1350*10**-4##cm**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "n=ni**2/p##per m**3\n",
+ "print \"Electron concentration = %0.3e per m**3 \"%n\n",
+ "sigma=q*(n*mu_n+p*mu_p)##(ohm-m)**-1\n",
+ "print \"Conductivity of Si = %0.f (ohm-m)**-1\"%sigma"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.12.1 Pg 1-39"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 107,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron concentration = 1e+17 per cm**3 \n",
+ "Holes = 2.25e+03 per cm**3 : \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "T=27+273##K\n",
+ "ND=10**17##per cm**3\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "n=ND##per m**3#ND>>n\n",
+ "print \"Electron concentration = %0.e per cm**3 \"%n\n",
+ "p=ni**2/n##per m**3\n",
+ "print \"Holes = %0.2e per cm**3 : \"%p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.12.2 Pg 1-39"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 108,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Total free electrons = 2.20e+17 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "Vol=4*50*1.5##mm**3\n",
+ "ni=2.4*10**19##per m**3\n",
+ "p=7.85*10**14##per m**3\n",
+ "n=ni**2/p##per m**3\n",
+ "Vol=Vol*10**-9##m**3\n",
+ "TotalElectron=n*Vol##no. of electrons\n",
+ "print \"Total free electrons = %0.2e per m**3\"%TotalElectron"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.13.1 Pg 1-41"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 109,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Total current density = 524.21 A/m**2 \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**14##per cm**3\n",
+ "NA=7*10**13##per cm**3\n",
+ "rho_i=60##ohm-cm\n",
+ "E=2##V/cm\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_p=1800##cm**2/V-s\n",
+ "mu_n=3800##cm**2/V-s\n",
+ "sigma_i=1/rho_i##(ohm-cm)**-1\n",
+ "ni=sigma_i/q/(mu_n+mu_p)##per cm**3\n",
+ "from sympy import symbols, solve\n",
+ "p = symbols('p')\n",
+ "n=p+(ND-NA)##per cm**3\n",
+ "#n*p=ni**2 \n",
+ "expr = n*p-ni**2 \n",
+ "#m=[1 (ND-NA) -ni**2]##polynomial\n",
+ "p=solve(expr,p)[1]##per m**3 #taking only +ve value\n",
+ "n=ni**2/p##per m**3\n",
+ "J=(n*mu_n+p*mu_p)*q*E/10**-4##A/m**2\n",
+ "print \"Total current density = %0.2f A/m**2 \"%J\n",
+ "#Answer in the textbook is not accurate."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.13.2 Pg 1-43"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 110,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Value of electrical field, E = 0.8150 V/cm \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**14##per cm**3\n",
+ "NA=7*10**3##per cm**3\n",
+ "rho_i=60##ohm-cm\n",
+ "J=52##mA/cm**2\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_p=1800##cm**2/V-s\n",
+ "mu_n=3800##cm**2/V-s\n",
+ "sigma_i=1/rho_i##(ohm-cm)**-1\n",
+ "ni=sigma_i/q/(mu_n+mu_p)##per cm**3\n",
+ "from sympy import symbols, solve\n",
+ "p = symbols('p')\n",
+ "n=p+(ND-NA)##per cm**3\n",
+ "#n*p=ni**2 \n",
+ "expr = n*p-ni**2 \n",
+ "#m=[1 (ND-NA) -ni**2]##polynomial\n",
+ "p=solve(expr,p)[1]##per m**3 #taking only +ve value\n",
+ "n=ni**2/p##per m**3\n",
+ "E=J*10**-3/q/(n*mu_n+p*mu_p)##V/m\n",
+ "print \"Value of electrical field, E = %0.4f V/cm \"%E"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.13.3 Pg 1-45"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 111,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Total current density = 382.75 A/m**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**14##per cm**3\n",
+ "NA=7*10**13##per cm**3\n",
+ "rho_i=60##ohm-cm\n",
+ "E=2##V/cm\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_p=500##cm**2/V-s\n",
+ "mu_n=1300##cm**2/V-s\n",
+ "sigma_i=1/rho_i##(ohm-cm)**-1\n",
+ "ni=sigma_i/q/(mu_n+mu_p)##per cm**3\n",
+ "from sympy import symbols, solve\n",
+ "p = symbols('p')\n",
+ "n=p+(ND-NA)##per cm**3\n",
+ "#n*p=ni**2 \n",
+ "expr = n*p-ni**2 \n",
+ "#m=[1 (ND-NA) -ni**2]##polynomial\n",
+ "p=solve(expr,p)[1]##per m**3 #taking only +ve value\n",
+ "n=ni**2/p##per m**3\n",
+ "J=(n*mu_n+p*mu_p)*q*E/10**-4##A/m**2\n",
+ "print \"Total current density = %0.2f A/m**2\"%J"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.13.4 Pg 1-46"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 112,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron mobility = 0.365 m**2/V-s \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=6*10**-2##m\n",
+ "V=12##volts\n",
+ "v=73##m/s\n",
+ "E=V/l##V/m\n",
+ "mu=v/E##m**2/V-s\n",
+ "print \"Electron mobility = %0.3f m**2/V-s \"%mu"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.15.1 Pg 1-54"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 113,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Magnitude of hall voltage = 65 mV \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**13##per cm**3\n",
+ "Bz=0.2##Wb/m**2\n",
+ "d=5##mm\n",
+ "E=5##V/cm\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=1300##cm**2/V-s\n",
+ "rho=ND*q##Coulomb/cm**3\n",
+ "J=rho*mu_n*E##A/cm**2\n",
+ "VH=Bz*10**-4*J*d*10**-1/rho##V\n",
+ "print \"Magnitude of hall voltage = %0.f mV \"%(VH*10**3)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.15.2 Pg 1-55"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 114,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Mobility = 0.050909 m**2/V-s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "rho=220*10**3*10**-2##ohm/m\n",
+ "d=2.2*10**-3##m\n",
+ "w=2*10**-3##m\n",
+ "B=0.1##Wb/m**2\n",
+ "I=5*10**-6##A\n",
+ "VH=28*10**-3##V\n",
+ "sigma=1/rho##(ohm-m)**-1\n",
+ "RH=VH*w/(B*I)##ohm\n",
+ "mu=sigma*RH##m**2/V-s\n",
+ "print \"Mobility = %0.6f m**2/V-s\"%mu"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.1 Pg 1-59"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 115,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Concentration of electron = 9.615e+20 per m**3\n",
+ "Electron velocity = 3.25 m/s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "l=4*10**-2##m\n",
+ "A=10*10**-6##m**2\n",
+ "V=1##Volt\n",
+ "I=5*10**-3##A\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu=1300##cm**2/V-s\n",
+ "J=I/A##A/m**2\n",
+ "E=V/l##V/m\n",
+ "n=J/(q*mu*10**-4*E)\n",
+ "v=mu*10**-4*E##m/s\n",
+ "print \"Concentration of electron = %0.3e per m**3\"%n\n",
+ "print \"Electron velocity = %0.2f m/s\"%v"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.2 Pg 1-59"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 116,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistivity of doped Ge = 3.738 ohm-cm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "mu_n=3800##cm**2/V-s\n",
+ "mu_p=1300##cm**2/V-s\n",
+ "ni=2.5*10**13##per cm**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "ND=4.4*10**22/10**8##per cm**3\n",
+ "sigma_n=ND*q*mu_n##(ohm-m)**-1\n",
+ "rho_n=1/sigma_n##ohm-cm\n",
+ "print \"Resistivity of doped Ge = %0.3f ohm-cm\"%rho_n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.3 Pg 1-60"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 117,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Minor carrier density = 4.5e+11 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.5*10**16##per m**3\n",
+ "n=5*10**20##per m**3\n",
+ "p=ni**2/n##per m**3\n",
+ "print \"Minor carrier density = %0.1e per m**3\"%p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.4 Pg 1-60"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 118,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron concentration, n = 8.00e+15 per cm**3\n",
+ "Hole concentration, p = 2.8125e+04 per cm**3 \n",
+ "Total electron concentration, nT = 8e+15 per cm**3\n",
+ "Total hole concentration, pT = 1.00e+14 per cm**3\n",
+ "Current, I = 36.45 mA\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.5*10**10##per cm**3\n",
+ "mu_n=1400##cm**2/V-s\n",
+ "mu_p=500##cm**2/V-s\n",
+ "l=1##cm\n",
+ "a=1##mm**2\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "del_n=10**14##per cm**3\n",
+ "del_p=10**14##per cm**3\n",
+ "Nd=8*10**15##per cm**3\n",
+ "n=Nd##per cm**3(Nd>>ni)\n",
+ "print \"Electron concentration, n = %0.2e per cm**3\"%n\n",
+ "p=ni**2/n##per m**3\n",
+ "print \"Hole concentration, p = %0.4e per cm**3 \"%p\n",
+ "nT=Nd+del_n##per cm**3\n",
+ "print \"Total electron concentration, nT = %0.e per cm**3\"%nT\n",
+ "pT=p+del_p##per cm**3\n",
+ "print \"Total hole concentration, pT = %0.2e per cm**3\"%pT\n",
+ "sigma=(nT*mu_n+pT*mu_p)*q##(ohm-cm)**-1\n",
+ "rho=1/sigma##ohm-cm\n",
+ "R=rho*l/(a*10**-2)##ohm\n",
+ "V=2##volt\n",
+ "I=V/R##A\n",
+ "print \"Current, I = %0.2f mA\"%(I*1000)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.5 Pg 1-61"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 119,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Current, I = 0.944 mA\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "A=2.3*10**-4##m**2\n",
+ "n=1.5*10**16##per m**3\n",
+ "l=1##mm\n",
+ "mu_n=1400##cm**2/V-s\n",
+ "mu_p=500##cm**2/V-s\n",
+ "p=n##per m**3\n",
+ "ni=n##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_i=ni*(mu_n*10**-4+mu_p*10**-4)*q##(ohm-m)**-1\n",
+ "rho_i=1/sigma_i##ohm-m\n",
+ "R=rho_i*l*10**-3/A##ohm\n",
+ "V=9##volt\n",
+ "I=V/R##A\n",
+ "print \"Current, I = %0.3f mA\"%(I*1000)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.6 Pg 1-62"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 120,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Concentration gradient, dn/dx = 1.624e+19\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**14##per m**3\n",
+ "Jn=10##mA/cm**2\n",
+ "E=3##V/cm\n",
+ "T=27+273##K\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=1500##cm**2/V-s\n",
+ "Dn=mu_n/39##Diffusion constant\n",
+ "n=ND##per m**3\n",
+ "dnBYdx=((Jn*10**-3/10**-4)-n*q*mu_n*E)/q/Dn##concentration gradient\n",
+ "print \"Concentration gradient, dn/dx = %0.3e\"%dnBYdx"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.7 Pg 1-63"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 121,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Total current density = 672.0 A/m**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**13##per m**3\n",
+ "NA=10**14##per m**3\n",
+ "rho_i=44##ohm-cm\n",
+ "E=3##V/cm\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=0.38##m**2/V-s\n",
+ "mu_p=0.18##m**2/V-s\n",
+ "ni=2.5*10**19##per m**3\n",
+ "from sympy import symbols, solve\n",
+ "p = symbols('p')\n",
+ "n=p+(ND-NA)##per cm**3\n",
+ "#n*p=ni**2 \n",
+ "expr = n*p-ni**2 \n",
+ "#m=[1 (ND-NA) -ni**2]##polynomial\n",
+ "p=solve(expr,p)[1]##per m**3 #taking only +ve value\n",
+ "n=ni**2/p##per m**3\n",
+ "J=(n*mu_n+p*mu_p)*q*(E/10**-2)##A/m**2\n",
+ "print \"Total current density = %0.1f A/m**2\"%J\n",
+ "#Ans in the textbook is not accurate."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.16.8 Pg 1-64"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 122,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Conductivity of intrinsic Ge = 2.24 (ohm-m)**-1\n",
+ "Conductivity after adding donor impurity = 267.52 (ohm-m)**-1\n",
+ "Conductivity after adding acceptor impurity = 126.72 (ohm-m)**-1 \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "T=300##K\n",
+ "ni=2.5*10**13##per cm**3\n",
+ "mu_n=3800##cm**2/V-s\n",
+ "mu_p=1800##cm**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_i=ni*(mu_n+mu_p)*q/10**-2##(ohm-m)**-1\n",
+ "print \"Conductivity of intrinsic Ge = %0.2f (ohm-m)**-1\"%sigma_i\n",
+ "ND=4.4*10**22/10**7##per cm**3\n",
+ "n=ND##per cm**3\n",
+ "sigma_n=n*mu_n*q/10**-2##(ohm-m)**-1\n",
+ "print \"Conductivity after adding donor impurity = %0.2f (ohm-m)**-1\"%sigma_n\n",
+ "NA=4.4*10**22/10**7##per cm**3\n",
+ "p=NA##per cm**3\n",
+ "sigma_p=p*mu_p*q/10**-2##(ohm-m)**-1\n",
+ "print \"Conductivity after adding acceptor impurity = %0.2f (ohm-m)**-1 \"%sigma_p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.1 Pg 1-102"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 123,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Equilibrium hole concentration = 2.25e+03 per cm**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=10**17##per cm**3\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "no=ND##per cm**3#/Nd>>ni\n",
+ "po=ni**2/no##per cm**3\n",
+ "print \"Equilibrium hole concentration = %0.2e per cm**3\"%po"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.3 Pg 1-103"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 124,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Fermi level, Ef = Ei + 0.407 eV\n"
+ ]
+ }
+ ],
+ "source": [
+ "from math import log\n",
+ "#Given : \n",
+ "ni=1.5*10**10##per cm**3\n",
+ "ND=10**17##per cm**3\n",
+ "no=ND##per cm**3#/Nd>>ni\n",
+ "po=ni**2/no##per cm**3\n",
+ "KT=0.0259##constant\n",
+ "delEf=KT*log(no/ni)##eV\n",
+ "print \"Fermi level, Ef = Ei +\",round(delEf,3),\"eV\""
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.4 Pg 1-104"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 125,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Diffusion coffiecients of electron = 4.40e-03 m**2/s\n",
+ "Diffusion coffiecients of holes = 6.47e-04 m**2/s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "K=1.38*10**-23##J/K\n",
+ "T=27+273##K\n",
+ "e=1.6*10**-19##constant\n",
+ "mu_n=0.17##m**2/V-s\n",
+ "mu_p=0.025##m**2/V-s\n",
+ "Dn=K*T/e*mu_n##m**2/s\n",
+ "print \"Diffusion coffiecients of electron = %0.2e m**2/s\"%Dn\n",
+ "Dp=K*T/e*mu_p##m**2/s\n",
+ "print \"Diffusion coffiecients of holes = %0.2e m**2/s\"%Dp"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.5 Pg 1-105"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 126,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Diffusion current density = 3.152e+03 A/m**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "K=1.38*10**-23##J/K\n",
+ "T=27+273##K\n",
+ "e=1.6*10**-19##constant\n",
+ "del_no=10**20##per.m**3\n",
+ "tau_n=10**-7##s\n",
+ "mu_n=0.15##m**2/V-s\n",
+ "Dn=K*T/e*mu_n##m**2/s\n",
+ "Ln=sqrt(Dn*tau_n)##m\n",
+ "Jn=e*Dn*del_no/Ln##A/m**2\n",
+ "print \"Diffusion current density = %0.3e A/m**2\"%Jn"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.6 Pg 1-105"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 127,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Concentration of holes = 4.680e+11 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "sigma_n=0.1##(ohm-cm)**-1\n",
+ "mu_n=1300##m**2/V-s\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "n_n=sigma_n/q/mu_n##per cm**3\n",
+ "p_n=ni**2/n_n##per cm**3\n",
+ "p_n=p_n*10**6##per m**3\n",
+ "print \"Concentration of holes = %0.3e per m**3\"%p_n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.7 Pg 1-106"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 128,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron transit time = 6.41e-09 s\n",
+ "Photoconductor gain : 216.0\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "L=100*10**-6##m\n",
+ "A=10**-7*10**-6##m**2\n",
+ "mu_e=0.13##m**2/V-s\n",
+ "mu_h=0.05##m**2/V-s\n",
+ "tau_h=10**-6##sec\n",
+ "V=12##volt\n",
+ "E=V/L##v/m\n",
+ "tn=L**2/(mu_e*V)##sec\n",
+ "print \"Electron transit time = %0.2e s\"%tn\n",
+ "Gain=tau_h/tn*(1+mu_h/mu_e)##\n",
+ "print \"Photoconductor gain :\",Gain"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.8 Pg 1-106"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 129,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistivity of intrinsic Ge at 300K = 45.00 ohm-cm\n",
+ "Resistivity of doped Ge = 3.74 ohm-cm \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "T=300##K\n",
+ "rho_i=45##ohm-cm\n",
+ "#part (i)\n",
+ "mu_n=3800##cm**2/V-s\n",
+ "mu_p=1800##cm**2/V-s\n",
+ "ni=2.5*10**13##per cm**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma=ni*q*(mu_n+mu_p)##(ohm-cm)**-1\n",
+ "rho=1/sigma##ohm-cm\n",
+ "print \"Resistivity of intrinsic Ge at 300K = %0.2f ohm-cm\"%round(rho)\n",
+ "#part (ii)\n",
+ "ND=4.4*10**22/10**8##per cm**3\n",
+ "sigma=ND*q*mu_n##(ohm-cm)**-1\n",
+ "rho=1/sigma##ohm-cm\n",
+ "print \"Resistivity of doped Ge = %0.2f ohm-cm \"%rho"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.9 Pg 1-107"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 130,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Electron concentration = 1e+22 per m**3\n",
+ "Electron concentration = 1e+10 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=10**16##per m**3\n",
+ "ND=10**22##per m**3\n",
+ "n=ND##per m**3#ND>>ni\n",
+ "print \"Electron concentration = %0.e per m**3\"%n\n",
+ "p=ni**2/n##per m**3\n",
+ "print \"Electron concentration = %0.e per m**3\"%p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.10 Pg 1-107"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 131,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Ratio of donor atom to Si atom : 1e-08\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "rho=9.6*10**-2##ohm-m\n",
+ "mu_n=1300##cm**2/V-s\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_n=1/rho##(ohm-m)**-1\n",
+ "ND=sigma_n/q/(mu_n*10**-4)##per m**3\n",
+ "ni=5*10**22*10**6##per m**3\n",
+ "print \"Ratio of donor atom to Si atom : %0.e\"%(ND/ni)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.11 Pg 1-108"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 132,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Equillibrium electron density = 2.25e+15 per cm**3\n",
+ "Equillibrium hole density = 1e+05 per cm**3 \n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.5*10**10##per cm**3\n",
+ "n_n=2.25*10**15##per cm**3\n",
+ "print \"Equillibrium electron density = %0.2e per cm**3\"%n_n\n",
+ "p_n=ni**2/n_n##per cm**3\n",
+ "print \"Equillibrium hole density = %0.e per cm**3 \"%p_n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.12 Pg 1-108"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 133,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Material is p-type & Carrier concentration = 1e+16 holes per cm**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "NA=2*10**16##per cm**3\n",
+ "ND=10**16##per cm**3\n",
+ "p=NA-ND##per cm**3\n",
+ "print \"Material is p-type & Carrier concentration = %0.e holes per cm**3\"%p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.13 Pg 1-108"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 134,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Rate of generation of minority carrier = 1e+20 electron hole pair/sec/cm**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "del_n=10**15##per cm**3\n",
+ "tau_p=10*10**-6##sec\n",
+ "rate=del_n/tau_p##rate of generation minority carrier\n",
+ "print \"Rate of generation of minority carrier = %0.e electron hole pair/sec/cm**3\"%rate"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.14 Pg 1-109"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 135,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Mobility = 5000 cm**2/V-s\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "E=10##V/cm\n",
+ "v=1/(20*10**-6)##m/s\n",
+ "mu=v/E##cm**2/V-s\n",
+ "print \"Mobility = %02.f cm**2/V-s\"%mu"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.15 Pg 1-109"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 136,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Hole diffucion current, Jp = -707.17 A/cm**2)\n",
+ "Electron diffucion current, Jp = 2552.08 A/cm**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ND=4.5*10**15 #per cm**3\n",
+ "A=1*10**-2 #cm**2\n",
+ "l=10 #cm\n",
+ "tau_p=1*10**-6 #sec\n",
+ "tau_n=1*10**-6 #sec\n",
+ "Dp=12 #cm**2/sec\n",
+ "Dn=30 #cm**2/sec\n",
+ "q=1.6*10**-19 #coulamb\n",
+ "del_p=10**21 #electron hole pair/cm**3/sec\n",
+ "x=34.6*10**-4 #cm\n",
+ "Kdash=26 #mV(Kdash is taken as K*T/q)\n",
+ "ni=1.5*10**10 #per cm**3\n",
+ "no=ND #per cm**3#ND<<ni\n",
+ "po=ni**2/no #per cm**3\n",
+ "ln=sqrt(Dn*tau_n) #cm\n",
+ "lp=sqrt(Dp*tau_p) #cm\n",
+ "dpBYdx=del_p*exp(-x/lp) #per cm**4\n",
+ "dnBYdx=del_p*exp(-x/ln) #per cm**4\n",
+ "Jp=-q*Dp*dpBYdx #A/cm**2\n",
+ "print \"Hole diffucion current, Jp = %0.2f A/cm**2)\"%Jp\n",
+ "Jn=q*Dn*dnBYdx #A/cm**2\n",
+ "print \"Electron diffucion current, Jp = %0.2f A/cm**2\"%Jn\n",
+ "#Answer is wrong in the book. Wrong calculation for dpBYdx and dnBYdx."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.16 Pg 1-111"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 137,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Energy band gap = 2.26 eV\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "h=6.626*10**-34##J-s\n",
+ "lamda=5490##Angstrum\n",
+ "c=3*10**8##m/s(speed of light)\n",
+ "f=c/(lamda*10**-10)##Hz\n",
+ "E=(h/1.6/10**-19)*f##eV\n",
+ "print \"Energy band gap = %0.2f eV\"%E"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.17 Pg 1-111"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 138,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEZCAYAAACEkhK6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xu0HGWd7vHvQwLI/SIKCAxxNDKAXKIQMM5gO3Bg60Fw\nvDEBQRxG9zoOXs4hmAWDi43juGQtdUSZmUO4CKNHwlGChovcMrREHQkbkhDYSQ5RNiZcVBQwiVwS\n8jt/VHXSdPburu7d1dfns1avdFW/VfV2r85++6n3fasUEZiZmdWyTbsrYGZm3cENhpmZZeIGw8zM\nMnGDYWZmmbjBMDOzTNxgmJlZJm4wzMwsEzcYZg2SdKykdZJ2GuO1xZI+1Y56meXFDYZZgyLiF8Aa\n4EPl6yW9FTgYuL4d9TLLixsMswqSRiWdJ2mppOckzZW0/TjFrwPOqlh3FnBrRDwr6TWSvivpGUnP\nSlok6fVjHPNNkn4vaVq6/AZJv5N0XHPfnVnjurrBkHSNpN9IWpah7HGSHpS0QdIHy9a/Oz19UHq8\nIOmUfGtuHS6ADwMnAW8EDgfOHqfsd4HjJO0PIGkbYCZJQwLwMWBXYH9gT2AQeGGrA0b8EpgNfFfS\nDsC3gW9HxL3NeUtmE9fVDQbJf6qBjGUfJ/nP+73ylRFxT0RMi4hpwF8DfwLubGotrRt9MyKejohn\ngZuBI8cqFBGrgSJwZrrqeGB74NZ0+WXgtcDUSCyOiLXj7OsqYBWwCNgb+McmvRezpujqBiMiFgLP\nlq9Lo/2PJQ1LulfSQWnZxyNiGbCpyi4/DNwWES/mV2vrEk+XPX8B2LlK2evY0mCcCVwfEa+ky98B\n7gDmSnpC0qWSJlfZ11XAocC3ImJDY1U3y0dXNxjjmAN8OiKOAs4H/q2Obf8Wd1Ta1mpd0vkmYH9J\n7wb+hi2no4iIjRHxxYg4FJgBnMzWfR4ASNoZ+AZJo3GJpD2aUXmzZqn2S6frpP/h3gF8X1Jp9XYZ\nt90XeCvJr0Gzcqr2YkSsl/QDklOkoxHx4OYNpQLwe2AEWAtsAF4Zaz/AZcCiiPikpCuA/w2cNvHq\nmzVHTzUYJInpubQ/opqxfjF+BJhXdirBrCSonTKuI+kYn12xfh+SP/z7A+uAuSSnqV5F0qnAicBh\n6ar/BSyRNDMinHqtIyjPGyhJugb478BvI+Kwccp8E3gPSWfz2RGxOF0/QBLPJwFXRcSl42w/Bbi5\ntH9JPwP+JSJ+oCRmHBYRD5WVvzYtf2PFfn4BzI6InzT+js3MelfefRhVRzFJei/w5oiYCnwS+Pd0\n/STg8nTbQ4CZkg4eY/vrgZ8DB0laLenjwBnAOZKWAA8Dp6Rlj5a0mmSS1RXlQ3HTRmc/NxZmZuPL\n9ZRURCxM/xiP5xTSDsKIuE/S7pL2IRn7vioiRgEkzQVOBZZX7H/mOPt9zxh1uR84YJx6jo73mpmZ\nJdo9Smo/YHXZ8pp03RvGWW9mZm3S7gYDaoxAMTOzztDuUVJP8OpTQfuTpIltK9YfkK5/FUn59dib\nmfWwiKj7x3q7E8Z80klMko4lGRL7G2AYmCppiqTtSMaizx9rB6OjwQknBEcfHTzySBDhR6OPiy++\nuO116KWHP09/lp36aFSuDcYYo5j+TtKgpEGAiLgN+JWkVcAVwKfS9RuBc0km0Y0AN0TE8rGOceCB\ncOedcM458K53wVe+Ahs35vmuzMz6U96jpMYbxVRe5txx1v8Y+HGW40gwOAgDA/D3fw/z5sG118Ih\nh9RXXzMzG1+7T0k1ldPGxBQKhXZXoaf482wef5adIdeZ3nmTFOPV//HHk7Tx/PNOG2Zm5SQRXdjp\nnRunDTOz5urZhFHOacPMbAsnjCqcNszMJq4vEkY5pw0z63dOGBk5bZiZNabvEkY5pw0z60dOGA1w\n2jAzy66vE0Y5pw0z6xdOGBPktGFmVp0TxhicNsyslzlhNJHThpnZ1pwwanDaMLNe44SRE6cNM7OE\nE0YdnDbMrBc4YbSA04aZ9TMnjAY5bZhZt3LCaDGnDTPrN7k2GJIGJK2Q9Kik2WO8voekmyQtlXSf\npEPLXhuV9JCkxZIW5VnPRpXuJT48DAsWwIwZMDLS7lqZmeUjtwZD0iTgcmAAOASYKengimIXAg9G\nxBHAWcBlZa8FUIiIaRExPa96NoPThpn1gzwTxnRgVUSMRsQGYC5wakWZg4F7ACJiJTBF0uvKXq/7\nHFu7OG2YWa/Ls8HYD1hdtrwmXVduKfABAEnTgQOB/dPXArhb0rCkT+RYz6Zy2jCzXjU5x31nGb70\nFeAySYuBZcBi4JX0tb+MiCfTxHGXpBURsbByB0NDQ5ufFwoFCoXCROs9YaW0MTCQjKSaN88jqcys\nfYrFIsViccL7yW1YraRjgaGIGEiXLwA2RcSlVbZ5DDgsItZVrL8YWBcRX6tY37ZhtVlFwJw5cNFF\ncN55MGsWTM6zmTYzq6ETh9UOA1MlTZG0HXAaML+8gKTd0tdITzv9JCLWSdpR0i7p+p2AE0kSSNdx\n34aZ9YrcGoyI2AicC9wBjAA3RMRySYOSBtNihwDLJK0ATgI+m67fG1goaQlwH3BLRNyZV11bwX0b\nZtbtPNO7DTxL3MzaqRNPSdk4nDbMrBs5YbSZ04aZtZoTRpdy2jCzbuGE0UGcNsysFZwweoDThpl1\nMieMDuW0YWZ5ccLoMU4bZtZpnDC6gNOGmTWTE0YPc9ows07ghNFlnDbMbKKcMPqE04aZtYsTRhdz\n2jCzRjhh9CGnDTNrJSeMHuG0YWZZOWH0OacNM8ubE0YPctows2qcMGwzpw0zy4MTRo9z2jCzSk4Y\nNianDTNrllwbDEkDklZIelTS7DFe30PSTZKWSrpP0qFZt7XsJBgchOFhWLAAZsyAkZF218rMuk1u\nDYakScDlwABwCDBT0sEVxS4EHoyII4CzgMvq2Nbq5LRhZhORZ8KYDqyKiNGI2ADMBU6tKHMwcA9A\nRKwEpkh6fcZtrQFOG2bWqDwbjP2A1WXLa9J15ZYCHwCQNB04ENg/47Y2AU4bZlavyTnuO8vwpa8A\nl0laDCwDFgOvZNwWgKGhoc3PC4UChUKhrkr2s1LaGBhIRlLNm+eRVGa9qFgsUiwWJ7yfmsNq047o\n44ApJH/IR4GFEfFIje2OBYYiYiBdvgDYFBGXVtnmMeAw4K1ZtvWw2uaJgDlz4KKL4LzzYNYsmJzn\nzwkza5umD6uVdKakRcBXgX2AX5E0FvsCX5V0v6SPVtn3MDBV0hRJ2wGnAfMrjrFb+hqSPgH8JCLW\nZdnWmst9G2ZWS7XfkHsAx0fE2rFelLQrcPZ4G0fERknnAncAk4CrI2K5pMH09StIRkBdKymAh4Fz\nqm1b75uz+pX6NubMSfo2nDbMrMQzvW1cniVu1pvyOCW1raSPSir1I3xM0uWSzpHkGeJ9wCOpzKzc\nuAlD0tXAbsB2wAvA9sCNwMnAryPi/FZVcjxOGK3jtGHWOxpNGNUajEci4lBJ2wK/AfaNiJckTSaZ\nnX34xKo8cW4wWssjqcx6Qx4XH9wAkM60vj8iXkqXN1LHPAnrHR5JZdbfqjUYT0vaBSAiTiqtlLQv\n8FLeFbPO5b4Ns/5U9ygpSTsBO0XEb/OpUl118SmpNnPfhln3yeV+GJI+l/77mdK6iFjfCY2FdQan\nDbP+UWt47FpJ5wPrW1EZ607u2zDrD9XmYVwMvAX4J+CgdNlsXE4bZr2tah+GpC8BdwInRsRFLatV\nRu7D6Fzu2zDrXHnd0/veiLgXKDZUK+tbThtmvadWg3FIxb9mmblvw6y3uNPbcue0YdYbql0a5GJg\nR+CzwDeB9RFxSQvrVpP7MLqP+zbM2q/pfRhp47ABOBF4udMaC+tOThtm3cud3tZy7tsw606ZLg0i\n6QiSe3qXrk0aETEvx3pl4lNS3c9XwDVrvaZf3rxsx98GDgMeATaV1kfEx+s9WLO5wegd7tswa508\nG4wR4NBO/MvsBqO3OG2YtUZeE/cA7qfBeRiSBiStkPSopNljvL6XpNslLZH0sKSzy14blfSQpMWS\nFjVyfOsu7tsw62xZEkYBmA88zZb7YEStO+5JmgSsBE4AniBpeGZGxPKyMkPA9hFxgaS90vJ7R8RG\nSY8Bb4+IP1Q5hhNGj3LaMMtPngnjauCjwADwvvRxSobtpgOrImI0vWvfXODUijJPAbumz3cFfp/e\n0a+k7jdkvcFpw6zzZGkwfhsR8yPiV+kf/9GIGM2w3X7A6rLlNem6clcCh0p6ElhKMkmwJIC7JQ1L\n+kSG41kP8rwNs86RJeQvlvQ94Gbg5XRdlmG1Wc4VXQgsiYiCpDcBd0k6IiLWAu+MiKckvS5dvyIi\nFlbuYGhoaPPzQqFAoVDIcFjrJqW0MTCQjKSaN88jqczqUSwWKRaLE95Plj6Maxnjj3+tYbWSjgWG\nImIgXb4A2BQRl5aVuQ3454j4Wbq8AJgdEcMV+7oYWBcRX6tY7z6MPuO+DbOJy21YbaMkTSbpxD4e\neBJYxNad3l8Hno+ISyTtDTwAHA68CEyKiLXpPcTvBC6JiDsrjuEGo0953oZZ43Lr9JZ0naTdy5b3\nkHRNre3SzutzgTuAEeCGiFguaVDSYFrsy8BRkpYCdwOfT0dF7QMslLQEuA+4pbKxsP7mvg2z1sty\nSmpJRBxZa107OGEYOG2Y1SvPYbWStGfZwp7ApHoPZJYXpw2z1siSMM4C/hH4vyTzIj5M0lH9H/lX\nrzonDKvktGFWW66d3pIOBf6aZLTUf0ZER0yhcoNhY/FIKrPqOm6UVCu4wbBqnDbMxpZnH4ZZV3Lf\nhllzOWFYX3DaMNsil4QhabKkexqvlllncNowm7iqDUY6+W5T+cQ9s27lK+CaTUyWsSPrgWWS7kqf\nQ3Lxwc/kVy2z/JTSxpw5SdrwSCqzbLLMwzg7fVoqKJIG47oc65WJ+zBsoty3Yf0o73kYOwJ/FhEr\nGqlcXtxgWDN43ob1mzwvPngKsBi4PV2eJml+/VU060zu2zDLJss8jCHgGOBZgIhYDPx5jnUyawuP\npDKrLkuDsSEinqtYtymPypi1m9OG2fiyNBiPSDoDmCxpqqRvAT/PuV5mbeW0Yba1LKOkdiK5Wu2J\n6ao7gH+KiBdzrltN7vS2VvBIKus1uV98UNJuJMNp/1jvQfLiBsNaxSOprJfk1mBIOhq4Btg1XfUc\ncE5EDNddyyZzg2Gt5rRhvSDPq9VeA3wqIg6MiAOBf0jXmfUd921YP8vSYGyMiIWlhYj4KZDpv4ik\nAUkrJD0qafYYr+8l6XZJSyQ9XDarvOa2Zu3ikVTWr7KckvoGsANwfbrqNOBF4DsAEfHgONtNAlYC\nJwBPAPcDMyNieVmZIWD7iLhA0l5p+b1JLkNSddt0e5+SsrZy34Z1ozxPSR0JvAW4OH38Rbrua+lj\nPNOBVRExGhEbgLnAqRVlnmJL38iuwO/TK+Rm2das7Zw2rJ/U/C0UEYUG970fsLpseQ3JjPFyVwL/\nKelJYBfgI3Vsa9YxfAVc6wd5fp2znCu6EFgSEQVJbwLuknREPQcZGhra/LxQKFAoFOrZ3KxpSmlj\nYCAZSTVvnkdSWWcoFosUi8UJ7ye3W7RKOhYYioiBdPkCYFNEXFpW5jbgnyPiZ+nyAmA2SUNWddt0\nvfswrCO5b8M6WZ59GI0aBqZKmiJpO5LO8sqr3K4g6dhG0t7AQcCvMm5r1rHct2G9KMvlzT8iadf0\n+Rck3STpbbW2SzuvzyW5lMgIcENELJc0KGkwLfZl4ChJS4G7gc9HxB/G27aRN2jWTp63Yb0ky7Da\nZRFxmKS/BL4EfBX4QkS0vRPap6Ssm3iWuHWKPE9JvZL+ezJwZUTcAmxX74HM+p3ThnW7LAnjVpLJ\nc/8NmEYyae++iKhrNFMenDCsWzltWDvlmTA+QtKXcGJ6I6U9gPPrPZCZbeG0Yd0o07BaSZOB11M2\nbyMifp1jvTJxwrBe4LRhrZZbwpD0aeBpklFMt5Y9zKwJnDasW2Tpw/glMD0ift+aKmXnhGG9xmnD\nWiHPPoxfAx1zlz2zXua0YZ0sS8K4huRqtbcCL6erIyK+nnPdanLCsF7mtGF5yTth3E0y92JnkqvK\n7lLvgcysPk4b1mkyX3xQ0i4AEbE21xrVwQnD+oXThjVTnqOkDpO0GHgEeETSA5Le2kglzawxThvW\nCbL0YfwXcGFE3JMuF4AvR8SM/KtXnROG9SOnDZuoPPswdiw1FgARUQR2qvdAZtYcThvWLlkSxg+B\nB4DvAALOAN4eEX+Tf/Wqc8Kwfue0YY3IM2F8nOSyIPOAG4HXAX9X74HMrPmcNqyVqiaM9BpSd0XE\nu1tXpeycMMy2cNqwrHJJGOmd7zZJ2r3hmplZSzhtWN6y9GHMJ7kPxl3A+nR1RMRncq5bTU4YZmNz\n2rBqGk0YWRqMj5F0dpcKiqTBuK7uWjaZGwyz8UXAnDlw0UVw3nkwaxZMnlx7O+t9eXZ67xER10bE\ndenjWpKbKGWp1ICkFZIelTR7jNdnSVqcPpZJ2lg6/SVpVNJD6WuL6npXZoYEg4MwPAwLFsCMGTAy\n0u5aWTfLkjAWR8S0inVLIuLIGttNAlYCJ5Dc4vV+YGZELB+n/MnA5yLihHT5MZLhu3+ocgwnDLMM\nnDasXNMThqSZkm4G3ijp5rJHEchyb4zpwKqIGI2IDcBc4NQq5U8Hrq+sRobjmFkNThvWDNV+Y/wc\neIpk3sVX2fLH+4/AQxn2vR+wumx5DXDMWAUl7QicBHyqbHUAd0t6BbgiIq7McEwzq6I0kmrOnGQk\nldOG1WPcr0lEPA48Lul04KmIeAFA0g7A/sBojX3Xc67ofcBPI+K5snXvjIinJL0OuEvSiohYWLnh\n0NDQ5ueFQoFCoVDHYc36TyltDAwkI6nmzfNIql5XLBYpFosT3k+WPoxhYEZEvJwubw/8LCKOqrHd\nscBQRAykyxcAmyLi0jHK3gTcEBFzx9nXxcC6iPhaxXr3YZhNgPs2+lOeo6QmlxoLgIh4Cdg2w3bD\nwFRJUyRtB5wGzK8sJGk34DjgR2Xrdizdf0PSTsCJwLIMxzSzOrhvw+qRpcF4RtLmzur0+TO1Nkpn\niZ8L3AGMkCSI5ZIGJQ2WFX0/cEfplFdqb2ChpCXAfcAtEXFnhrqaWQM8S9yyyHJK6s3A/wHekK5a\nA5wZEatyrltNPiVl1nyeJd77cpvpXXaAndPyvkWrWY9z30Zvy/MWrftIuhr4QUSslXSIpHMaqqWZ\ndQX3bdhYsvRhXAvcyZZTUo8C/zOvCplZ53DfhpXL0mDsFRE3AK8ApLO2/ZUx6xNOG1aSpcFYJ+m1\npYV0fsXz+VXJzDqR04ZlGSX1duBbwKHAIySXCvlQRCzNv3rVudPbrD08kqq75TpKStK2wEHp4sr0\ntFTbucEwax+PpOpeTW8wJH2Q5HpQKvuX9DkRMa+xqjaPGwyz9nPa6D55NBjXUuUCghHx8XoP1mxu\nMMw6g9NGd8l94l4ncoNh1lmcNrpD7hP3JN2eLnvinpmNySOpepsn7plZU3neRu/yxD0zy4XTRu/x\nxD0zy43TRm/J0mCcB9wM/LmknwPfAT6Ta63MrKc4bfSGeifuiWTi3ss1NmkJj5Iy6z4eSdV+TR8l\nJalQeh4RGyLi4YhYVt5YSHp33TU1s77mtNG9qk3c+yrJvbbvJrk/91MkDcw+wFHACcA9EfH51lR1\nzDo6YZh1MaeN9shl4p6kXYBTgXcCB6arHwd+CvwoItY1UNemcYNh1v08S7z1OnKmt6QB4BvAJOCq\niLi04vVZwBnp4mTgYJJhvM/V2jbd3g2GWY9w2mid3BoMSa8BPghMIfnjLSAi4os1tpsErCQ5dfUE\ncD8wMyKWj1P+ZOBzEXFC1m3dYJj1FqeN1sjt0iDAj4BTgA3AemBd+m8t04FVETGaTvabS3J6azyn\nA9c3uK2Z9QDP2+hsWdru/SLipAb2vR+wumx5DXDMWAUl7QicBHyq3m3NrPeURlLNmZOMpHLa6AxZ\nPv6fSzo8Ih6qc9/1nCt6H/DTiHiu3m2HhoY2Py8UChQKhToOa2adqpQ2BgaSvo1589y30ahisUix\nWJzwfrL0YSwH3gw8BryUro6IOLzGdscCQxExkC5fAGwap/P6JuCGiJhbz7buwzDrD+7baK48O72n\njLU+IkZrbDeZpOP6eOBJYBFjd1zvBvwK2D8iXqhzWzcYZn3EI6maI7dO77TjeatHhu02AucCdwAj\nJAliuaRBSYNlRd8P3FFqLKptW9c7M7Oe41ni7eU77plZV3LaaFyew2rNzDqO00brOWGYWddz2qiP\nE4aZ9S2njdZwwjCznuK0UZsThpkZTht5csIws57ltDE2JwwzswpOG83lhGFmfcFpYwsnDDOzKpw2\nJs4Jw8z6Tr+nDScMM7OMnDYa44RhZn2tH9OGE4aZWQOcNrJzwjAzS/VL2nDCMDObIKeN6pwwzMzG\n0MtpwwnDzKyJnDa25oRhZlZDr6WNjkwYkgYkrZD0qKTZ45QpSFos6WFJxbL1o5IeSl9blGc9zcyq\ncdpI5JYwJE0CVgInAE8A9wMzI2J5WZndgZ8BJ0XEGkl7RcQz6WuPAW+PiD9UOYYThpm1VC+kjU5M\nGNOBVRExGhEbgLnAqRVlTgdujIg1AKXGokzdb8jMLE/9nDbybDD2A1aXLa9J15WbCuwp6R5Jw5LO\nLHstgLvT9Z/IsZ5mZnWRYHAQhodhwQKYMQNGRtpdq/zl2WBkOVe0LfA24L3AScAXJE1NX/vLiJgG\nvAf4B0l/lU81zcwa029pY3KO+34COKBs+QCSlFFuNfBMRLwAvCDpXuAI4NGIeBIgIn4n6SaSU1wL\nKw8yNDS0+XmhUKBQKDTxLZiZVVdKGwMDSd/GvHmd17dRLBYpFosT3k+end6TSTq9jweeBBaxdaf3\nXwCXk6SL7YH7gNOAUWBSRKyVtBNwJ3BJRNxZcQx3eptZx4iAOXPgoovgvPNg1iyYnOfP8gZ1XKd3\nRGwEzgXuAEaAGyJiuaRBSYNpmRXA7cBDJI3FlRExAuwDLJS0JF1/S2VjYWbWaXq9b8MT98zMctDJ\naaPRhOEGw8wsR504b6PjTkmZmVlvjaRywjAza5FOSRtOGGZmHa7b04YThplZG7QzbThhmJl1kW5M\nG04YZmZt1uq04YRhZtaluiVtOGGYmXWQVqQNJwwzsx7QyWnDCcPMrEPllTacMMzMekynpQ0nDDOz\nLtDMtOGEYWbWwzohbThhmJl1mYmmDScMM7M+0a604YRhZtbFGkkbThhmZn2olWnDCcPMrEdkTRsd\nmTAkDUhaIelRSbPHKVOQtFjSw5KK9WxrZmZb5J02cmswJE0CLgcGgEOAmZIOriizO/CvwPsi4q3A\nh7Jua81XLBbbXYWe4s+zefxZZifB4CAMD8OCBTBjBoyMNGffeSaM6cCqiBiNiA3AXODUijKnAzdG\nxBqAiHimjm2tyfyfsrn8eTaPP8v65ZE28mww9gNWly2vSdeVmwrsKekeScOSzqxjWzMzq6LZaWNy\n86q2lSy90dsCbwOOB3YE/kvSLzJua2ZmGZTSxpw5SdpoVG6jpCQdCwxFxEC6fAGwKSIuLSszG9gh\nIobS5auA20kSRdVt0/VuWMzMGtDIKKk8E8YwMFXSFOBJ4DRgZkWZHwGXp53c2wPHAF8H/l+GbRt6\nw2Zm1pjcGoyI2CjpXOAOYBJwdUQslzSYvn5FRKyQdDvwELAJuDIiRgDG2javupqZWW1dPXHPzMxa\npysuDZJxAuA309eXSprW6jp2k1qfZzqZ8vl0QuViSRe1o57dQNI1kn4jaVmVMv5uZlDrs/T3sj6S\nDkhHoD6SToz+zDjlsn8/I6KjHySnpFYBU0hGVS0BDq4o817gtvT5McAv2l3vTn1k/DwLwPx217Ub\nHsBfAdOAZeO87u9m8z5Lfy/r+zz3AY5Mn+8MrJzo385uSBhZJvGdAlwHEBH3AbtL2ru11ewaWSdF\nekBBBhGxEHi2ShF/NzPK8FmCv5eZRcTTEbEkfb4OWA68oaJYXd/PbmgwskziG6vM/jnXq1tl+TwD\nmJFG1NskNenW833J383m8feyQemI02nAfRUv1fX9zHNYbbNk7ZWv/OXh3vyxZflcHgQOiIg/SXoP\n8EPgLflWq6f5u9kc/l42QNLOwA+Az6ZJY6siFcvjfj+7IWE8ARxQtnwASStYrcz+6TrbWs3PMyLW\nRsSf0uc/BraVtGfrqthT/N1sEn8v6ydpW+BG4LsR8cMxitT1/eyGBmPzBEBJ25FM4ptfUWY+cBZs\nnmH+XET8prXV7Bo1P09Je0tS+nw6yfDrP7S+qj3B380m8feyPulndTUwEhHfGKdYXd/Pjj8lFdkm\nAN4m6b2SVgHrgY+3scodLcvnSXKZ+f8haSPwJ+Bv21bhDifpeuBdwF6SVgMXk4w+83ezTrU+S/y9\nrNc7gY8CD0lanK67EPgzaOz76Yl7ZmaWSTeckjIzsw7gBsPMzDJxg2FmZpm4wTAzs0zcYJiZWSZu\nMMzMLBM3GNbzJL1e0q11bnOJpOPzqlMjJJ0q6eAGtjtF0hfyqJP1F8/DsJ4n6Yskl8z+fk77nxwR\nG/PYd8VxrgVujogb69hmEsndLBcDR6dXKDZriBsM6wmSjgauIrl8+2SSq3J+JCJGJI0AR6UXrTsb\neD+wIzAV+BrwGuB04CXgvRHxbPkf53Tf3wB2Al4ETiCZdfyBdN026fNvA28kmYX8yYh41Y2A6jj2\nm4DLgdel+/oE8FrgZuD59PGB9LivKhcRK9O6vwgcCfw0ImZJ+nfgloioK2mZlev4S4OYZRER90ua\nD3wJ2AH4TtpY7AO8UrpoXepQkj+mOwC/BM6PiLdJ+jrJdXUuI7liZ6TX25pL0vg8kF7584V0P9OA\nwyLiOUnfAh6IiPdLejfwH+nrlbIcew4wGBGrJB0D/FtEHJ++v5sjYh6ApAWV5YDSabQ3AO+ILb8I\nFwHHAW4wrGFuMKyXfJHk4oovAJ9O1x0IPFVWJoB7ImI9sF7ScyS/3AGWAYeXlRVwEPBURDwAm29E\ng6QA7opHUIWzAAABlUlEQVSI59Ky7yT51U9E3CPptZJ2rricdM1jS9oJmAF8P73OHsB2FXUqXbL6\nHeOUC+D7ZY0FwJPAwNYfmVl2bjCsl+xFcopoEskv+FKqqLze/0tlzzeVLW9i6/8T1c7Zrq9YznI3\nuFrH3gZ4NiLGu7dyqT7bkFxZdLxyf6pY3gbfh8MmyKOkrJdcAVwEfA+4NF33OMm9jUuq/VEf60Yy\nK4F9JR0FIGmXtCO5suxC4Iy0TAH43Rg3q6l57IhYCzwm6UPpviSplHrWArum5f5YpdxY9iX5LMwa\n5gbDeoKks4CXImIu8BXgaEmFiHgamCxpx7Ro8Opf2pXPX/UrPB1VdBrwLUlLSC4L/5oxyg4Bb5e0\nFPgy8LExqpn12GcA56THe5jkvsuQ9KWcL+kBSW+sUq5y35AMBrh3jDqZZeZRUtbzJA0ByyPihnbX\npR0kbUNye9OjWjH813qXE4b1g39l7F/8/eJk4AduLGyinDDMzCwTJwwzM8vEDYaZmWXiBsPMzDJx\ng2FmZpm4wTAzs0zcYJiZWSb/H7crH3T5904sAAAAAElFTkSuQmCC\n",
+ "text/plain": [
+ "<matplotlib.figure.Figure at 0x7f5da1c9e290>"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Current density = 1120.00 A/cm**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "q=1.6*10**-19##Coulomb\n",
+ "Dn=35##cm**2/s\n",
+ "x=[0, 2]##micro meter\n",
+ "n=[10**17 ,6*10**16]##per cm**3\n",
+ "%matplotlib inline\n",
+ "from matplotlib.pyplot import plot, title, show, xlabel, ylabel\n",
+ "plot(x,n)#\n",
+ "title('n Vs x')#\n",
+ "xlabel('x(micro meter)')#\n",
+ "ylabel('n(electrons per cm**3)')#\n",
+ "show()\n",
+ "dnBYdx=(n[1]-n[0])/(x[0]-x[1])/10**-4##gradient\n",
+ "Jn=q*Dn*dnBYdx##A/cm**2\n",
+ "print \"Current density = %0.2f A/cm**2\"%Jn"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.18 Pg 1-112"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 139,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Resistance of the bar = 1 Mohm\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "q=1.6*10**-19##Coulomb\n",
+ "l=0.1##cm\n",
+ "A=100*10**-8##cm**2\n",
+ "n_n=5*10**20*10**-6##per cm**3\n",
+ "mu_n=0.13*10**4##cm**2/V-s\n",
+ "sigma_n=q*n_n*mu_n##(ohm-cm)**-1\n",
+ "rho=1/sigma_n##ohm-cm\n",
+ "R=rho*l/A##ohm\n",
+ "print \"Resistance of the bar = %0.f Mohm\"%round(R/10**6)"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.19 Pg 1-113"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 140,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Answer is (B). Depletion width on p-side = 0.33 micro meter\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "NA=9*10**16##per cm**3\n",
+ "ND=1*10**16##per cm**3\n",
+ "w_total=3##micro meter\n",
+ "w_p=w_total*ND/NA##micro meter\n",
+ "print \"Answer is (B). Depletion width on p-side = %0.2f micro meter\"%w_p"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.20 Pg 1-113"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 141,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Majority carrier density = 4.50e+11 per m**3\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "ni=1.5*10**16##per m**3\n",
+ "n_n=5*10**20##per m**3\n",
+ "p_n=ni**2/n_n##per m**3\n",
+ "print \"Majority carrier density = %0.2e per m**3\"%p_n"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.21 Pg 1-113"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 142,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHXdJREFUeJzt3XuUXWWd5vHvQwIqAVSkBQUURKQBwUYBQUcsGwcDo2Ab\nWyZilMsaGR1vjM2gqEO8tbKW9qAwAiogoyMRgXaC4r0psdWWiwHCtYmKEvDGTSEgBvLMH3tXcurk\nVNWuU2ef6/NZ66ycves9+/xqr0q99Zx3v++WbSIiIiZs0usCIiKiv6RjiIiISdIxRETEJOkYIiJi\nknQMERExSTqGiIiYJB1DRERMko4hYgaSDpD0oKQFLb62QtJbe1FXRF3SMUTMwPa/AauB1zbul/Rc\nYHfggl7UFVGXdAwxsiTdLundkq6TdL+kZZIeN0Xz84E3Nu17I/AN2/dJerykL0m6W9J9kq6U9NQW\n77mLpHsk7VNuP13SHyQd1NnvLqJ9A9ExSDpX0u8krazQ9iBJP5O0VtKiFl/fStJqSafXU20MEAN/\nD7wC2BnYGzh6irZfAg6StAOApE2AxRQdBsCbgK2AHYCtgeOBhzd6Q/vnwEnAlyQ9ATgPOM/2FZ35\nliLmbiA6Bor/PAsrtv0VxX/SL0/x9Q8DP+hEUTEUPm37t7bvAy4F/qZVI9t3AOPAknLXwcDjgG+U\n238BngLs6sIK2w9McazPA6uAK4Ftgfd16HuJ6IiB6Bhs/xC4r3FfGcm/KelqSVdI2q1s+yvbK4F1\nzceR9ALgqcB3ulF3DITfNjx/GNhimrbns6FjWAJcYPuxcvuLwLeBZZLulHSqpPnTHOvzwJ7A6bbX\ntld6RD0GomOYwmeBt9veFzgR+Mx0jcvo/wng3V2oLQbTTEsN/zOwg6SXAX/Hho+RsP2o7Q/Z3hN4\nEfBKNh6TAEDSFsBpFJ3DByU9uRPFR3TKdH/R9K3yP9aBwFclTezebIaXvRW4zPZdanhRRINpfy5s\nr5F0EcVHm7fb/tn6F0pjwD3ATcADwFrgsVbHAT4FXGn7zZLOBs4Cjpx7+RGdMZAdA0XSud/2PjO0\na/wL8ADgJeU151sAm0l6wPbJdRUZA8fMnBrOpxigPqlp/3YUv+B3AB4EllF8vDSJpCOAQ4C9yl3/\nHbhW0mLbuew1+oLqvFGPpHOB/wT83vZeU7T5NHAo8BBwtO0VU7TbCbh04jiSfgT8L9sXlQlgL9vX\nN7T/Qtn+4hbHehOwr+23z+Hbi4gYSnWPMUx7NZGkw4Bn294VeDNw5hTtLgB+DOwm6Q5JxwBHAcdJ\nuha4ATi8bLufpDsoJiOdPc0lrrl1XUREC7UmBtj4L/2mr50FXG77K+X2LcBLbf+u1qIiImJKvb4q\naXvgjobt1RSf0UZERI/0umOAja8EyUc8ERE91Ourku4EdmzY3qHcN4mkdBYREW2wPevL83udGJZT\nTgKSdADFJagtxxd+8AOzyy7mDW8w99xj7NF8nHLKKT2voV8eORc5FzkX0z/aVWvH0OJqomMlHS/p\neADblwG/kLQKOJtiElpLBx0E110HW28Ne+0Fl15aZ+UREaOr1o+SbC+u0OZtVY+3YAF86lOwaBEc\neyxceGGxvfXWc6szIiI26PVHSW0Z5fQwNjbW6xL6Rs7FBjkXG+RczF3t8xg6QZKnqvOKK4r0cOCB\nSQ8REY0k4QEcfJ6zUU4PERF1GPjE0CjpISJig5FNDI2SHiIi5m6oEkOjpIeIGHVJDE2SHiIi2jO0\niaFR0kNEjKIkhmkkPUREVDcSiaFR0kNEjIokhoqSHiIipjdyiaFR0kNEDLMkhjYkPUREbGykE0Oj\npIeIGDZJDHOU9BARUUhiaCHpISKGQRJDByU9RMQoS2KYQdJDRAyqJIaaJD1ExKhJYpiFpIeIGCRJ\nDF2Q9BARoyCJoU1JDxHR75IYuizpISKGVRJDByQ9REQ/SmLooaSHiBgmSQwdlvQQEf0iiaFPJD1E\nxKBLYqhR0kNE9FISQx9KeoiIQZTE0CVJDxHRbUkMfS7pISIGRRJDDyQ9REQ3JDEMkKSHiOhnSQw9\nlvQQEXVJYhhQSQ8R0W+SGPpI0kNEdFJfJgZJCyXdIuk2SSe1+Po2kr4l6VpJN0g6us56+l3SQ0T0\ng9oSg6R5wK3Ay4E7gauAxbZvbmizFHic7fdK2qZsv63tR5uONRKJoVHSQ0TMVT8mhv2BVbZvt70W\nWAYc0dTmN8BW5fOtgHuaO4VRlfQQEb1SZ8ewPXBHw/bqcl+jzwF7SroLuA54Z431DJwFC4q0cMEF\ncMIJsGQJ3Htvr6uKiGE3v8ZjV/ns52TgWttjknYBvivpebYfaG64dOnS9c/HxsYYGxvrVJ19byI9\nnHxykR7OOgte9apeVxUR/WZ8fJzx8fE5H6fOMYYDgKW2F5bb7wXW2T61oc1lwEdt/6jc/j5wku2r\nm441cmMMU8nYQ0RU1Y9jDFcDu0raSdJmwJHA8qY2t1AMTiNpW2A34Bc11jTwMvYQEXWrdR6DpEOB\n04B5wDm2PybpeADbZ5dXIp0HPIOik/qY7S+3OE4SQwtJDxExnXYTQya4Dbg1a4qxh4suythDREyW\njmHEJT1ERLN+HGOILsrYQ0R0ShLDEEp6iAhIYogGSQ8RMRdJDEMu6SFidCUxREtJDxExW0kMIyTp\nIWK01JYYJO0p6S2STpX0cUn/VdKe7ZUZvZT0EBFVTJkYJC0B3g7cA1wJ3AUIeBrFktrbAJ+y/aXa\ni0xi6Likh4jh125imG511ScDB7da6bR8w62Ao2f7htEfsmJrREwlYwyR9BAxpDo+xiBpU0lvkDSx\nbPabJJ0h6ThJuZppiGTsISIaTTfGcA7wRGAz4GHgccDFwCuBX9s+sWtFJjF0TdJDxPDo+CJ6km60\nvaekTYHfAU+z/Yik+cDPbO89t5JnUWQ6hq7Kiq0Rw6GOy1XXAtheC1xl+5Fy+1Gq3bYzBlTuNR0x\n2qbrGH4raUsA26+Y2CnpacAjdRcWvZexh4jRNOurkiQtABbY/n09JbV8z3yU1GMZe4gYPLXMfJb0\nrvLfd0zss72mm51C9Iekh4jRMW1ikHQcsDVwr+1zulbVxnUkMfSRpIeIwVDHPIZTgOcAHwZ2K7cj\nkh4ihtxMieEjwHeAQ2y/v2tVbVxHEkOfSnqI6F91ra56he0rgPG2qoqhl/QQMXxm6hj2aPo3YiOZ\n9xAxXGbqGB6QdCKwphvFxGBLeogYDtMtiXEKsDnwTuDTwBrbH+xibY21ZIxhwGTsIaL3Oj7GUHYC\na4FDgL/0qlOIwZT0EDG4Zroq6RDb35H0ctvf62JdzXUkMQywpIeI3qjlqiTb3ymf/kHSEZIWlY/X\ntFVljKSkh4jBMuNaSZLOA/YCbgTWTey3fUy9pU2qIYlhSCQ9RHRPx+/H0HDgm4A9e/mbOR3DcMn9\nHiK6o64JbgBXkXkM0UGZ9xDR36p0DOcBP5H075JWlo/r6y4shl/GHiL6U5WPkn4OnADcwOQxhttr\nrWxyDfkoachl7CGi8+r8KOn3tpfb/oXt2ycesy8xYmpJDxH9o0pi+AzwJOBS4C/lbtu+pObaGmtI\nYhghSQ8RnVFnYtic4h7PhwCvLB+VriORtFDSLZJuk3TSFG3GJK2QdIOk8Yp1xxBLeojorVnf87ny\ngaV5wK3Ay4E7Ka5uWmz75oY2TwJ+BLzC9mpJ29i+u8WxkhhGVNJDRPtqSwySzi9/gU9sP1nSuRWO\nvT+wqhyTWAssA45oavN64GLbqwFadQox2pIeIrqvykdJz7N9/8SG7fuA51d43fbAHQ3bq8t9jXYF\ntpZ0uaSrJS2pcNwYMZn3ENFdVToGSdq6YWNrYF6F11X57GdTik7mMOAVwAck7VrhdTGCkh4iumN+\nhTafpJjgdiEg4O+Bj1Z43Z3Ajg3bO1KkhkZ3AHfbfhh4WNIVwPOA25oPtnTp0vXPx8bGGBsbq1BC\nDJuJ9LBoUTH2cOGFGXuImDA+Ps74+Picj1Np8FnSnsDfUqSAf7F9U4XXzKcYfD4YuAu4ko0Hn/8a\nOIMiLTwO+ClwZPPxM/gcrWTNpYjp1baI3lxIOhQ4jeKjp3Nsf0zS8QC2zy7b/ANwDMWs6s/Z/nSL\n46RjiCnlyqWI1vqyY+iUdAwxk6SHiI2lY4gg6SGiUS3zGCTNl3R5+2VFdFeuXIqYuyprJX0fWNQ4\nl6HbkhiiHUkPMerqXCtpDbBS0rmSTi8fGw0QR/SbpIeI9lRJDEeXTycaimJ11fNrrKu5hiSGmJOk\nhxhFtQ4+S9oceIbtW9opbq7SMUQn5MqlGDV1LqJ3OLAC+Fa5vY+k5bMvMaK3suZSRDVVxhiWAi8E\n7gOwvQJ4Vo01RdQqYw8R06vSMaxtcUXSupYtIwZE0kPE1Kp0DDdKOgqYL2lXSacDP665roiuSHqI\n2FiVq5IWAO+juLUnwLeBD9v+c821NdaQweeoXa5cimFT+5IYkp5IcZnqn2b7JnOVjiG6JVcuxTCp\nrWOQtB9wLrBVuet+4DjbV8+6yjalY4huS3qIYVDnzOdzgbfafqbtZwL/rdwXMbQy9hCjrEpiWGF7\nn6Z9P7Nd5b7PHZHEEL2U9BCDqs7E8ANJZ0saKx9nlvueL6lrnUNEryQ9xKipkhjG2bBO0kZsv6zD\nNbWqIYkh+kLSQwyS3Kgnokty5VIMinQMEV2W9BD9rs4xhohoIWMPMaySGCI6IOkh+lGdy26/TtJW\n5fMPSPrnXI0UMVnSQwyTKlclrbS9l6T/AHwE+ATwAdsv7EaBZQ1JDDEwkh6iX9Q5xvBY+e8rgc/Z\n/jqw2WzfKGJUJD3EoKuSGL4B3An8R2Af4M/AT20/r/7y1teQxBADKekheqnOxPA6iqW2Dylv2PNk\n4MTZvlHEKEp6iEFU6aokSfOBpwLzJ/bZ/nWNdTW/fxJDDLykh+i2Oq9KejvwW+B7wDcaHhExC0kP\nMSiqjDH8HNjf9j3dKallDUkMMVSSHqIb6hxj+DXQ9bu2RQyzpIfoZ1USw7nAcyg+PvpLudu2/6nm\n2hprSGKIoZX0EHWpOzF8j2LuwhbAluUjIjog6SH6TeW1kiRtCWD7gVorav3eSQwxEpIeopPqvCpp\nL0krgBuBGyVdI+m57RQZEdNLeoh+UGWM4SfAybYvL7fHgH+0/aL6y1tfQxJDjJykh5irOscYNp/o\nFABsjwMLZvtGETE7SQ/RK1U6hl+Wy23vJGlnSe8HflHl4JIWSrpF0m2STpqm3X6SHpX0mqqFR4yC\nBQuKtHDBBXDCCbBkCdx7b6+rimFXpWM4hmI5jEuAi4G/Ao6d6UWS5gFnAAuBPYDFknafot2pwLeA\nWUeeiFGQ9BDdNO0YQ7lG0ndtv2zWB5YOBE6xvbDcfg+A7Y83tXsXxfyI/YCv2764xbEyxhBRythD\nVFXLGIPtR4F1kp7URk3bA3c0bK8u960naXvgCODMibds430iRkrSQ9Rt/sxNWAOslPTd8jkUM5/f\nMcPrqvySPw14j21LEvkoKaKSibGHRYuK9HDhhUkP0TlVOoaLKcYXJn7Ri2q/9O8EdmzY3pEiNTR6\nAbCs6BPYBjhU0lrby5sPtnTp0vXPx8bGGBsbq1BCxHCbSA8nn1ykh7POgle9qtdVRa+Mj48zPj4+\n5+NUmcfwLtunzbSvxevmA7cCBwN3AVcCi23fPEX784BLbV/S4msZY4iYQcYeolmd8xje1GLf0TO9\nqByfeBvF3d9uAr5i+2ZJx0s6flZVRsSMMvYQnTJlYpC0GHg98BLghw1f2hJ4zPbB9Ze3vpYkhohZ\nSHoIaD8xTNcxPBPYGfg4cBIbBob/BFxfJoKuSMcQMXtr1hRjDxddlLGHUdXxjqHhwM8CfmP74XL7\nCcC2tm9vp9B2pGOIaF/Sw+iqc4zhQuCxhu11wEWzfaOI6I2MPcRsVekY5tueuHMbth8BNq2vpIjo\ntKy5FLNRpWO4W9IRExvl87vrKyki6pL0EFVUGWN4NvB/gaeXu1YDS2yvqrm2xhoyxhDRYRl7GH61\njTHYXmX7hcDuwB62D+xmpxAR9Uh6iKlUSQzbAR8Ftre9UNIewIG2z+lGgWUNSQwRNUp6GE51XpX0\nBeA7bPgo6TbghNm+UUT0r6SHaFSlY9jG9lcoL1m1vRbo2uS2iOiOXLkUE6p0DA9KesrEhqQDgD/W\nV1JE9FLSQ1QZY3gBcDqwJ3Ajxa09X2v7uvrLW19DxhgieiBjD4OtzquSrgFeCrwYeDPFlUld6xQi\noneSHkbTdIvoLaK4IY8a/qV8Tqv7JtQliSGi95IeBk8dq6t+gWnu1Gb7mNm+WbvSMUT0h6zYOlhq\nW121H6RjiOgvSQ+DobYxBknbSTpH0rfK7T0kHddOkRExHDL2MNyqXJX0LeA84H2295a0KbDC9nO7\nUWBZQxJDRJ9Keuhfdc58zgS3iJhS0sPwyQS3iJizzJoeLlU6hncDlwLPkvRj4IvAO2qtKiIGUtLD\ncKh0VVI5rrAbxVyGWxvv6NYNGWOIGDwZe+i9jo8xSBqbeG57re0bbK9s7BQkvWzWlUbESEh6GFzT\nTXD7BHAQ8D3gauA3FB3JdsC+wMuBy23/j9qLTGKIGGhJD71RywQ3SVsCR1Csk/TMcvevgH8F/p/t\nB9uoddbSMUQMvsya7r7MfI6IgZD00D21dQySHg8sAnYC5lEuqmf7Q23U2ZZ0DBHDJemhO+rsGL4N\n3A9cQznJDcD2J2f7Zu1KxxAxnJIe6lVnx3BDN5e/mKKGdAwRQyrpoT51dgyfBc6wfX27xc1VOoaI\n4Zf00Hl1rpX0EuAaSf8uaWX56FknERHDKfMe+keVxLBTq/22b+98OVPWkMQQMUKSHjqjzns+397q\n0VaVEREVJD30VuYxRERfS3poX51jDBERPZP00H1JDBExMJIeZqdvE4OkhZJukXSbpJNafP0oSddJ\nul7SjyTtXXdNETGYkh66o9bEIGkecCvFSqx3AlcBi23f3NDmQOAm23+UtBBYavuApuMkMUTEJEkP\nM+vXxLA/sKq8kmktsIxitdb1bP/E9sStQn8K7FBzTRExBJIe6lN3x7A9cEfD9upy31SOAy6rtaKI\nGBq513Q95td8/Mqf/5R3gzuW4t4PG1m6dOn652NjY4yNjc2xtIgYFhPp4eSTi/QwqmsujY+PMz4+\nPufj1D3GcADFmMHCcvu9wDrbpza12xu4BFhoe1WL42SMISIqydjDBv06xnA1sKuknSRtBhwJLG9s\nIOkZFJ3CG1p1ChERs5Gxh7mrfR6DpEOB0yhu8nOO7Y9JOh7A9tmSPg/8HfDr8iVrbe/fdIwkhoiY\ntVFPD7m1Z0REC6N8v4d0DBER0xjF9NCvYwwREX0hYw/VJTFExMgZlfSQxBARUVHSw/SSGCJipA1z\nekhiiIhoQ9LDxpIYIiJKw5YekhgiIuYo6aGQxBAR0cIwpIckhoiIDhrl9JDEEBExg0FND0kMERE1\nGbX0kMQQETELg5QekhgiIrpgFNJDEkNERJv6PT0kMUREdNmwpockhoiIDujH9JDEEBHRQ8OUHpIY\nIiI6rF/SQxJDRESfGPT0kMQQEVGjXqaHJIaIiD40iOkhiSEioku6nR6SGCIi+tygpIckhoiIHuhG\nekhiiIgYIP2cHpIYIiJ6rK70kMQQETGg+i09JDFERPSRTqaHJIaIiCHQD+khiSEiok/NNT0kMURE\nDJlepYckhoiIAdBOekhiiIgYYt1MD0kMEREDpmp66MvEIGmhpFsk3SbppCnafLr8+nWS9qmznoiI\nYVB3eqitY5A0DzgDWAjsASyWtHtTm8OAZ9veFXgzcGZd9QyL8fHxXpfQN3IuNsi52GBUzsWCBUVa\nuOACOOEEWLIE7r23M8euMzHsD6yyfbvttcAy4IimNocD5wPY/inwJEnb1ljTwBuVH/oqci42yLnY\nYNTORR3poc6OYXvgjobt1eW+mdrsUGNNERFDp9Ppoc6OoepocfPASEaZIyLa0Jwe2lXbVUmSDgCW\n2l5Ybr8XWGf71IY2ZwHjtpeV27cAL7X9u6ZjpbOIiGhDO1clza+jkNLVwK6SdgLuAo4EFje1WQ68\nDVhWdiT3N3cK0N43FhER7amtY7D9qKS3Ad8G5gHn2L5Z0vHl18+2fZmkwyStAtYAx9RVT0REVDMQ\nE9wiIqJ7+mpJjEyI22CmcyHpryX9RNKfJb27FzV2S4VzcVT583C9pB9J2rsXdXZDhXNxRHkuVki6\nRtLf9qLOulX5XVG220/So5Je0836uqnCz8SYpD+WPxMrJL1/xoPa7osHxcdNq4CdgE2Ba4Hdm9oc\nBlxWPn8h8G+9rruH5+KvgH2BjwDv7nXNPT4XBwJPLJ8vHPGfiwUNz/eimEvU89q7fR4a2v0L8HVg\nUa/r7uHPxBiwfDbH7afEkAlxG8x4Lmz/wfbVwNpeFNhFVc7FT2z/sdz8KcM7F6bKuVjTsLkFcHcX\n6+uWKr8rAN4OXAT8oZvFdVnVczGrC3j6qWPIhLgNqpyLUTHbc3EccFmtFfVOpXMh6dWSbga+Cbyj\nS7V104znQdL2FL8gJ5bZGdbB1Co/EwZeVH7EeJmkPWY6aJ2Xq85WJsRtMIzfU7sqnwtJLwOOBV5c\nXzk9Velc2P4a8DVJLwG+COxWa1XdV+U8nAa8x7YliVn+xTxAqpyLnwE72n5I0qHA14DnTPeCfkoM\ndwI7NmzvSNH7Tddmh3LfsKlyLkZFpXNRDjh/Djjc9n1dqq3bZvVzYfuHwHxJT6m7sC6rch5eQDE/\n6pfAIuAzkg7vUn3dNOO5sP2A7YfK598ENpU07W1++qljWD8hTtJmFBPilje1WQ68EdbPrG45IW4I\nVDkXE4b1L6EJM54LSc8ALgHeYHtVD2rslirnYpfyL2QkPR/A9j1dr7ReM54H28+yvbPtnSnGGd5i\ne6r/Q4Osys/Etg0/E/tTTFOYdiWlvvkoyZkQt16VcyFpO+AqYCtgnaR3AnvYfrBnhdegyrkA/ifw\nZODM8ud/re39e1VzXSqei0XAGyWtBR4E/nPPCq5JxfMwEiqei9cCb5H0KPAQFX4mMsEtIiIm6aeP\nkiIiog+kY4iIiEnSMURExCTpGCIiYpJ0DBERMUk6hoiImCQdQwwFSU+V9I1ZvuaDkg6uq6Z2lMtm\n797G6w6X9IE6aorRk3kMMRQkfQhYafurNR1/vu1H6zh20/t8AbjU9sWzeM08YB2wAtivXGUzom3p\nGGJgSNoP+DzFUsPzKZbYfp3tmyTdBOxbLhR2NPBqYHNgV+CTwOOB1wOPAIfZvq/xl3B57NOABcCf\ngZdTzBh9Tblvk/L5ecDOFDNI32x7ZVONVd97F+AMivtqPAT8F+ApwKXAH8vHa8r3ndTO9q1l7X8G\n/gb4V9v/IOlM4Ou2Z5WcIpr1zZIYETOxfZWk5RQ3J3oC8MWyU9gOeGxiobDSnhS/NJ8A/Bw40fbz\nJf0TxXpbn6JYmdLlGjPLKDqZayRtATxcHmcfYC/b90s6HbjG9qvLlVz/T/n1ZlXe+7PA8bZXSXoh\n8BnbB5ff36W2LwGQ9P3mdsDEx19PBw70hr/urgQOAtIxxJykY4hB8yGKhcMeprgRC8Azgd80tDFw\neXnTmjWS7qf4SxxgJdB4609RLEv9G9vXAEysNyXJwHdt31+2fTHFX/HYvlzSUyRt0bQ+1YzvLWkB\n8CLgq+XaTgCbNdVE2UEdOEU7A19t6BQA7qK4g13EnKRjiEGzDcVHO/Mo/iKfSAnNq8w+0vB8XcP2\nOjb+uZ/u89Q1TdtVVrOd6b03Ae6zPdU9yyfq2YRiBeGp2j3UtL0Js7h/RcRUclVSDJqzgfcDXwZO\nLff9Ctiuoc10v7xb3ejpVuBpkvYFkLRlOaDb3PaHwFFlmzHgDy1Ws53xvW0/APxS0mvLY6m8nwTA\nAxQr5mL7T9O0a+VpFOciYk7SMcTAkPRG4BHby4CPA/tJGrP9W4ob0mxeNjWT/3Jufj7pr+ryKp4j\ngdMlXUuxhPHjW7RdCrxA0nXAPwJvalFm1fc+CjiufL8bKO5nDsVYx4mSrpG08zTtmo8NxaD8FS1q\nipiVXJUUQ0HSUuBm21/pdS29IGkTils47tuNy2pjuCUxxLD437T+C35UvBK4KJ1CdEISQ0RETJLE\nEBERk6RjiIiISdIxRETEJOkYIiJiknQMERExSTqGiIiY5P8DtduvqIC+TFYAAAAASUVORK5CYII=\n",
+ "text/plain": [
+ "<matplotlib.figure.Figure at 0x7f5d8376ebd0>"
+ ]
+ },
+ "metadata": {},
+ "output_type": "display_data"
+ },
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Current density = 8.00 A/cm**2\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "q=1.6*10**-19##Coulomb\n",
+ "Dn=25##cm**2/s\n",
+ "x=[0 ,0.5]##micro meter(base width)\n",
+ "n=[10**14 ,0]##per cm**3\n",
+ "from matplotlib.pyplot import plot, title, show, xlabel, ylabel\n",
+ "plot(x,n)#\n",
+ "title('n Vs x')#\n",
+ "xlabel('x(micro meter)')#\n",
+ "ylabel('n(electrons per cm**3)')#\n",
+ "show()\n",
+ "dnBYdx=(n[1]-n[0])/(x[0]-x[1])/10**-4##gradient\n",
+ "Jn=q*Dn*dnBYdx##A/cm**2\n",
+ "print \"Current density = %0.2f A/cm**2\"%Jn"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.22 Pg 1-114"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 143,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Band gap = 1.431 eV\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "h=6.64*10**-34##planks constant\n",
+ "c=3*10**8##m/s(speed of light)\n",
+ "lamda=0.87*10**-6##m\n",
+ "Eg=h*c/lamda/(1.6*10**-19)##eV\n",
+ "print \"Band gap = %0.3f eV\"%Eg"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.23 Pg 1-114"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 144,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "(a) Absorbed power = 9.00 mW\n",
+ "(b) Rate of excess thermal energy = 2.564e-03 J/s\n",
+ "(c) No. of photons per sec : 2.81e+16\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "t=0.46*10**-4##cm\n",
+ "E=2##eV\n",
+ "alfa=5*10**4##cm**-1\n",
+ "Io=10##mW\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "It=Io*exp(-alfa*t)##mW\n",
+ "Pabs=Io-It##mW\n",
+ "print \"(a) Absorbed power = %0.2f mW\"%round(Pabs)\n",
+ "Eg=1.43##eV(Band gap)\n",
+ "heat_fraction=(E-Eg)/E#\n",
+ "E_heat=heat_fraction*Pabs*10**-3##J/s(energy converted to heat)\n",
+ "print \"(b) Rate of excess thermal energy = %0.3e J/s\"%E_heat\n",
+ "photons=Pabs*10**-3/q/E##no. of photons per sec\n",
+ "print \"(c) No. of photons per sec : %0.2e\"%photons"
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.24 Pg 1-115"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 145,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Steady state separation between Fp & Ec = 0.70 eV\n",
+ "Hole current = 1900.55 A\n",
+ "Excess stored hole charge = 1.44e-07 Coulomb\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "Kdash=0.0259##constant(taken as K*T/q)\n",
+ "A=0.5##cm**2\n",
+ "Na=10**17##per cm**3\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "delta_p=5*10**16##per cm**3\n",
+ "x=1000##Angstrum\n",
+ "mu_p=500##cm**2/V-s\n",
+ "tau_p=10**-10##sec\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "\n",
+ "Dp=Kdash*mu_p##cm/s\n",
+ "Lp=sqrt(Dp*tau_p)##cm\n",
+ "p0=Na##per cm**3\n",
+ "p=p0+delta_p*exp(x*10**-8/Lp)##per cm**3\n",
+ "delE1=log(p/ni)*Kdash##eV(taken as Ei-Fp)\n",
+ "Eg=1.12##eV(Band gap)\n",
+ "delE2=Eg-delE1##eV(taken as Ec-Fp)\n",
+ "print \"Steady state separation between Fp & Ec = %0.2f eV\"%delE2\n",
+ "Ip=q*A*Dp/Lp*delta_p*exp(x*10**-8/Lp)##A\n",
+ "print \"Hole current = %0.2f A\"%Ip\n",
+ "Qp=q*A*delta_p*Lp##C\n",
+ "print \"Excess stored hole charge = %0.2e Coulomb\"%Qp\n",
+ "#Answer in the book is wrong beacause of calculation mistake in the value of p & Ip."
+ ]
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "## Ex 1.40.25 Pg 1-116"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 146,
+ "metadata": {
+ "collapsed": false
+ },
+ "outputs": [
+ {
+ "name": "stdout",
+ "output_type": "stream",
+ "text": [
+ "Steady state separation between Fp & Ec = 0.70 eV\n",
+ "Hole current = 1900.55 A \n",
+ "Excess stored hole charge = 1.44e-07 Coulomb\n"
+ ]
+ }
+ ],
+ "source": [
+ "#Given : \n",
+ "Kdash=0.0259##constant(taken as K*T/q)\n",
+ "A=0.5##cm**2\n",
+ "Na=10**17##per cm**3\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "delta_p=5*10**16##per cm**3\n",
+ "x=1000##Angstrum\n",
+ "mu_p=500##cm**2/V-s\n",
+ "tau_p=10**-10##sec\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "\n",
+ "Dp=Kdash*mu_p##cm/s\n",
+ "Lp=sqrt(Dp*tau_p)##cm\n",
+ "p0=Na##per cm**3\n",
+ "p=p0+delta_p*exp(x*10**-8/Lp)##per cm**3\n",
+ "delE1=log(p/ni)*Kdash##eV(taken as Ei-Fp)\n",
+ "Eg=1.12##eV(Band gap)\n",
+ "delE2=Eg-delE1##eV(taken as Ec-Fp)\n",
+ "print \"Steady state separation between Fp & Ec = %0.2f eV\"%delE2\n",
+ "Ip=q*A*Dp/Lp*delta_p*exp(x*10**-8/Lp)##A\n",
+ "print \"Hole current = %0.2f A \"%Ip\n",
+ "Qp=q*A*delta_p*Lp##C\n",
+ "print \"Excess stored hole charge = %0.2e Coulomb\"%Qp\n",
+ "#Answer in the book is wrong beacause of calculation mistake in the value of p & Ip."
+ ]
+ }
+ ],
+ "metadata": {
+ "kernelspec": {
+ "display_name": "Python 2",
+ "language": "python",
+ "name": "python2"
+ },
+ "language_info": {
+ "codemirror_mode": {
+ "name": "ipython",
+ "version": 2
+ },
+ "file_extension": ".py",
+ "mimetype": "text/x-python",
+ "name": "python",
+ "nbconvert_exporter": "python",
+ "pygments_lexer": "ipython2",
+ "version": "2.7.9"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter2.ipynb b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter2.ipynb
new file mode 100644
index 00000000..fff383e3
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter2.ipynb
@@ -0,0 +1,1262 @@
+{
+ "metadata": {
+ "name": "",
+ "signature": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter-2 Junctions and interfaces"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.6.1 Pg 2-21"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from __future__ import division\n",
+ "from math import log\n",
+ "#Given : \n",
+ "Ge=4.4*10**22##atoms/cm**3\n",
+ "NA=Ge/10**8##per cm**3\n",
+ "NA=NA*10**6##per m**3\n",
+ "ND=NA*10**3##per m**3\n",
+ "ni=2.5*10**13##per cm**3\n",
+ "ni=ni*10**6##per m**3\n",
+ "VT=26##mV\n",
+ "Vj=VT*log(NA*ND/ni**2)##mV\n",
+ "print \"Junction potential = %0.1f mV\"%Vj"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Junction potential = 328.7 mV\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.6.2 Pg 2-22"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "ni=2.5*10**15##per cm**3\n",
+ "Ge=4.4*10**22##atoms/cm**3\n",
+ "NA=Ge/10**8##per cm**3\n",
+ "NA=NA*10**6##per m**3\n",
+ "ND=NA*10**3##per m**3\n",
+ "ni=ni*10**6##per m**3\n",
+ "T=27+273##K\n",
+ "VT=T/11600##V\n",
+ "Vo=VT*log(NA*ND/ni**2)##V\n",
+ "print \"Contact potential = %0.4f V\"%Vo"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Contact potential = 0.0888 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.6.3 Pg 2-23"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "mu_n=1500*10**-4##m**2/V-s\n",
+ "mu_p=475*10**-4##m**2/V-s\n",
+ "ni=1.45*10**10*10**6##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "rho_p=10##ohm-cm\n",
+ "rho_p=rho_p*10**-2##ohm-m\n",
+ "rho_n=3.5##ohm-cm\n",
+ "rho_n=rho_n*10**-2##ohm-m\n",
+ "sigma_p=1/rho_p##(ohm-m)**-1\n",
+ "NA=sigma_p/q/mu_p##m**3\n",
+ "sigma_n=1/rho_n##(ohm-m)**-1\n",
+ "ND=sigma_p/q/mu_n##m**3\n",
+ "VT=26*10**-3##V\n",
+ "Vj=VT*log(NA*ND/ni**2)##V\n",
+ "print \"Height of potential barrier = %0.3f V\"%Vj\n",
+ "#Answer in the book is wrong."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Height of potential barrier = 0.564 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 6
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.6.4 Pg 2-24"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "rho_p=2##ohm-cm\n",
+ "rho_p=rho_p*10**-2##ohm-m\n",
+ "rho_n=1##ohm-cm\n",
+ "rho_n=rho_n*10**-2##ohm-m\n",
+ "mu_n=1500*10**-4##m**2/V-s\n",
+ "mu_p=2100*10**-4##m**2/V-s\n",
+ "ni=2.5*10**13##per m**3\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "sigma_p=1/rho_p##(ohm-m)**-1\n",
+ "NA=sigma_p/q/mu_p##m**3\n",
+ "sigma_n=1/rho_n##(ohm-m)**-1\n",
+ "ND=sigma_p/q/mu_n##m**3\n",
+ "T=27+273##K\n",
+ "VT=T/11600##V\n",
+ "Vj=VT*log(NA*ND/ni**2)##V\n",
+ "print \"Height of potential barrier = %0.4f V\"%Vj\n",
+ "#Anser in the book is wrong."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Height of potential barrier = 0.9347 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 8
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.7.1 Pg 2-27"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Vgamma=0.6##Volt\n",
+ "rf=12##ohm\n",
+ "V=5##Volts\n",
+ "R=1##kohm\n",
+ "IF=(V-Vgamma)/(R*1000+rf)##A\n",
+ "print \"Diode current = %0.1f mA\"%(IF*1000)\n",
+ "VF=Vgamma+IF*rf##volts\n",
+ "print \"Diode voltage = %0.2f Volts\"%VF"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Diode current = 4.3 mA\n",
+ "Diode voltage = 0.65 Volts\n"
+ ]
+ }
+ ],
+ "prompt_number": 10
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.7.2 Pg 2-35"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Vgamma=0.6##Volt\n",
+ "Rf=10##ohm\n",
+ "Eta=2#\n",
+ "Vm=0.2##Volts\n",
+ "Vdc=10##Volts\n",
+ "RL=1##kohm\n",
+ "IDQ=(Vdc-Vgamma)/(RL*1000+Rf)##A\n",
+ "VT=25*10**-3##Volts\n",
+ "rd=Eta*VT/IDQ##ohm\n",
+ "print \"Alternating component of voltage across RL, Vo(ac) = \",round((RL*1000/(RL*1000+rd)*Vm),4),\"*sin(omega*t)\"\n",
+ "Vo_DC=IDQ*RL*1000##Volts\n",
+ "print \"Total load voltage = \",round(Vo_DC,1),\"+\",round((RL*1000/(RL*1000+rd)*Vm),4),\"*sin(omega*t)\""
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Alternating component of voltage across RL, Vo(ac) = 0.1989 *sin(omega*t)\n",
+ "Total load voltage = 9.3 + 0.1989 *sin(omega*t)\n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.7.3 Pg 2-37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import exp\n",
+ "#Given : \n",
+ "Eta=2##for Si diode\n",
+ "T=300##K\n",
+ "VT=T/11600##V\n",
+ "IbyIo=90/100#\n",
+ "#I=Io*(exp(V/Eta/VT)-1)\n",
+ "V=log(IbyIo+1)*Eta*VT##V\n",
+ "print \"Saturation value of voltage = %0.2f mV\"%(V*1000)\n",
+ "VF=0.5##Volts\n",
+ "VR=-0.5##Volts\n",
+ "IFbyIR=(exp(VF/Eta/VT)-1)/(exp(VR/Eta/VT)-1)##ratio\n",
+ "print \"Ratio of forward to reverse current = %0.2f\"%IFbyIR\n",
+ "#Answer in the book is wrong."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Saturation value of voltage = 33.20 mV\n",
+ "Ratio of forward to reverse current = -15782.65\n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.7.4 Pg 2-37"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Eta=2##for Si diode\n",
+ "T=300##K\n",
+ "VT=T/11600##V\n",
+ "IbyIo=90/100#\n",
+ "#I=Io*(exp(V/Eta/VT)-1)\n",
+ "V=log(IbyIo+1)*Eta*VT##V\n",
+ "print \"Saturation value of voltage = %0.2f mV\"%(V*1000)\n",
+ "VF=0.2##Volts\n",
+ "VR=-0.2##Volts\n",
+ "IFbyIR=(exp(VF/Eta/VT)-1)/(exp(VR/Eta/VT)-1)##ratio\n",
+ "print \"Ratio of forward to reverse current : %0.2f \"%IFbyIR\n",
+ "#Answer in the book is wrong."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Saturation value of voltage = 33.20 mV\n",
+ "Ratio of forward to reverse current : -47.78 \n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.1 Pg 2-61"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "IF=10##mA\n",
+ "VF=0.75##volts\n",
+ "T=27+273##K\n",
+ "Eta=2##for Si diode\n",
+ "VT=T/11600##V\n",
+ "Io=IF/(exp(VF/Eta/VT)-1)##mA\n",
+ "print \"Reverse saturation current = %0.3f nA\"%(Io*10**6)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Reverse saturation current = 5.043 nA\n"
+ ]
+ }
+ ],
+ "prompt_number": 18
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.2 Pg -61"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "IF=10##mA\n",
+ "VF=0.3##Volts\n",
+ "T=27+273##K\n",
+ "Eta=1##for Ge diode\n",
+ "VT=T/11600##V\n",
+ "Io=IF/(exp(VF/Eta/VT)-1)##mA\n",
+ "print \"Reverse saturation current = %0.2f nA\"%(Io*10**6)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Reverse saturation current = 91.66 nA\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.3 Pg 2-61"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io=1*10**-9##A\n",
+ "T=27+273##K\n",
+ "VT=T/11600##V\n",
+ "VF=0.3##Volts\n",
+ "Eta=1##for Ge diode\n",
+ "IF=Io*(exp(VF/Eta/VT)-1)##mA\n",
+ "print \"Forwad current = %0.4f mA\"%(IF*10**3)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Forwad current = 0.1091 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 21
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.4 Pg 2-62"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "T=27+273##K\n",
+ "V1=0.4##V\n",
+ "V2=0.42##V\n",
+ "I1=10##mA\n",
+ "I2=20##mA\n",
+ "VT=T/11600##V\n",
+ "Eta=1/log(I1/I2)*(V1-V2)/VT\n",
+ "print \"Value of Eta : %0.2f\"%Eta\n",
+ "Io=I1/(exp(V1/Eta/VT)-1)*10**-3##A\n",
+ "print \"Current, Io = %0.2f nA\"%(Io*10**9)\n",
+ "#Ans in the book is not accurate."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Value of Eta : 1.12\n",
+ "Current, Io = 9.54 nA\n"
+ ]
+ }
+ ],
+ "prompt_number": 22
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.5 Pg 2-63"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io1=10**-12##A\n",
+ "Io2=10**-10##A\n",
+ "I=2##mA\n",
+ "Eta=1##constant\n",
+ "T=27+273##K\n",
+ "VT=26/1000##V\n",
+ "#I=I1+I2\n",
+ "V=(log(I*10**-3/(Io1+Io2))+1)*Eta*VT##V\n",
+ "print \"Voltage across the diodes = %0.4f V\"%V\n",
+ "#Ans in the book is not accurate."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Voltage across the diodes = 0.4628 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.6 Pg 2-64"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io1=10*10**-9##A\n",
+ "Io2=10*10**-9##A\n",
+ "Eta=1.1##constant\n",
+ "T=25+273##K\n",
+ "V=0.2##V(assumed)\n",
+ "VT=T/11600##V\n",
+ "I1=Io1*(exp(V/Eta/VT)-1)##A\n",
+ "I2=Io2*(exp(V/Eta/VT)-1)##A\n",
+ "I=I1+I2##A\n",
+ "print \"Source current = %0.2f micro Ampere\"%(I*10**6)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Source current = 23.68 micro Ampere\n"
+ ]
+ }
+ ],
+ "prompt_number": 25
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.7 Pg 2-65"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io=10**-13##A\n",
+ "T=27+273##K\n",
+ "Eta=1##constant\n",
+ "V=0.6##V\n",
+ "VT=26/1000##V\n",
+ "I3=Io*(exp(V/Eta/VT)-1)##A\n",
+ "R=1*1000##ohm\n",
+ "Ir=V/R##A\n",
+ "Itotal=I3+Ir##A\n",
+ "VD1=log(Itotal/Io)*Eta*VT##V\n",
+ "VD2=VD1##V\n",
+ "Vin=VD1+VD2+V##V\n",
+ "print \"Voltage Vin = %0.3f V\"%Vin"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Voltage Vin = 1.823 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 27
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.8 Pg 2-66"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Vs=10##V\n",
+ "print \"Case(i) : Vb=9.8V\"\n",
+ "Vb=9.8##V\n",
+ "#D1 forward & D2 reverse biased: Breakdown D2\n",
+ "VD2=Vb##V\n",
+ "VD1=Vs-Vb##V\n",
+ "print \"VD1 = %0.3f V\"%VD1\n",
+ "print \"VD2 = %0.3f V\"%VD2\n",
+ "print \"Case(ii) : Vb=10.2V\"\n",
+ "Vb=10.2##V\n",
+ "#D1 forward & D2 reverse biased: none will be breakdown\n",
+ "VD2=Vb##V\n",
+ "#I=I0 so exp(V1/Eta/VT)-1=1\n",
+ "Eta=1##constant\n",
+ "VT=26/1000##V\n",
+ "VD1=log(1+1)*Eta*VT##V\n",
+ "VD2=Vs-VD1##V\n",
+ "print \"VD1 = %0.3f V\"%VD1\n",
+ "print \"VD2 = %0.3f V\"%VD2"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Case(i) : Vb=9.8V\n",
+ "VD1 = 0.200 V\n",
+ "VD2 = 9.800 V\n",
+ "Case(ii) : Vb=10.2V\n",
+ "VD1 = 0.018 V\n",
+ "VD2 = 9.982 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 30
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.9.9 Pg 2-67"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Vs=5##Volt\n",
+ "Eta=1##constant\n",
+ "VT=26/1000##V\n",
+ "#I=I0 so exp(V1/Eta/VT)-1=1\n",
+ "V1=log(1+1)*Eta*VT##Volt\n",
+ "V2=Vs-V1##Volt\n",
+ "print \"Voltage across diode D1 = %0.3f V\"%V1\n",
+ "print \"Voltage across diode D2 = %0.3f V\"%V2"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Voltage across diode D1 = 0.018 V\n",
+ "Voltage across diode D2 = 4.982 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 32
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.10.2 Pg 2-70"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "rho_n=10##ohm-cm\n",
+ "rho_p=3.5##ohm-cm\n",
+ "ni=1.5*10**10##per cm**3\n",
+ "Vj=0.56##volt\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "mu_n=1500##cm**2/V-s\n",
+ "mu_p=500##cm**2/V-s\n",
+ "sigma_p=1/rho_p##(ohm-cm)**-1\n",
+ "NA=sigma_p/q/mu_p##per cm**3\n",
+ "sigma_n=1/rho_n##(ohm-cm)**-1\n",
+ "ND=sigma_n/q/mu_n##per cm**3\n",
+ "VT=Vj/log(NA*ND/ni**2)##V\n",
+ "T=11600*VT##K\n",
+ "print \"Temperature of junction = %0.2f degree K\"%T\n",
+ "t=T-273##degree C\n",
+ "print \"Temperature of junction = %0.2f degree C\"%t"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Temperature of junction = 287.28 degree K\n",
+ "Temperature of junction = 14.28 degree C\n"
+ ]
+ }
+ ],
+ "prompt_number": 34
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.11.1 Pg 2-75"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io=10##nA\n",
+ "T1=27+273##K\n",
+ "T2=87+273##K\n",
+ "VT=T1/11600##V\n",
+ "Eta=2##for Si\n",
+ "m=1.5##for Si\n",
+ "VGO=-1.21##volt\n",
+ "K=Io*10**-9/T1**m/exp(VGO/Eta/VT)##constant\n",
+ "VT=T2/11600##V\n",
+ "Io2=K*T2**m*exp(VGO/Eta/VT)##A\n",
+ "print \"Reverse saturation current at 87 degree C = %0.2f nA\"%(Io2*10**9)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Reverse saturation current at 87 degree C = 648.69 nA\n"
+ ]
+ }
+ ],
+ "prompt_number": 35
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.11.2 Pg 2-76"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "V=0.45##volt\n",
+ "Eta=2##for Si\n",
+ "T1=27+273##K\n",
+ "T2=125+273##K\n",
+ "VT1=T1/11600##V\n",
+ "VT2=T2/11600##V\n",
+ "I1BYIo1=exp(V/Eta/VT1)#\n",
+ "I2BYIo2=exp(V/Eta/VT2)#\n",
+ "m=1.5##for Si\n",
+ "VGO=1.21##volt\n",
+ "Io1BYIo2=(T1/T2)**m*exp(-VGO/Eta/VT1+VGO/Eta/VT2)##constant\n",
+ "I2BYI1=I2BYIo2/I1BYIo1/Io1BYIo2#\n",
+ "print \"Factor by which current increases : %0.2f \"%I2BYI1\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Factor by which current increases : 56.94 \n"
+ ]
+ }
+ ],
+ "prompt_number": 36
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.11.3 Pg 2-78"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Io1=2##nA\n",
+ "T1=10+273##K\n",
+ "V=0.4##volt\n",
+ "VT1=T1/11600##V\n",
+ "m=1.5##for Si\n",
+ "Eta=2##for Si\n",
+ "VGO=-1.21##volt\n",
+ "K=Io1*10**-9/T1**m/exp(VGO/Eta/VT1)##constant\n",
+ "I1=Io1*10**-9*(exp(V/Eta/VT1)-1)##nA\n",
+ "T2=70+273##K\n",
+ "VT2=T2/11600##V\n",
+ "Io2=K*T2**m*(exp(VGO/Eta/VT2))##A\n",
+ "I2=Io2*(exp(V/Eta/VT2)-1)##nA\n",
+ "change=(I2-I1)/I1*100##%\n",
+ "print \"%% change = %0.2f diode current\"%change\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "% change = 2332.39 diode current\n"
+ ]
+ }
+ ],
+ "prompt_number": 37
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.11.4 Pg 2-79"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "T=300##K\n",
+ "m_Si=1.5##for Si\n",
+ "m_Ge=1.5##for Ge\n",
+ "EGO_Si=1.21##Volt\n",
+ "EGO_Ge=0.785##Volt\n",
+ "Eta_Si=2#\n",
+ "Eta_Ge=1#\n",
+ "VT=26/1000##V\n",
+ "print \"Part(i)\"\n",
+ "d_logIoBYdt_Ge=m_Ge/T+EGO_Ge/(Eta_Ge*T*VT)##per degree C\n",
+ "print \"d(log(Io))/dt for Ge = %0.2f per degree C\"%d_logIoBYdt_Ge\n",
+ "d_logIoBYdt_Si=m_Si/T+EGO_Si/(Eta_Si*T*VT)##per degree C\n",
+ "print \"d(log(Io))/dt for Si = %0.2f per degree C \"%d_logIoBYdt_Si\n",
+ "print \"Part(ii)\"\n",
+ "V=0.2##Volt\n",
+ "dVBYdt_Ge=V/T-Eta_Ge*VT*d_logIoBYdt_Ge#\n",
+ "print \"dV/dt for Si = %0.2f mV per degree C \"%(dVBYdt_Ge*1000)\n",
+ "V=0.6##Volt\n",
+ "dVBYdt_Si=V/T-Eta_Si*VT*d_logIoBYdt_Si\n",
+ "print \"dV/dt for Si = %0.2f mV per degree C \"%(dVBYdt_Si*1000)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Part(i)\n",
+ "d(log(Io))/dt for Ge = 0.11 per degree C\n",
+ "d(log(Io))/dt for Si = 0.08 per degree C \n",
+ "Part(ii)\n",
+ "dV/dt for Si = -2.08 mV per degree C \n",
+ "dV/dt for Si = -2.29 mV per degree C \n"
+ ]
+ }
+ ],
+ "prompt_number": 38
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.1 Pg 2-85"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt\n",
+ "#Given : \n",
+ "NA=4*10**20##per m**3\n",
+ "Vj=0.2##Volt\n",
+ "V1=-1##Volts\n",
+ "V2=-5##Volts\n",
+ "epsilon_r=16##for Ge\n",
+ "epsilon_o=8.85*10**-12##permitivity\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "W1=sqrt(2*epsilon_r*epsilon_o*(Vj-V1)/q/NA)##m\n",
+ "print \"Width of depletion region = %0.2f micro meter \"%(W1*10**6)\n",
+ "W2=sqrt(2*epsilon_r*epsilon_o*(Vj-V2)/q/NA)##m\n",
+ "print \"New value of Width of depletion region = %0.2f micro meter \"%(W2*10**6)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Width of depletion region = 2.30 micro meter \n",
+ "New value of Width of depletion region = 4.80 micro meter \n"
+ ]
+ }
+ ],
+ "prompt_number": 39
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.2 Pg 2-86"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "NA=4*10**20##per m**3\n",
+ "Vj=0.2##Volt\n",
+ "V1=-1##Volts\n",
+ "V2=-5##Volts\n",
+ "A=0.8*10**-6##m**2\n",
+ "epsilon_r=16##for Ge\n",
+ "epsilon_o=8.85*10**-12##permitivity\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "W1=sqrt(2*epsilon_r*epsilon_o*(Vj-V1)/q/NA)##m\n",
+ "CT1=epsilon_r*epsilon_o*A/W1##\n",
+ "print \"Transition capacitance = %0.2f pF \"%(CT1*10**12)\n",
+ "W2=sqrt(2*epsilon_r*epsilon_o*(Vj-V2)/q/NA)##m\n",
+ "CT2=epsilon_r*epsilon_o*A/W2##\n",
+ "print \"New value of Transition capacitance = %0.2f pF \"%( CT2*10**12)"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Transition capacitance = 49.16 pF \n",
+ "New value of Transition capacitance = 23.62 pF \n"
+ ]
+ }
+ ],
+ "prompt_number": 40
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.3 Pg 2-87"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "NA=3*10**20##per m**3\n",
+ "Vj=0.2##Volt\n",
+ "V=-10##Volts\n",
+ "A=1*10**-6##m**2\n",
+ "epsilon_r=16##for Ge\n",
+ "epsilon_o=8.854*10**-12##permitivity\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "W=sqrt(2*epsilon_r*epsilon_o*(Vj-V)/q/NA)##m\n",
+ "print \"Width of depletion region = %0.2f micro meter\"%(W*10**6)\n",
+ "CT=epsilon_r*epsilon_o*A/W##\n",
+ "print \"Transition capacitance = %0.2f pF\"%(CT*10**12)\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Width of depletion region = 7.76 micro meter\n",
+ "Transition capacitance = 18.26 pF\n"
+ ]
+ }
+ ],
+ "prompt_number": 41
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.4 Pg 2-88"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "W=2*10**-4*10**-2##m\n",
+ "A=1*10**-6##m**2\n",
+ "epsilon_r=16##for Ge\n",
+ "epsilon_o=8.854*10**-12##permitivity\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "CT=epsilon_r*epsilon_o*A/W##\n",
+ "print \"Barrier capacitance = %0.2f pF \"%(CT*10**12)\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Barrier capacitance = 70.83 pF \n"
+ ]
+ }
+ ],
+ "prompt_number": 42
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.5 Pg 2-88"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "#Given : \n",
+ "Vj=0.5##Volt\n",
+ "V=-4.5##Volt\n",
+ "rho_p=5*10**-2##ohm-m\n",
+ "epsilon_r=12##for Si\n",
+ "epsilon_o=8.854*10**-12##permitivity\n",
+ "q=1.6*10**-19##Coulomb\n",
+ "CT=100*10**-12##F\n",
+ "mu_p=500*10**-4##m**2/V-s\n",
+ "sigma_p=1/rho_p##(ohm-m)**-1\n",
+ "NA=sigma_p/q/mu_p##per m**3\n",
+ "W=sqrt(2*epsilon_r*epsilon_o*(Vj-V)/q/NA)##m\n",
+ "A=CT*W/(epsilon_r*epsilon_o)##\n",
+ "r=sqrt(A/pi)##m\n",
+ "D=2*r##m\n",
+ "print \"Diameter = %0.2f micro meter\"%(D*10**6)\n",
+ "#Answer is wrong = %0.2f the textbook. Sqrt is not taken while calculatng W value and also other mistakes."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Diameter = 1397.53 micro meter\n"
+ ]
+ }
+ ],
+ "prompt_number": 43
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.6 Pg 2-90"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Eta=2##for Si\n",
+ "T=300##K\n",
+ "VT=26/1000##V\n",
+ "IbyIo=0.9#\n",
+ "#part (i)\n",
+ "V=log(IbyIo+1)*Eta*VT##Volt\n",
+ "print \"Value of reverse voltage = %0.2f mV\"%(V*1000)\n",
+ "#part (ii)\n",
+ "VF=0.2##Volt\n",
+ "VR=-0.2##Volt\n",
+ "IFbyIR=(exp(VF/Eta/VT)-1)/(exp(VR/Eta/VT)-1)#\n",
+ "print \"Ratio of forward bias current to reverse saturation current = %0.2f \"%IFbyIR\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Value of reverse voltage = 33.38 mV\n",
+ "Ratio of forward bias current to reverse saturation current = -46.81 \n"
+ ]
+ }
+ ],
+ "prompt_number": 44
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.7 Pg 2-91"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Vs=100##V\n",
+ "Rf1=20##ohm\n",
+ "Vgamma1=0.2##Volts\n",
+ "Rf2=15##ohm\n",
+ "Vgamma2=0.6##Volts\n",
+ "Vb_Ge=0.2##Volts\n",
+ "Vb_Si=0.6##Volts\n",
+ "R1=10*10**3##ohm\n",
+ "R2=1*10**3##ohm\n",
+ "#Case(i)\n",
+ "Imax=Vs/R1##A\n",
+ "#D1 ON & D2 off\n",
+ "V=Vb_Ge+Rf1*Imax##Volt\n",
+ "#D2 off as V<Vb_Si\n",
+ "I2=0##A\n",
+ "I1=(Vs-V)/(R1+Rf1)##A\n",
+ "print \"For R=10 kohm\"\n",
+ "print \"I1 = %0.2f mA\"%(I1*1000)\n",
+ "print \"I2 = %0.2f mA\"%I2\n",
+ "#Case(ii)\n",
+ "R=R2##ohm#D1 & D2 ON \n",
+ "#V=Vb_Ge+Rf1*I1#V=Vb_Si+Rf2*I2\n",
+ "#V=Vs-I*R#V=Vs-(I1+I2)*R\n",
+ "#20*I1-15*I2=Vb_Si-Vb_Ge\n",
+ "#1020*I1+1000*I2=99.8\n",
+ "from numpy import mat, linalg\n",
+ "A=mat([[20, 1020],[-Rf2, R]])#\n",
+ "B=mat([[Vb_Ge-Vb_Ge],[Vs-Vb_Ge]])#\n",
+ "X = linalg.solve(A,B)\n",
+ "I1=X[0]*1000##mA\n",
+ "I2=X[1]*1000##mA\n",
+ "print \"For R=1 kohm\"\n",
+ "print \"I1 = %0.2f mA\"%I1\n",
+ "print \"I2 = %0.2f mA\"%I2\n",
+ "#Answer for 2nd part is not accurate in the book."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "For R=10 kohm\n",
+ "I1 = 9.94 mA\n",
+ "I2 = 0.00 mA\n",
+ "For R=1 kohm"
+ ]
+ },
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "\n",
+ "I1 = -2883.74 mA\n",
+ "I2 = 56.54 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 45
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 2.12.8 Pg 2-93"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "Rf=10##ohm\n",
+ "Vgamma=0.5##Volt\n",
+ "RL=20##ohm\n",
+ "V=3##Volt\n",
+ "#Loop 1: 75*I1-50*I=V-Vgamma\n",
+ "#Loop 2: -50*I1+80*I=-Vgamma\n",
+ "import numpy as np\n",
+ "A=np.mat([[75 ,-50],[-50, 80]])#\n",
+ "B=np.mat([[V-Vgamma], [-Vgamma]])#\n",
+ "X = linalg.solve(A,B)\n",
+ "I1=X[0]*1000##mA\n",
+ "I2=X[1]*1000##mA\n",
+ "print \"For R=1 kohm\"\n",
+ "Vx=-Vgamma+50*I1##Volt\n",
+ "print \"DC source = %0.2f Volts\"%Vx[0,0]\n",
+ "#Answer is wrong in the textbook."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "For R=1 kohm\n",
+ "DC source = 2499.50 Volts\n"
+ ]
+ }
+ ],
+ "prompt_number": 46
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter5.ipynb b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter5.ipynb
new file mode 100644
index 00000000..725c9736
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/chapter5.ipynb
@@ -0,0 +1,64 @@
+{
+ "metadata": {
+ "name": "",
+ "signature": ""
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "Chapter-5 Metal semiconductor field effect transistors"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Ex 5.6.1 Pg 5-22"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Given : \n",
+ "VTN=0.7##V\n",
+ "W=45##micro m\n",
+ "L=4##micro m\n",
+ "mu_n=700##cm**2/V-s\n",
+ "t_ox=450##Angstrum\n",
+ "epsilon_ox=3.9*8.85*10**-14##F/cm\n",
+ "VGS=2*VTN##V\n",
+ "Kn=(W*10**-4)*mu_n*epsilon_ox/(2*(L*10**-4)*(t_ox*10**-8))##A/V**2\n",
+ "Kn=Kn*10**3##mA/V**2\n",
+ "print \"Kn = %0.3f mA/V**2\"%Kn\n",
+ "ID=Kn*(VGS-VTN)**2##A\n",
+ "print \"Current = %0.2f mA\"%ID\n",
+ "#Answer is wrong in the book. Calculation mistake whle calculating value for Kn."
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Kn = 0.302 mA/V**2\n",
+ "Current = 0.15 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+}
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ1_chapter1.png b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ1_chapter1.png
new file mode 100644
index 00000000..63a8c57c
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ1_chapter1.png
Binary files differ
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ_ch1.png b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ_ch1.png
new file mode 100644
index 00000000..5cd8e98d
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_econ_ch1.png
Binary files differ
diff --git a/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_percentChangeinDiodeCurrent_chapter2.png b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_percentChangeinDiodeCurrent_chapter2.png
new file mode 100644
index 00000000..a378a6f1
--- /dev/null
+++ b/Advance_Semiconductor_Devices_by_K._C._Nandi/screenshots/KC_percentChangeinDiodeCurrent_chapter2.png
Binary files differ
diff --git a/sample_notebooks/SUMITPRADHAN/chapter6.ipynb b/sample_notebooks/SUMITPRADHAN/chapter6.ipynb
new file mode 100644
index 00000000..f9a54cd4
--- /dev/null
+++ b/sample_notebooks/SUMITPRADHAN/chapter6.ipynb
@@ -0,0 +1,1632 @@
+{
+ "metadata": {
+ "name": "",
+ "signature": "sha256:24c423084544eb74b3903e47cf4df50df61e5b0825c0f2f723c85cbcdae27d10"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+ {
+ "cells": [
+ {
+ "cell_type": "heading",
+ "level": 1,
+ "metadata": {},
+ "source": [
+ "CHAPTER 6: SEMICONDUCTOR DIODE"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.2, Page number 81"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration \n",
+ "Vf =20; #Peak Input Voltage in V\n",
+ "rf=10; #Forward Resistance in ohms\n",
+ "RL=500.0; #Load Resistance in ohms\n",
+ "V0=0.7; #Potential Barrier Voltage of the diodes in V\n",
+ "\n",
+ "#Calculation\n",
+ "#(1)\n",
+ "If_peak=(Vf-V0)/(rf+RL); #Peak current through the diode in A\n",
+ "If_peak=If_peak*1000; #Peak current through the diode in mA\n",
+ "#(2)\n",
+ "V_out_peak =If_peak * RL/1000 ; #Peak output voltage in V\n",
+ "\n",
+ "#For an Ideal diode\n",
+ "If_peak_ideal=Vf/RL; #Peak current through the ideal diode in A\n",
+ "If_peak_ideal=If_peak_ideal*1000; #Peak current through the ideal diode in mA\n",
+ "\n",
+ "V_out_peak_ideal=If_peak_ideal * RL/1000; # Peak output voltage in case of the ideal diode in V\n",
+ "\n",
+ "#Result\n",
+ "print '(i) Peak current through the diode = %.1f mA '%If_peak;\n",
+ "print '(ii) Peak output voltage = %.1f V'%V_out_peak;\n",
+ "print '(iii) Peak current through the ideal diode = %d mA '%If_peak_ideal;\n",
+ "print '(iv) Peak output voltage in case of the ideal diode = %d V'%V_out_peak_ideal;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "(i) Peak current through the diode = 37.8 mA \n",
+ "(ii) Peak output voltage = 18.9 V\n",
+ "(iii) Peak current through the ideal diode = 40 mA \n",
+ "(iv) Peak output voltage in case of the ideal diode = 20 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 23
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.3, Page number 82"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "V =10.0; #Battery voltage in V\n",
+ "R1=50.0; #Resistor 1's resistance in ohms\n",
+ "R2=5.0; #Resistor 2's resistance in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "#Using Thevenin's Theorem to find current in the diode\n",
+ "E0=(R2/(R1+R2))*V; #Thevenin's Voltage in V\n",
+ "R0=(R1*R2)/(R1+R2); #Thevenin's Resistance in ohms\n",
+ "\n",
+ "I0=E0/R0; #Current through the diode in A\n",
+ "I0=I0*1000; #Current through the diode in mA\n",
+ "\n",
+ "#Result\n",
+ "print 'Current through the diode = %d mA '%Io;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Current through the diode = 200 mA \n"
+ ]
+ }
+ ],
+ "prompt_number": 14
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.4, Page number 82-83 "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "V =10.0; #Battery voltage in V\n",
+ "R0=48.0; #Resistance of the resistor in ohms\n",
+ "Rd=1.0; #Forward resistance of the diodes in ohms\n",
+ "Vd=0.7; #Potential barrier of the diodes in V\n",
+ "#Calculation\n",
+ "V_net=V-Vd-Vd; #Net voltage in the circuit in V\n",
+ "R_net=R0+Rd+Rd #Net resistance of the circuit in ohms\n",
+ "I_net=V_net/R_net; #Net current in the circuit in A\n",
+ "I_net=I_net*1000; #Net current in mA\n",
+ "\n",
+ "#Result\n",
+ "print 'Net current in the circuit = %d mA '%I_net;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Net current in the circuit = 172 mA \n"
+ ]
+ }
+ ],
+ "prompt_number": 15
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.5, Page number 83"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "E1=24; #Voltage of first source in V\n",
+ "E2=4; #Voltage of second source in V\n",
+ "V0=0.7; #Potential barrier of diodes in V\n",
+ "R=2000; #Resistance of the given resistor in ohms\n",
+ "Rd=0; #Forward resistance of the diodes in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "I=(E1-E2-V0)/(R+Rd); #Current in the circuit in A\n",
+ "I=I*1000; #Current in the circuit in mA \n",
+ "\n",
+ "#Result\n",
+ "print 'Current in the circuit = %.2f mA '%I;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Current in the circuit = 9.65 mA \n"
+ ]
+ }
+ ],
+ "prompt_number": 18
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.6, Page number 83-84"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "V=20; #Voltage of source in V\n",
+ "V0=0.3; #Potential barrier of Germanium diode in V\n",
+ "V0_Si=0.7; #Potetial barrier of Silicon diode in V \n",
+ "\n",
+ "#Calculation\n",
+ "#As only Ge diode is turned on due to less potential barrier,\n",
+ "VA=V-V0; #Voltage VA acroos resistor of 3k ohms\n",
+ "\n",
+ "#Result\n",
+ "print 'Voltage VA = %.1f mA '%VA;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Voltage VA = 19.7 mA \n"
+ ]
+ }
+ ],
+ "prompt_number": 20
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.7, Page number 84"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "V=10; #Voltage of source in V\n",
+ "V0=0.7; #Potetial barrier of Silicon diode in V \n",
+ "# Resistance of all resistors in ohms\n",
+ "R1=2000;\n",
+ "R2=2000;\n",
+ "R3=2000;\n",
+ "\n",
+ "#Calculation\n",
+ "Id=(V-V0)/(R2+2*R3); #Current through the diodes in A\n",
+ "VQ=2*Id*R3; #Voltage VQ across the grounded 2k ohm resistor in V\n",
+ "Id=Id*1000; #Current through the diodes in mA\n",
+ "\n",
+ "#Result\n",
+ "print 'Voltage VQ = %.1f V '%VQ;\n",
+ "print 'Current through the diodes, Id = %.2f mA '%Id;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Voltage VQ = 6.2 V \n",
+ "Current through the diodes, Id = 1.55 mA \n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.8, Page number 84"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "V=15; #Voltage of source in V\n",
+ "V0=0.7; #Potetial barrier of Silicon diode in V \n",
+ "R=500 # Resistance of all resistors in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "I1=(V-V0)/R; #total current in the circuit in A\n",
+ "Id1=I1/2; #current in first diode in A\n",
+ "Id1=Id1*1000; #current in first diode in mA\n",
+ "Id2=Id1 #current in second diode in mA\n",
+ "\n",
+ "#Result\n",
+ "print ('Current in first diode = %.1f mA'%Id1);\n",
+ "print ('Current in second diode = %.1f mA'%Id2);\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Current in first diode = 14.3 mA\n",
+ "Current in second diode = 14.3 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.9, Page number 85"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "E=20; #Voltage of source in V\n",
+ "V0_d1=0.7; #Potetial barrier of first Silicon diode in V\n",
+ "V0_d2=0.7; #Potetial barrier of second Silicon diode in V\n",
+ "R1=5600; # Resistance of first resistor in ohms\n",
+ "R2=3300; # Resistance of second resistor in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "I2=V0_d2/R2; #Current I2 through resistor R2 in A\n",
+ "I2=round((I2*1000),3); #Current I2 through resistor R2 in mA\n",
+ "I1=(E-V0_d1-V0_d2)/R1; #Current I1 through resistor R1 in A\n",
+ "I1=round((I1*1000),2); #Current I1 through resistor R1 in mA\n",
+ "I3=I1-I2; #Current I3 through diode D2 in mA\n",
+ "\n",
+ "#Result\n",
+ "print 'Current I1= %.2f mA'%I1;\n",
+ "print 'Current I1= %.3f mA'%I2;\n",
+ "print 'Current I1= %.3f mA'%I3;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Current I1= 3.32 mA\n",
+ "Current I1= 0.212 mA\n",
+ "Current I1= 3.108 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.10, Page number 85-86"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable Declaration\n",
+ "E=10.0; #Voltage of source in V\n",
+ "V0=0.7; #Potetial barrier of Silicon diode in V\n",
+ "R1=2000; # Resistance of first resistor in ohms\n",
+ "R2=8000; # Resistance of second resistor in ohms\n",
+ "R3=4000; #Resistance of third resistor in ohms\n",
+ "R4=6000; #Resistance of fourth resistor in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "#Assuming the given diode to be reverse bised and calculating voltage across it's terminals\n",
+ "V1=(E/(R1+R2))*R2; #voltage at the P side of the diode, i.e, voltage across R2 resistor,according to voltage divider rule, in V\n",
+ "V2=(E/(R3+R4))*R4; #voltage at the N side of the diode, i.e, voltage across R4 resistor,according to voltage divider rule, in V\n",
+ "\n",
+ "#Result\n",
+ "if((V1-V2)>=V0):\n",
+ " print 'Our assumption was wrong and, the diode is forward biased';\n",
+ "else:\n",
+ " print 'The diode is reverse biased';\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Our assumption was wrong and, the diode is forward biased\n"
+ ]
+ }
+ ],
+ "prompt_number": 3
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.11, Page number 86"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "V=2; #Supply voltage in V\n",
+ "V0=0.7; #Potential barrier voltage of the diode in V \n",
+ "R1=4000.0; #Resistance of first resistor in \u03a9\n",
+ "R2=1000.0; ##Resistance of second resistor in \u03a9\n",
+ "\n",
+ "#Calculation\n",
+ "#Assuming the diode to be in ON state\n",
+ "I1=((V-V0)/R1)*1000; #Current through resistor R1, in mA\n",
+ "I2=(V0/R2)*1000; #Current through resistor R2, in mA\n",
+ "ID=I1-I2; #Diode current, in mA\n",
+ "\n",
+ "if(ID<0):\n",
+ " #Since the diode current is negative, the diode must be OFF \n",
+ " ID=0; #True value of diode current, mA\n",
+ " \n",
+ "#As the diode is in OFF state it can be replaced by an open ciruit equivalent \n",
+ "VD=V*R2/(R1 +R2); #Voltage across the diode, in V\n",
+ "\n",
+ "#Result\n",
+ "print 'ID =%d mA'%ID;\n",
+ "print 'VD =%.1f V'%VD;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "ID =0 mA\n",
+ "VD =0.4 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 9
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.12, Page number 89-90"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "AC_Input_Power=100.0; #Input AC Power in watts\n",
+ "AC_Output_Power=40.0; #Output AC Power in watts\n",
+ "Accepted_Power=50.0; #Power accepted by the half-wave rectifier in watt\n",
+ "\n",
+ "#Calculation\n",
+ "R_eff=(AC_Output_Power/AC_Input_Power)*100; #Rectification efficiency of the half-wave rectifier\n",
+ "Unused_power=AC_Input_Power-Accepted_Power; #Power not used by the half_wave rectifier due to open circuited condition of the diode in watt\n",
+ "Power_dissipated=Accepted_Power-AC_Output_Power; #Power dissipated by the diode watt\n",
+ "\n",
+ "#Result\n",
+ "print 'The rectification efficiency of the half-wave rectifier= %d%% '%R_eff;\n",
+ "\n",
+ "print 'Rest 60%% of the power is the unused power and power dissipated by the diode = %d watts and %d watts' %(Unused_power ,Power_dissipated);\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The rectification efficiency of the half-wave rectifier= 40% \n",
+ "Rest 60% of the power is the unused power and power dissipated by the diode = 50 watts and 10 watts\n"
+ ]
+ }
+ ],
+ "prompt_number": 24
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.13, Page number 90"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "from math import sqrt\n",
+ "#Variable declaration\n",
+ "Vrms=230.0; #AC supply RMS voltage in V\n",
+ "Turns_Ratio=10/1; #turn ratio of the transformer \n",
+ "\n",
+ "#Calculation\n",
+ "Vpm=sqrt(2)*Vrms; #Maximum primary voltage in V\n",
+ "Vsm=Vpm/Turns_Ratio; #Maximum secondary voltage in V\n",
+ "#Case 1\n",
+ "Vdc=Vsm/(round(pi,2)); #Output D.C voltage, which is the average voltage in V\n",
+ "Vdc=round(Vdc,2);\n",
+ "#Case 2\n",
+ "PIV=Vsm; #Peak Inverse Voltage in V\n",
+ "\n",
+ "#Result\n",
+ "print 'The output d.c voltage= %.2f V'%Vdc;\n",
+ "print 'The peak inverse voltage= %.2f V'%PIV;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The output d.c voltage= 10.36 V\n",
+ "The peak inverse voltage= 32.53 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 27
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.14, Page number 90-91"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "#Variable declaration\n",
+ "rf=20.0; #Internal resistance of the crystal diode in ohms\n",
+ "Vm=50.0; #Maximum applied voltage in V\n",
+ "RL=800.0; #Load Resistance in ohms\n",
+ "\n",
+ "#Calculation\n",
+ "# 1\n",
+ "Im=Vm/(rf+RL); #Maximum current in A\n",
+ "Im=Im*1000; #Maximum current in \n",
+ "Im=round(Im,0);\n",
+ "Idc=Im/pi; #Average voltage in mA\n",
+ "Idc=round(Idc,1);\n",
+ "Irms=Im/2; #RMS value of the current in mA\n",
+ "Irms=round(Irms,1)\n",
+ "\n",
+ "# 2\n",
+ "AC_Input_Power=pow(Irms/1000,2)*(rf+RL); #Input a.c power in watt\n",
+ "\n",
+ "DC_Output_Power=pow(Idc/1000,2)*RL; #Output d.c power in watt\n",
+ "\n",
+ "# 3\n",
+ "DC_Output_Voltage=(Idc/1000)*RL; #Output d.c voltage in V\n",
+ "\n",
+ "# 4\n",
+ "Rectifier_efficiency=(DC_Output_Power/AC_Input_Power)*100; # Efficiency of rectification of the half-wave rectifier\n",
+ "\n",
+ "#Result\n",
+ "print ' i:';\n",
+ "print ' Im = %d mA'%Im;\n",
+ "print ' Idc = %.1f mA'%Idc;\n",
+ "print ' Irms = %.1f mA'%Irms;\n",
+ "print ' ii: ';\n",
+ "print ' a.c input power= %.3f watt'%AC_Input_Power;\n",
+ "print ' d.c output power= %.3f watt'%DC_Output_Power;\n",
+ "print ' iii: ';\n",
+ "print ' d.c output voltage = %.2f volts'%DC_Output_Voltage;\n",
+ "print ' iv: '\n",
+ "print ' Efficiency of rectification = %.1f%%'%Rectifier_efficiency;\n",
+ "\n",
+ " "
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " i:\n",
+ " Im = 61 mA\n",
+ " Idc = 19.4 mA\n",
+ " Irms = 30.5 mA\n",
+ " ii: \n",
+ " a.c input power= 0.763 watt\n",
+ " d.c output power= 0.301 watt\n",
+ " iii: \n",
+ " d.c output voltage = 15.52 volts\n",
+ " iv: \n",
+ " Efficiency of rectification = 39.5%\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.15, Page number 91"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "#Variable declaration\n",
+ "Vdc=50.0; #Output d.c voltage in V\n",
+ "rf=25; #Diode resistance in ohm\n",
+ "RL=800; #Load resistance in ohm\n",
+ "\n",
+ "\n",
+ "#Calculation\n",
+ "Vm=(pi*(rf+RL)*Vdc)/RL; #[ Vdc=Vm*RL/(pi*(rf+RL)) ]Maximum value of a.c voltage required to get a volatge of Vdc from the half-wave rectifier, in V\n",
+ "Vm=round(Vm,0); \n",
+ "#Result\n",
+ "print 'The a.c voltage required should have maximum value of = %d V' %Vm;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The a.c voltage required should have maximum value of = 162 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 18
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.16, Page number 95"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt \n",
+ "from math import pi\n",
+ "#Variable declaration\n",
+ "rf=20; #Internal resistance of the diodes in ohm\n",
+ "Vrms=50; #RMS value of transformer's secondary voltage from centre tap to each end of secondary\n",
+ "RL=980; #Load resistance in ohm\n",
+ "\n",
+ "#Calculation\n",
+ "Vm=Vrms*sqrt(2); #Maximum a.c voltage in V\n",
+ "Im=Vm/(rf+RL); #Maximum load current in A\n",
+ "Im=Im*1000; #Maximum load current in mA\n",
+ " \n",
+ "# 1:\n",
+ "Idc=2*Im/pi; #Mean load current\n",
+ "\n",
+ "# 2:\n",
+ "Irms=Im/sqrt(2); #RMS value of load current in A\n",
+ "\n",
+ "#Result\n",
+ "print 'i:';\n",
+ "print' The mean load current= %d mA'%Idc;\n",
+ "print 'ii:';\n",
+ "print ' The r.m.s value of the load current = %d mA'%Irms; "
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i:\n",
+ " The mean load current= 45 mA\n",
+ "ii:\n",
+ " The r.m.s value of the load current = 50 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 20
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.17, Page number 95-96"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt \n",
+ "#Variable declaration\n",
+ "RL=100; #Load resistance in ohm \n",
+ "rf=0; #Internal resistance of the diodes in ohm\n",
+ "Turns_ratio=5/1; #Primary to secondary turns ratio of transformer \n",
+ "P_Vrms=230; #R.M.S value of voltage in primary winding in V\n",
+ "S_Vrms=P_Vrms/Turns_ratio; #R.M.S value of voltage in secondary winding in V\n",
+ "S_Vm=S_Vrms*sqrt(2); #Maximum voltage across secondary winding in V\n",
+ "Vm=S_Vm/2; #Maximum voltage across half seconfdary winding in V\n",
+ "\n",
+ "\n",
+ "#Calculation\n",
+ "# 1:\n",
+ "Idc=2*Vm/(pi*RL); #Average current in A\n",
+ "Vdc=Idc*RL; #d.c output voltage in V\n",
+ "\n",
+ "# 2:\n",
+ "PIV=S_Vm; #Peak Invers Voltage(= Maximum secondary voltage) in V\n",
+ "\n",
+ "# 3:\n",
+ "Pac=pow(Vm/(RL*sqrt(2)),2)*(rf+RL); #a.c input power in watt\n",
+ "Pdc=(pow(Idc,2)*RL); #d.c output power in watt\n",
+ "R_eff=(Pdc/Pac)*100; #Rectification efficiency\n",
+ "R_eff=round(R_eff,1);\n",
+ "\n",
+ "#Result\n",
+ "print 'i:';\n",
+ "print ' The d.c output voltage= %.1f V'%Vdc;\n",
+ "print 'ii:';\n",
+ "print ' The peak inverse voltage= %d V'%PIV;\n",
+ "print 'iii:';\n",
+ "print ' Rectification efficiency= %.1f%%'%R_eff;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i:\n",
+ " The d.c output voltage= 20.7 V\n",
+ "ii:\n",
+ " The peak inverse voltage= 65 V\n",
+ "iii:\n",
+ " Rectification efficiency= 81.1%\n"
+ ]
+ }
+ ],
+ "prompt_number": 27
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "NOTE: The value of rectification efficiency is calculated as 81.2% in the textbook using the formula 0.812/(1 + (rf/RL)), but by calculating using the correct values in the formula we get 81.1%."
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.18, Page number 96-97"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt \n",
+ "#Variable declaration\n",
+ "fin=50; #frequency of input ac source in Hz\n",
+ "RL=200; #Load resistance in ohm\n",
+ "Turns_ratio=4/1; #Transformers turns ratio, primary to secondary.\n",
+ "P_Vrms=230.0; #R.M.S value of voltage in primary winding in V\n",
+ "S_Vrms=P_Vrms/Turns_ratio #R.M.S value of voltage in secondary winding in V\n",
+ "Vm=S_Vrms*sqrt(2); #Maximum voltage across secondary winding in V\n",
+ "\n",
+ "#Calculation\n",
+ "# 1:\n",
+ "Idc=2*Vm/(pi*RL); # Average current in A\n",
+ "Vdc=Idc*RL; #Output d.c voltage in V\n",
+ "Vdc=round(Vdc,0);\n",
+ "# 2:\n",
+ "PIV= Vm; #Peak Inverse Voltage(= Maximum volutage across secondary winding) in V\n",
+ "\n",
+ "# 3:\n",
+ "fout=2*fin; #Output frequency in Hz\n",
+ "\n",
+ "#Result\n",
+ "print 'i:';\n",
+ "print ' The d.c output voltage = %d V' %Vdc;\n",
+ "print 'ii:';\n",
+ "print ' The peak inverse voltage = %.1f V'%PIV;\n",
+ "print 'iii:';\n",
+ "print ' The output frequency = %d Hz'%fout;\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i:\n",
+ " The d.c output voltage = 52 V\n",
+ "ii:\n",
+ " The peak inverse voltage = 81.3 V\n",
+ "iii:\n",
+ " The output frequency = 100 Hz\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.19, Page number 97"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "from math import sqrt\n",
+ "\n",
+ "#Variable declaration\n",
+ "RL=100.0; #Load Resistance in ohm\n",
+ "Turns_ratio=5/1; #Primary to secondary turns ratio of the transformer\n",
+ "Vin=230.0; #R.M.S value of input voltage in V\n",
+ "fin=50; #Input frequency in Hz\n",
+ "\n",
+ "#Calculation\n",
+ "Vs_rms=Vin/Turns_ratio; #R.M.S value of the voltage in secondary winding, in v\n",
+ "Vs_max=Vs_rms*sqrt(2); #Maximum voltage across secondary, in V\n",
+ "\n",
+ "# (i)\n",
+ "#Case i: Centre-tap circuit\n",
+ "Vm=Vs_max/2; #Maximum voltage across half secondary winding, in V \n",
+ "Vdc=2*Vm*RL/(pi*RL); #DC output voltage, in V \n",
+ "print 'The d.c output voltage for the centre-tap circuit = %.1f V'%Vdc;\n",
+ "\n",
+ "#Case ii:\n",
+ "Vm=Vs_max; #Maximum voltage across secondary, in V\n",
+ "Vdc=2*Vm*RL/(pi*RL); #DC output voltage, in V \n",
+ "print 'The d.c output voltage for the bridge circuit = %.1f V'%Vdc; \n",
+ "\n",
+ "# ii:\n",
+ "#Case i: Centre-tap circuit\n",
+ "Turns_ratio=5/1;\n",
+ "Vs_rms=Vin/Turns_ratio;\n",
+ "Vs_max=Vs_rms*sqrt(2);\n",
+ "Vm=Vs_max/2;\n",
+ "PIV=2*Vm;\n",
+ "print 'PIV in case of centre-tap circuit = %d V'%PIV;\n",
+ "\n",
+ "#Case ii: Bridge circuit\n",
+ "Turns_ratio=10/1;\n",
+ "Vs_rms=Vin/Turns_ratio;\n",
+ "Vs_max=Vs_rms*sqrt(2);\n",
+ "PIV=Vm;\n",
+ "print 'PIV in case of bridge circuit = %.1f V'%PIV;\n",
+ "\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The d.c output voltage for the centre-tap circuit = 20.7 V\n",
+ "The d.c output voltage for the bridge circuit = 41.4 V\n",
+ "PIV in case of centre-tap circuit = 65 V\n",
+ "PIV in case of bridge circuit = 32.5 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 44
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.20, Page number 98"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "from math import sqrt\n",
+ "#Variable declaration\n",
+ "rf=1; #forward resistance of diodes of the rectifier in ohm\n",
+ "RL=480; #Load resistance in ohm\n",
+ "Vrms=240.0; #a.c supply voltage in V\n",
+ "Vm=Vrms*sqrt(2); #Maximum a.c voltage in V \n",
+ "\n",
+ "#Calculation\n",
+ "# 1:\n",
+ "Rt=2*rf+RL; #Total circuit resistance at any instance in ohm\n",
+ "Im=Vm/Rt; #Maximum load current in A\n",
+ "Idc=2*Im/pi; #Mean load current in A\n",
+ "\n",
+ "# 2:\n",
+ "Irms=Im/2; #R.M.S value of current in A\n",
+ "P=pow(Irms,2)*rf; #Power dissipated in each diode in watt\n",
+ "\n",
+ "\n",
+ "#Result\n",
+ "print 'i:';\n",
+ "print ' Mean load current = %.2f A'%Idc;\n",
+ "print 'ii:';\n",
+ "print ' Power dissipated in each diode= %.3f W'%P;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i:\n",
+ " Mean load current = 0.45 A\n",
+ "ii:\n",
+ " Power dissipated in each diode= 0.124 W\n"
+ ]
+ }
+ ],
+ "prompt_number": 39
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "NOTE: The value of power dissipated is approximately 0.124 W , but in the textbook it is approximated as 0.123W."
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.21, Page number 98-99"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt,pi\n",
+ "#Variable declaration\n",
+ "RL=12000; #Load resistance in ohm\n",
+ "V0=0.7; #Potential barrier voltage of diodes in V\n",
+ "Vrms=12; #R.M.S value of input a.c voltage in V\n",
+ "Vs_pk=Vrms*sqrt(2); #Peak secondary voltage in V\n",
+ "\n",
+ "#Calculation\n",
+ "# 1:\n",
+ "Vout_pk=Vs_pk-(2*V0); #Peak output voltage in V\n",
+ "Vav=2*Vout_pk/pi; #Average output voltage in V\n",
+ "Vav=round(Vav,2);\n",
+ "\n",
+ "# 2:\n",
+ "Iav=Vav/RL; #Average output current in A\n",
+ "Iav=Iav*pow(10,6); #Average output current in \u03bcA\n",
+ "\n",
+ "\n",
+ "#Result\n",
+ "print 'i:';\n",
+ "print ' Average output voltage=%.2f V'%Vav;\n",
+ "print 'ii:';\n",
+ "print ' Average output current=%.1f \u03bcA'%Iav;\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i:\n",
+ " Average output voltage=9.91 V\n",
+ "ii:\n",
+ " Average output current=825.8 \u03bcA\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.22, Page number 102"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Vdc_A=10; #Supply voltage of A in V\n",
+ "Vdc_B=25; #Supply voltage of B in V\n",
+ "Vac_rms_a=0.5; #Ripples in power supply A in V\n",
+ "Vac_rms_b=0.001; #Ripples in power supply B in V\n",
+ "\n",
+ "#Calculation\n",
+ "#For power supply A\n",
+ "ripple_factor_A=Vac_rms_a/Vdc_A; #Ripple factor of power supply A\n",
+ "\n",
+ "#For power supply B\n",
+ "ripple_factor_B=Vac_rms_b/Vdc_B; #Ripple factor of power supply B\n",
+ "\n",
+ "#Result\n",
+ "if(ripple_factor_A<ripple_factor_B):\n",
+ " print 'Power supply A is better';\n",
+ "else :\n",
+ " print 'Power supply B is better';"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Power supply B is better\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.23, Page number 105-106"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import sqrt\n",
+ "#Variable declaration\n",
+ "RL=2200; #Load resistance in ohm\n",
+ "C=50*pow(10,-6); #Capacitance of the capacitor used in filter circuit in F\n",
+ "V0=0.7; #Potential barrier voltage of the diodes of the rectifier in V\n",
+ "Vrms=115.0; #R.M.S value of input a.c voltage in V \n",
+ "fin=60; #Frequency of input a.c voltage in Hz\n",
+ "Turns_ratio=10/1; #Primary to secondary, turns ratio of the transformer \n",
+ "\n",
+ "#Calculation\n",
+ "Vp_prim=Vrms*sqrt(2); #Peak primary voltage in V\n",
+ "Vp_sec=Vp_prim/Turns_ratio; #Peak secondary voltage in V\n",
+ "Vp_in= Vp_sec - 2*V0; #Peak full wave rectified voltage at the filter input in V\n",
+ "f=2*fin; #Output frequency in Hz\n",
+ "Vdc=Vp_in*(1-(1/(2*f*RL*C))); #Output d.c voltage in V\n",
+ "\n",
+ "#Result\n",
+ "print 'The output d.c voltage is = %.1f V'%Vdc;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The output d.c voltage is = 14.3 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.24, Page number 106"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "from math import pi\n",
+ "#Variable declaration\n",
+ "R=25; #d.c resistance of the choke in ohm\n",
+ "RL=750; #Load resistance in ohm\n",
+ "Vm=25.7; #Maximum value of the pulsating output from the rectifier in V\n",
+ "\n",
+ "#Calculation\n",
+ "V_dc=2*Vm/pi; #d.c component of the pulsating output in V\n",
+ "V_dc=round(V_dc,1);\n",
+ "V_dc_out=(V_dc*RL)/(R+RL); #Output d.c voltage in V\n",
+ "V_dc_out=round(V_dc_out,1);\n",
+ "\n",
+ "#Result\n",
+ "print ' The output d.c voltage accross the load resistance is = %.1f V'%V_dc_out;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ " The output d.c voltage accross the load resistance is = 15.9 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 28
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.25, Page number 113-114"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Ei=120.0; #Input Voltage in V\n",
+ "Vz=50.0; #Zener Voltage in V\n",
+ "R=5000.0; #Resistance of the series resistor in ohm\n",
+ "RL=10000.0; #Load resistance in ohm\n",
+ "\n",
+ "#Calculation\n",
+ "V=Ei*RL/(R+RL); #Voltage across the open circuit if the zener diode is removed\n",
+ "if(V>Vz):\n",
+ " #Zener diode is in ON state\n",
+ " # i:\n",
+ " Output_voltage=Vz; #Voltage across load resistance, in V\n",
+ " #ii:\n",
+ " Voltage_R=Ei-Vz; #Voltage across the series resistance R, in V\n",
+ " #iii:\n",
+ " IL=Vz/RL; #Load current through RL in A\n",
+ " IL=IL*1000; #Load current through RL in mA\n",
+ " I=Voltage_R/R; #Current through the series resistance in A\n",
+ " I=I*1000; #Current through the series resistance in mA\n",
+ " Iz=I-IL; #Applying Kirchhoff's first law, Zener current in mA\n",
+ " \n",
+ " #Result\n",
+ " print 'i) The output voltage across the load resistance RL = %d V'%Output_voltage;\n",
+ " print 'ii) The voltage drop across the series resistance R = %d V'%Voltage_R;\n",
+ " print 'iii) The current through the zener diode = %d mA'%Iz;\n",
+ " "
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i) The output voltage across the load resistance RL = 50 V\n",
+ "ii) The voltage drop across the series resistance R = 70 V\n",
+ "iii) The current through the zener diode = 9 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.26, Page number 114-115"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Max_V=120.0; #Maximum input voltage in V\n",
+ "Min_V=80.0; #Minimum input voltage in V\n",
+ "R=5000.0; #Series resistance in ohm\n",
+ "RL=10000.0; #Load resistance in ohm\n",
+ "Vz=50.0; #Zener voltage in V\n",
+ "\n",
+ "\n",
+ "#Calculation\n",
+ "#Case i: Maximum zener current\n",
+ "#Zener current will be maximum when the input voltage is maximum\n",
+ "V_R_max=Max_V-Vz; #Voltage across series resistance R, in V\n",
+ "I_max=V_R_max/R; #Current through series resistance R, in A\n",
+ "I_max=I_max*1000; #Current through series resistance R, in mA\n",
+ "IL_max=Vz/RL; #Load current in A\n",
+ "IL_max=IL_max*1000; #Load current in mA\n",
+ "Iz_max=I_max-IL_max; #Applying Kirchhoff's first law, Zener current in mA;\n",
+ "\n",
+ "#Case ii: Minimum zener current\n",
+ "#The zener will conduct minimum current when the input voltage is minimum\n",
+ "V_R_min=Min_V-Vz; #Voltage across series resistance R, in V\n",
+ "I_min=V_R_min/R; #Current through series resistance R, in A\n",
+ "I_min=I_min*1000; #Current through series resistance R, in mA\n",
+ "IL_min=Vz/RL; #Load current in A\n",
+ "IL_min=IL_min*1000; #Load current in mA\n",
+ "Iz_min=I_min-IL_min; #Applying Kirchhoff's first law, Zener current in mA\n",
+ "\n",
+ "#Result\n",
+ "print 'Case i: ';\n",
+ "print 'Maximum zener current = %d mA'%Iz_max;\n",
+ "print 'Case ii: ';\n",
+ "print 'Minimum zener current = %d mA'%Iz_min;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Case i: \n",
+ "Maximum zener current = 9 mA\n",
+ "Case ii: \n",
+ "Minimum zener current = 1 mA\n"
+ ]
+ }
+ ],
+ "prompt_number": 11
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.27, Page number 115"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Ei=12; #Input voltage in V\n",
+ "Vz=7.2; #Zener voltage in V\n",
+ "E0=Vz; #Voltage to be maintained across the load in V\n",
+ "IL_max=0.1; #Maximum load current in A\n",
+ "IL_min=0.012; #Minimum load current in A\n",
+ "Iz_min=0.01; #Minimum zener current in A\n",
+ "\n",
+ "#Calculation\n",
+ "#When the load current is maximum at minimum value of RL, the zener current is minimum and, as the load current decreases due to increase in value of RL\n",
+ "R=(Ei-E0)/(Iz_min+IL_max); #The value of series resistance R to maintain a voltage=E0 across load, in ohm\n",
+ "\n",
+ "#Result\n",
+ "print 'The minimum value of series resistance R to maintain a constant value of 7.2 V is = %.1f \u03a9'%R;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The minimum value of series resistance R to maintain a constant value of 7.2 V is = 43.6 \u03a9\n"
+ ]
+ }
+ ],
+ "prompt_number": 16
+ },
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "NOTE: The actual value of R is 43.636363 (recurring) but, in the textbook the value of R is wrongly approximated 43.5 \u03a9"
+ ]
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.28, Page number 115"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Ei_min=22; #Minimum input voltage in V\n",
+ "Ei_max=28; #Maximum input voltage in V\n",
+ "Vz=18; #Zener voltage in V\n",
+ "E0=Vz; #Constant voltage maintained across the load resistance in V\n",
+ "Iz_min=0.2; #Minimum zener current in A\n",
+ "Iz_max=2; #Maximum zener current in A\n",
+ "RL=18; #Load resistance in \u03a9\n",
+ "\n",
+ "#Calculation\n",
+ "IL=Vz/RL; #Constant value of load current in A\n",
+ "#When the input voltage is minimum, the zener current will be minimum\n",
+ "R=(Ei_min-E0)/(Iz_min+IL) #The value of series resistance so that the voltage E0 across RL remains constant\n",
+ "\n",
+ "print 'The value of series resistance R, to maintain constant voltage E0 across RL = %.2f \u03a9.'%R;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The value of series resistance R, to maintain constant voltage E0 across RL = 3.33 \u03a9.\n"
+ ]
+ }
+ ],
+ "prompt_number": 19
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.29, Page number 116 "
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Vz=10 #Zener voltage in V\n",
+ "Ei_min=13; #Minimum input voltage in V\n",
+ "Ei_max=16; #Maximum input voltage in V\n",
+ "Iz_min=0.015; #Minimum zener current in A\n",
+ "IL_min=0.01; #Minimum load current in A \n",
+ "IL_max=0.085; #Maximum load curremt in A\n",
+ "E0=Vz; #Constant voltage to be maintained in V \n",
+ "\n",
+ "#Calculation\n",
+ "#The zener current will be minimum when the input voltage will be minimum and at that time the load current will be maximum\n",
+ "R=(Ei_min-E0)/(Iz_min+IL_max); #The value of series resistance R to maintain a constant voltage across load\n",
+ "\n",
+ "\n",
+ "#Result\n",
+ "print 'The value of series resistance to maintain a constant voltage across the load resistance is = %d \u03a9'%R;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The value of series resistance to maintain a constant voltage across the load resistance is = 30 \u03a9\n"
+ ]
+ }
+ ],
+ "prompt_number": 23
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.30, Page number 116"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Iz=0.2; #Current rating of each zener in A\n",
+ "Vz=15; #Voltage rating of each zener in V\n",
+ "Ei=45; #Input voltage in V\n",
+ "\n",
+ "#Calculation\n",
+ "# i: Regulated output voltage across the two zener diodes \n",
+ "E0=2*Vz; # V\n",
+ "\n",
+ "# ii: Value of series resistance \n",
+ "R=(Ei-E0)/Iz; # \u03a9\n",
+ "\n",
+ "#Result\n",
+ "print 'i) The regulated output voltage = %d V'%E0;\n",
+ "print 'ii) The value of the series resistance = %d \u03a9'%R;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "i) The regulated output voltage = 30 V\n",
+ "ii) The value of the series resistance = 75 \u03a9\n"
+ ]
+ }
+ ],
+ "prompt_number": 28
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.31, Page number 116-117"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Vz=10; #Voltage rating of each zener in V\n",
+ "Iz=1; #Current rating of each zener in A\n",
+ "Ei=45; #Input unregulated voltage in V\n",
+ "\n",
+ "#Calculation\n",
+ "#Regulated output voltage across the three zener diodes\n",
+ "E0=3*Vz; # V\n",
+ "\n",
+ "#Value of series resistance to obtain a 30V regulated output voltage\n",
+ "R=(Ei-E0)/Iz; # \u03a9\n",
+ "\n",
+ "#Result\n",
+ "print 'Value of series resistance to obtain a 30V regulated output voltage = %d \u03a9'%R;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Value of series resistance to obtain a 30V regulated output voltage = 15 \u03a9\n"
+ ]
+ }
+ ],
+ "prompt_number": 29
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.32, Page number 117"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#variable declaration\n",
+ "RL=2000.0; #Load resistance in \u03a9\n",
+ "R=200.0; #Series resistance in \u03a9\n",
+ "Iz=0.025; #Zener current rating in A\n",
+ "E0=30.0; #Output regulated voltage in V \n",
+ "\n",
+ "#Calculation\n",
+ "#Minimum input voltage will be required when Iz=0 A, and at this condition\n",
+ "IL=E0/RL; #Load current during Iz=0, in A\n",
+ "I=IL; #According to Kirchhoff's law, total current, in A\n",
+ "Ei_min=E0+(I*R); #Minimum input voltage in V\n",
+ "\n",
+ "#The maximum input voltage required will be when Iz=0.025 A, and at that condition \n",
+ "I=IL+Iz; #According to Kirchhoff's law, total current, in A\n",
+ "Ei_max=E0+(I*R); #maximum input voltage in V\n",
+ "\n",
+ "\n",
+ "#Result\n",
+ "print 'The required range of input voltage is from %d V to %d V'%(Ei_min,Ei_max); \n",
+ "\n",
+ "\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The required range of input voltage is from 33 V to 38 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 2
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.33, Page number 117-118"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "Ei=16; #Unregulated input voltage in V\n",
+ "E0=12; #Output regulated voltage in V\n",
+ "IL_min=0; #Minimum load current in A\n",
+ "IL_max=0.2; #Maximum load current in A\n",
+ "Iz_min=0; #Minimum zener current in A\n",
+ "Iz_max=0.2; #Maximum zener current in A\n",
+ "\n",
+ "#Calculation\n",
+ "#As the regulated voltage required across the load is 12V\n",
+ "Vz=E0; #Voltage rating of zener diode in V\n",
+ "V_R=Ei-E0; #Constant Voltage that should remain across series resistance \n",
+ "#The minimum zener current will occur when the curent in the load in maximum\n",
+ "R=V_R/(Iz_min+IL_max); #Series resistance in \u03a9\n",
+ "\n",
+ "Max_power_rating=Vz*Iz_max; #Maximum power rating of zener diode in W\n",
+ "\n",
+ "#Result\n",
+ "print 'The regulator is designed using a Seris resistance of %d \u03a9 and a zener diode of zener voltage %d V'%(R,Vz);\n",
+ "print 'The maximum power rating of the zener diode is = %.1f W '%Max_power_rating;"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "The regulator is designed using a Seris resistance of 20 \u03a9 and a zener diode of zener voltage 12 V\n",
+ "The maximum power rating of the zener diode is = 2.4 W \n"
+ ]
+ }
+ ],
+ "prompt_number": 1
+ },
+ {
+ "cell_type": "heading",
+ "level": 2,
+ "metadata": {},
+ "source": [
+ "Example 6.34, Page number 118"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [
+ "#Variable declaration\n",
+ "V=12; #Source voltage in V\n",
+ "R=1000; #Series resistance in \u03a9\n",
+ "RL=5000; #Load resistance in \u03a9\n",
+ "Vz=6; #Voltage rating of zener in V\n",
+ "\n",
+ "#Calculation\n",
+ "#Case i: zener is working properly\n",
+ "#The output voltage across the load will be equal to the zener voltage.\n",
+ "V0=Vz; # V\n",
+ "\n",
+ "#Result\n",
+ "print 'Case i: Output voltage when zener is working properly is %d V'%V0;\n",
+ "\n",
+ "#Case ii: zener is shorted\n",
+ "#As the zener is shorted, the potential difference across the load will be zero\n",
+ "V0=0; #V\n",
+ "\n",
+ "#Result\n",
+ "print 'Case ii: Output voltage when zener is short circuited is %d V'%V0;\n",
+ " \n",
+ "#Case iii: zener is open circuited\n",
+ "#If the zener is open circuited, the total voltage will drop across R and RL according to the voltage divider rule\n",
+ "V0=V*RL/(R+RL); #V\n",
+ "\n",
+ "#Result\n",
+ "print 'Case iii: Output voltage when zener is open circuited is %d V'%V0;\n"
+ ],
+ "language": "python",
+ "metadata": {},
+ "outputs": [
+ {
+ "output_type": "stream",
+ "stream": "stdout",
+ "text": [
+ "Case i: Output voltage when zener is working properly is 6 V\n",
+ "Case ii: Output voltage when zener is short circuited is 0 V\n",
+ "Case iii: Output voltage when zener is open circuited is 10 V\n"
+ ]
+ }
+ ],
+ "prompt_number": 4
+ },
+ {
+ "cell_type": "code",
+ "collapsed": false,
+ "input": [],
+ "language": "python",
+ "metadata": {},
+ "outputs": []
+ }
+ ],
+ "metadata": {}
+ }
+ ]
+} \ No newline at end of file