From e461df7573c2b7b7e26c965d8cf2d8e175d67378 Mon Sep 17 00:00:00 2001
From: Nguyễn Gia Phong <vn.mcsinyx@gmail.com>
Date: Mon, 16 Dec 2019 21:31:18 +0700
Subject: [usth/MATH2.2] Numerical Methods

The future starts now.
---
 usth/MATH2.2/project/script | 22 ++++++++++++++++++++++
 1 file changed, 22 insertions(+)
 create mode 100755 usth/MATH2.2/project/script

(limited to 'usth/MATH2.2/project/script')

diff --git a/usth/MATH2.2/project/script b/usth/MATH2.2/project/script
new file mode 100755
index 0000000..99b61f7
--- /dev/null
+++ b/usth/MATH2.2/project/script
@@ -0,0 +1,22 @@
+#!/usr/bin/env octave
+[cp lambda rho] = deal (900, 237, 2700);
+[Tg Td T0] = deal (100, 50, 0);
+[L dx dt N] = deal (0.4, 0.01, 0.5, 841);
+T = heatrans (cp, lambda, rho, Tg, Td, T0, L, dx, dt, N);
+[X Y Z] = deal (0:dx:L, 0:dt:420, T');
+
+figure (1);   # mesh
+mesh (X, Y, Z);
+title ("Temperature of the bar during the first 420 seconds");
+xlabel ("x (m)");
+ylabel ("t (s)");
+zlabel ("T (°C)");
+
+figure (2);   # contour
+contourf (X, Y, Z, 0:7:100, "showtext", "on");
+title ("Temperature of the bar during the first 420 seconds");
+xlabel ("x (m)");
+ylabel ("t (s)");
+
+disp ("Press any key to coninue...");
+kbhit;
-- 
cgit 1.4.1