From fe6e959ae32a7e0b7db574d18b3df40de82b16bc Mon Sep 17 00:00:00 2001
From: Mavroudo <mavroudo@csd.auth.gr>
Date: Sun, 11 Feb 2018 12:18:37 +0200
Subject: [PATCH 1/2] Setting the files in the same folder

---
 Bisection.py => ArithmeticAnalysis/Bisection.py       | 0
 Intersection.py => ArithmeticAnalysis/Intersection.py | 0
 NeutonMethod.py => ArithmeticAnalysis/NeutonMethod.py | 0
 3 files changed, 0 insertions(+), 0 deletions(-)
 rename Bisection.py => ArithmeticAnalysis/Bisection.py (100%)
 rename Intersection.py => ArithmeticAnalysis/Intersection.py (100%)
 rename NeutonMethod.py => ArithmeticAnalysis/NeutonMethod.py (100%)

diff --git a/Bisection.py b/ArithmeticAnalysis/Bisection.py
similarity index 100%
rename from Bisection.py
rename to ArithmeticAnalysis/Bisection.py
diff --git a/Intersection.py b/ArithmeticAnalysis/Intersection.py
similarity index 100%
rename from Intersection.py
rename to ArithmeticAnalysis/Intersection.py
diff --git a/NeutonMethod.py b/ArithmeticAnalysis/NeutonMethod.py
similarity index 100%
rename from NeutonMethod.py
rename to ArithmeticAnalysis/NeutonMethod.py

From a48836998b46131ec7698995e0603fbf31eefd55 Mon Sep 17 00:00:00 2001
From: Mavroudo <mavroudo@csd.auth.gr>
Date: Sun, 11 Feb 2018 21:43:38 +0200
Subject: [PATCH 2/2] LU decomposition

---
 ArithmeticAnalysis/LUdecomposition.py | 34 +++++++++++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100644 ArithmeticAnalysis/LUdecomposition.py

diff --git a/ArithmeticAnalysis/LUdecomposition.py b/ArithmeticAnalysis/LUdecomposition.py
new file mode 100644
index 00000000..bb0827bb
--- /dev/null
+++ b/ArithmeticAnalysis/LUdecomposition.py
@@ -0,0 +1,34 @@
+import math
+import numpy
+
+def LUDecompose (table): #table that contains our data
+    #table has to be a square array so we need to check first
+    rows,columns=numpy.shape(table)
+    L=numpy.zeros((rows,columns))
+    U=numpy.zeros((rows,columns))
+    if rows!=columns:
+        return
+    for i in range (columns):
+        for j in range(i-1):
+            sum=0
+            for k in range (j-1):
+                sum+=L[i][k]*U[k][j]
+            L[i][j]=(table[i][j]-sum)/U[j][j]
+        L[i][i]=1
+        for j in range(i-1,columns):
+            sum1=0
+            for k in range(i-1):
+                sum1+=L[i][k]*U[k][j]
+            U[i][j]=table[i][j]-sum1
+    return L,U
+
+
+
+
+
+
+
+matrix =numpy.array([[2,-2,1],[0,1,2],[5,3,1]])
+L,U = LUDecompose(matrix)
+print(L)
+print(U)
\ No newline at end of file