1 #!/usr/bin/env python3 2 3 import sys 4 import numpy as np 5 import scipy.linalg as la 6 7 n = 3 8 9 if sys.argv.__len__()!=n*n+1: 10 print( "bad # args" ) 11 sys.exit() 12 13 14 A = np.empty([n,n], dtype='f8') 15 16 for j in range(n): 17 for i in range(n): 18 A[j][i] = sys.argv[1+j*n+i] 19 20 print( "A:" ) 21 print( A ) 22 23 Q,R,P = la.qr(A, pivoting=True) 24 25 print( "Q:" ) 26 print( Q ) 27 28 print( "R:" ) 29 print( R ) 30 31 print( "P:" ) 32 print( P ) 33
