summaryrefslogtreecommitdiff
path: root/1970/CH13/EX13.6/CH13Exa6.sce
diff options
context:
space:
mode:
Diffstat (limited to '1970/CH13/EX13.6/CH13Exa6.sce')
-rwxr-xr-x1970/CH13/EX13.6/CH13Exa6.sce18
1 files changed, 18 insertions, 0 deletions
diff --git a/1970/CH13/EX13.6/CH13Exa6.sce b/1970/CH13/EX13.6/CH13Exa6.sce
new file mode 100755
index 000000000..2047de159
--- /dev/null
+++ b/1970/CH13/EX13.6/CH13Exa6.sce
@@ -0,0 +1,18 @@
+// Scilab code Exa13.6 : : Page-601 (2011)
+clc; clear;
+r = 1/2; // Radius of the tube, metre
+a = %pi*r^2; // Area of the torus, square metre
+V = 3*%pi*a; // Volume of the torus, cubic metre
+P = 10^-5*13.6e+3*9.81; // Pressure of the gas, newton per square metre
+C = 1200e-6; // Capacitance, farad
+v = 4e+4; // potential, volts
+T_room = 293; // Room temperature, kelvin
+N_k = P*V/T_room; // From gas equation
+E = 1/2*C*v^2; // Energy stored, joules
+T_k = 1/6*E/(N_k*10); // Temperature attained by thermonuclear device, kelvin
+printf("\nThe temperature attained by thermonuclear device = %4.2e K", T_k);
+
+// Result
+// The temperature attained by thermonuclear device = 4.75e+005 K
+
+ \ No newline at end of file