hw06

week06/hw06.py

+
+59110440204 นางสาวนฤมล   โมคภา
+
+def Ab2Uc2(A, b):
+	
+	import numpy as np
+
+	A = np.array(A, float)
+	b = np.array(b, float)
+	n = len(b)
+
+	for k in range(0,n-1):
+    
+    		for i in range(k+1,n):
+        	lam = A[i,k]/A[k,k]
+        	A[i,k] = 0
+        	A[i,k+1:n] = A[i,k+1:n] - lam*A[k,k+1:n]
+        	b[i] -= lam*b[k]
+        
+        	U = A.copy()
+        	c = b.copy()
+        
+	print("U = ", U)
+	print("c = ", c)
+