1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
// Scilab code Ex17.4 : Pg:890 (2011)
clc;clear;
nucleus = cell(4,3);
// For Li nuclides
nucleus(1,1).entries = 'Li';
nucleus(1,2).entries = 3;
nucleus(1,3).entries = 7;
nucleus(2,1).entries = 'Li';
nucleus(2,2).entries = 3;
nucleus(2,3).entries = 8;
// For Be nuclides
nucleus(3,1).entries = 'Be';
nucleus(3,2).entries = 4;
nucleus(3,3).entries = 9;
nucleus(4,1).entries = 'Be';
nucleus(4,2).entries = 4;
nucleus(4,3).entries = 10;
a_c = 0.7053; // Asymmetry energy constant, MeV
a_a = 23.702; // Coulomb energy constant, MeV
for i = 1:1:4
Z = nucleus(i,3).entries/(2+a_c/(2*a_a)*nucleus(i,3).entries^(2/3));
if abs(Z-int(Z)) < 0.5 then
printf("\n%s(%d,%d) is more stable", nucleus(i,1).entries, nucleus(i,2).entries, nucleus(i,3).entries);
end
end
// Result
// Li(3,7) is more stable
// Be(4,9) is more stable
|
