initial commit / add all books

8 files changed, 110 insertions, 0 deletions

diff --git a/1970/CH5/EX5.1/Ch05Exa1.sce b/1970/CH5/EX5.1/Ch05Exa1.sce
new file mode 100755
index 000000000..3be2fadc6
--- /dev/null
+++ b/1970/CH5/EX5.1/Ch05Exa1.sce
@@ -0,0 +1,16 @@
+// Scilab code Exa5.1 : : Page 203 (2011)

+clc; clear;

+E_a = 8.766;    // Energy of the alpha particle, MeV

+A = 212;          // Atomic mass of Po-212, amu

+M_a = 4;         // Atomic mass of alpha particle, amu

+e = 1.6e-019;        // Charge of an electron, coulomb

+Z = 82;                // Atomic number of Po-212

+R_0 = 1.4e-015;        // Distance of closest approach,metre

+K = 8.99e+09;            // Coulomb constant

+E = E_a*A/(A-M_a);    // Disintegration energy, mega electron volts

+B_H = 2*Z*e^2*K/(R_0*A^(1/3)*1.6*10^-13);    // Barrier height for an alpha particle within the nucleus, MeV

+printf("\nDisintegration energy : %5.3f MeV \nBarrier height for alpha-particle: %5.2f MeV", E,B_H);

+

+// Result

+// Disintegration energy : 8.935 MeV 

+// Barrier height for alpha-particle: 28.26 MeV 
\ No newline at end of file
diff --git a/1970/CH5/EX5.10/Ch05Exa10.sce b/1970/CH5/EX5.10/Ch05Exa10.sce
new file mode 100755
index 000000000..2238f4624
--- /dev/null
+++ b/1970/CH5/EX5.10/Ch05Exa10.sce
@@ -0,0 +1,14 @@
+// Scilab code Exa5.10 : : Page 206 (2011)

+clc; clear;

+h_kt = 1.05457e-34;        // Reduced Planck's constant, joule sec

+e = 1.60218e-19;        // Charge of an electron, coulomb

+l = 2;                // Orbital angular momentum

+eps_0 = 8.5542e-12;        // Absolute permittivity of free space, coulomb square per newton per metre square

+Z_D = 90;            // Atomic number of daughter nucleus

+m = 6.644e-27;        // Mass of alpha particle, Kg

+R = 8.627e-15;        // Radius of daughter nucleus, metre

+T1_by_T0 = exp(2*l*(l+1)*h_kt/e*sqrt(%pi*eps_0/(Z_D*m*R)));    // Hindrance factor

+printf("\nThe hindrance factor for alpha particle = %5.3f" ,T1_by_T0);

+

+// Result

+// The hindrance factor for alpha particle = 1.768 
\ No newline at end of file
diff --git a/1970/CH5/EX5.2/Ch05Exa2.sce b/1970/CH5/EX5.2/Ch05Exa2.sce
new file mode 100755
index 000000000..ab17c351e
--- /dev/null
+++ b/1970/CH5/EX5.2/Ch05Exa2.sce
@@ -0,0 +1,18 @@
+// Scilab code Exa5.2 : : Page 203 (2011)

+// We have to make calculation for alpha particle and for  proton

+clc; clear;

+E_a = 8.766;     // Energy of the alpha particle, mega electron volts

+A_Bi = 209;      // Atomic mass of Bi-209, atomic mass unit

+A_a = 4;         // Atomic mass of alpha particle, atomic mass unit

+A_p = 1;         // Atomic mass of proton, atomic mass unit

+e = 1.6e-019;    // Charge of an electron, coulomb

+Z = 83;          // Atomic number of bismuth

+R_0 = 1.4e-015;  // Distance of closest approach,metre

+K = 8.99e+09;    // Coulomb constant

+B_H_a = 2*Z*e^2*K/(R_0*1.6e-013*(A_Bi^(1/3)+A_a^(1/3)));    // Barrier height for an alpha particle, mega electron volts

+B_H_p = 1*Z*e^2*K/(R_0*1.6e-013*(A_Bi^(1/3)+A_p^(1/3)));    // Barrier height for proton, mega electron volts

+printf("\nBarrier height for the alpha particle = %5.2f MeV \nBarrier height for the proton = %5.2f MeV", B_H_a,B_H_p);

+

+// Result

+// Barrier height for the alpha particle = 22.67 MeV 

+// Barrier height for the proton = 12.30 MeV 
\ No newline at end of file
diff --git a/1970/CH5/EX5.3/Ch05Exa3.sce b/1970/CH5/EX5.3/Ch05Exa3.sce
new file mode 100755
index 000000000..f5c1d36cc
--- /dev/null
+++ b/1970/CH5/EX5.3/Ch05Exa3.sce
@@ -0,0 +1,14 @@
+// Scilab code Exa5.3 : : Page 203 (2011)

+// We have also calculate the value of magnetic field in a particular orbit. 

+clc; clear;

+C = 3e+08;                // Velocity of light, m/S

+M_0 = 6.644e-027*(C)^2/(1.60218e-013);        // Rest mass of alpha particle, MeV

+T = 5.998;                // Kinetic energy of alpha particle emitted by Po-218

+q = 2*1.60218e-019;        // Charge of alpha particle, C

+V = sqrt(C^2*T*(T+2*M_0)/(T+M_0)^2);            // Velocity of alpha particle,metre per sec

+B_r = V*M_0*(1.60218e-013)/(C^2*q*sqrt(1-V^2/C^2));                // magnetic field in a particular orbit, Web per mtere

+printf("\nThe velocity of alpha particle : %5.3e m/s\nThe magnetic field in a particular orbit : %6.4f Wb/m", V , B_r);

+

+// Result

+// The velocity of alpha particle : 1.699e+007 m/s

+// The magnetic field in a particular orbit : 0.3528 Wb/m 
\ No newline at end of file
diff --git a/1970/CH5/EX5.4/Ch05Exa4.sce b/1970/CH5/EX5.4/Ch05Exa4.sce
new file mode 100755
index 000000000..da53b07b4
--- /dev/null
+++ b/1970/CH5/EX5.4/Ch05Exa4.sce
@@ -0,0 +1,12 @@
+// Scilab code Exa5.4: : Page 204 (2011)

+clc; clear;

+a = 10^-14;        // Width of the  potential barrier, m

+E = 5*1.60218e-013;        // Energy of the alpha particle, joule

+V = 10*1.60218e-013;        // Potential height, joule

+M_0 = 6.644e-027;        // Rest mass of the alpha particle, joule

+h_red = 1.05457e-034;        // Reduced value of Planck's constant,joule sec 

+T = 4*exp(-2*a*sqrt(2*M_0*(V-E)/h_red^2));   // Probability of leakage through through potential barrier

+printf("\nThe probability of leakage of alpha-particle through potential barrier = %5.3e ",T);

+

+// Result

+// The probability of leakage of alpha-particle through potential barrier = 1.271e-008  
\ No newline at end of file
diff --git a/1970/CH5/EX5.6/Ch05Exa6.sce b/1970/CH5/EX5.6/Ch05Exa6.sce
new file mode 100755
index 000000000..75fc456b4
--- /dev/null
+++ b/1970/CH5/EX5.6/Ch05Exa6.sce
@@ -0,0 +1,12 @@
+// Scilab code Exa5.6: : Page 204 (2011)

+clc; clear; 

+Z_D = 82;    // Atomic number of Po

+E_Po210 = 5.3;    // Alpha-source for Po210, MeV

+E_Po214 = 7.7;    // Alpha-source for Po214, MeV

+log_lambda_Po210 = -1*1.72*Z_D*E_Po210^(-1/2);    

+log_lambda_Po214 = -1*1.72*Z_D*E_Po214^(-1/2);    

+delta_OM_t = log_lambda_Po214 - log_lambda_Po210;    // Difference in order of magnitude of life times of Po214 and Po210

+printf("\nThe disintegration constant increases by a factor of some 10^%2d", delta_OM_t);

+

+// Result

+// The disintegration constant increases by a factor of some 10^10  
\ No newline at end of file
diff --git a/1970/CH5/EX5.8/Ch05Exa8.sce b/1970/CH5/EX5.8/Ch05Exa8.sce
new file mode 100755
index 000000000..7ccf0e867
--- /dev/null
+++ b/1970/CH5/EX5.8/Ch05Exa8.sce
@@ -0,0 +1,11 @@
+// Scilab code Exa5.8:  : Page 205 (2011)

+clc; clear;

+N = 120.1*6.023e+023/239;    // Number of Pu nuclei

+P_rel = 0.231;        // Power released, watt

+E_rel = 5.323*1.6026e-13;        // Energy released, joule

+decay_rate = P_rel/E_rel;        // Decay rate of Pu239, per hour

+t_half = N*log(2)/(decay_rate*365*86400);    // Half life of Po239, sec

+printf("\nThe half life of Pu = %4.2e yr", t_half);

+

+// Result

+// The half life of Pu = 2.46e+004 yr 
\ No newline at end of file
diff --git a/1970/CH5/EX5.9/Ch05Exa9.sce b/1970/CH5/EX5.9/Ch05Exa9.sce
new file mode 100755
index 000000000..f640baa38
--- /dev/null
+++ b/1970/CH5/EX5.9/Ch05Exa9.sce
@@ -0,0 +1,13 @@
+// Scilab code Exa5.9 : : Page 205(2011)

+clc; clear;

+a_v = 14;        // Volume energy constant, MeV

+a_s = 13;        // Surface energy constant, MeV

+a_c = 0.60;      // Coulomb energy constant, MeV

+a_a = 19;        // Asymmetric energy constant, MeV

+A = 202;        // Mass number

+Z = 82;         //  Atomic number 

+dE_by_dN = -8/9*a_s/A^(4/3)-4/3*a_c*Z/A^(4/3)*(1-4*Z/(3*A))-16*a_a*Z/A^2*(1-2*Z/A);        // Slope, mega electron volts per nucleon

+printf("\nThe slope of alpha decay energy versus atomic number = %7.5f MeV/nucleon", dE_by_dN);

+

+// Result

+// The slope of alpha decay energy versus atomic number = -0.15007 MeV/nucleon 
\ No newline at end of file