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 Ex1.11 : : Page 55 (2011)
clc; clear;
M_n = 1.008665; // Mass of neutron, amu
M_p = 1.007825; // Mass of proton, amu
N_Ni = 36; // Number of neutron in Ni-64
Z_Ni = 28; // Atomic number of Ni-64
N_Cu = 35; // Number of neutron in Cu-64
Z_Cu = 29; // Atomic number of Cu-64
A = 64; // Mass number, amu
M_Ni = 63.927958; // Mass of Ni-64
M_Cu = 63.929759; // Mass of Cu-64
m_e = 0.511; // Mass of an electron, MeV
d_M_Ni = N_Ni*M_n+Z_Ni*M_p-M_Ni; // Mass defect, amu
d_M_Cu = N_Cu*M_n+Z_Cu*M_p-M_Cu; // Mass defect, amu
B_E_Ni = d_M_Ni*931.49; // Binding energy of Ni-64, MeV
B_E_Cu = d_M_Cu*931.49; // Binding energy of Cu-64, MeV
Av_B_E_Ni = B_E_Ni/A; // Average binding energy of Ni-64, MeV
Av_B_E_Cu = B_E_Cu/A; // Average binding energy of Cu-64, MeV
printf("\nBinding energy of Ni-64 : %7.3f MeV \nBinding energy of CU-64 : %7.3f MeV \nAverage binding energy of Ni-64 : %5.3f MeV \nAverage binding energy of Cu-64 : %5.3f MeV ", B_E_Ni, B_E_Cu, Av_B_E_Ni, Av_B_E_Cu);
if (M_Cu - M_Ni)*931.48 > 2*m_e then
printf("\nNi-64 is not a beta emitter but Cu-64 is a beta emitter");
end
// Result
// Binding energy of Ni-64 : 561.765 MeV
// Binding energy of CU-64 : 559.305 MeV
// Average binding energy of Ni-64 : 8.778 MeV
// Average binding energy of Cu-64 : 8.739 MeV
// Ni-64 is not a beta emitter but Cu-64 is a beta emitter
|
