import numpy as np

A = np.matrix( '-2,  2, -5,     ?,  ?,  ?;\
                 2, -3,  7,     ?,  ?,  ?;\
                -4,  3, -7,     ?,  ?,  ?' )

print( 'A = ' )
print( A )

G0 = np.matrix( ' 1,  0,  0;\
                  ?,  1,  0;\
                  ?,  0,  1' );

A0 = G0 * A

print( 'A0 = ' )
print( A0 )

print( 'A0 =' )
print( A0 )

G1 = np.matrix( ' 1,  ?,  0;\
                  0,  1,  0;\
                  0,  ?,  1' );

A1 = G1 * A0

print( 'A1 = ' )
print( A1 )

print( 'A1 =' )
print( A1 )

G2 = np.matrix( ' 1,  0,  ?;\
                  0,  1,  ?;\
                  0,  0,  1' );

A2 = G2 * A1

print( 'A2 = ' )
print( A2 )

print( 'A2 =' )
print( A2 )

D = np.matrix( '  ?,  0,  0;\
                  0,  ?,  0;\
                  0,  0,  ?' );

A3 = D * A2

print( 'A3 = ' )
print( A3 )

B = A3[ 0:3, 3:6 ]

print( 'B = ' )
print( B )

A = np.matrix( '-2,  2, -5;\
                 2, -3,  7;\
                -4,  3, -7' )

print( 'A = ' )
print( A )

print( 'B = ' )
print( B )

print( ' A * B = ' )
print( A * B )