blob: a6fc9dd8dd427e0640e8c4fbcdacba2faf27f4bc (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
// Scilab Code Ex9.1.3:Page-411 (2008)
clc; clear;
u = 931.508; // Energy equivalent of 1 amu, MeV
Z = 8; // Atomic number of O-16
A = 16; // Mass number of O-16
m_p = 1.008142; // Mass of a proton, u
m_n = 1.008982; // Mass of a neutron, u
M_O = 15.994915; // Atomic mass of O-16 nucleus, u
delta_m = Z*m_p + (A-Z)*m_n - M_O; // Mass difference, u
BE = delta_m*u; // Binding energy of O-16 nucleus, MeV
BE_bar = BE/A; // Binding energy per nucleon of O-16 nucleus, MeV
delta_m = abs(M_O - A); // Mass difference, u
PF = delta_m/A; // Packing fraction for O-16 nucleus, u
printf("\nThe binding energy per nucleon for O-16 nucleus = %4.2f MeV/nucleon", BE_bar);
printf("\nThe packing fraction for O-16 nucleus = %5.3e u", PF);
// Result
// The binding energy per nucleon for O-16 nucleus = 8.27 MeV/nucleon
// The packing fraction for O-16 nucleus = 3.178e-004 u
|
