add

week05/back.py

-import numpy as np
-
-def upperTriSol(A, b):
-    n = np.size(b)
-    x = np.zeros_like(b)
-
-    x[-1] = 1. / A[-1, -1] * b[-1]
-    for i in xrange(n-2, -1, -1):
-        x[i] = 1. / A[i, i] * (b[i] - np.sum(A[i, i+1:] * x[i+1:]))
-
-    return x

week05/ge.py

-def ge(A, b):
-    pass
-
-def solve(A, b):
-    import numpy.linalg
-    result = numpy.linalg.solve(A,b)
-    for i,x in zip(1,range(len(result)+1),result):
-        print('$x_{} = {}'.format(i,x))
-    
-m,n,_ = map(int,input().split('x'))
-print(m, n)
-A= []
-for row in range(m):
-    A.append(list(map(int,input().split(','))))
-#solve(A,b)
-print(n,_)
-b=list(map(int,input().split(',')))
-solve(A,b)
-

week05/readme.md

+เขียนฟังก์ชันเพื่อหา upper-triangular matrix (U) จากสมการ Ax = b 
+โดยเขียนฟังก์ชัน Ux = c เป็นสมการที่ลดรูปแล้ว