5 files changed, 72 insertions, 0 deletions
diff --git a/3535/CH4/EX4.1/Ex4_1.sce b/3535/CH4/EX4.1/Ex4_1.sce
new file mode 100644
index 000000000..069bd8e25
--- /dev/null
+++ b/3535/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,12 @@
+//Chapter 4, Example 4.1, Page 86
+clc
+clear
+//Binding energy
+mn = 1.0078250
+mp = 1.0086649
+M = 4.0026032 // mass of He
+MD = 2*mn+2*mp-M //Mass defect
+BE = MD*931.5
+printf("\n Mass defect = %f u",MD);
+printf("\n Nuclear binding energy = %f MeV",BE); // answer provided in the textbook is wrong
+//Answer may vary due to round off error
diff --git a/3535/CH4/EX4.2/Ex4_2.sce b/3535/CH4/EX4.2/Ex4_2.sce
new file mode 100644
index 000000000..bbeb86040
--- /dev/null
+++ b/3535/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,11 @@
+//Chapter 4, Example 4.2, Page 89
+clc
+clear
+//Binding energy
+O15 = 15.0030654 // atomic mass of O15 isotope
+mn = 1.00866492
+O16 = 15.9949146 // atomic mass of O16 isotope
+c2 = 931.5 // C^2 in MeV
+S = (O15+mn-O16)*c2
+printf("\n Binding energy = %f MeV",S); 
+//Answer may vary due to round off error
diff --git a/3535/CH4/EX4.3/Ex4_3.sce b/3535/CH4/EX4.3/Ex4_3.sce
new file mode 100644
index 000000000..c5a1ed738
--- /dev/null
+++ b/3535/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,26 @@
+//Chapter 4, Example 4.3, Page 94
+clc
+clear
+// Q value of an endothermic and exothermic reaction
+//Exothermic reaction
+Be = 9.012182 //Reactants
+He = 4.002603 //Reactants
+C12 = 12 //Product
+n = 1.008664 //Product
+C2 = 931.5 // C^2 in MeV
+Exo1 = Be+He
+Exo2 = C12+n
+Dif1 = Exo1-Exo2
+Q1 = Dif1*C2
+printf("\n Q of the exothermic reaction = %f MeV",Q1); 
+//Endothermic reaction
+O = 15.994915 //Reactants
+n = 1.008664 //Reactant
+C13 = 13.003354 //Product
+He = 4.002603 //product
+End1 = O+n
+End2 = C13+He
+Dif2 = End1-End2
+Q2 = Dif2*C2
+printf("\n Q of the exothermic reaction = %f MeV",Q2); 
+//Answer may vary due to round off error 
diff --git a/3535/CH4/EX4.4/Ex4_4.sce b/3535/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..4e714fb75
--- /dev/null
+++ b/3535/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,11 @@
+//Chapter 4, Example 4.4, Page 95
+clc
+clear
+// Q value in a reaction
+MH = 1.00782503
+MD = 2.01410178
+me = 0.00054858
+C2 = 931.5
+Q = (2*MH-MD-2*me)*C2
+printf("\n Q of the reaction = %f MeV",Q);// Answer provided in the text is wrong 
+//Answer may vary due to round off error 
diff --git a/3535/CH4/EX4.5/Ex4_5.sce b/3535/CH4/EX4.5/Ex4_5.sce
new file mode 100644
index 000000000..f1f26d0ee
--- /dev/null
+++ b/3535/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,12 @@
+//Chapter 4, Example 4.5, Page 96
+clc
+clear
+// Q value of the reaction
+mn = 1.0086649
+MB = 10.0129370
+MHe = 4.0026032
+MLi = 7.0160040
+C2 = 931.5
+Q = (mn+MB-MHe-MLi)*C2 -0.48
+printf("\n Q of the reaction = %f MeV",Q);
+//Answer may vary due to round off error