initial commit / add all books

1 files changed, 22 insertions, 0 deletions

diff --git a/1970/CH1/EX1.13/Ch01Exa13.sce b/1970/CH1/EX1.13/Ch01Exa13.sce
new file mode 100755
index 000000000..f98991ada
--- /dev/null
+++ b/1970/CH1/EX1.13/Ch01Exa13.sce
@@ -0,0 +1,22 @@
+// Scilab code Ex1.13 : : P.No.55 (2011)

+// We have to determine for mass numbers 80 and 97.

+clc; clear;

+A = [80, 97];   // Matrix of Mass numbers

+Element = ["Br","Mo"];  // Matrix of elements

+M_n = 939.6;    // Mass of neutron, MeV

+M_H = 938.8;    // Mass of proton, MeV

+a_v = 14.0;     // Volume energy, MeV

+a_s = 13.0;     // Surface energy, MeV

+a_c = 0.583;    // Coulomb energy, MeV

+a_a = 19.3;     // Asymmetry energy, MeV

+a_p = 33.5;     // Pairing energy, MeV

+for i = 1:1:2

+Z = poly(0,'Z');        // Declare the polynomial variable

+M_AZ = M_n*(A(i)-Z)+M_H*Z-a_v*A(i)+a_s*A(i)^(2/3)+a_c*Z*(Z-1)*A(i)^(-1/3)+a_a*(A(i)-2*Z)^2/A(i)+a_p*A(i)^(-3/4); // Mass of the nuclide, MeV/c^2

+Z = roots(derivat(M_AZ));

+printf("\nFor A = %d, the most stable isobar is %s(%d,%d)", A(i), Element(i), Z, A(i));    

+end

+

+// Result

+// For A = 80, the most stable isobar is Br(35,80)

+// For A = 97, the most stable isobar is Mo(42,97)