1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
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)
|
