blob: edb0407845f49032e7e62f6c7b0c7e2e301cfcec (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
// Scilab code Exa10.13 : : Page-459 (2011)
clc; clear;
h = 6.625e-34; // Planck's constant, joule sec
m_n = 1.67e-27; // Mass of neutron, Kg
E = 4.906; // Energy, joule
w_y = 0.124; // radiation width, eV
w_n = 0.007*E^(1/2); // Probability of elastic emission of neutron, eV
I = 3; // Total angular momentum
I_c = 2; // Total angular momentum in the compound state
sigma = ((h^2)*(2*I_c+1)*w_y*w_n)*10^28/(2*%pi*m_n*E*1.602e-019*(2*I+1)*(w_y+w_n)^2); // Cross section, barns
printf("\nThe cross section of neutron capture = %5.3e barns", sigma);
// Result
// The cross section of neutron capture = 3.755e+004 barns
|
