diff options
Diffstat (limited to '1970/CH1/EX1.11')
| -rwxr-xr-x | 1970/CH1/EX1.11/Ch01Exa11.sce | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/1970/CH1/EX1.11/Ch01Exa11.sce b/1970/CH1/EX1.11/Ch01Exa11.sce new file mode 100755 index 000000000..5185a4f1b --- /dev/null +++ b/1970/CH1/EX1.11/Ch01Exa11.sce @@ -0,0 +1,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
|