diff options
Diffstat (limited to '3014')
222 files changed, 1856 insertions, 0 deletions
diff --git a/3014/CH1/EX1.1/Ex1_1.sce b/3014/CH1/EX1.1/Ex1_1.sce new file mode 100755 index 000000000..6ad049d82 --- /dev/null +++ b/3014/CH1/EX1.1/Ex1_1.sce @@ -0,0 +1,10 @@ + clc
+//given that
+M = 6e24 // Mass of earth in Kg
+v = 3e4 // Orbital velocity of earth in m/s
+h = 6.625e-34 // Plank constant
+
+printf("Example 1.1")
+lambda = h/(M*v) // calculation of de Broglie wavelength
+
+printf("\n de Broglie wavelength of earth is %e m.\n\n\n",lambda)
diff --git a/3014/CH1/EX1.1/Ex1_1.txt b/3014/CH1/EX1.1/Ex1_1.txt new file mode 100755 index 000000000..318aed484 --- /dev/null +++ b/3014/CH1/EX1.1/Ex1_1.txt @@ -0,0 +1,2 @@ + Example 1.1
+ de Broglie wavelength of earth is 3.680556e-63 m.
\ No newline at end of file diff --git a/3014/CH1/EX1.10/Ex1_10.sce b/3014/CH1/EX1.10/Ex1_10.sce new file mode 100755 index 000000000..9103c9fe8 --- /dev/null +++ b/3014/CH1/EX1.10/Ex1_10.sce @@ -0,0 +1,14 @@ +clc
+//Given that
+E = 10 // Energy of electron in KeV
+m_e = 9.1e-31 // Mass of electron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.10")
+v = sqrt(2*E*1.6e-16/m_e) // Calculation of velocity of moving electron
+p = m_e*v //Calculation of momentum of moving electron
+lambda = h/p // calculation of de Broglie wavelength
+printf("\n velocity of electron is %0.2e m/s.",v)
+printf("\n momentum of electron is %.3e Kgm/s.",p)
+printf("\n de Broglie wavelength of electron is %.2f angstrom.\n\n\n", lambda*1e10)
+ // Answers in book are v = 5.93e6 m/s, p = 5.397e-24 kgm/s, lambda = 1.23 angstrom
+// Which is due to wrong calculation
diff --git a/3014/CH1/EX1.10/Ex1_10.txt b/3014/CH1/EX1.10/Ex1_10.txt new file mode 100755 index 000000000..9376b4733 --- /dev/null +++ b/3014/CH1/EX1.10/Ex1_10.txt @@ -0,0 +1,5 @@ +
+Example 1.10
+ velocity of electron is 5.929995e+07 m/s.
+ momentum of electron is 5.396295e-23 Kgm/s.
+ de Broglie wavelength of electron is 0.122862 angstrom.
diff --git a/3014/CH1/EX1.11/Ex1_11.sce b/3014/CH1/EX1.11/Ex1_11.sce new file mode 100755 index 000000000..b62895426 --- /dev/null +++ b/3014/CH1/EX1.11/Ex1_11.sce @@ -0,0 +1,12 @@ +
+clc
+//Given that
+lambda = 1 // de Broglie wavelength of neutron in angstrom
+m = 1.67e-27 // Mass of electron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.11")
+
+v = h/(m*lambda*1e-10) // Calculation of velocity of moving neutron
+E = 1/2*m*v^2 // Calculation of kinetic energy of moving neutron
+printf("\n velocity of neutron is %e m/s.",v)
+printf("\n Kinetic energy of neutron is %f eV.\n\n\n",E/1.6e-19)
diff --git a/3014/CH1/EX1.11/Ex1_11.txt b/3014/CH1/EX1.11/Ex1_11.txt new file mode 100755 index 000000000..5c7763352 --- /dev/null +++ b/3014/CH1/EX1.11/Ex1_11.txt @@ -0,0 +1,4 @@ +
+Example 1.11
+ velocity of neutron is 3.964072e+03 m/s.
+ Kinetic energy of neutron is 0.082007 eV.
\ No newline at end of file diff --git a/3014/CH1/EX1.12/Ex1_12.sce b/3014/CH1/EX1.12/Ex1_12.sce new file mode 100755 index 000000000..f6351694f --- /dev/null +++ b/3014/CH1/EX1.12/Ex1_12.sce @@ -0,0 +1,11 @@ +
+
+clc
+//Given that
+E = 2 // Energy of accelerated electron in KeV
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.12")
+lambda = h/sqrt(2*m*E*1e3*1.6e-19) // Calculation of velocity of moving electron
+printf("\n Wavelength of electron is %e m.\n\n\n",lambda)
+// Answer in book is 2.74e-12m
diff --git a/3014/CH1/EX1.12/Ex1_12.txt b/3014/CH1/EX1.12/Ex1_12.txt new file mode 100755 index 000000000..9c3962280 --- /dev/null +++ b/3014/CH1/EX1.12/Ex1_12.txt @@ -0,0 +1,3 @@ +
+Example 1.12
+ Wavelength of electron is 2.743136e-11 m.
diff --git a/3014/CH1/EX1.13/Ex1_13.sce b/3014/CH1/EX1.13/Ex1_13.sce new file mode 100755 index 000000000..8707b1f46 --- /dev/null +++ b/3014/CH1/EX1.13/Ex1_13.sce @@ -0,0 +1,11 @@ +
+
+clc
+//Given that
+v = 2e8 // speed of moving proton in m/s
+c = 3e8 // speed of light in m/s
+m = 1.67e-27 // Mass of proton in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.13")
+lambda = h/(m*v/sqrt(1-(v/c)^2)) // Calculation of velocity of moving electron
+printf("\n Wavelength of electron is %e angstrom.\n\n\n", lambda*1e10)
diff --git a/3014/CH1/EX1.13/Ex1_13.txt b/3014/CH1/EX1.13/Ex1_13.txt new file mode 100755 index 000000000..cb78b0f85 --- /dev/null +++ b/3014/CH1/EX1.13/Ex1_13.txt @@ -0,0 +1,3 @@ +
+Example 1.13
+ Wavelength of electron is 1.477322e-05 angstrom.
diff --git a/3014/CH1/EX1.14/Ex1_14.sce b/3014/CH1/EX1.14/Ex1_14.sce new file mode 100755 index 000000000..427d8587b --- /dev/null +++ b/3014/CH1/EX1.14/Ex1_14.sce @@ -0,0 +1,18 @@ +clc
+//given that
+lambda = 1// wavelength in m/s
+m_e = 9.1e-31 // Mass of electron in Kg
+m_p = 1.67e-27 // Mass of proton in kg
+c = 3e8 // speed of light in m/s
+h = 6.63e-34 // Plank constant
+printf("Example 1.14")
+p_p = h/(lambda*1e-10) // Momentum of photon
+p_e = h/(lambda*1e-10) // Momentum of electron
+E_e = p_e^2/(2*m_e) +m_e*c^2 // Total energy of electron
+E_p = h*c/(lambda*1e-10) // Total energy of photon
+K_e = p_e^2/(2*m_e) // Kinetic energy of electron
+K_p = h*c/(lambda*1e-10)// Kinetic energy of photon
+r_K = K_e/K_p // Ratio of kinetic energies
+printf("\n Momentum of photon is %e Kgm/s while Momentum of electron is %e Kgm/s \n which are equal.",p_p,p_e)
+printf("\n Total Energy of photon is %f KeV while Total Energy of electron is %f MeV ",E_p/(1.6e-19*1e3),E_e/(1.6e-19*1e6))
+printf("\n Ratio of kinetic energies is %e \n\n\n",r_K)
diff --git a/3014/CH1/EX1.14/Ex1_14.txt b/3014/CH1/EX1.14/Ex1_14.txt new file mode 100755 index 000000000..0b57d98b1 --- /dev/null +++ b/3014/CH1/EX1.14/Ex1_14.txt @@ -0,0 +1,6 @@ +
+Example 1.14
+ Momentum of photon is 6.630000e-24 Kgm/s while Momentum of electron is 6.630000e-24 Kgm/s
+ which are equal
+ Total Energy of photon is 12.431250 KeV while Total Energy of electron is 0.512026 MeV
+ Ratio of kinetic energies is 1.214286e-02
diff --git a/3014/CH1/EX1.15/Ex1_15.sce b/3014/CH1/EX1.15/Ex1_15.sce new file mode 100755 index 000000000..7d800a1bf --- /dev/null +++ b/3014/CH1/EX1.15/Ex1_15.sce @@ -0,0 +1,16 @@ +
+clc
+//Given that
+e = 25 // Energy of neutron in eV
+c = 3e8 // speed of light in m/s
+m = 1.67e-27 // Mass of neutron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.15")
+rest_e = m*c^2/(1e6*1.6e-19)// rest mass energy of neutron in MeV
+if e/rest_e < 0.015 then
+ E = e;
+ else E = rest_e +e;
+end
+lambda = h/(sqrt(2*m*e*1.6e-19)) // calculation of de Broglie wavelength
+printf("\n de Broglie wavelength of neutron is %f angstrom.\n\n\n", lambda*1e10)
+// Answer in book is 8.04e-5 angstrom
diff --git a/3014/CH1/EX1.15/Ex1_15.txt b/3014/CH1/EX1.15/Ex1_15.txt new file mode 100755 index 000000000..9a7a37809 --- /dev/null +++ b/3014/CH1/EX1.15/Ex1_15.txt @@ -0,0 +1,3 @@ +
+Example 1.15
+ de Broglie wavelength of neutron is 0.057274 angstrom.
diff --git a/3014/CH1/EX1.16/Ex1_16.sce b/3014/CH1/EX1.16/Ex1_16.sce new file mode 100755 index 000000000..19cbe59cb --- /dev/null +++ b/3014/CH1/EX1.16/Ex1_16.sce @@ -0,0 +1,11 @@ +clc
+//Given that
+e = 2*1.6e-19 // charge on alpha particle in coulomb
+V = 200 // Applied voltage in volts
+m = 4*1.67e-27 // Mass of alpha particle in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.16")
+
+lambda = h/(sqrt(2*e*V*m)) // calculation of de Broglie wavelength
+printf("\n de Broglie wavelength of neutron is %f angstrom.\n\n\n", lambda*1e10)
+// while answer in book is 0.00715 angstrom
diff --git a/3014/CH1/EX1.16/Ex1_16.txt b/3014/CH1/EX1.16/Ex1_16.txt new file mode 100755 index 000000000..312921424 --- /dev/null +++ b/3014/CH1/EX1.16/Ex1_16.txt @@ -0,0 +1,4 @@ +
+
+Example 1.16
+ de Broglie wavelength of neutron is 0.007170 angstrom.
diff --git a/3014/CH1/EX1.17/Ex1_17.sce b/3014/CH1/EX1.17/Ex1_17.sce new file mode 100755 index 000000000..16aafac29 --- /dev/null +++ b/3014/CH1/EX1.17/Ex1_17.sce @@ -0,0 +1,15 @@ +clc
+//Given that
+M = 20 // Mass of ball in Kg
+V = 5 // velocity of of ball in m/s
+m = 9.1e-31 //Mass of electron in Kg
+v = 1e6 // velocity of of electron in m/s
+h = 6.62e-34 // Plank constant
+
+printf("Example 1.17")
+lambda_b = h/(M*V) // calculation of de Broglie wavelength for ball
+lambda_e = h/(m*v) // calculation of de Broglie wavelength electron
+printf("\n de Broglie wavelength of ball is %e angstrom.",lambda_b*1e10)
+printf("\n de Broglie wavelength of electron is %f angstrom.\n\n\n",lambda_e*1e10)
+// answer in book is 6.62e-22 angstrom for ball
+
diff --git a/3014/CH1/EX1.17/Ex1_17.txt b/3014/CH1/EX1.17/Ex1_17.txt new file mode 100755 index 000000000..083e94630 --- /dev/null +++ b/3014/CH1/EX1.17/Ex1_17.txt @@ -0,0 +1,5 @@ +
+
+Example 1.17
+ de Broglie wavelength of ball is 6.620000e-26 angstrom.
+ de Broglie wavelength of electron is 7.274725 angstrom.
diff --git a/3014/CH1/EX1.18/Ex1_18.sce b/3014/CH1/EX1.18/Ex1_18.sce new file mode 100755 index 000000000..e51069e67 --- /dev/null +++ b/3014/CH1/EX1.18/Ex1_18.sce @@ -0,0 +1,10 @@ +
+clc
+//Given that
+E = 1 // Energy of neutron in eV
+m = 1.67e-27 // Mass of neutron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.18")
+lambda = h/sqrt(2*m*E*1.6e-19) // Calculation of velocity of moving electron
+printf("\n Wavelength of electron is %f angstrom.\n\n\n",lambda*1e10)
+// Answer in book is 6.62e-22 angstrom
diff --git a/3014/CH1/EX1.18/Ex1_18.txt b/3014/CH1/EX1.18/Ex1_18.txt new file mode 100755 index 000000000..37a460734 --- /dev/null +++ b/3014/CH1/EX1.18/Ex1_18.txt @@ -0,0 +1,4 @@ +
+
+Example 1.18
+ Wavelength of electron is 0.286368 angstrom.
diff --git a/3014/CH1/EX1.19/Ex1_19.sce b/3014/CH1/EX1.19/Ex1_19.sce new file mode 100755 index 000000000..410928a4e --- /dev/null +++ b/3014/CH1/EX1.19/Ex1_19.sce @@ -0,0 +1,11 @@ +
+clc
+//Given that
+lambda = 0.5// wavelength of electron in angstrom
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.62e-34 // Plank constant
+q = 1.6e-19 // charge on electron in coulomb
+printf("Example 1.19")
+V = h^2/(2*m*q*(lambda*1e-10)^2) // Calculation of velocity of moving electron
+printf("\n Applied voltage on electron is %f V.\n\n\n",V)
+// Answer in book is 601.6 Volt
diff --git a/3014/CH1/EX1.19/Ex1_19.txt b/3014/CH1/EX1.19/Ex1_19.txt new file mode 100755 index 000000000..5a5004b3d --- /dev/null +++ b/3014/CH1/EX1.19/Ex1_19.txt @@ -0,0 +1,3 @@ +
+Example 1.19
+ Applied voltage on electron is 601.983516 V.
diff --git a/3014/CH1/EX1.2/Ex1_2.sce b/3014/CH1/EX1.2/Ex1_2.sce new file mode 100755 index 000000000..a5cb088f4 --- /dev/null +++ b/3014/CH1/EX1.2/Ex1_2.sce @@ -0,0 +1,10 @@ +clc
+//given that
+M = 1 // Mass of object in Kg
+v = 10 // velocity of object in m/s
+h = 6.625e-34 // Plank constant
+
+printf("Example 1.2")
+lambda = h/(M*v) // calculation of de Broglie wavelength
+
+printf("\n de Broglie wavelength of body is %e m.\n\n\n",lambda)
diff --git a/3014/CH1/EX1.2/Ex1_2.txt b/3014/CH1/EX1.2/Ex1_2.txt new file mode 100755 index 000000000..2f8187055 --- /dev/null +++ b/3014/CH1/EX1.2/Ex1_2.txt @@ -0,0 +1,3 @@ +
+Example 1.2
+ de Broglie wavelength of body is 6.625000e-35 m.
diff --git a/3014/CH1/EX1.21/Ex1_21.sce b/3014/CH1/EX1.21/Ex1_21.sce new file mode 100755 index 000000000..cc1e81b08 --- /dev/null +++ b/3014/CH1/EX1.21/Ex1_21.sce @@ -0,0 +1,10 @@ +clc
+//Given that
+k = 8.6e-5 // Boltzmann constant
+t = 37 // Temperature in degree Celsius
+h = 6.62e-34 // Plank constant
+m = 1.67e-27 // Mass of neutron
+printf("Example 1.21")
+lambda = h/sqrt(3*m*(k*1.6e-19)*(t+273))// Calculation of wavelength
+printf("\n Wavelength of neutron at %d degree Celsius is %f angstrom.\n\n\n",t,lambda*1e10)
+
diff --git a/3014/CH1/EX1.21/Ex1_21.txt b/3014/CH1/EX1.21/Ex1_21.txt new file mode 100755 index 000000000..3c63dd681 --- /dev/null +++ b/3014/CH1/EX1.21/Ex1_21.txt @@ -0,0 +1,3 @@ +
+Example 1.21
+ Wavelength of neutron at 37 degree Celsius is 1.432020 angstrom.
diff --git a/3014/CH1/EX1.22/Ex1_22.sce b/3014/CH1/EX1.22/Ex1_22.sce new file mode 100755 index 000000000..5c13cea55 --- /dev/null +++ b/3014/CH1/EX1.22/Ex1_22.sce @@ -0,0 +1,9 @@ +clc
+//Given that
+k = 8.6e-5 // Boltzmann constant
+t = 27 // Temperature in degree Celsius
+h = 6.62e-34 // Plank constant
+m = 6.7e-27 // Mass of helium atom
+printf("Example 1.22")
+lambda = h/sqrt(3*m*(k*1.6e-19)*(t+273))// Calculation of wavelength
+printf("\n Wavelength of helium at %d degree Celsius is %f angstrom.\n\n\n",t,lambda*1e10)
diff --git a/3014/CH1/EX1.22/Ex1_22.txt b/3014/CH1/EX1.22/Ex1_22.txt new file mode 100755 index 000000000..443148772 --- /dev/null +++ b/3014/CH1/EX1.22/Ex1_22.txt @@ -0,0 +1,3 @@ +
+Example 1.22
+ Wavelength of helium at 27 degree Celsius is 0.726758 angstrom.
diff --git a/3014/CH1/EX1.23/Ex1_23.sce b/3014/CH1/EX1.23/Ex1_23.sce new file mode 100755 index 000000000..781d59fb9 --- /dev/null +++ b/3014/CH1/EX1.23/Ex1_23.sce @@ -0,0 +1,12 @@ +clc
+//Given that
+E = 200 // energy of electrons in eV
+x = 20 // distance of screen in cm
+D = 2 // diameter of ring in cm
+h = 6.62e-34 // Plank constant
+m = 9.1e-31 // Mass of electron in kg
+printf("Example 1.23")
+lambda = h/sqrt(2*m*E*1.6e-19) // Calculation of wavelength
+theta = atan(D/(2*x))
+d = lambda/(2*sin(theta))// calculation of interatomic spacing of crystal
+printf("\n Interatomic spacing of crystal is %f angstrom.\n\n\n",d*1e10)
diff --git a/3014/CH1/EX1.23/Ex1_23.txt b/3014/CH1/EX1.23/Ex1_23.txt new file mode 100755 index 000000000..c026cc08d --- /dev/null +++ b/3014/CH1/EX1.23/Ex1_23.txt @@ -0,0 +1,4 @@ +
+Example 1.23
+ Interatomic spacing of crystal is 8.685393 angstrom.
+
diff --git a/3014/CH1/EX1.24/Ex1_24.sce b/3014/CH1/EX1.24/Ex1_24.sce new file mode 100755 index 000000000..51f83352d --- /dev/null +++ b/3014/CH1/EX1.24/Ex1_24.sce @@ -0,0 +1,10 @@ +
+clc
+//Given that
+r = 0.5 // Bohr radius of hydrogen in angstrom
+m = 9.1e-31 // Mass of neutron in Kg
+h = 6.6e-34 // Plank constant
+printf("Example 1.24")
+v = h/(2*%pi*r*1e-10*m) // velocity of electron in ground state
+printf("\n Velocity of electron in ground state is %e m/s.\n\n\n",v)
+// Answer in book is 2.31e6 m/s
diff --git a/3014/CH1/EX1.24/Ex1_24.txt b/3014/CH1/EX1.24/Ex1_24.txt new file mode 100755 index 000000000..cbcf78c5a --- /dev/null +++ b/3014/CH1/EX1.24/Ex1_24.txt @@ -0,0 +1,4 @@ +
+
+Example 1.24
+ Velocity of electron in ground state is 2.308621e+06 m/s.
diff --git a/3014/CH1/EX1.25/Ex1_25.sce b/3014/CH1/EX1.25/Ex1_25.sce new file mode 100755 index 000000000..80898b401 --- /dev/null +++ b/3014/CH1/EX1.25/Ex1_25.sce @@ -0,0 +1,10 @@ +
+clc
+//Given that
+lambda = 5890 // wavelength of yellow radiation in angstrom
+m = 9.1e-31 // Mass of neutron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.25")
+v = h/(lambda*1e-10*m) // velocity of electron in ground state
+printf("\n Velocity of electron in ground state is %e m/s.\n\n\n",v)
+// Answer in book is 1.24e3 m/s
diff --git a/3014/CH1/EX1.25/Ex1_25.txt b/3014/CH1/EX1.25/Ex1_25.txt new file mode 100755 index 000000000..6135995bf --- /dev/null +++ b/3014/CH1/EX1.25/Ex1_25.txt @@ -0,0 +1,3 @@ +
+Example 1.25
+ Velocity of electron in ground state is 1.236963e+03 m/s.
diff --git a/3014/CH1/EX1.26/Ex1_26.sce b/3014/CH1/EX1.26/Ex1_26.sce new file mode 100755 index 000000000..146f8e67a --- /dev/null +++ b/3014/CH1/EX1.26/Ex1_26.sce @@ -0,0 +1,12 @@ +
+clc
+//Given that
+lambda = 2 // wavelength of neutron in angstrom
+m = 1.67e-27 // Mass of neutron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.26")
+v = h/(lambda*1e-10*m) // velocity of neutron
+k = 0.5*m*v^2 // Kinetic energy of neutron
+printf("\n Velocity of neutron is %e m/s.",v)
+printf("\n Kinetic energy of neutron is %f eV.\n\n\n",k/1.6e-19)
+// Answer in book is 0.021eV
diff --git a/3014/CH1/EX1.26/Ex1_26.txt b/3014/CH1/EX1.26/Ex1_26.txt new file mode 100755 index 000000000..f419c6812 --- /dev/null +++ b/3014/CH1/EX1.26/Ex1_26.txt @@ -0,0 +1,5 @@ +
+
+Example 1.26
+ Velocity of neutron is 1.985030e+03 m/s.
+ Kinetic energy of neutron is 0.020564 eV.
diff --git a/3014/CH1/EX1.29/Ex1_29.sce b/3014/CH1/EX1.29/Ex1_29.sce new file mode 100755 index 000000000..994e1f041 --- /dev/null +++ b/3014/CH1/EX1.29/Ex1_29.sce @@ -0,0 +1,16 @@ +clc
+//given that
+v1 = 50 // Previous applied voltage
+v2 = 65 // final applied voltage
+k = 12.28
+d = 0.91 // Spacing in a crystal in angstrom
+printf("Example 1.29")
+
+lambda = k/sqrt(v1)
+theta= asin(lambda/(2*d))// Angel for initial applied voltage
+lambda1 = k/sqrt(v2)// wavelength for final applied voltage
+theta1 = asin(lambda1/(2*d))// Angel for final applied voltage
+printf("\n For first order, sin(theta) is %f \n For second order sin(theta) must be %f \n which is not possible for any value of angle. \n So no maxima occur for higher orders \n\n\n",sin(theta),2*sin(theta))
+printf("\n Angle of diffraction for first order of beam \n is %f degree at %d Volts\n\n\n",theta1*180/%pi,v2)
+// Answer in book is 57.14 degree
+
diff --git a/3014/CH1/EX1.29/Ex1_29.txt b/3014/CH1/EX1.29/Ex1_29.txt new file mode 100755 index 000000000..3fc5d4bce --- /dev/null +++ b/3014/CH1/EX1.29/Ex1_29.txt @@ -0,0 +1,12 @@ +
+Example 1.29
+ For first order, sin(theta) is 0.954206
+ For second order sin(theta) must be 1.908411
+ which is not possible for any value of angle.
+ So no maxima occur for higher orders
+
+
+
+ Angle of diffraction for first order of beam
+ is 56.813542 degree at 65 Volts
+
diff --git a/3014/CH1/EX1.3/Ex1_3.sce b/3014/CH1/EX1.3/Ex1_3.sce new file mode 100755 index 000000000..f0af4190c --- /dev/null +++ b/3014/CH1/EX1.3/Ex1_3.sce @@ -0,0 +1,10 @@ +clc
+// Given that
+m = 1e-30 // Mass of any object in Kg
+v = 1e5 // velocity of object in m/s
+h = 6.625e-34 // Plank constant
+
+printf("Example 1.3")
+lambda = h/(m*v) // calculation of de Broglie wavelength
+
+printf("\n de Broglie wavelength of body is %e m.\n\n\n",lambda)
diff --git a/3014/CH1/EX1.3/Ex1_3.txt b/3014/CH1/EX1.3/Ex1_3.txt new file mode 100755 index 000000000..3b4cf032a --- /dev/null +++ b/3014/CH1/EX1.3/Ex1_3.txt @@ -0,0 +1,3 @@ +
+Example 1.3
+ de Broglie wavelength of body is 6.625000e-09 m.
\ No newline at end of file diff --git a/3014/CH1/EX1.30/Ex1_30.sce b/3014/CH1/EX1.30/Ex1_30.sce new file mode 100755 index 000000000..d1d38596d --- /dev/null +++ b/3014/CH1/EX1.30/Ex1_30.sce @@ -0,0 +1,8 @@ +clc
+//Given that
+lambda = 680 // Wavelength in m
+g = 9.8 //Acceleration due to gravity
+printf("Example 1.30")
+v_g = 1/2*sqrt(g*lambda/(2*%pi)) // Calculation of group velocity
+printf("\n Group velocity of seawater waves is %f m/s.\n\n\n",v_g)
+// Answer in book is 16.29 m/s
diff --git a/3014/CH1/EX1.30/Ex1_30.txt b/3014/CH1/EX1.30/Ex1_30.txt new file mode 100755 index 000000000..6f7a09ead --- /dev/null +++ b/3014/CH1/EX1.30/Ex1_30.txt @@ -0,0 +1,4 @@ +
+
+Example 1.30
+ Group velocity of seawater waves is 16.283493 m/s.
diff --git a/3014/CH1/EX1.32/Ex1_32.sce b/3014/CH1/EX1.32/Ex1_32.sce new file mode 100755 index 000000000..98a97d20c --- /dev/null +++ b/3014/CH1/EX1.32/Ex1_32.sce @@ -0,0 +1,13 @@ +clc
+//Given that
+lambda = 2e-13 // de Broglie wavelength of an electron in m
+c = 3e8 // Speed of light in m/s
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.32")
+E = h*c/(lambda*1.6e-19)
+E_rest = m*c^2/(1.6e-19) // Calculation of rest mass energy
+E_total = sqrt(E^2+E_rest^2) // Total energy in eV
+v_g = c*sqrt(1-(E_rest/E_total)^2) // Group velocity
+v_p = c^2/v_g // Phase velocity
+printf("\n Group velocity of de Broglie waves is %fc and\n phase velocity is %fc .\n\n\n",v_g/c,v_p/c)
diff --git a/3014/CH1/EX1.32/Ex1_32.txt b/3014/CH1/EX1.32/Ex1_32.txt new file mode 100755 index 000000000..d7dd425cc --- /dev/null +++ b/3014/CH1/EX1.32/Ex1_32.txt @@ -0,0 +1,4 @@ +
+Example 1.32
+ Group velocity of de Broglie waves is 0.996626c and
+ phase velocity is 1.003385c .
diff --git a/3014/CH1/EX1.33/Ex1_33.sce b/3014/CH1/EX1.33/Ex1_33.sce new file mode 100755 index 000000000..6f3fb21b8 --- /dev/null +++ b/3014/CH1/EX1.33/Ex1_33.sce @@ -0,0 +1,18 @@ +
+clc
+//Given that
+lambda = 2e-12 // de Broglie wavelength of an electron in m
+c = 3e8 // Speed of light in m/s
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.33")
+E = h*c/(lambda*1.6e-19) // Energy due to momentum
+E_rest = m*c^2/(1.6e-19) // Calculation of rest mass energy
+E_total = sqrt(E^2+E_rest^2) // Total energy in eV
+KE = E_total - E_rest // Kinetic energy
+v_g = c*sqrt(1-(E_rest/E_total)^2) // Group velocity
+v_p = c^2/v_g // Phase velocity
+
+printf("\n Kinetic energy of electron is %f KeV.",KE/1000)
+printf("\n Group velocity of de Broglie waves is %fc m/s and\n phase velocity is %fc m/s.\n\n\n",v_g/c,v_p/c)
+// Answer in book is v_g = 0.6035c & v_p = 1.657c
diff --git a/3014/CH1/EX1.33/Ex1_33.txt b/3014/CH1/EX1.33/Ex1_33.txt new file mode 100755 index 000000000..f66d0dc70 --- /dev/null +++ b/3014/CH1/EX1.33/Ex1_33.txt @@ -0,0 +1,7 @@ +
+
+Example 1.33
+ Kinetic energy of electron is 293.330537 KeV.
+ Group velocity of de Broglie waves is 0.771930c m/s and
+ phase velocity is 1.295454c m/s.
+
diff --git a/3014/CH1/EX1.4/Ex1_4.sce b/3014/CH1/EX1.4/Ex1_4.sce new file mode 100755 index 000000000..0a9fddd53 --- /dev/null +++ b/3014/CH1/EX1.4/Ex1_4.sce @@ -0,0 +1,13 @@ +clc
+// Given that
+KE = 4.55e-25 // Kinetic energy of an electron in Joule
+m = 9.1e-31 // Mass of any object in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.4")
+v = sqrt(2*KE/m) // Calculation of velocity of moving electron
+p = m*v //Calculation of momentum of moving electron
+lambda = h/p // calculation of de Broglie wavelength
+printf("\n velocity of electron is %e m/s.",v)
+printf("\n momentum of electron is %e Kgm/s.",p)
+printf("\n de Broglie wavelength of electron is %e m.\n\n\n",lambda)
+
diff --git a/3014/CH1/EX1.4/Ex1_4.txt b/3014/CH1/EX1.4/Ex1_4.txt new file mode 100755 index 000000000..a943301f0 --- /dev/null +++ b/3014/CH1/EX1.4/Ex1_4.txt @@ -0,0 +1,6 @@ +
+
+Example 1.4
+ velocity of electron is 1.000000e+03 m/s.
+ momentum of electron is 9.100000e-28 Kgm/s.
+ de Broglie wavelength of electron is 7.274725e-07 m.
diff --git a/3014/CH1/EX1.5/Ex1_5.sce b/3014/CH1/EX1.5/Ex1_5.sce new file mode 100755 index 000000000..f96337272 --- /dev/null +++ b/3014/CH1/EX1.5/Ex1_5.sce @@ -0,0 +1,11 @@ +
+clc
+//Given that
+c = 3e8 // speed of light in m/s
+v = c/20 // Speed of proton in m/s
+m = 1.67e-27 // Mass of proton in Kg
+h = 6.625e-34 // Plank constant
+printf("Example 1.5")
+lambda = h/(m*v) // calculation of de Broglie wavelength
+printf("\n de Broglie wavelength of proton is %e m.\n\n\n",lambda)
+// Answer in book is 6.645e-14m which is a calculation mistake
diff --git a/3014/CH1/EX1.5/Ex1_5.txt b/3014/CH1/EX1.5/Ex1_5.txt new file mode 100755 index 000000000..e75a8a5d5 --- /dev/null +++ b/3014/CH1/EX1.5/Ex1_5.txt @@ -0,0 +1,4 @@ +
+
+Example 1.5
+ de Broglie wavelength of proton is 2.644711e-14 m.
diff --git a/3014/CH1/EX1.6/Ex1_6.sce b/3014/CH1/EX1.6/Ex1_6.sce new file mode 100755 index 000000000..cc61e20d5 --- /dev/null +++ b/3014/CH1/EX1.6/Ex1_6.sce @@ -0,0 +1,16 @@ +clc
+//Given that
+e = 12.8 // Energy of neutron in MeV
+c = 3e8 // speed of light in m/s
+m = 1.675e-27 // Mass of neutron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.6")
+rest_e = m*c^2/(1e6*1.6e-19)// rest mass energy of neutron in MeV
+if e/rest_e < 0.015 then
+ E = e;
+ else E = rest_e +e;
+end
+lambda = h/(sqrt(2*m*e*1e6*1.6e-19)) // calculation of de Broglie wavelength
+
+printf("\n de Broglie wavelength of neutron is %e angstrom.\n\n\n", lambda*1e10)
+// Answer in book is 8.04e-5 angstrom which is misprinted
diff --git a/3014/CH1/EX1.6/Ex1_6.txt b/3014/CH1/EX1.6/Ex1_6.txt new file mode 100755 index 000000000..c020af70a --- /dev/null +++ b/3014/CH1/EX1.6/Ex1_6.txt @@ -0,0 +1,4 @@ +
+Example 1.6
+ de Broglie wavelength of neutron is 7.992279e-05 angstrom.
+
diff --git a/3014/CH1/EX1.7/Ex1_7.sce b/3014/CH1/EX1.7/Ex1_7.sce new file mode 100755 index 000000000..0ce5bd9d2 --- /dev/null +++ b/3014/CH1/EX1.7/Ex1_7.sce @@ -0,0 +1,11 @@ +clc
+//Given that
+e = 1.632e-19 // charge on electron in coulomb
+V = 50 // Applied voltage in volts
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.62e-34 // Plank constant
+printf("Example 1.7")
+
+lambda = h/(sqrt(2*e*V*m)) // calculation of de Broglie wavelength
+printf("\n de Broglie wavelength of neutron is %f angstrom.\n\n\n", lambda*1e10)
+// Answer in book is 1.735 angstrom which is misprinted
diff --git a/3014/CH1/EX1.7/Ex1_7.txt b/3014/CH1/EX1.7/Ex1_7.txt new file mode 100755 index 000000000..5cc9c93d1 --- /dev/null +++ b/3014/CH1/EX1.7/Ex1_7.txt @@ -0,0 +1,3 @@ +
+Example 1.7
+ de Broglie wavelength of neutron is 1.717818 angstrom.
diff --git a/3014/CH1/EX1.9/Ex1_9.sce b/3014/CH1/EX1.9/Ex1_9.sce new file mode 100755 index 000000000..841e599e2 --- /dev/null +++ b/3014/CH1/EX1.9/Ex1_9.sce @@ -0,0 +1,10 @@ +clc
+//Given that
+e = 1.6e-19 // charge on electron in coulomb
+V = 54 // Applied voltage in volts
+m = 9.1e-31 // Mass of electron in Kg
+h = 6.63e-34 // Plank constant
+printf("Example 1.9")
+
+lambda = h/(sqrt(2*e*V*m)) // calculation of de Broglie wavelength
+printf("\n de Broglie wavelength of neutron is %f angstrom.\n\n\n", lambda*1e10)
diff --git a/3014/CH1/EX1.9/Ex1_9.txt b/3014/CH1/EX1.9/Ex1_9.txt new file mode 100755 index 000000000..c350ab0bf --- /dev/null +++ b/3014/CH1/EX1.9/Ex1_9.txt @@ -0,0 +1,3 @@ +
+Example 1.9
+ de Broglie wavelength of neutron is 1.671941 angstrom.
diff --git a/3014/CH2/EX2.1/Ex2_1.sce b/3014/CH2/EX2.1/Ex2_1.sce new file mode 100755 index 000000000..e6a294d73 --- /dev/null +++ b/3014/CH2/EX2.1/Ex2_1.sce @@ -0,0 +1,9 @@ + clc
+//given that
+del_x = 0.2 // Uncertainty in position in angstrom
+h = 6.63e-34 // Plank constant
+
+printf("Example 2.1")
+h_bar = h / (2*%pi) // constant
+del_p = h_bar/(2*del_x*1e-10) // Calculation of uncertainty in momentum
+printf("\n Uncertainty in momentum of particle is %e kgm/sec \n\n\n",del_p)
diff --git a/3014/CH2/EX2.1/Ex2_1.txt b/3014/CH2/EX2.1/Ex2_1.txt new file mode 100755 index 000000000..ba18b7e14 --- /dev/null +++ b/3014/CH2/EX2.1/Ex2_1.txt @@ -0,0 +1,2 @@ + Example 2.1
+ Uncertainty in momentum of particle is 2.637993e-24 kgm/sec.
\ No newline at end of file diff --git a/3014/CH2/EX2.10/Ex2_10.sce b/3014/CH2/EX2.10/Ex2_10.sce new file mode 100755 index 000000000..3cf24dcc9 --- /dev/null +++ b/3014/CH2/EX2.10/Ex2_10.sce @@ -0,0 +1,15 @@ +clc
+
+//given that
+r = 0.53 // radius of hydrogen atom in angstrom
+m_e = 9.1e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.10")
+h_bar = h / (2*%pi) // constant
+del_x = 2*r // calculation of uncertainty in position
+del_p = h_bar/(2*del_x*1e-10) // calculation of uncertainty in momentum
+p = del_p
+E = p^2/(2*m_e*1.6e-19)// Calculation of energy in eV
+printf("\n Kinetic energy needed by an electron to be \n confined in electron is %f eV.\n\n\n",E)
+// When problem is solved by del_x*del_p = h_bar, then minimum value of kinetic energy will become 13.6eV
+
diff --git a/3014/CH2/EX2.10/Ex2_10.txt b/3014/CH2/EX2.10/Ex2_10.txt new file mode 100755 index 000000000..8a1fb3267 --- /dev/null +++ b/3014/CH2/EX2.10/Ex2_10.txt @@ -0,0 +1,4 @@ +
+Example 2.10
+ Kinetic energy needed by an electron to be
+ confined in electron is 0.850754 eV.
diff --git a/3014/CH2/EX2.11/Ex2_11.sce b/3014/CH2/EX2.11/Ex2_11.sce new file mode 100755 index 000000000..633fdfe1b --- /dev/null +++ b/3014/CH2/EX2.11/Ex2_11.sce @@ -0,0 +1,10 @@ +clc
+
+//given that
+del_t = 2.5e-14 // lifetime in exited state in micro sec
+h = 6.63e-34 // Plank constant
+printf("Example 2.11")
+h_bar = h / (2*%pi) // constant
+del_E = h_bar/(1.6e-19*del_t*1e-6) // calculation of uncertainty in momentum
+printf("\n Minimum error in measurement of energy of this state is %e eV.\n\n\n",del_E)
+
diff --git a/3014/CH2/EX2.11/Ex2_11.txt b/3014/CH2/EX2.11/Ex2_11.txt new file mode 100755 index 000000000..6a0d2541b --- /dev/null +++ b/3014/CH2/EX2.11/Ex2_11.txt @@ -0,0 +1,3 @@ +
+Example 2.11
+ Minimum error in measurement of energy of this state is 2.637993e+04 eV.
diff --git a/3014/CH2/EX2.12/Ex2_12.sce b/3014/CH2/EX2.12/Ex2_12.sce new file mode 100755 index 000000000..7b9aa959e --- /dev/null +++ b/3014/CH2/EX2.12/Ex2_12.sce @@ -0,0 +1,15 @@ +clc
+
+//given that
+E_eV = 0.5// kinetic energy of electron in KeV
+del_x = 0.4 // Uncertainty in position in nm
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in kg
+printf("Example 2.12")
+h_bar = h / (2*%pi) // constant
+E_J = E_eV*1e3*1.6e-19
+p = sqrt(2*m*E_J) // Calculation of momentum in kgm/s
+del_p = h_bar/(2*del_x*1e-9) // Calculation of uncertainty in momentum
+per_error = del_p*100 / p // calculation of percentage error in momentum
+printf("\n Percentage error in momentum is %f percent.\n\n\n",per_error)
+// Answer in book is 1.08 percentage
diff --git a/3014/CH2/EX2.12/Ex2_12.txt b/3014/CH2/EX2.12/Ex2_12.txt new file mode 100755 index 000000000..5168c739a --- /dev/null +++ b/3014/CH2/EX2.12/Ex2_12.txt @@ -0,0 +1,3 @@ +
+Example 2.12
+ Percentage error in momentum is 1.093108 percent.
\ No newline at end of file diff --git a/3014/CH2/EX2.13/Ex2_13.sce b/3014/CH2/EX2.13/Ex2_13.sce new file mode 100755 index 000000000..e7e5df6d9 --- /dev/null +++ b/3014/CH2/EX2.13/Ex2_13.sce @@ -0,0 +1,11 @@ +clc
+
+//given that
+del_x = 2e-9 // Uncertainty in position in m
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in Kg
+printf("Example 2.13")
+h_bar = h / (2*%pi) // constant
+del_p = h_bar/(2*del_x) // Calculation of uncertainty in momentum
+del_v = del_p/m
+printf("\n Uncertainty in velocity of particle is %e m/s.\n\n\n",del_v)
diff --git a/3014/CH2/EX2.13/Ex2_13.txt b/3014/CH2/EX2.13/Ex2_13.txt new file mode 100755 index 000000000..2f12c1e89 --- /dev/null +++ b/3014/CH2/EX2.13/Ex2_13.txt @@ -0,0 +1,3 @@ +
+Example 2.13
+ Uncertainty in velocity of particle is 2.898894e+04 m/s.
\ No newline at end of file diff --git a/3014/CH2/EX2.15/Ex2_15.sce b/3014/CH2/EX2.15/Ex2_15.sce new file mode 100755 index 000000000..ce8c8d299 --- /dev/null +++ b/3014/CH2/EX2.15/Ex2_15.sce @@ -0,0 +1,15 @@ +clc
+
+//given that
+del_x = 5000 // Uncertainty in position in angstrom
+h = 6.63e-34 // Plank constant
+m = 200 // mass of ball in gram
+v = 6 // velocity of moving ball in m/s
+printf("Example 2.15")
+h_bar = h / (2*%pi) // constant
+del_p = h_bar/(2*del_x*1e-10) // Calculation of uncertainty in momentum
+p = m*v/1000 // Calculation of momentum
+per_error = del_p*100/p // Calculation of percentage error in calculation of momentum
+printf("\n Uncertainty in momentum of ball is %e kgm/s.",del_p)
+printf("\n Percentage error in calculation of momentum is %e.\n\n\n",per_error)
+
diff --git a/3014/CH2/EX2.15/Ex2_15.txt b/3014/CH2/EX2.15/Ex2_15.txt new file mode 100755 index 000000000..a0a2d5693 --- /dev/null +++ b/3014/CH2/EX2.15/Ex2_15.txt @@ -0,0 +1,4 @@ +
+Example 2.15
+ Uncertainty in momentum of ball is 1.055197e-28 kgm/s.
+ Percentage error in calculation of momentum is 8.793311e-27.
diff --git a/3014/CH2/EX2.16/Ex2_16.sce b/3014/CH2/EX2.16/Ex2_16.sce new file mode 100755 index 000000000..2c7e4a73e --- /dev/null +++ b/3014/CH2/EX2.16/Ex2_16.sce @@ -0,0 +1,15 @@ +clc
+
+//given that
+c = 3e8 // speed of light in m/s
+v = c/10 // Velocity of moving proton in m/s
+v_error = 1 // Percentage error in measurement of velocity
+m = 1.67e-27 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+
+printf("Example 2.16")
+h_bar = h / (2*%pi) // constant
+del_v = v*v_error/100// calculation of uncertainty in position
+del_x = h_bar/(2*m*del_v) // calculation of uncertainty in momentum
+printf("\n Uncertainty in position of particle is %e m.\n\n\n",del_x)
+// Answer in book is 1.04e-13 m
diff --git a/3014/CH2/EX2.16/Ex2_16.txt b/3014/CH2/EX2.16/Ex2_16.txt new file mode 100755 index 000000000..e0b7aecbb --- /dev/null +++ b/3014/CH2/EX2.16/Ex2_16.txt @@ -0,0 +1,4 @@ +
+
+Example 2.16
+ Uncertainty in position of particle is 1.053091e-13 m.
diff --git a/3014/CH2/EX2.17/Ex2_17.sce b/3014/CH2/EX2.17/Ex2_17.sce new file mode 100755 index 000000000..3fa55fa59 --- /dev/null +++ b/3014/CH2/EX2.17/Ex2_17.sce @@ -0,0 +1,11 @@ +clc
+
+//given that
+del_x = 1e-9 // Uncertainty in position in m
+h = 6.63e-34 // Plank constant
+m = 200 // mass of ball in gram
+printf("Example 2.17")
+h_bar = h / (2*%pi) // constant
+del_v = h_bar/(2*del_x*m/1000) // Calculation of uncertainty in momentum
+printf("\n Uncertainty in velocity of ball is %e m/s.\n\n\n",del_v)
+// Answer in book is 2.64e-25 m/s
diff --git a/3014/CH2/EX2.17/Ex2_17.txt b/3014/CH2/EX2.17/Ex2_17.txt new file mode 100755 index 000000000..81919635e --- /dev/null +++ b/3014/CH2/EX2.17/Ex2_17.txt @@ -0,0 +1,3 @@ +
+Example 2.17
+ Uncertainty in velocity of ball is 2.637993e-25 m/s.
diff --git a/3014/CH2/EX2.18/Ex2_18.sce b/3014/CH2/EX2.18/Ex2_18.sce new file mode 100755 index 000000000..a8a499c0d --- /dev/null +++ b/3014/CH2/EX2.18/Ex2_18.sce @@ -0,0 +1,12 @@ +clc
+//given that
+del_t = 2e-12 // lifetime of exited state in sec
+h = 6.63e-34 // Plank constant
+printf("Example 2.18")
+h_bar = h / (2*%pi) // constant
+del_E = h_bar/(1.6e-19*2*del_t) // calculation of uncertainty in momentum
+printf("\n Minimum error in measurement of energy of this state is %e eV.\n\n\n",del_E)
+// Answer in book is 1.65e-4 eV
+
+
+
diff --git a/3014/CH2/EX2.18/Ex2_18.txt b/3014/CH2/EX2.18/Ex2_18.txt new file mode 100755 index 000000000..5f8e25a3f --- /dev/null +++ b/3014/CH2/EX2.18/Ex2_18.txt @@ -0,0 +1,3 @@ +
+Example 2.18
+ Minimum error in measurement of energy of this state is 1.648746e-04 eV.
diff --git a/3014/CH2/EX2.19/Ex2_19.sce b/3014/CH2/EX2.19/Ex2_19.sce new file mode 100755 index 000000000..f7fa6c4cc --- /dev/null +++ b/3014/CH2/EX2.19/Ex2_19.sce @@ -0,0 +1,9 @@ +clc
+//given that
+del_t = 1e-8 // lifetime of exited state in sec
+h = 6.63e-34 // Plank constant
+printf("Example 2.19")
+h_bar = h / (2*%pi) // constant
+del_nu = h_bar/(2*del_t*h) // calculation of uncertainty in frequency
+printf("\n Minimum error in measurement of \n frequency of photon is %e per second.\n\n\n",del_nu)
+// Answer in book is 8e6 per second
diff --git a/3014/CH2/EX2.19/Ex2_19.txt b/3014/CH2/EX2.19/Ex2_19.txt new file mode 100755 index 000000000..c9a8f72ef --- /dev/null +++ b/3014/CH2/EX2.19/Ex2_19.txt @@ -0,0 +1,4 @@ +
+Example 2.19
+ Minimum error in measurement of
+ frequency of photon is 7.957747e+06 per second.
diff --git a/3014/CH2/EX2.2/Ex2_2.sce b/3014/CH2/EX2.2/Ex2_2.sce new file mode 100755 index 000000000..fc9241b1b --- /dev/null +++ b/3014/CH2/EX2.2/Ex2_2.sce @@ -0,0 +1,11 @@ +
+clc
+//given that
+del_x = 4e-10 // Uncertainty in position in m
+h = 6.63e-34 // Plank constant
+
+printf("Example 2.2")
+h_bar = h / (2*%pi) // constant
+del_p = h_bar/(2*del_x) // Calculation of uncertainty in momentum
+printf("\n Uncertainty in momentum of particle is %e kgm/sec.\n\n\n",del_p)
+// Answer in book is given as 1.32e-23 kgm/sec
diff --git a/3014/CH2/EX2.2/Ex2_2.txt b/3014/CH2/EX2.2/Ex2_2.txt new file mode 100755 index 000000000..dba81e11b --- /dev/null +++ b/3014/CH2/EX2.2/Ex2_2.txt @@ -0,0 +1,3 @@ +
+Example 2.2
+ Uncertainty in momentum of particle is 1.318997e-25 kgm/sec.
\ No newline at end of file diff --git a/3014/CH2/EX2.20/Ex2_20.sce b/3014/CH2/EX2.20/Ex2_20.sce new file mode 100755 index 000000000..ef85f594b --- /dev/null +++ b/3014/CH2/EX2.20/Ex2_20.sce @@ -0,0 +1,11 @@ +clc
+//given that
+del_v = 5.5e-20 // Uncertainty in velocity in m/s
+h = 6.63e-34 // Plank constant
+m = 1 // mass of dust particle in mg
+printf("Example 2.20")
+h_bar = h / (2*%pi) // constant
+del_x = h_bar/(2*del_v*m*1e-6) // Calculation of uncertainty in momentum
+printf("\n Uncertainty in position of ball is %f angstrom.\n\n\n",del_x*1e10)
+// Answer in book is 9.6 angstrom
+
diff --git a/3014/CH2/EX2.20/Ex2_20.txt b/3014/CH2/EX2.20/Ex2_20.txt new file mode 100755 index 000000000..ef76d8497 --- /dev/null +++ b/3014/CH2/EX2.20/Ex2_20.txt @@ -0,0 +1,4 @@ +
+
+Example 2.20
+ Uncertainty in position of ball is 9.592702 angstrom.
diff --git a/3014/CH2/EX2.21/Ex2_21.sce b/3014/CH2/EX2.21/Ex2_21.sce new file mode 100755 index 000000000..6a4ef10f9 --- /dev/null +++ b/3014/CH2/EX2.21/Ex2_21.sce @@ -0,0 +1,14 @@ +clc
+
+
+//given that
+l = 1 // width of potential well in angstrom
+n = 1 // order corresponding to ground state
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in Kg
+printf("Example 2.21")
+E = n^2*h^2/(8*m*(l*1e-10)^2) // Calculation of energy in Joule
+E_eV = E/1.6e-19 // Calculation of energy in eV
+
+printf("\n Energy of electron is %f eV.\n\n\n",E_eV)
+// Answer in book is 37.74 eV angstrom
diff --git a/3014/CH2/EX2.21/Ex2_21.txt b/3014/CH2/EX2.21/Ex2_21.txt new file mode 100755 index 000000000..932aff438 --- /dev/null +++ b/3014/CH2/EX2.21/Ex2_21.txt @@ -0,0 +1,6 @@ +
+
+
+Example 2.21
+ Energy of electron is 37.737723 eV.
+
diff --git a/3014/CH2/EX2.24/Ex2_24.sce b/3014/CH2/EX2.24/Ex2_24.sce new file mode 100755 index 000000000..3cc4f425a --- /dev/null +++ b/3014/CH2/EX2.24/Ex2_24.sce @@ -0,0 +1,13 @@ +clc
+
+//given that
+l = 2.5e-10 // width of potential well in m
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in Kg
+printf("Example 2.24")
+for n = 1:2
+ E = n^2*h^2/(8*m*l^2) // Calculation of energy in Joule
+E_eV = E/1.6e-19 // Calculation of energy in eV
+
+printf("\n Energy of electron for state %d is %f eV.",n,E_eV);
+end
diff --git a/3014/CH2/EX2.24/Ex2_24.txt b/3014/CH2/EX2.24/Ex2_24.txt new file mode 100755 index 000000000..f9b6ee4c7 --- /dev/null +++ b/3014/CH2/EX2.24/Ex2_24.txt @@ -0,0 +1,4 @@ +
+Example 2.24
+ Energy of electron for state 1 is 6.038036 eV.
+ Energy of electron for state 2 is 24.152143 eV.
\ No newline at end of file diff --git a/3014/CH2/EX2.26/Ex2_26.sce b/3014/CH2/EX2.26/Ex2_26.sce new file mode 100755 index 000000000..9fff3fadf --- /dev/null +++ b/3014/CH2/EX2.26/Ex2_26.sce @@ -0,0 +1,11 @@ +clc
+// given that
+L = 1// let unit length
+l1 = 0.45*L // initial point
+l2 = 0.55*L // Final point
+
+
+printf("Example 2.26 \n")
+p = (1/L)*((l2-(L/(2*%pi) *sin(2*l2*%pi/L)))- (l1-(L/(2*%pi) *sin(2*l1*%pi/L)))) // Calculation of probability of finding particle
+p_per = p*100 // probability of finding particle in percentage
+printf("\n Probability of finding electron between \n %fL and %fL is %f percent.",l2,l1,p_per)
diff --git a/3014/CH2/EX2.26/Ex2_26.txt b/3014/CH2/EX2.26/Ex2_26.txt new file mode 100755 index 000000000..c37c26fec --- /dev/null +++ b/3014/CH2/EX2.26/Ex2_26.txt @@ -0,0 +1,4 @@ +Example 2.26
+
+ Probability of finding electron between
+ 0.550000L and 0.450000L is 19.836316 percent.
\ No newline at end of file diff --git a/3014/CH2/EX2.27/Ex2_27.sce b/3014/CH2/EX2.27/Ex2_27.sce new file mode 100755 index 000000000..0f562e720 --- /dev/null +++ b/3014/CH2/EX2.27/Ex2_27.sce @@ -0,0 +1,16 @@ +clc
+
+
+//given that
+l = 1e-8 // width of potential well in cm
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in Kg
+printf("\nExample 2.27")
+E_1 = (h)^2/(8*m*(l*1e-2)^2) // Calculation of energy of ground state in Joule
+E_1_eV = E_1/1.6e-19 // Calculation of energy in eV
+E_2 = (2)^2*h^2/(8*m*(l*1e-2)^2) // Calculation of energy of first state in Joule
+E_2_eV = E_2/1.6e-19 // Calculation of energy in eV
+del_E = E_2_eV - E_1_eV // calculation of difference between first state and ground state
+printf("\n Difference between first state \n and ground state energies is %f eV.\n\n\n",del_E);
+// Answer in book is 113.04 eV
+
diff --git a/3014/CH2/EX2.27/Ex2_27.txt b/3014/CH2/EX2.27/Ex2_27.txt new file mode 100755 index 000000000..883aa5c59 --- /dev/null +++ b/3014/CH2/EX2.27/Ex2_27.txt @@ -0,0 +1,4 @@ +
+Example 2.27
+ Difference between first state
+ and ground state energies is 113.213170 eV.
diff --git a/3014/CH2/EX2.28/Ex2_28.sce b/3014/CH2/EX2.28/Ex2_28.sce new file mode 100755 index 000000000..bb8cb9c91 --- /dev/null +++ b/3014/CH2/EX2.28/Ex2_28.sce @@ -0,0 +1,13 @@ +clc
+
+//given that
+l = 1 // width of potential well in angstrom
+h = 6.63e-34 // Plank constant
+m = 9.1e-31 // mass of electron in Kg
+printf("Example 2.28")
+for n = 1:3
+ lambda = 2*l/n // Calculation of wavelength
+ E = n^2*h^2/(8*m*(l*1e-10)^2) // Calculation of energy in Joule
+E_eV = E/1.6e-19 // Calculation of energy in eV
+printf("\n For state:%d Energy is %f eV & wavelength is %f angstrom ",n,E_eV,lambda);
+end
diff --git a/3014/CH2/EX2.28/Ex2_28.txt b/3014/CH2/EX2.28/Ex2_28.txt new file mode 100755 index 000000000..de22674c0 --- /dev/null +++ b/3014/CH2/EX2.28/Ex2_28.txt @@ -0,0 +1,6 @@ +
+
+Example 2.28
+ For state:1 Energy is 37.737723 eV & wavelength is 2.000000 angstrom
+ For state:2 Energy is 150.950893 eV & wavelength is 1.000000 angstrom
+ For state:3 Energy is 339.639509 eV & wavelength is 0.666667 angstrom
\ No newline at end of file diff --git a/3014/CH2/EX2.29/Ex2_29.sce b/3014/CH2/EX2.29/Ex2_29.sce new file mode 100755 index 000000000..15b6f45a8 --- /dev/null +++ b/3014/CH2/EX2.29/Ex2_29.sce @@ -0,0 +1,12 @@ +clc
+
+//given that
+m = 100 //mass of ball in gram
+l = 1 // length of box in m
+h = 6.63e-34 // Plank constant
+printf("\nExample 2.29")
+for n = 1:3
+ E = (n^2*h^2)/(8*m*1e-3*l^2*1.6e-19)
+ printf("\n Energy state E%d of ball is %e eV",n,E)
+end
+printf("\n As energy difference is very small so we cannot see energy states.\n")
diff --git a/3014/CH2/EX2.29/Ex2_29.txt b/3014/CH2/EX2.29/Ex2_29.txt new file mode 100755 index 000000000..4ddda5ec8 --- /dev/null +++ b/3014/CH2/EX2.29/Ex2_29.txt @@ -0,0 +1,6 @@ +
+Example 2.29
+ Energy state E1 of ball is 3.434133e-48 eV
+ Energy state E2 of ball is 1.373653e-47 eV
+ Energy state E3 of ball is 3.090720e-47 eV
+ As energy difference is very small so we cannot see energy states.
diff --git a/3014/CH2/EX2.3/Ex2_3.sce b/3014/CH2/EX2.3/Ex2_3.sce new file mode 100755 index 000000000..2afc4914a --- /dev/null +++ b/3014/CH2/EX2.3/Ex2_3.sce @@ -0,0 +1,14 @@ +clc
+
+//given that
+v = 3e7 // Velocity of moving electron in m/s
+m = 9.1e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+c = 3e8 // speed of light in m/s
+printf("Example 2.3")
+h_bar = h / (2*%pi) // constant
+del_p = m*v/(sqrt(1-(v/c)^2)) // calculation of uncertainty in momentum
+del_x = h_bar/(2*del_p) // Calculation of uncertainty in position
+printf("\n Uncertainty in position of particle is %f angstrom.\n\n\n",del_x*1e10)
+//Answer in book is 0.0194 angstrom which is due to using approximate values at intermediate steps
+
diff --git a/3014/CH2/EX2.3/Ex2_3.txt b/3014/CH2/EX2.3/Ex2_3.txt new file mode 100755 index 000000000..600c1ab98 --- /dev/null +++ b/3014/CH2/EX2.3/Ex2_3.txt @@ -0,0 +1,3 @@ +
+Example 2.3
+ Uncertainty in position of particle is 0.019229 angstrom.
diff --git a/3014/CH2/EX2.30/Ex2_30.sce b/3014/CH2/EX2.30/Ex2_30.sce new file mode 100755 index 000000000..154543a6d --- /dev/null +++ b/3014/CH2/EX2.30/Ex2_30.sce @@ -0,0 +1,13 @@ +clc
+
+//given that
+l = 30 // width of potential well in angstrom
+x = l/2
+del_x = 2 // interval of length at centre in angstrom
+h = 6.63e-34 // Plank constant
+n = 1 // ground state
+printf("\nExample 2.30")
+phi_x = ((sqrt(2/l))*sin(n*%pi*x/l))^2
+p = phi_x*del_x // Calculation of probability at centre
+printf("\n Probability of finding particle at centre is %d percent.\n\n\n",p*100)
+// Answer given in book is 16 percent. It is due to wrong calculation
diff --git a/3014/CH2/EX2.30/Ex2_30.txt b/3014/CH2/EX2.30/Ex2_30.txt new file mode 100755 index 000000000..25a24c29a --- /dev/null +++ b/3014/CH2/EX2.30/Ex2_30.txt @@ -0,0 +1,3 @@ +
+Example 2.30
+ Probability of finding particle at centre is 13 percent.
diff --git a/3014/CH2/EX2.5/Ex2_5.sce b/3014/CH2/EX2.5/Ex2_5.sce new file mode 100755 index 000000000..aeae45abe --- /dev/null +++ b/3014/CH2/EX2.5/Ex2_5.sce @@ -0,0 +1,14 @@ +clc
+//given that
+v = 1.05e4 // Velocity of moving electron in m/s
+v_error = 0.02 //Percentage error in measurement of velocity
+
+m = 9e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.5")
+h_bar = h / (2*%pi) // constant
+p = m*v
+del_p = v_error*p/100 // calculation of uncertainty in momentum
+del_x = h_bar/del_p
+printf("\n Uncertainty in position of particle is %e m.\n\n\n",del_x)
+// Answer in book is given as 5.58e-3 m
diff --git a/3014/CH2/EX2.5/Ex2_5.txt b/3014/CH2/EX2.5/Ex2_5.txt new file mode 100755 index 000000000..57c3dff07 --- /dev/null +++ b/3014/CH2/EX2.5/Ex2_5.txt @@ -0,0 +1,3 @@ +
+Example 2.5
+ Uncertainty in position of particle is 1.922908e-12 m.
\ No newline at end of file diff --git a/3014/CH2/EX2.6/Ex2_6.sce b/3014/CH2/EX2.6/Ex2_6.sce new file mode 100755 index 000000000..aeceb712a --- /dev/null +++ b/3014/CH2/EX2.6/Ex2_6.sce @@ -0,0 +1,13 @@ +clc
+//given that
+v = 600 // Velocity of moving electron in m/s
+v_error = 0.005 //Percentage error in measurement of velocity
+m = 9.1e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.6")
+h_bar = h / (2*%pi) // constant
+p = m*v
+del_p = v_error*p/100 // calculation of uncertainty in momentum
+del_x = h_bar/(del_p) // Calculation of uncertainty in position
+printf("\n Uncertainty in position of particle is %e m.\n\n\n",del_x)
+// Answer in book is 0.39e-2 m
diff --git a/3014/CH2/EX2.6/Ex2_6.txt b/3014/CH2/EX2.6/Ex2_6.txt new file mode 100755 index 000000000..51aaf1205 --- /dev/null +++ b/3014/CH2/EX2.6/Ex2_6.txt @@ -0,0 +1,3 @@ +
+Example 2.6
+ Uncertainty in position of particle is 3.865191e-03 m.
\ No newline at end of file diff --git a/3014/CH2/EX2.7/Ex2_7.sce b/3014/CH2/EX2.7/Ex2_7.sce new file mode 100755 index 000000000..269df61e0 --- /dev/null +++ b/3014/CH2/EX2.7/Ex2_7.sce @@ -0,0 +1,10 @@ +clc
+//given that
+del_x = 1 // let uncertainty in position is unity
+m_e = 9.1e-31 // mass of electron in kg
+m_p = 1.67e-27 // mass of proton in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.7")
+h_bar = h / (2*%pi) // constant
+del_v_ratio = m_p/m_e // calculation in uncertainties in the velocity of electron and proton
+printf("\n Ratio of uncertainties in the velocity of electron to proton is %d.\n\n\n",del_v_ratio)
diff --git a/3014/CH2/EX2.7/Ex2_7.txt b/3014/CH2/EX2.7/Ex2_7.txt new file mode 100755 index 000000000..75c3fb7a4 --- /dev/null +++ b/3014/CH2/EX2.7/Ex2_7.txt @@ -0,0 +1,3 @@ +
+Example 2.7
+ Ratio of uncertainties in the velocity of electron to proton is 1835.
diff --git a/3014/CH2/EX2.8/Ex2_8.sce b/3014/CH2/EX2.8/Ex2_8.sce new file mode 100755 index 000000000..d42f02a64 --- /dev/null +++ b/3014/CH2/EX2.8/Ex2_8.sce @@ -0,0 +1,13 @@ +clc
+//given that
+r = 0.5 // radius of hydrogen atom in angstrom
+m_e = 9.1e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.8")
+h_bar = h / (2*%pi) // constant
+del_x = 2*r // calculation of uncertainty in position
+del_p = h_bar/(2*del_x*1e-10) // calculation of uncertainty in momentum
+p = del_p
+E = p^2/(2*m_e*1.6e-19)// Calculation of energy in eV
+printf("\n Kinetic energy needed by an electron to be \n confined in electron is %f eV.\n\n\n",ceil(E*100)/100)
+
diff --git a/3014/CH2/EX2.8/Ex2_8.txt b/3014/CH2/EX2.8/Ex2_8.txt new file mode 100755 index 000000000..c66325d48 --- /dev/null +++ b/3014/CH2/EX2.8/Ex2_8.txt @@ -0,0 +1,4 @@ +
+Example 2.8
+ Kinetic energy needed by an electron to be
+ confined in electron is 0.955908 eV.
diff --git a/3014/CH2/EX2.9/Ex2_9.sce b/3014/CH2/EX2.9/Ex2_9.sce new file mode 100755 index 000000000..7bed2b03a --- /dev/null +++ b/3014/CH2/EX2.9/Ex2_9.sce @@ -0,0 +1,13 @@ +clc
+//given that
+v = 5e3 // Velocity of moving electron in m/s
+v_error = 0.003 //Percentage error in measurement of velocity
+
+m = 9.1e-31 // mass of electron in kg
+h = 6.63e-34 // Plank constant
+printf("Example 2.9")
+h_bar = h / (2*%pi) // constant
+p = m*v
+del_p = v_error*p/100 // calculation of uncertainty in momentum
+del_x = h_bar/(2*del_p) // Calculation of uncertainty in position
+printf("\n Uncertainty in position of particle is %e m.\n\n\n",del_x)
diff --git a/3014/CH2/EX2.9/Ex2_9.txt b/3014/CH2/EX2.9/Ex2_9.txt new file mode 100755 index 000000000..5b94a08ae --- /dev/null +++ b/3014/CH2/EX2.9/Ex2_9.txt @@ -0,0 +1,3 @@ +
+Example 2.9
+ Uncertainty in position of particle is 3.865191e-04 m.
diff --git a/3014/CH3/EX3.1/Ex3_1.sce b/3014/CH3/EX3.1/Ex3_1.sce new file mode 100755 index 000000000..198388a30 --- /dev/null +++ b/3014/CH3/EX3.1/Ex3_1.sce @@ -0,0 +1,9 @@ +clc
+//given that
+d = 2.82 // crystal spacing in angstrom
+n = 2 // order for longest passing wavelength
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.1")
+lambda = 2*d*sin(theta*%pi/180)/n // Calculation of longest wavelength
+
+printf("\n Longest wavelength is %f angstrom. \n\n\n",lambda)
diff --git a/3014/CH3/EX3.1/Ex3_1.txt b/3014/CH3/EX3.1/Ex3_1.txt new file mode 100755 index 000000000..0d2e37b77 --- /dev/null +++ b/3014/CH3/EX3.1/Ex3_1.txt @@ -0,0 +1,2 @@ + Example 3.1
+ Longest wavelength is 2.820000 angstrom.
\ No newline at end of file diff --git a/3014/CH3/EX3.10/Ex3_10.sce b/3014/CH3/EX3.10/Ex3_10.sce new file mode 100755 index 000000000..fee4c2953 --- /dev/null +++ b/3014/CH3/EX3.10/Ex3_10.sce @@ -0,0 +1,11 @@ + clc
+//given that
+E = 100 // Energy of X ray beam in KeV
+theta = 30 // Scattering angle in degree
+m = 9.1e-31 // mass of electron in kg
+c = 3e8 // Speed of light in m/s
+printf("Example 3.10")
+E_rest = m*c^2/(1.6e-19*1e3) // Rest mass energy in KeV
+k = 1/E + (1-cos(theta*%pi/180))/(E_rest)
+del_e = E - 1/k // Energy of recoiled electron
+ printf("\n Energy of recoiled electron is %f KeV\n\n\n",del_e)
diff --git a/3014/CH3/EX3.10/Ex3_10.txt b/3014/CH3/EX3.10/Ex3_10.txt new file mode 100755 index 000000000..332227ac7 --- /dev/null +++ b/3014/CH3/EX3.10/Ex3_10.txt @@ -0,0 +1,3 @@ +
+Example 3.10
+ Energy of recoiled electron is 2.550573 KeV
\ No newline at end of file diff --git a/3014/CH3/EX3.11/Ex3_11.sce b/3014/CH3/EX3.11/Ex3_11.sce new file mode 100755 index 000000000..97855a53c --- /dev/null +++ b/3014/CH3/EX3.11/Ex3_11.sce @@ -0,0 +1,12 @@ + clc
+//given that
+lambda = 1 //wavelength in angstrom
+h = 6.62e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.11")
+d_lambda= h*(1-cos(theta*%pi/180))/(m_e*c) // calculation of wavelength shift
+
+printf("\nWavelength shift is %f angstrom. \n\n\n",d_lambda*1e10)
+
diff --git a/3014/CH3/EX3.11/Ex3_11.txt b/3014/CH3/EX3.11/Ex3_11.txt new file mode 100755 index 000000000..dbc11e9fc --- /dev/null +++ b/3014/CH3/EX3.11/Ex3_11.txt @@ -0,0 +1,3 @@ +
+Example 3.11
+Wavelength shift is 0.024249 angstrom.
diff --git a/3014/CH3/EX3.12/Ex3_12.sce b/3014/CH3/EX3.12/Ex3_12.sce new file mode 100755 index 000000000..6aa91d8c9 --- /dev/null +++ b/3014/CH3/EX3.12/Ex3_12.sce @@ -0,0 +1,12 @@ + clc
+//given that
+lambda = 0.015 //wavelength in angstrom
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 60 // angle for longest passing wavelength
+printf("Example 3.12")
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda
+
+printf("\n Wavelength shift is %f angstrom. \n\n\n",lambda_n)
diff --git a/3014/CH3/EX3.12/Ex3_12.txt b/3014/CH3/EX3.12/Ex3_12.txt new file mode 100755 index 000000000..4b29bc5ed --- /dev/null +++ b/3014/CH3/EX3.12/Ex3_12.txt @@ -0,0 +1,3 @@ +
+Example 3.12
+ Wavelength shift is 0.027143 angstrom.
diff --git a/3014/CH3/EX3.13/Ex3_13.sce b/3014/CH3/EX3.13/Ex3_13.sce new file mode 100755 index 000000000..62561d2d8 --- /dev/null +++ b/3014/CH3/EX3.13/Ex3_13.sce @@ -0,0 +1,14 @@ + clc
+//given that
+lambda = 1 //wavelength in angstrom
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.13")
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda // Calculation of recoiled electron wavelength
+d_E = h*c*(lambda_n-lambda)*1e10/(1.6e-19*lambda_n*lambda)// Calculation of recoiled electron energy in eV
+printf("\nWavelength shift is %f angstrom.",lambda_n)
+printf("\nEnergy of recoiled electron is %deV. \n\n\n",ceil (d_E))
+
diff --git a/3014/CH3/EX3.13/Ex3_13.txt b/3014/CH3/EX3.13/Ex3_13.txt new file mode 100755 index 000000000..785500b08 --- /dev/null +++ b/3014/CH3/EX3.13/Ex3_13.txt @@ -0,0 +1,5 @@ +
+
+Example 3.13
+Wavelength shift is 1.024286 angstrom.
+Energy of recoiled electron is 295eV.
diff --git a/3014/CH3/EX3.14/Ex3_14.sce b/3014/CH3/EX3.14/Ex3_14.sce new file mode 100755 index 000000000..eeb0c7a07 --- /dev/null +++ b/3014/CH3/EX3.14/Ex3_14.sce @@ -0,0 +1,14 @@ + clc
+//given that
+lambda = 1 //let wavelength in angstrom
+lambda_n = 2*lambda // recoiled electron wavelength
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.14")
+lambda = h*1e10/(m_e*c) // calculation of wavelength in angstrom
+E = h*c*1e10/(lambda*1.6e-19) // calculation of energy of electron
+
+printf("\nWavelength shift is %f angstrom. ",lambda)
+printf("\nEnergy of recoiled electron is %f KeV. \n\n\n",E/1e3)
diff --git a/3014/CH3/EX3.14/Ex3_14.txt b/3014/CH3/EX3.14/Ex3_14.txt new file mode 100755 index 000000000..71ce2dbf4 --- /dev/null +++ b/3014/CH3/EX3.14/Ex3_14.txt @@ -0,0 +1,4 @@ +
+Example 3.14
+Wavelength shift is 0.024286 angstrom.
+Energy of recoiled electron is 511.875000eV.
\ No newline at end of file diff --git a/3014/CH3/EX3.15/Ex3_15.sce b/3014/CH3/EX3.15/Ex3_15.sce new file mode 100755 index 000000000..22f7c2704 --- /dev/null +++ b/3014/CH3/EX3.15/Ex3_15.sce @@ -0,0 +1,14 @@ + clc
+//given that
+lambda = 2 //wavelength in angstrom
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 45 // scattering angle
+printf("Example 3.15")
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda // Calculation of recoiled electron wavelength
+
+f = d_lambda/lambda // Calculation of fraction of energy lost by photon
+
+printf("\nFraction of energy lost by photon is %f\n\n\n",f)
diff --git a/3014/CH3/EX3.15/Ex3_15.txt b/3014/CH3/EX3.15/Ex3_15.txt new file mode 100755 index 000000000..67a0a2c3f --- /dev/null +++ b/3014/CH3/EX3.15/Ex3_15.txt @@ -0,0 +1,4 @@ +
+
+Example 3.15
+Fraction of energy lost by photon is 0.003557
diff --git a/3014/CH3/EX3.16/Ex3_16.sce b/3014/CH3/EX3.16/Ex3_16.sce new file mode 100755 index 000000000..8187d02cc --- /dev/null +++ b/3014/CH3/EX3.16/Ex3_16.sce @@ -0,0 +1,17 @@ +
+ clc
+//given that
+E_eV = 510 // Energy of gamma ray in keV
+lambda = 2 //wavelength in angstrom
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // scattering angle in degree
+printf("Example 3.16")
+E_j = E_eV*1e3*1.6e-19 // Energy of gamma ray in Joule
+lambda = h*c*1e10/E_j // Calculation of wavelength in angstrom
+
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda // Calculation of recoiled electron wavelength
+printf("\nWavelength of scattered radiation is %f Angstrom \n\n\n",lambda_n)
+
diff --git a/3014/CH3/EX3.16/Ex3_16.txt b/3014/CH3/EX3.16/Ex3_16.txt new file mode 100755 index 000000000..a4a91a0d8 --- /dev/null +++ b/3014/CH3/EX3.16/Ex3_16.txt @@ -0,0 +1,3 @@ +
+Example 3.16
+Wavelength of scattered radiation is 0.048661 Angstrom.
diff --git a/3014/CH3/EX3.17/Ex3_17.sce b/3014/CH3/EX3.17/Ex3_17.sce new file mode 100755 index 000000000..e08f1a852 --- /dev/null +++ b/3014/CH3/EX3.17/Ex3_17.sce @@ -0,0 +1,14 @@ +
+ clc
+//given that
+lambda = 2 //wavelength in angstrom
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.17")
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda // Calculation of recoiled electron wavelength
+d_E = h*c*(lambda_n-lambda)*1e10/(1.6e-19*lambda_n*lambda)// Calculation of recoiled electron energy in eV
+printf("\n Scattered wavelength is %f angstrom.",lambda_n)
+printf("\n Energy of recoiled electron is %feV. \n\n\n",d_E)
diff --git a/3014/CH3/EX3.17/Ex3_17.txt b/3014/CH3/EX3.17/Ex3_17.txt new file mode 100755 index 000000000..71b50f1ea --- /dev/null +++ b/3014/CH3/EX3.17/Ex3_17.txt @@ -0,0 +1,4 @@ +
+Example 3.17
+ Scattered wavelength is 2.024286 angstrom.
+ Energy of recoiled electron is 74.569954eV.
diff --git a/3014/CH3/EX3.18/Ex3_18.sce b/3014/CH3/EX3.18/Ex3_18.sce new file mode 100755 index 000000000..d07677c97 --- /dev/null +++ b/3014/CH3/EX3.18/Ex3_18.sce @@ -0,0 +1,21 @@ +
+ clc
+//given that
+E_eV = 510 // Energy of gamma ray in keV
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // scattering angle in degree
+printf("Example 3.18")
+E_j = E_eV*1e3*1.6e-19 // Energy of gamma ray in Joule
+lambda = h*c/E_j // Calculation of wavelength in meter
+
+d_lambda= h*(1-cos(theta*%pi/180))*1e10/(m_e*c) // calculation of wavelength shift in angstrom
+lambda_n = lambda+d_lambda/1e10 // Calculation of recoiled electron wavelength
+d_E = h*c*(d_lambda/1e10)/(1.6e-19*lambda_n*lambda)// Calculation of recoiled electron energy in eV
+psi= atan(1/(tan((theta*%pi/180)/2)/(1+(h/(lambda*m_e*c)))))
+phi_deg = 90 - psi*180/%pi // Calculation of degree part of angle of recoiled electron
+phi_min = 60*(phi_deg - floor(phi_deg))// Calculation of minute part of angle of recoiled electron
+printf("\nWavelength of scattered radiation is %e m ",lambda_n)
+printf("\nEnergy of recoiled electron is %f MeV.",d_E/1e6)
+printf("\nRecoiled electron angle is %d degree%d minute \n\n\n",phi_deg,phi_min)
diff --git a/3014/CH3/EX3.18/Ex3_18.txt b/3014/CH3/EX3.18/Ex3_18.txt new file mode 100755 index 000000000..6af012803 --- /dev/null +++ b/3014/CH3/EX3.18/Ex3_18.txt @@ -0,0 +1,5 @@ +
+Example 3.18
+Wavelength of scattered radiation is 4.866071e-12 m
+Energy of recoiled electron is 0.254532 MeV.
+Recoiled electron angle is 26 degree36 minute
diff --git a/3014/CH3/EX3.19/Ex3_19.sce b/3014/CH3/EX3.19/Ex3_19.sce new file mode 100755 index 000000000..19d4d0f79 --- /dev/null +++ b/3014/CH3/EX3.19/Ex3_19.sce @@ -0,0 +1,13 @@ + clc
+//given that
+nu = 2e19 // initial frequency of X ray photon
+h = 6.63e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+c = 3e8 // speed of light in m/sec
+theta = 90 // scattering angle in degree
+printf("Example 3.19")
+d_lambda = h/(m_e*c) // calculation of wavelength shift
+k = 1/nu + d_lambda/c
+nu_1 = 1/k // Frequency after collision
+nu_1 = floor(nu_1/1e18)*1e18 // rounding off
+printf("\nFrequency after collision is %e Hz \n\n\n",nu_1)
diff --git a/3014/CH3/EX3.19/Ex3_19.txt b/3014/CH3/EX3.19/Ex3_19.txt new file mode 100755 index 000000000..436a0845f --- /dev/null +++ b/3014/CH3/EX3.19/Ex3_19.txt @@ -0,0 +1,4 @@ +
+
+Example 3.19
+Frequency after collision is 1.721311e+19 Hz
\ No newline at end of file diff --git a/3014/CH3/EX3.2/Ex3_2.sce b/3014/CH3/EX3.2/Ex3_2.sce new file mode 100755 index 000000000..94ec1ffef --- /dev/null +++ b/3014/CH3/EX3.2/Ex3_2.sce @@ -0,0 +1,13 @@ + clc
+//given that
+lambda = 0.3 // Wavelength in angstrom
+d = 0.5 // crystal spacing in angstrom
+n = 2 // order
+m = 3 // order
+printf("Example 3.2")
+theta_n = asin(n*lambda/(2*d))*180/%pi // Calculation of angle for order n
+theta_m = asin(m*lambda/(2*d))*180/%pi // Calculation of angle for order m
+
+printf("\nAngle for %dnd order maxima is %f degree. ",n,theta_n)
+printf("\nAngle for %drd order maxima is %f degree. \n\n\n",m,theta_m)
+// Answers in book are 40.97 degree and 72.29 degree which are due to wrong calculation
diff --git a/3014/CH3/EX3.2/Ex3_2.txt b/3014/CH3/EX3.2/Ex3_2.txt new file mode 100755 index 000000000..2e43f4bf0 --- /dev/null +++ b/3014/CH3/EX3.2/Ex3_2.txt @@ -0,0 +1,4 @@ +
+Example 3.2
+Angle for 2nd order maxima is 36.869898 degree.
+Angle for 3rd order maxima is 64.158067 degree.
\ No newline at end of file diff --git a/3014/CH3/EX3.3/Ex3_3.sce b/3014/CH3/EX3.3/Ex3_3.sce new file mode 100755 index 000000000..ff24f73bf --- /dev/null +++ b/3014/CH3/EX3.3/Ex3_3.sce @@ -0,0 +1,12 @@ +
+
+ clc
+//given that
+d = 1.87 // crystal spacing in angstrom
+n = 2 // order for longest passing wavelength
+theta = 30 // angle for longest passing wavelength
+printf("Example 3.3")
+lambda = 2*d*sin(theta*%pi/180)/n // Calculation of longest wavelength
+
+printf("\n Longest wavelength is %f angstrom. \n\n\n",lambda)
+
diff --git a/3014/CH3/EX3.3/Ex3_3.txt b/3014/CH3/EX3.3/Ex3_3.txt new file mode 100755 index 000000000..5efa2ef06 --- /dev/null +++ b/3014/CH3/EX3.3/Ex3_3.txt @@ -0,0 +1,3 @@ +
+Example 3.3
+ Longest wavelength is 0.935000 angstrom.
\ No newline at end of file diff --git a/3014/CH3/EX3.4/Ex3_4.sce b/3014/CH3/EX3.4/Ex3_4.sce new file mode 100755 index 000000000..3b29a034a --- /dev/null +++ b/3014/CH3/EX3.4/Ex3_4.sce @@ -0,0 +1,11 @@ + clc
+//given that
+lambda = 3.6e-9 // Wavelength in cm
+theta = 4.8 // glancing angle in degree
+n = 1 // order
+
+printf("Example 3.4")
+d = n*lambda/(2*sin(theta*%pi/180)) // calculation of crystal spacing in angstrom
+
+printf("\n Crystal spacing in angstrom is %e cm. \n\n\n",d)
+
diff --git a/3014/CH3/EX3.4/Ex3_4.txt b/3014/CH3/EX3.4/Ex3_4.txt new file mode 100755 index 000000000..60ae3c83f --- /dev/null +++ b/3014/CH3/EX3.4/Ex3_4.txt @@ -0,0 +1,3 @@ +
+Example 3.4
+ Crystal spacing in angstrom is 2.151107e-08 cm.
\ No newline at end of file diff --git a/3014/CH3/EX3.5/Ex3_5.sce b/3014/CH3/EX3.5/Ex3_5.sce new file mode 100755 index 000000000..c30bc12b7 --- /dev/null +++ b/3014/CH3/EX3.5/Ex3_5.sce @@ -0,0 +1,11 @@ + clc
+//given that
+d = 2.5 // crystal spacing in angstrom
+n = 1 // order for longest passing wavelength
+theta = 20 // angle for longest passing wavelength
+printf("Example 3.5")
+lambda = 2*d*sin(theta*%pi/180)/n // Calculation of longest wavelength
+
+printf("\nLongest wavelength is %f angstrom. \n\n\n",lambda)
+
+
diff --git a/3014/CH3/EX3.5/Ex3_5.txt b/3014/CH3/EX3.5/Ex3_5.txt new file mode 100755 index 000000000..ed2b4dde9 --- /dev/null +++ b/3014/CH3/EX3.5/Ex3_5.txt @@ -0,0 +1,4 @@ +
+
+Example 3.5
+Longest wavelength is 1.710101 angstrom
\ No newline at end of file diff --git a/3014/CH3/EX3.6/Ex3_6.sce b/3014/CH3/EX3.6/Ex3_6.sce new file mode 100755 index 000000000..37961fd77 --- /dev/null +++ b/3014/CH3/EX3.6/Ex3_6.sce @@ -0,0 +1,9 @@ + clc
+//given that
+d = 2.5 // crystal spacing in angstrom
+n = 1 // order for longest passing wavelength
+theta = 90 // angle for longest passing wavelength
+printf("Example 3.6")
+lambda = 2*d*sin(theta*%pi/180)/n // Calculation of longest wavelength
+
+printf("\nLongest wavelength is of %d angstrom. \n\n\n",lambda)
diff --git a/3014/CH3/EX3.6/Ex3_6.txt b/3014/CH3/EX3.6/Ex3_6.txt new file mode 100755 index 000000000..c8c9bd0f7 --- /dev/null +++ b/3014/CH3/EX3.6/Ex3_6.txt @@ -0,0 +1,3 @@ +
+Example 3.6
+Longest wavelength is 5 angstrom.
diff --git a/3014/CH3/EX3.7/Ex3_7.sce b/3014/CH3/EX3.7/Ex3_7.sce new file mode 100755 index 000000000..a2e097f82 --- /dev/null +++ b/3014/CH3/EX3.7/Ex3_7.sce @@ -0,0 +1,18 @@ +clc
+// given that
+theta1_deg = 5 // Absolut degree part of angle for first angle
+theta1_min = 23//remainder minute part of angle for first angle
+theta2_deg = 7 // Absolut degree part of angle for second angle
+theta2_min = 37//remainder minute part of angle for second angle
+theta3_deg = 9 // Absolut degree part of angle for third angle
+theta3_min = 25//remainder minute part of angle for third angle
+
+printf("Example 3.7 \n")
+val1 = sin((theta1_deg+ theta1_min/60)*%pi/180)// Sin value for first angle
+val2 = sin((theta2_deg+ theta2_min/60)*%pi/180) //Sin value for second angle
+val3 = sin((theta3_deg+ theta3_min/60)*%pi/180)//Sin value for third angle
+ratio_21 = val2/val1
+ratio_31 = val3/val1
+printf("\n Interatomic layer separation ratios in crystal are as\n 1 : %f : %f",ratio_21,ratio_31)
+printf("\n Above relation shows that crystal is simple cubic crystal structure.")
+
diff --git a/3014/CH3/EX3.7/Ex3_7.txt b/3014/CH3/EX3.7/Ex3_7.txt new file mode 100755 index 000000000..767c63626 --- /dev/null +++ b/3014/CH3/EX3.7/Ex3_7.txt @@ -0,0 +1,5 @@ +Example 3.7
+
+ Interatomic layer seperation ratios in crystal are as
+ 1 : 1.412775 : 1.743926
+ Above relation shows that cristal is simple cubic crystal structure.
\ No newline at end of file diff --git a/3014/CH3/EX3.8/Ex3_8.sce b/3014/CH3/EX3.8/Ex3_8.sce new file mode 100755 index 000000000..16e967980 --- /dev/null +++ b/3014/CH3/EX3.8/Ex3_8.sce @@ -0,0 +1,13 @@ + clc
+//given that
+lambda = 1.2 // wavelength in angstrom
+theta_deg = 9 // angle fraction in degree
+theta_min = 30 // Angle fraction in minute
+printf("Example 3.8\n")
+theta = theta_deg+theta_min/60 // Total angel
+for n = 1:4
+ d = lambda/(n*2*sin(theta*%pi/180)) // Inter layer spacing
+ printf(" If order is %d then spacing is %f angstrom.\n",n,d)
+end
+
+
diff --git a/3014/CH3/EX3.8/Ex3_8.txt b/3014/CH3/EX3.8/Ex3_8.txt new file mode 100755 index 000000000..5d2b34f98 --- /dev/null +++ b/3014/CH3/EX3.8/Ex3_8.txt @@ -0,0 +1,6 @@ +
+Example 3.8
+ If order is 1 then spacing is 3.635315 angstrom.
+ If order is 2 then spacing is 1.817657 angstrom.
+ If order is 3 then spacing is 1.211772 angstrom.
+ If order is 4 then spacing is 0.908829 angstrom.
\ No newline at end of file diff --git a/3014/CH3/EX3.9/Ex3_9.sce b/3014/CH3/EX3.9/Ex3_9.sce new file mode 100755 index 000000000..da3aba8ef --- /dev/null +++ b/3014/CH3/EX3.9/Ex3_9.sce @@ -0,0 +1,13 @@ + clc
+//given that
+h = 6.62e-34 // Planks constant
+m_e = 9.1e-31 // mass of electron in kg
+e = 1.6e-19 // charge on electron in coulomb
+v = 340 // Applied voltage in volt
+n = 1 // order for longest passing wavelength
+theta = 60 // angle for longest passing wavelength
+printf("Example 3.9")
+lambda= h/sqrt(2*m_e*e*v) // calculation of wavelength
+d = n*lambda/(2*sin(theta*%pi/180))// calculation of spacing of crystal
+
+printf("\nSpacing of crystal is %f angstrom. \n\n\n",d*1e10)
diff --git a/3014/CH3/EX3.9/Ex3_9.txt b/3014/CH3/EX3.9/Ex3_9.txt new file mode 100755 index 000000000..f02627cb2 --- /dev/null +++ b/3014/CH3/EX3.9/Ex3_9.txt @@ -0,0 +1,3 @@ +
+Example 3.9
+Spacing of crystal is 0.384116 angstrom.
diff --git a/3014/CH4/EX4.10/Ex4_10.sce b/3014/CH4/EX4.10/Ex4_10.sce new file mode 100755 index 000000000..03b5a08c0 --- /dev/null +++ b/3014/CH4/EX4.10/Ex4_10.sce @@ -0,0 +1,12 @@ +clc
+//Given that
+epsilon = 2.4e-10 // permitivity of a dielectric material in C^2/N?m^2
+epsilon_0 = 8.854e-12 // Permittivity of free space
+
+printf("Example 4.10")
+K = epsilon/epsilon_0 // Calculation of dielectric constant
+zai_e = epsilon_0*(K-1) // Calculation of electrical susceptibility
+
+printf("\n Relative permittivity is %f",K)
+printf("\n Electrical susceptibility is %e C^2/Nm^2\n\n\n",zai_e)
+
diff --git a/3014/CH4/EX4.10/Ex4_10.txt b/3014/CH4/EX4.10/Ex4_10.txt new file mode 100755 index 000000000..d642b7f3b --- /dev/null +++ b/3014/CH4/EX4.10/Ex4_10.txt @@ -0,0 +1,4 @@ +
+Example 4.10
+ Relative permittivity is 27.106393
+ Electrical susceptibility is 2.311460e-10 C^2/N/m^2
\ No newline at end of file diff --git a/3014/CH4/EX4.11/Ex4_11.sce b/3014/CH4/EX4.11/Ex4_11.sce new file mode 100755 index 000000000..f496914e8 --- /dev/null +++ b/3014/CH4/EX4.11/Ex4_11.sce @@ -0,0 +1,20 @@ +
+clc
+//Given that
+V = 100 // Applied potential in Volt
+d = 1 // Separation between plates in cm
+k1 = 8 // Dielectric constant
+k2 = 9 //dielectric constant
+epsilon_0 = 8.854e-12 // Permittivity of free space
+
+printf("Example 4.11")
+E_0 = V/(d*1e-2) // Calculation of electric field
+E = E_0/k1*k2 // Calculation of electric field
+D = k1*epsilon_0*E // Calculation of electrical displacement vector
+P = (k1-1)*epsilon_0*E // Calculation of electrical polarization
+
+printf("\n Magnitude of Electrical vector is %e Volt/meter",E) // Answer in book is 1.125e3 Volt/meter
+
+printf("\n Magnitude of Electrical Displacement vector is %e C/m^2",D)// Answer in book is 8.85e-8C/m^2
+
+printf("\n Magnitude of Electric polarization vector is %e C/m^2\n\n\n",P)// Answer in book is 7.774e-8C/m^2
diff --git a/3014/CH4/EX4.11/Ex4_11.txt b/3014/CH4/EX4.11/Ex4_11.txt new file mode 100755 index 000000000..71acc4ddd --- /dev/null +++ b/3014/CH4/EX4.11/Ex4_11.txt @@ -0,0 +1,6 @@ +
+Example 4.11
+ Magnitude of Electrical vector is 1.125000e+04 Volt/meter.
+ Magnitude of Electrical Displacement vector is 7.968600e-07 C/m^2.
+ Magnitude of Electric polarization vector is 6.972525e-07 C/m^2.
+
diff --git a/3014/CH4/EX4.12/Ex4_12.sce b/3014/CH4/EX4.12/Ex4_12.sce new file mode 100755 index 000000000..af8cce616 --- /dev/null +++ b/3014/CH4/EX4.12/Ex4_12.sce @@ -0,0 +1,15 @@ +
+clc
+//Given that
+alpha_300 = 2.5e-39 // total polarisability in C^2m/N at 300 K
+alpha_600 = 1.75e-39 // total polarisability in C^2m/N at 600 K
+T1 = 300 // Initial temperature in Kelvin
+T2 = 600 // Final Temperature in Kelvin
+printf("Example 4.12\n")
+b = (alpha_300-alpha_600)*T2
+al_def_300 = alpha_300 - b/300
+al_oriant_300 = b/300
+al_oriant_600 = b/600
+printf("\n Deformational Polarizability is %e C^2mN^-1",al_def_300)
+printf("\n Orientational Polarizability at %d degree Celcius is %e C^2mN^-1",T1,al_oriant_300)
+printf("\n Orientational Polarizability at %d degree Celcius is %e C^2mN^-1",T2,al_oriant_600)
diff --git a/3014/CH4/EX4.12/Ex4_12.txt b/3014/CH4/EX4.12/Ex4_12.txt new file mode 100755 index 000000000..6b72216c9 --- /dev/null +++ b/3014/CH4/EX4.12/Ex4_12.txt @@ -0,0 +1,5 @@ +Example 4.12
+
+ Deformational Polarizability is 1.000000e-39 C^2mN^-1
+ Orientational Polarizability at 300 degree Celcius is 1.500000e-39 C^2mN^-1
+ Orientational Polarizability at 600 degree Celcius is 7.500000e-40 C^2mN^-1
\ No newline at end of file diff --git a/3014/CH4/EX4.13/Ex4_13.sce b/3014/CH4/EX4.13/Ex4_13.sce new file mode 100755 index 000000000..3f86677a1 --- /dev/null +++ b/3014/CH4/EX4.13/Ex4_13.sce @@ -0,0 +1,12 @@ +
+clc
+//Given that
+alpha_e = 1.5e-40 // Electronic polarizability in Fm^2
+N = 4e28 // density in atoms/m^3
+epsilon_0 = 8.85e-12 // Permittivity of free space
+
+printf("Example 4.13")
+k = N*alpha_e/(3*epsilon_0)
+epsilon_r = (1+ k*2)/(1-k)// Calculation of relative permittivity
+printf("\n Relative permittivity is %f\n\n\n",epsilon_r)
+
diff --git a/3014/CH4/EX4.13/Ex4_13.txt b/3014/CH4/EX4.13/Ex4_13.txt new file mode 100755 index 000000000..93c1c0480 --- /dev/null +++ b/3014/CH4/EX4.13/Ex4_13.txt @@ -0,0 +1,2 @@ +Example 4.13
+ Relative permittivity is 1.875912
\ No newline at end of file diff --git a/3014/CH4/EX4.14/Ex4_14.sce b/3014/CH4/EX4.14/Ex4_14.sce new file mode 100755 index 000000000..73264f4cc --- /dev/null +++ b/3014/CH4/EX4.14/Ex4_14.sce @@ -0,0 +1,16 @@ +
+clc
+//Given that
+m = 32 // Atomic weight of sulphur
+d = 2.08 // Density in g/cm^3
+alpha_e = 3.5e-40 // Electronic polarizability in Fm^2
+N_a = 6.022e23 // Avogadro Number
+epsilon_0 = 8.85e-12 // Permittivity of free space
+
+printf("Example 4.14")
+N = N_a*d*1e6/m // Calculation of Atoms per unit
+k = N*alpha_e/(3*epsilon_0)
+
+epsilon_r = (1+ k*2)/(1-k)// Calculation of relative permittivity
+printf("\n Relative permittivity is %f\n\n\n",epsilon_r)
+// Answer in book is 4.17
diff --git a/3014/CH4/EX4.14/Ex4_14.txt b/3014/CH4/EX4.14/Ex4_14.txt new file mode 100755 index 000000000..31ae3df6a --- /dev/null +++ b/3014/CH4/EX4.14/Ex4_14.txt @@ -0,0 +1,4 @@ +
+Example 4.14
+ Relative permittivity is 4.198468
+
diff --git a/3014/CH4/EX4.15/Ex4_15.sce b/3014/CH4/EX4.15/Ex4_15.sce new file mode 100755 index 000000000..633f55bbc --- /dev/null +++ b/3014/CH4/EX4.15/Ex4_15.sce @@ -0,0 +1,9 @@ +
+clc
+//Given that
+n = 1.5 // Refractive index
+epsilon = 5.6 // Static dielectric constant
+printf("Example 4.15")
+per = (1-((n^2-1)/(n^2+2))*(epsilon+2)/(epsilon-1))*100 // Pecentage ionic polarisability
+printf("\n Percentage ionic polarizability is %f pecent\n\n\n",per)
+// Answer in book is 5.14 %
diff --git a/3014/CH4/EX4.15/Ex4_15.txt b/3014/CH4/EX4.15/Ex4_15.txt new file mode 100755 index 000000000..35704226e --- /dev/null +++ b/3014/CH4/EX4.15/Ex4_15.txt @@ -0,0 +1,3 @@ +
+Example 4.15
+ Percentage ionic polarizability is 51.406650 pecent
diff --git a/3014/CH4/EX4.16/Ex4_16.sce b/3014/CH4/EX4.16/Ex4_16.sce new file mode 100755 index 000000000..d92b510e5 --- /dev/null +++ b/3014/CH4/EX4.16/Ex4_16.sce @@ -0,0 +1,16 @@ +
+
+clc
+//Given that
+m = 32 // Atomic weight of sulphur
+d = 2050 // Density in Kg/m^3
+N_a = 6.022e23 // Avogadro Number
+epsilon_0 = 8.85e-12 // Permittivity of free space
+epsilon_r = 3.75 // Dielectric constant of sulphur
+
+printf("Example 4.16")
+N = N_a*d*1e3/m // Calculation of Atoms per unit
+alpha_e = 3*epsilon_0*((epsilon_r-1)/(epsilon_r+2)) / N
+
+
+printf("\n Electronic polarizability is %e Fm^2\n\n\n",alpha_e)
diff --git a/3014/CH4/EX4.16/Ex4_16.txt b/3014/CH4/EX4.16/Ex4_16.txt new file mode 100755 index 000000000..3fe386f75 --- /dev/null +++ b/3014/CH4/EX4.16/Ex4_16.txt @@ -0,0 +1,3 @@ +
+Example 4.16
+ Electronic polarizability is 3.291431e-40 Fm^2
diff --git a/3014/CH4/EX4.17/Ex4_17.sce b/3014/CH4/EX4.17/Ex4_17.sce new file mode 100755 index 000000000..345d13344 --- /dev/null +++ b/3014/CH4/EX4.17/Ex4_17.sce @@ -0,0 +1,12 @@ +
+clc
+//Given that
+n = 1.5 // Refractive index
+epsilon = 4 // Static dielectric constant
+epsilon_0 = 8.85e-12 // permittivity of free space
+printf("Example 4.17")
+k1 = (epsilon-1)/(epsilon+2)
+k2 = (n^2-1)/(n^2+2)
+ratio = 1/((k1/k2)-1)
+printf("\n Ratio of electronic to ionic polarizability is %f .\n\n\n",ratio)
+// Answer in book is 1.43
diff --git a/3014/CH4/EX4.17/Ex4_17.txt b/3014/CH4/EX4.17/Ex4_17.txt new file mode 100755 index 000000000..ca5b0709a --- /dev/null +++ b/3014/CH4/EX4.17/Ex4_17.txt @@ -0,0 +1,3 @@ +
+Example 4.17
+ Ratio of electronic to ionic polarizability is 1.428571 .
diff --git a/3014/CH4/EX4.18/Ex4_18.sce b/3014/CH4/EX4.18/Ex4_18.sce new file mode 100755 index 000000000..a894f4969 --- /dev/null +++ b/3014/CH4/EX4.18/Ex4_18.sce @@ -0,0 +1,13 @@ +
+clc
+//Given that
+t = 1.8e-5 // Relaxation time in second
+epsilon_r = 1 // let
+printf("Example 4.18")
+f = 1/(2*%pi*t) // Calculation of frequency
+delta = atan(epsilon_r/epsilon_r)
+phi = 90 - delta*180/%pi // Calculation of phase difference
+printf("\n Frequency is %f KHz\n",f/1e3)
+printf(" Phase difference between current and voltage is %d degree.",phi)
+
+
diff --git a/3014/CH4/EX4.18/Ex4_18.txt b/3014/CH4/EX4.18/Ex4_18.txt new file mode 100755 index 000000000..273aff63b --- /dev/null +++ b/3014/CH4/EX4.18/Ex4_18.txt @@ -0,0 +1,4 @@ +
+Example 4.18
+ Frequency is 8.841941 KHz
+ Phase difference between current and voltage is 45 degree.
\ No newline at end of file diff --git a/3014/CH4/EX4.4/Ex4_4.sce b/3014/CH4/EX4.4/Ex4_4.sce new file mode 100755 index 000000000..bd458f3d5 --- /dev/null +++ b/3014/CH4/EX4.4/Ex4_4.sce @@ -0,0 +1,16 @@ +
+clc
+//Given that
+epsilon_r = 1.000074 // Dielectric constant of He at 0C and 1atm
+epsilon_0 = 8.854e-12 // Permittivity of free space
+E = 100 // Electric field in V/m
+n = 2.68e27 // Electron density in no,/m^
+N_a = 6e23 // Avogadro number
+V = 22.4 // Volume at STP in litter
+printf("Example 4.4")
+P = epsilon_0*(epsilon_r-1)*E // Calculation of polarization
+
+N = N_a/(V*1e-3)// Calculation of total number of atoms
+p = P/N // dipole moment per atom
+printf("\n Dipole moment per atom is %e Coulomb-meter \n\n\n",p)
+// Answer in book is in different form and as 24.45e-40 coulomb-meter
diff --git a/3014/CH4/EX4.4/Ex4_4.txt b/3014/CH4/EX4.4/Ex4_4.txt new file mode 100755 index 000000000..af03c5d31 --- /dev/null +++ b/3014/CH4/EX4.4/Ex4_4.txt @@ -0,0 +1,2 @@ +Example 4.4
+ Dipole moment per atom is 2.446065e-39 Coulomb-meter
diff --git a/3014/CH4/EX4.6/Ex4_6.sce b/3014/CH4/EX4.6/Ex4_6.sce new file mode 100755 index 000000000..2fdd4a31a --- /dev/null +++ b/3014/CH4/EX4.6/Ex4_6.sce @@ -0,0 +1,16 @@ +
+clc
+//Given that
+r = 0.055 // Radius of hydrogen atom in nm
+n = 9.8e26 // Number of atoms/cc
+
+epsilon_0 = 8.854e-12 // Permittivity of free space
+
+printf("Example 4.6")
+alpha_e = 4*%pi*epsilon_0*(r*1e-9)^3 // Calculation of electronic polarisability
+epsilon_r = 1+n*alpha_e/epsilon_0 // Calculation of relative permeability
+
+printf("\n Electronic polarisability is %eFm^2 \n Relative permeability is %f \n\n\n",alpha_e,epsilon_r)
+
+
+
diff --git a/3014/CH4/EX4.6/Ex4_6.txt b/3014/CH4/EX4.6/Ex4_6.txt new file mode 100755 index 000000000..1d24f515a --- /dev/null +++ b/3014/CH4/EX4.6/Ex4_6.txt @@ -0,0 +1,5 @@ +
+Example 4.6
+ Electronic polarisability is 1.851132e-41Fm^2
+ Relative permeability is 1.002049
+
diff --git a/3014/CH4/EX4.8/Ex4_8.sce b/3014/CH4/EX4.8/Ex4_8.sce new file mode 100755 index 000000000..abae3cbe3 --- /dev/null +++ b/3014/CH4/EX4.8/Ex4_8.sce @@ -0,0 +1,12 @@ +
+clc
+//Given that
+epsilon_0 = 8.854e-12 // Permittivity of free space
+E = 2000 // Electric field in V/m
+P = 6.4e-8 // Polarization in C/m^2
+printf("Example 4.8")
+epsilon_r = 1+ P/(epsilon_0*E) // Calculation of relative permittivity
+
+printf("\n Relative permittivity is %f\n\n\n",epsilon_r)
+
+
diff --git a/3014/CH4/EX4.8/Ex4_8.txt b/3014/CH4/EX4.8/Ex4_8.txt new file mode 100755 index 000000000..a27dbe60f --- /dev/null +++ b/3014/CH4/EX4.8/Ex4_8.txt @@ -0,0 +1,3 @@ +
+Example 4.8
+ Relative permittivity is 4.614186
\ No newline at end of file diff --git a/3014/CH4/EX4.9/Ex4_9.sce b/3014/CH4/EX4.9/Ex4_9.sce new file mode 100755 index 000000000..103e6bd00 --- /dev/null +++ b/3014/CH4/EX4.9/Ex4_9.sce @@ -0,0 +1,13 @@ +
+clc
+//Given that
+alpha_e = 2e-40 // Electronic polarisability in Fm^2
+N = 4e28 // density in atoms/m^3
+epsilon_0 = 8.85e-12 // Permittivity of free space
+
+printf("Example 4.9")
+epsilon_r = 1+ N*alpha_e/(epsilon_0) // Calculation of relative permittivity
+printf("\n Relative permittivity is %f\n\n\n",epsilon_r)
+
+
+
diff --git a/3014/CH4/EX4.9/Ex4_9.txt b/3014/CH4/EX4.9/Ex4_9.txt new file mode 100755 index 000000000..b2f1b8b79 --- /dev/null +++ b/3014/CH4/EX4.9/Ex4_9.txt @@ -0,0 +1,3 @@ +
+Example 4.9
+ Relative permittivity is 1.903955
diff --git a/3014/CH6/EX6.1/Ex6_1.sce b/3014/CH6/EX6.1/Ex6_1.sce new file mode 100755 index 000000000..2b42ec99a --- /dev/null +++ b/3014/CH6/EX6.1/Ex6_1.sce @@ -0,0 +1,11 @@ +
+clc
+//Given that
+E = 7.9e10 // Young’s modulus in N/m^2
+rho = 2650 // Density in Kg/m^3
+t = 0.003 // Thickness of quartz crystal in m
+printf("Example 6.1\n")
+v = sqrt(E/rho)// Calculation of velocity
+lambda = 2*t // Calculation of fundamental wavelength
+nu = v/lambda // Calculation of fundamental frequency
+printf("Fundamental frequency is %e Hz.\n\n\n",nu)
diff --git a/3014/CH6/EX6.1/Ex6_1.txt b/3014/CH6/EX6.1/Ex6_1.txt new file mode 100755 index 000000000..8480c8df7 --- /dev/null +++ b/3014/CH6/EX6.1/Ex6_1.txt @@ -0,0 +1,2 @@ + Example 6.1
+Fundamental frequency is 9.099957e+05 Hz.
\ No newline at end of file diff --git a/3014/CH6/EX6.2/Ex6_2.sce b/3014/CH6/EX6.2/Ex6_2.sce new file mode 100755 index 000000000..875b9ca11 --- /dev/null +++ b/3014/CH6/EX6.2/Ex6_2.sce @@ -0,0 +1,9 @@ +
+clc
+//Given that
+v = 5760 // Velocity in m/s
+T = 1.6 // Thickness of quartz crystal in mm
+printf("Example 6.2\n")
+nu = v/(2*T*1e-3)// Calculation of fundamental frequency
+printf("Fundamental frequency of crystal is %f MHz.\n\n\n",nu/1e6)
+
diff --git a/3014/CH6/EX6.2/Ex6_2.txt b/3014/CH6/EX6.2/Ex6_2.txt new file mode 100755 index 000000000..794c8f6ed --- /dev/null +++ b/3014/CH6/EX6.2/Ex6_2.txt @@ -0,0 +1,3 @@ +
+Example 6.2
+Fundamental frequency of crystal is 1.800000 MHz.
diff --git a/3014/CH6/EX6.3/Ex6_3.sce b/3014/CH6/EX6.3/Ex6_3.sce new file mode 100755 index 000000000..51313c93d --- /dev/null +++ b/3014/CH6/EX6.3/Ex6_3.sce @@ -0,0 +1,9 @@ +
+clc
+//Given that
+T =40 // Thickness of steel bar in cm
+t1 = 40 // Time in ms
+t2 = 80 // Time in ms
+printf("Example 6.3\n")
+X = T*t1/t2 // Calculation of depth of defect
+printf("Depth of defect is %d cm.\n\n\n",X)
diff --git a/3014/CH6/EX6.3/Ex6_3.txt b/3014/CH6/EX6.3/Ex6_3.txt new file mode 100755 index 000000000..0644bda81 --- /dev/null +++ b/3014/CH6/EX6.3/Ex6_3.txt @@ -0,0 +1,3 @@ +
+Example 6.3
+Depth of defect is 20 cm.
diff --git a/3014/CH6/EX6.4/Ex6_4.sce b/3014/CH6/EX6.4/Ex6_4.sce new file mode 100755 index 000000000..ce81ea668 --- /dev/null +++ b/3014/CH6/EX6.4/Ex6_4.sce @@ -0,0 +1,11 @@ +clc
+//Given that
+E = 7.9e10 // Young’s modulus in N/m^2
+rho = 2650 // Density in Kg/m^3
+t = 0.006 // Thickness of quartz crystal in m
+printf("Example 6.4\n")
+v = sqrt(E/rho)// Calculation of velocity
+lambda = 2*t // Calculation of fundamental wavelength
+nu = v/lambda // Calculation of fundamental frequency
+printf("Fundamental frequency is %e Hz.\n\n\n",nu)
+
diff --git a/3014/CH6/EX6.4/Ex6_4.txt b/3014/CH6/EX6.4/Ex6_4.txt new file mode 100755 index 000000000..1b30d8157 --- /dev/null +++ b/3014/CH6/EX6.4/Ex6_4.txt @@ -0,0 +1,3 @@ +
+Example 6.4
+Fundamental frequency is 4.549979e+05 Hz.
diff --git a/3014/CH6/EX6.5/Ex6_5.sce b/3014/CH6/EX6.5/Ex6_5.sce new file mode 100755 index 000000000..3155fcd7d --- /dev/null +++ b/3014/CH6/EX6.5/Ex6_5.sce @@ -0,0 +1,9 @@ +
+clc
+//Given that
+L = 1 // Inductance in Hanery
+nu = 2e6 // Frequency in Hz
+printf("Example 6.5\n")
+C= 1/(4*((%pi)^2)*nu^2*L) // Calculation of capacitance
+printf("Capacitance is %e microfarad.\n\n\n",C*1e6)
+// Answer in book is 0.00634 micro Farad
diff --git a/3014/CH6/EX6.5/Ex6_5.txt b/3014/CH6/EX6.5/Ex6_5.txt new file mode 100755 index 000000000..76984269a --- /dev/null +++ b/3014/CH6/EX6.5/Ex6_5.txt @@ -0,0 +1,3 @@ +
+Example 6.5
+Capacitance is 6.332574e-09microfarad.
diff --git a/3014/CH7/EX7.1/Ex7_1.sce b/3014/CH7/EX7.1/Ex7_1.sce new file mode 100755 index 000000000..0209dc224 --- /dev/null +++ b/3014/CH7/EX7.1/Ex7_1.sce @@ -0,0 +1,13 @@ +
+clc
+//Given that
+p = 1000 // power in watt
+d = 2 // Distance from lamp in m
+epsilon_0 = 8.854e-12 // Permittivity of free space
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+printf("Example 7.1")
+s = p/(4*%pi*d^2)// Calculation of pointing vector
+E_H_ratio = sqrt(mu_0/epsilon_0) // Calculation of ratio of Electric field and magnetic field
+E= sqrt(E_H_ratio*s) // Calculation of Electric field
+printf("\n Average value of electric field at distance %d m is %f Volt/m \n\n\n",d,E)
+// Answer in book is 48.87 volt/m which is due to wrong calculation at intermediate steps
diff --git a/3014/CH7/EX7.1/Ex7_1.txt b/3014/CH7/EX7.1/Ex7_1.txt new file mode 100755 index 000000000..a7be705dc --- /dev/null +++ b/3014/CH7/EX7.1/Ex7_1.txt @@ -0,0 +1,3 @@ + Example 7.1
+ Average value of electric field at distance 2 m is 86.573038 Volt/m
+
diff --git a/3014/CH7/EX7.10/Ex7_10.sce b/3014/CH7/EX7.10/Ex7_10.sce new file mode 100755 index 000000000..d2bbf802e --- /dev/null +++ b/3014/CH7/EX7.10/Ex7_10.sce @@ -0,0 +1,15 @@ +
+clc
+//Given that
+p = 500 // power in watt
+d = 1 // Distance from lamp in m
+epsilon_0 = 8.854e-12 // Permittivity of free space
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+printf("Example 7.10")
+s = p/(4*%pi*d^2)// Calculation of pointing vector
+E_H_ratio = sqrt(mu_0/epsilon_0) // Calculation of ratio of Electric field and magnetic field
+H = s/E_H_ratio // Calculation of Electric field
+h = ceil(H*100)/100 // rounding off for 2 decimal places
+E= p/(4*%pi*h) // Calculation of Electric field
+printf("\n Average value of electric field at distance %d m is %f Volt/m ",d,E)
+printf("\n Average value of magnetic field at distance %d m is %f Amp-turn/m \n\n\n",d,h)
diff --git a/3014/CH7/EX7.10/Ex7_10.txt b/3014/CH7/EX7.10/Ex7_10.txt new file mode 100755 index 000000000..762603d79 --- /dev/null +++ b/3014/CH7/EX7.10/Ex7_10.txt @@ -0,0 +1,4 @@ +
+Example 7.10
+ Average value of electric field at distance 1 m is 122.432765 Volt/m
+ Average value of magnetic field at distance 1 m is 0.324984 Amp-turn/m
diff --git a/3014/CH7/EX7.11/Ex7_11.sce b/3014/CH7/EX7.11/Ex7_11.sce new file mode 100755 index 000000000..c262d1919 --- /dev/null +++ b/3014/CH7/EX7.11/Ex7_11.sce @@ -0,0 +1,13 @@ +
+clc
+//Given that
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+mu = mu_0 //Permeability of silver
+sigma = 3.5e7 // conductivity in simens /m
+delta = 0.03 // Skin depth penetration in mm
+
+printf("Example 7.11")
+
+f = 2/((delta*1e-3)^2*sigma*mu*2*%pi) // Calculation of skin depth penetration
+printf("\n Required frequency is %d MHz.",f/1e6)
+printf("\n The incident electromagnetic wave is the radio part of spectrum")
diff --git a/3014/CH7/EX7.11/Ex7_11.txt b/3014/CH7/EX7.11/Ex7_11.txt new file mode 100755 index 000000000..5050c1461 --- /dev/null +++ b/3014/CH7/EX7.11/Ex7_11.txt @@ -0,0 +1,4 @@ +
+
+Example 7.11
+ Required frequency is 8 MHz.
diff --git a/3014/CH7/EX7.12/Ex7_12.sce b/3014/CH7/EX7.12/Ex7_12.sce new file mode 100755 index 000000000..af83bfbb0 --- /dev/null +++ b/3014/CH7/EX7.12/Ex7_12.sce @@ -0,0 +1,14 @@ +
+clc
+//Given that
+p = 3.8e26 // power radiated by moon in watt
+d_sun = 1.44e11 // Distance between sun and earth in meter
+d_moon = 3e8 //Distance between moon and earth in meter
+epsilon_0 = 8.854e-12 // Permittivity of free space
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+printf("Example 7.12")
+s = p/(4*%pi*d_sun^2)// Calculation of solar energy received during solar eclipse in watt /m^2
+S = s*60/(4.2*1e4) // Unit conversion
+
+printf("\n Solar energy received during solar eclipse is %f Cal per min per m^2 \n\n\n",S)
+// Ansewr in book is 2.1 cal per min per m^2
diff --git a/3014/CH7/EX7.12/Ex7_12.txt b/3014/CH7/EX7.12/Ex7_12.txt new file mode 100755 index 000000000..37fd989bc --- /dev/null +++ b/3014/CH7/EX7.12/Ex7_12.txt @@ -0,0 +1,3 @@ +
+Example 7.12
+ Solar energy received during solar eclipse is 1.999130 Cal per min per m^2
diff --git a/3014/CH7/EX7.13/Ex7_13.sce b/3014/CH7/EX7.13/Ex7_13.sce new file mode 100755 index 000000000..89f689d15 --- /dev/null +++ b/3014/CH7/EX7.13/Ex7_13.sce @@ -0,0 +1,17 @@ +
+
+clc
+//Given that
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+mu = mu_0 //Permeability of silver
+sigma = 3.5e7 // conductivity in simens /m
+lambda = 6328 // Wavelength in angstrom
+c = 3e8// Speed of light in m/sec
+
+printf("Example 7.13")
+f = c/(lambda*1e-10)
+omega = 2*%pi/f // Calculation of time period
+f = c/(lambda*1e-10) // Calculation of frequency in Hz
+delta = sqrt(1/(%pi*f*sigma*mu)) // Calculation of skin depth penetration
+printf("\n Skin depth penetration is %f nm. \n\n\n",delta*1e9)
+// Answer in book is 3.9 mm, unit used in book is wrong
diff --git a/3014/CH7/EX7.13/Ex7_13.txt b/3014/CH7/EX7.13/Ex7_13.txt new file mode 100755 index 000000000..7d9f55b4c --- /dev/null +++ b/3014/CH7/EX7.13/Ex7_13.txt @@ -0,0 +1,3 @@ +
+Example 7.13
+ Skin depth penetration is 3.907138 nm.
diff --git a/3014/CH7/EX7.2/Ex7_2.sce b/3014/CH7/EX7.2/Ex7_2.sce new file mode 100755 index 000000000..41b1cd005 --- /dev/null +++ b/3014/CH7/EX7.2/Ex7_2.sce @@ -0,0 +1,15 @@ +
+clc
+//Given that
+p = 2 // power in cal/min/cm^2
+
+epsilon_0 = 8.854e-12 // Permittivity of free space
+mu_0 = 4*%pi*1e-7 // permeability of free space
+printf("Example 7.2")
+s = p*4.2e4/60 // Calculation of pointing vector
+E_H_ratio = sqrt(mu_0/epsilon_0) // Calculation of ratio of Electric field and magnetic field
+E= sqrt(E_H_ratio*s) // Calculation of Electric field
+H = s/E // Calculation of Electric field
+
+printf("\n Average value of electric field is %f Volt/m ",E*sqrt(2))
+printf(" \nAverage value of magnetic field is %f Amp turn/m \n\n\n",H*sqrt(2))
diff --git a/3014/CH7/EX7.2/Ex7_2.txt b/3014/CH7/EX7.2/Ex7_2.txt new file mode 100755 index 000000000..2959e8640 --- /dev/null +++ b/3014/CH7/EX7.2/Ex7_2.txt @@ -0,0 +1,4 @@ +
+Example 7.2
+ Average value of electric field is 1027.061861 Volt/m.
+Average value of magnetic field is 2.726223 Amp turn/m.
\ No newline at end of file diff --git a/3014/CH7/EX7.3/Ex7_3.sce b/3014/CH7/EX7.3/Ex7_3.sce new file mode 100755 index 000000000..167f73ce2 --- /dev/null +++ b/3014/CH7/EX7.3/Ex7_3.sce @@ -0,0 +1,13 @@ +
+clc
+//Given that
+mu_0 = 4*%pi*1e-7 // permeability of free space
+mu = mu_0 //permeability of silver
+sigma = 3e7 // conductivity in mhos/m
+f = 1e8 // frequency in Hz
+printf("Example 7.3")
+omega = 2*%pi/f // Calculation of time period
+delta = sqrt(2/(omega*sigma*mu)) // Calculation of skin depth penetration
+Delta = floor (delta/100)*100 // Rounding off
+printf("\n Skin depth penetration is %e cm. \n\n\n",Delta*1e-6)
+
diff --git a/3014/CH7/EX7.3/Ex7_3.txt b/3014/CH7/EX7.3/Ex7_3.txt new file mode 100755 index 000000000..6cc9f7539 --- /dev/null +++ b/3014/CH7/EX7.3/Ex7_3.txt @@ -0,0 +1,3 @@ +
+Example 7.3
+ Skin depth penetration is 9.188815e-04 cm.
diff --git a/3014/CH7/EX7.5/Ex7_5.sce b/3014/CH7/EX7.5/Ex7_5.sce new file mode 100755 index 000000000..5d50688ed --- /dev/null +++ b/3014/CH7/EX7.5/Ex7_5.sce @@ -0,0 +1,19 @@ +
+clc
+//Given that
+k = 80 // relative Dielectric constant of sea water
+epsilon_0 = 1/9e9 // Permittivity of free space
+epsilon = 80*epsilon_0 // Permittivity of free space
+sigma = 4.3 // conductivity in mho/m
+delta = 10 // penetration depth in cm
+mu_0 = 4*%pi*1e-7 // permeability f free space
+F = 1e8 // Given frequency in Hz
+printf("Example 7.5")
+f = (1/(%pi*mu_0*sigma))/(delta*1e-2)^2 // Calculation of frequency
+f1= ceil(f/1e8)*1e8 // Rounding off
+printf("\nFrequency required for penetration of depth %d cm is %e Hz",delta,f1)
+omega = 2*%pi*F
+x = 2*sigma/(epsilon*omega)
+if x>1 then
+ printf("\n Sea water is good conductor at frequency lesser than 1e8 Hz\n\n ")
+end
diff --git a/3014/CH7/EX7.5/Ex7_5.txt b/3014/CH7/EX7.5/Ex7_5.txt new file mode 100755 index 000000000..136f4d936 --- /dev/null +++ b/3014/CH7/EX7.5/Ex7_5.txt @@ -0,0 +1,4 @@ +
+Example 7.5
+Frequency required for penetration of depth 10 cm is 1.716212e+03 Hz.
+ Sea water is good conductor at frequency lesser than 1e8 Hz.
diff --git a/3014/CH7/EX7.7/Ex7_7.sce b/3014/CH7/EX7.7/Ex7_7.sce new file mode 100755 index 000000000..5ff36e6e0 --- /dev/null +++ b/3014/CH7/EX7.7/Ex7_7.sce @@ -0,0 +1,18 @@ +
+clc
+//Given that
+k = 12 // relative Dielectric constant of sea water
+epsilon_0 = 1/9e9 // Permittivity of free space
+sigma = 2 // conductivity in mho/cm
+mu_0 = 4*%pi*1e-7 // permeability f free space
+f= 1e9 // Given frequency in Hz
+F = 1e6 // Given frequency in Hz
+printf("\nExample 7.7")
+delta = sqrt(2/(2*%pi*F*mu_0*sigma*100)) // Calculation of frequency
+printf("\n For %eHz frequency, Penetration depth is %f cm",F,delta*100)
+omega = 2*%pi*f
+x = 2*sigma*100/(k*epsilon_0*omega)
+if x>1 then
+ printf("\n Silicon is good conductor at frequency lesser than 1e9 Hz \n\n\n")
+end
+// Answer in book is 3.6 cm
diff --git a/3014/CH7/EX7.7/Ex7_7.txt b/3014/CH7/EX7.7/Ex7_7.txt new file mode 100755 index 000000000..854cd5f3b --- /dev/null +++ b/3014/CH7/EX7.7/Ex7_7.txt @@ -0,0 +1,4 @@ +
+Example 7.7
+ For frequency 1.000000e+06 Penetration depth is 3.558813 cm
+ Silicon is good conductor at frequency lesser than 1e9 Hz
diff --git a/3014/CH7/EX7.8/Ex7_8.sce b/3014/CH7/EX7.8/Ex7_8.sce new file mode 100755 index 000000000..198827463 --- /dev/null +++ b/3014/CH7/EX7.8/Ex7_8.sce @@ -0,0 +1,13 @@ +clc
+//Given that
+mu_0 = 4*%pi*1e-7 // permeability of free space
+mu = mu_0 //permeability of silver
+sigma = 5.8e7 // conductivity in simens /m
+delta = 0.1 // Skin depth penetration in mm
+
+printf("Example 7.8")
+f = 2/((delta*1e-3)^2*sigma*mu*2*%pi) // Calculation of skin depth penetration
+printf("\n Required frequency is %.2e Hz",f)
+printf("\n The incident electromagnetic wave is the radio part of spectrum.")
+// Answer in book is 3.36e5 Hz. Difference is due to approximation at intermediate stages
+
diff --git a/3014/CH7/EX7.8/Ex7_8.txt b/3014/CH7/EX7.8/Ex7_8.txt new file mode 100755 index 000000000..175c8f6f8 --- /dev/null +++ b/3014/CH7/EX7.8/Ex7_8.txt @@ -0,0 +1,3 @@ +
+Example 7.8
+ Required frequency is 4.367292e+05 Hz
diff --git a/3014/CH7/EX7.9/Ex7_9.sce b/3014/CH7/EX7.9/Ex7_9.sce new file mode 100755 index 000000000..04c29812e --- /dev/null +++ b/3014/CH7/EX7.9/Ex7_9.sce @@ -0,0 +1,12 @@ +
+
+clc
+//Given that
+mu_0 = 4*%pi*1e-7 // Permeability of free space
+mu = mu_0 //Permeability of silver
+sigma = 3e7 // conductivity in mhos/m
+f = 1e10 // frequency in Hz
+printf("Example 7.9")
+delta = sqrt(1/(%pi*sigma*f*mu)) // Calculation of skin depth penetration
+printf("\n Skin depth penetration is %f micrometre. \n\n\n",delta*1e6)
+// Answer in book is 0.93 micrometer
diff --git a/3014/CH7/EX7.9/Ex7_9.txt b/3014/CH7/EX7.9/Ex7_9.txt new file mode 100755 index 000000000..7eec8e2be --- /dev/null +++ b/3014/CH7/EX7.9/Ex7_9.txt @@ -0,0 +1,3 @@ +
+Example 7.9
+ Skin depth penetration is 0.918881 micrometre.
diff --git a/3014/CH8/EX8.1/Ex8_1.sce b/3014/CH8/EX8.1/Ex8_1.sce new file mode 100755 index 000000000..2d0d7ec5a --- /dev/null +++ b/3014/CH8/EX8.1/Ex8_1.sce @@ -0,0 +1,10 @@ + clc
+// Given that
+H_c_0= 0.0306// Critical Field in tesla
+T_c = 3.7 // Critical temperature in kelvin
+T = 2 // Temperature in kelvin
+printf("Example 8.1\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+H_c = H_c_0*(1-(T/T_c)^2) // Calculation of critical field
+
+printf("Magnetic Field at %d K is %f tesla.\n\n\n",T,H_c)
diff --git a/3014/CH8/EX8.1/Ex8_1.txt b/3014/CH8/EX8.1/Ex8_1.txt new file mode 100755 index 000000000..ef3f5de33 --- /dev/null +++ b/3014/CH8/EX8.1/Ex8_1.txt @@ -0,0 +1,3 @@ + Example 8.1
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Magnetic Field at 2 K is 0.021659 tesla.
\ No newline at end of file diff --git a/3014/CH8/EX8.10/Ex8_10.sce b/3014/CH8/EX8.10/Ex8_10.sce new file mode 100755 index 000000000..1fdbe6fd0 --- /dev/null +++ b/3014/CH8/EX8.10/Ex8_10.sce @@ -0,0 +1,13 @@ +
+ clc
+// Given that
+T_1 = 3.5 // Temperature in kelvin
+T_c = 4.153 // Critical temp in kelvin
+lambda_t = 750 // Penetration depth at T_1 in angstrom
+printf("Example 8.10\n")
+printf("Standard formula used \nlambda_0 = lambda_t*sqrt(1-(T_1/T_c)^4) \n")
+
+lambda_0 = lambda_t*sqrt(1-(T_1/T_c)^4) // Calculation of penetration depth at 3.5K
+printf("\n Penetration depth at 0 K is %f angstrom.\n\n\n",lambda_0)
+
+
diff --git a/3014/CH8/EX8.10/Ex8_10.txt b/3014/CH8/EX8.10/Ex8_10.txt new file mode 100755 index 000000000..43033cec0 --- /dev/null +++ b/3014/CH8/EX8.10/Ex8_10.txt @@ -0,0 +1,5 @@ + Example 8.10
+Standard formula used
+lambda_0 = lambda_t*sqrt(1-(T_1/T_c)^4)
+
+ Penetration depth at 0 K is 527.960928 angstrom.
\ No newline at end of file diff --git a/3014/CH8/EX8.2/Ex8_2.sce b/3014/CH8/EX8.2/Ex8_2.sce new file mode 100755 index 000000000..277e8fa44 --- /dev/null +++ b/3014/CH8/EX8.2/Ex8_2.sce @@ -0,0 +1,10 @@ +
+ clc
+// Given that
+H_c= 3.3e4 // // Magnetic field in A/m
+T_c = 7.2 // Critical temperature in kelvin
+T = 5 // Temperature in kelvin
+printf("Example 8.2\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+H_c_0 = H_c*(1-(T/T_c)^2)^(-1) // Calculation of critical field
+printf("Magnetic Field at %d K is %e A/m\n\n\n",T,H_c_0)
diff --git a/3014/CH8/EX8.2/Ex8_2.txt b/3014/CH8/EX8.2/Ex8_2.txt new file mode 100755 index 000000000..58fa2ce1a --- /dev/null +++ b/3014/CH8/EX8.2/Ex8_2.txt @@ -0,0 +1,4 @@ +
+Example 8.2
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Magnetic Field at 5 K is 6.373770e+04 A/m
diff --git a/3014/CH8/EX8.3/Ex8_3.sce b/3014/CH8/EX8.3/Ex8_3.sce new file mode 100755 index 000000000..37f10beb5 --- /dev/null +++ b/3014/CH8/EX8.3/Ex8_3.sce @@ -0,0 +1,12 @@ +
+ clc
+// Given that
+H_c_0= 1 // Let
+H_c= 0.1 * H_c_0 // Magnetic field in A/m
+T_c = 7.2 // Critical temperature in kelvin
+
+printf("Example 8.3\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+T = T_c*sqrt(1- (H_c/H_c_0)) // Calculation of Temperature
+
+printf("Required temperature is %f K.\n\n\n",T)
diff --git a/3014/CH8/EX8.3/Ex8_3.txt b/3014/CH8/EX8.3/Ex8_3.txt new file mode 100755 index 000000000..e8452ca5f --- /dev/null +++ b/3014/CH8/EX8.3/Ex8_3.txt @@ -0,0 +1,6 @@ +
+
+Example 8.3
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Required temperature is 6.830520 K.
+
diff --git a/3014/CH8/EX8.4/Ex8_4.sce b/3014/CH8/EX8.4/Ex8_4.sce new file mode 100755 index 000000000..1d3da19bf --- /dev/null +++ b/3014/CH8/EX8.4/Ex8_4.sce @@ -0,0 +1,13 @@ +
+
+ clc
+// Given that
+H_c_0= 0.0803// Critical Field in tesla
+T_c = 7.2 // Critical temperature in kelvin
+T = 4.2 // Temperature in kelvin
+printf("Example 8.4\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+H_c = H_c_0*(1-(T/T_c)^2) // Calculation of critical field
+
+printf("Magnetic Field at %d K is %f tesla.\n\n\n",T,H_c)
+// Answer in book is 0.0548 tesla
diff --git a/3014/CH8/EX8.4/Ex8_4.txt b/3014/CH8/EX8.4/Ex8_4.txt new file mode 100755 index 000000000..6095e71e6 --- /dev/null +++ b/3014/CH8/EX8.4/Ex8_4.txt @@ -0,0 +1,6 @@ +
+
+
+Example 8.4
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Magnetic Field at 4 K is 0.052976 tesla.
\ No newline at end of file diff --git a/3014/CH8/EX8.5/Ex8_5.sce b/3014/CH8/EX8.5/Ex8_5.sce new file mode 100755 index 000000000..3d818e98a --- /dev/null +++ b/3014/CH8/EX8.5/Ex8_5.sce @@ -0,0 +1,13 @@ +
+
+ clc
+// Given that
+H_c_0= 1.5e5// Critical field in A/m
+H_c= 1.05e5 // Magnetic field in A/m
+T_c = 9.2 // Critical temperature in kelvin
+
+printf ("Example 8.5\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+T = T_c*sqrt(1- (H_c/H_c_0)) // Calculation of Temperature
+
+printf("Required temperature is %f K.\n\n\n",T)
diff --git a/3014/CH8/EX8.5/Ex8_5.txt b/3014/CH8/EX8.5/Ex8_5.txt new file mode 100755 index 000000000..9a39daf9f --- /dev/null +++ b/3014/CH8/EX8.5/Ex8_5.txt @@ -0,0 +1,4 @@ +
+Example 8.5
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Required temperature is 5.039048 K.
diff --git a/3014/CH8/EX8.6/Ex8_6.sce b/3014/CH8/EX8.6/Ex8_6.sce new file mode 100755 index 000000000..318f8e158 --- /dev/null +++ b/3014/CH8/EX8.6/Ex8_6.sce @@ -0,0 +1,13 @@ +
+
+ clc
+// Given that
+H_c_0= 2e5// Critical field in A/m
+H_c= 1e5 // Magnetic field in A/m
+T_c = 8 // Critical temperature in kelvin
+
+printf("Example 8.6\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+T = T_c/sqrt(1- (H_c/H_c_0)) // Calculation of Temperature
+
+printf("Required temperature is %f K.\n\n\n",T)
diff --git a/3014/CH8/EX8.6/Ex8_6.txt b/3014/CH8/EX8.6/Ex8_6.txt new file mode 100755 index 000000000..bb539ca1e --- /dev/null +++ b/3014/CH8/EX8.6/Ex8_6.txt @@ -0,0 +1,4 @@ +
+Example 8.6
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Required temperature is 11.313708 K.
diff --git a/3014/CH8/EX8.7/Ex8_7.sce b/3014/CH8/EX8.7/Ex8_7.sce new file mode 100755 index 000000000..7160962a5 --- /dev/null +++ b/3014/CH8/EX8.7/Ex8_7.sce @@ -0,0 +1,12 @@ +
+ clc
+// Given that
+H_c_0= 8e5// Critical field in A/m
+H_c= 4e4 // Magnetic field in A/m
+T_c = 7.26 // Critical temperature in kelvin
+
+printf("Example 8.7\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+T = T_c*sqrt(1- (H_c/H_c_0)) // Calculation of Temperature
+
+printf("Required temperature is %f K.\n\n\n",T)
diff --git a/3014/CH8/EX8.7/Ex8_7.txt b/3014/CH8/EX8.7/Ex8_7.txt new file mode 100755 index 000000000..9b6086059 --- /dev/null +++ b/3014/CH8/EX8.7/Ex8_7.txt @@ -0,0 +1,4 @@ +
+Example 8.7
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Required temperature is 7.076173 K.
diff --git a/3014/CH8/EX8.8/Ex8_8.sce b/3014/CH8/EX8.8/Ex8_8.sce new file mode 100755 index 000000000..321eeaa14 --- /dev/null +++ b/3014/CH8/EX8.8/Ex8_8.sce @@ -0,0 +1,20 @@ +
+ clc
+// Given that
+T1 = 14 // Temp in K
+T2 = 13 // Temp in K
+T = 4.2 // Temp in K
+Hc_T1 = 0.176 // Critical field at Temp T1
+Hc_T2 = 0.528 // Critical field at Temp T2
+
+printf("Example 8.8\n")
+printf("Standard formula used \tH_c = H_c_0*(1-(T/T_c)^2) \n")
+T_c = sqrt((T1^2*(Hc_T2/Hc_T1)- T2^2) /(Hc_T2/Hc_T1 - 1)) // Calculation of transition temperature
+t_c = ceil(T_c*10)/10 // Rounding off two two decimal places
+Hc_0 = Hc_T1/(1-(T1/t_c)^2) // Calculation of critical field
+Hc_T = Hc_0*(1-(T/t_c)^2) // Calculation of critical field
+
+printf("\n Transition temperature is %f K.\n",t_c)
+printf("Critical field at %f K is %fT.\n",T,Hc_0)
+printf("Critical field at 0 K is %fT.\n\n\n",Hc_T)
+// Answer in book is 2.588 T for 0 K and 2.37 for 4.2 K
diff --git a/3014/CH8/EX8.8/Ex8_8.txt b/3014/CH8/EX8.8/Ex8_8.txt new file mode 100755 index 000000000..c32d57dc5 --- /dev/null +++ b/3014/CH8/EX8.8/Ex8_8.txt @@ -0,0 +1,8 @@ +
+
+
+Example 8.8
+Standard formula used H_c = H_c_0*(1-(T/T_c)^2)
+Transition temperature is 14.474115 K.
+Critical field at 4.200000 K is 2.731259T.
+Critical field at 0 K is 2.982377T.
diff --git a/3014/CH8/EX8.9/Ex8_9.sce b/3014/CH8/EX8.9/Ex8_9.sce new file mode 100755 index 000000000..4415a0c05 --- /dev/null +++ b/3014/CH8/EX8.9/Ex8_9.sce @@ -0,0 +1,18 @@ +
+ clc
+// Given that
+m_0 = 9.1e-31 // Mass of electron in kg
+mu_0 = 1.256e-6// SI
+e = 1.6e-19 // Charge on electron in coulomb
+eta_s = 1e28 // superelectron density in no. per cube
+T_1 = 0 // First temp in kelvin
+T_2 = 1 // Second temp in kelvin
+T_c = 3 // Critical temp in kelvin
+
+printf("Example 8.9\n")
+printf("Standard formula used \tlambda_0 = sqrt(m_0/(mu_0*eta_s*e^2))\n")
+lambda_0 = sqrt(m_0/(mu_0*eta_s*e^2))// Calculation of penetration depth at 0K
+lambda_t = lambda_0/sqrt(1-(T_2/T_c)^4) // Calculation of penetration depth at 2K
+
+printf("Penetration depth at %d K is %d angestrom.",T_1,lambda_0*1e10)
+printf("\nPenetration depth at %d K is %f angestrom.\n\n\n",T_2,lambda_t*1e10)
diff --git a/3014/CH8/EX8.9/Ex8_9.txt b/3014/CH8/EX8.9/Ex8_9.txt new file mode 100755 index 000000000..99f099fff --- /dev/null +++ b/3014/CH8/EX8.9/Ex8_9.txt @@ -0,0 +1,4 @@ + Example 8.9
+Standard formula used lambda_0 = sqrt(m_0/(mu_0*eta_s*e^2))
+Penetration depth at 0 K is 531 angestrom.
+Penetration depth at 1 K is 535.307601 angestrom.
\ No newline at end of file |