1 files changed, 25 insertions, 0 deletions
diff --git a/767/CH4/EX4.8.3/Ch04Exa4_8_3.sci b/767/CH4/EX4.8.3/Ch04Exa4_8_3.sci
new file mode 100755
index 000000000..6ed557a61
--- /dev/null
+++ b/767/CH4/EX4.8.3/Ch04Exa4_8_3.sci
@@ -0,0 +1,25 @@
+// Scilab code Exa4.8.3: To calculate the energy released and temperature required for fusion of given gases : Page 194 (2011)
+// Firstly calculate for B-10
+Z_B = 5; // Atomic number of B-10
+r_B = 5.17; //  Seperation of two nuclei, fm
+K = 1.38e-023; // Boltzmann's constant
+F = 1/137; // Fine structure constant
+E = 197.5*1.6e-013; // Energy, J
+V_c_B = F*Z_B^2*E/r_B; // Coulomb barrier for B-10, J
+T_B = 2/3*V_c_B/K; // Temperature required to overcome the barrier for B-10, K
+// Now calculate for Mg-24
+Z_Mg = 12; // Atomic number of Mg-24
+r_Mg = 6.92; //  Seperation of two nuclei, fm
+K = 1.38e-023; // Boltzmann's constant
+F = 1/137; // Fine structure constant
+E = 197.5*1.6e-013; // Energy, J
+V_c_Mg = F*Z_Mg^2*E/r_Mg; // Coulomb barrier for Mg-24, J
+T_Mg = 2/3*V_c_Mg/K; // Temperature required to overcome the barrier for Mg-24, K
+printf("\nFor B-10 \n Energy released   = %4.2e J \n Temperature required   = %4.1e K   \nFor Mg-24  \n Energy released    = %4.2e J \n Temperature required  =  %4.2e K", V_c_B,T_B, V_c_Mg,  T_Mg)
+// Result  
+ //      For B-10 
+ // Energy released   = 1.12e-012 J 
+ // Temperature required   = 5.4e+010 K   
+// For Mg-24  
+//  Energy released    = 4.80e-012 J 
+//  Temperature required  =  2.32e+011 K