blob: c633371aa7b8ba7d1612201a6cb71dacfd2b8a35 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
// Scilab code Exa10.8 : : Page-457 (2011)
clc; clear;
N_0 = 6.02252e+026; // Avogadro's constant
rho = 8.9*10^3; // Nuclear density of Co-59, Kg per cubic metre
M = 59; // Mass number
sigma = 30e-028; // Cross section, per square metre
phi = 10^16; // Neutron flux, neutrons per square metre per sec
d = 0.04e-02; // Thickness of Co-59 sheet, metre
t = 3*60*60; // Total reaction time, sec
t_half = 5.2*365*86400; // Half life of Co-60, sec
lambda = 0.693/t_half; // Disintegration constant, per sec
N_nuclei = round(N_0*rho/M*sigma*phi*d*t); // Number of nuclei of Co-60 produced
Init_activity = lambda*N_nuclei; // Initial activity, decays per sec
printf("\nThe number of nuclei of Co60 produced = %5.2e \nThe initial activity per Sq. metre = %1.0g decays per sec", N_nuclei, Init_activity);
// Result
// The number of nuclei of Co60 produced = 1.18e+019
// The initial activity per Sq. metre = 5e+010 decays per sec
|
