blob: 19ebee71767343876fd9db4cf2b41460c169ff01 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
|
// Scilab code: Ex11.2: Threshold energy for the reaction: Pg:229 (2008)
mu = 931.5; // Energy equivalent of 1 amu, MeV
mx = 1.008665; // Mass of neutron, amu
Mx = 13.003355; // Mass of carbon atom, amu
M_alpha = 4.002603; // Mass of alpha particle, amu
M_Be = 10.013534; // Mass of beryllium, amu
MD = (Mx + mx - M_Be - M_alpha); // Mass defect of the reaction, amu
Q = MD*mu; // Q-value of the nuclear reaction, MeV
E_th = -Q*(1 + mx/Mx); // Threshold energy for the reaction in the laboratory, MeV
printf("\nThe threshold energy of the reaction is = %4.2f MeV", E_th);
// Result
// The threshold energy of the reaction is = 4.13 MeV
|