From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 1970/CH5/EX5.1/Ch05Exa1.sce | 16 ++++++++++++++++ 1970/CH5/EX5.10/Ch05Exa10.sce | 14 ++++++++++++++ 1970/CH5/EX5.2/Ch05Exa2.sce | 18 ++++++++++++++++++ 1970/CH5/EX5.3/Ch05Exa3.sce | 14 ++++++++++++++ 1970/CH5/EX5.4/Ch05Exa4.sce | 12 ++++++++++++ 1970/CH5/EX5.6/Ch05Exa6.sce | 12 ++++++++++++ 1970/CH5/EX5.8/Ch05Exa8.sce | 11 +++++++++++ 1970/CH5/EX5.9/Ch05Exa9.sce | 13 +++++++++++++ 8 files changed, 110 insertions(+) create mode 100755 1970/CH5/EX5.1/Ch05Exa1.sce create mode 100755 1970/CH5/EX5.10/Ch05Exa10.sce create mode 100755 1970/CH5/EX5.2/Ch05Exa2.sce create mode 100755 1970/CH5/EX5.3/Ch05Exa3.sce create mode 100755 1970/CH5/EX5.4/Ch05Exa4.sce create mode 100755 1970/CH5/EX5.6/Ch05Exa6.sce create mode 100755 1970/CH5/EX5.8/Ch05Exa8.sce create mode 100755 1970/CH5/EX5.9/Ch05Exa9.sce (limited to '1970/CH5') 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 -- cgit