From 82e6cf7d1046d6cee16f7e8b044ec33e7ec6c4b7 Mon Sep 17 00:00:00 2001
From: Nguyễn Gia Phong <vn.mcsinyx@gmail.com>
Date: Sun, 16 Feb 2020 14:26:55 +0700
Subject: [usth] Numerical Method is MATH2.4

---
 usth/MATH2.4/labwork/3/LU.m | 10 ++++++++++
 1 file changed, 10 insertions(+)
 create mode 100644 usth/MATH2.4/labwork/3/LU.m

(limited to 'usth/MATH2.4/labwork/3/LU.m')

diff --git a/usth/MATH2.4/labwork/3/LU.m b/usth/MATH2.4/labwork/3/LU.m
new file mode 100644
index 0000000..860486d
--- /dev/null
+++ b/usth/MATH2.4/labwork/3/LU.m
@@ -0,0 +1,10 @@
+function [L U] = LU (A)
+  [n n1] = size (A);
+  [L U] = deal (eye (n), A);
+  for k = 1:n
+    for i = k + 1 : n
+      L(i, k) = U(i, k) / U(k, k);
+      U(i, :) = U(i, :) - L(i, k)*U(k, :);
+    end
+  end
+end
-- 
cgit 1.4.1