summaryrefslogtreecommitdiff
path: root/2411/CH9/EX9.2.13/Ex9_2_13.sce
blob: 80763a6e41f4179f81e2bb116a2873a476191a47 (plain)
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
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
// Scilab Code Ex9.2.13: Page-418(2008)
clc; clear;
function p = Find(Z, A)
    if Z == 2 & A == 4 then
        p = 'alpha';
    end
    if Z == -1 & A == 0 then
        p = 'beta-';
    end
    if  Z == 1 & A == 0 then
        p = 'beta+';
    end
endfunction   
R1 = cell(4,3);
R2 = cell(4,3);
// Enter data for first cell (Reaction)
R1(1,1).entries = 'Li'; // Element
R1(1,2).entries = 3;   // Atomic number
R1(1,3).entries = 6;    // Mass number
R1(2,1).entries = 'd';
R1(2,2).entries = 1;
R1(2,3).entries = 2;
R1(3,1).entries = 'X';
R1(3,2).entries = 0;
R1(3,3).entries = 0;
R1(4,1).entries = 'He';
R1(4,2).entries = 2;
R1(4,3).entries = 4;
// Enter data for second cell (Reaction)
R2(1,1).entries = "Te";
R2(1,2).entries = 52;
R2(1,3).entries = 122;
R2(2,1).entries = 'X';
R2(2,2).entries = 0;
R2(2,3).entries = 0;
R2(3,1).entries = 'I';
R2(3,2).entries = 53;
R2(3,3).entries = 124;
R2(4,1).entries = 'd';
R2(4,2).entries = 1;
R2(4,3).entries = 2;
R1(3,2).entries = R1(1,2).entries+R1(2,2).entries-R1(4,2).entries
R1(3,3).entries = R1(1,3).entries+R1(2,3).entries-R1(4,3).entries
particle = Find(R1(3,2).entries, R1(3,3).entries);    // Find the unknown particle
printf("\nFor the reaction\n")
            printf("\t%s(%d) + %s(%d) --> %s + %s(%d)\n X must be an %s particle", R1(1,1).entries, R1(1,3).entries, R1(2,1).entries, R1(2,3).entries, R1(3,1).entries, R1(4,1).entries, R1(4,3).entries, particle);
R2(2,2).entries = R2(3,2).entries+R2(4,2).entries-R2(1,2).entries
R2(2,3).entries = R2(3,3).entries+R2(4,3).entries-R2(1,3).entries
particle = Find(R2(2,2).entries, R2(2,3).entries);    // Find the unknown particle
printf("\n\nFor the reaction\n")
            printf("\t%s(%d) + %s --> %s(%d)+%s(%d)\n X must be an %s particle", R2(1,1).entries, R2(1,3).entries, R2(2,1).entries, R2(3,1).entries, R2(3,3).entries, R2(4,1).entries, R2(4,3).entries, particle);
            
// Result
// For the reaction
//	Li(6) + d(2) --> X + He(4)
//  X must be an alpha particle

// For the reaction
// 	Te(122) + X --> I(124)+d(2)
//  X must be an alpha particle