From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 767/CH6/EX6.3.2/Ch06Exa6_3_2.sci | 52 ++++++++++++++++++++++++++++++++++++++++ 1 file changed, 52 insertions(+) create mode 100755 767/CH6/EX6.3.2/Ch06Exa6_3_2.sci (limited to '767/CH6/EX6.3.2') diff --git a/767/CH6/EX6.3.2/Ch06Exa6_3_2.sci b/767/CH6/EX6.3.2/Ch06Exa6_3_2.sci new file mode 100755 index 000000000..13419dcba --- /dev/null +++ b/767/CH6/EX6.3.2/Ch06Exa6_3_2.sci @@ -0,0 +1,52 @@ +// Scilab code Exa6.3.2 : To calculate the kinetic energy of protons and no. of possibile reactions: Page 265 (2011) +V = 5; // Voltage of accelerator, MV +// Declare three cells (for three reactions): Page no. : 133(2011) + R1 = cell(3,2) + R2 = cell(10,2) +// Enter data for first cell (Reaction) +R1(1,1).entries = "p"; +R1(1,2).entries = 1; +R1(2,1).entries = 'd'; +R1(2,2).entries = 1; +R1(3,1).entries = 'He'; +R1(3,2).entries = 2; +E_p = (R1(1,2).entries)*V +E_d = (R1(2,2).entries)*V +E_He = (R1(3,2).entries)*V + // Enter data for second cell (Reaction) + R2(1,1).entries = "p" + R2(1,2).entries = 1 + R2(2,1).entries = "N" + R2(2,2).entries = 14 + R2(3,1).entries = "O" + R2(3,2).entries = 15 + R2(4,1).entries = "y" + R2(4,2).entries = 0 + R2(5,1).entries = "d" + R2(5,2).entries = 1 + R2(6,1).entries = "n" + R2(6,2).entries = 0 + R2(7,1).entries = "He" + R2(7,2).entries = 3 + R2(8,1).entries = "C" + R2(8,2).entries = 13 + R2(9,1).entries = "He" + R2(9,2).entries = 4 + R2(10,1).entries = "C" + R2(10,2).entries = 12 + printf("\nProtons energy = -%d MeV \n Deuterons energy = -%d MeV \n Double charged He-3 = -%d MeV", E_p, E_d, E_He) + printf("\n Possible reaction at these energies are") + printf("\n %s + %s(%d) ---> %s(%d)+ %s", R2(1,1).entries,R2(2,1).entries,R2(2,2).entries,R2(3,1).entries,R2(3,2).entries,R2(4,1).entries) +printf("\n %s + %s(%d) ---> %s(%d) + %s ", R2(5,1).entries,R2(2,1).entries,R2(2,2).entries,R2(3,1).entries,R2(3,2).entries,R2(6,1).entries) +printf("\n %s(%d) +%s(%d) ---> %s(%d)+ %s", R2(7,1).entries,R2(7,2).entries,R2(8,1).entries,R2(8,2).entries,R2(3,1).entries,R2(3,2).entries,R2(6,1).entries) + printf("\n %s(%d) + %s(%d) ---> %s(%d) +%s", R2(9,1).entries,R2(9,2).entries,R2(10,1).entries,R2(10,2).entries,R2(3,1).entries,R2(3,2).entries,R2(6,1).entries) + +// Result +// Protons energy = -5 MeV +// Deuterons energy = -5 MeV +// Double charged He-3 = -10 MeV +// Possible reaction at these energies are +// p + N(14) ---> O(15)+ y +// d + N(14) ---> O(15) + n +// He(3) +C(13) ---> O(15)+ n +// He(4) + C(12) ---> O(15) +n -- cgit