import numpy as np
import laff
import flame

L = np.matrix( ' 1, 0, 0. 0;\
                -2, 1, 0, 0;\
                 1,-3, 1, 0;\
                 2, 3,-1, 1' )

U = np.matrix( ' 2,-1, 3,-2;\
                 0,-2, 1,-1;\
                 0, 0, 1, 2;\
                 0, 0, 0, 3' )

A = L * U
Aold = np.matrix( np.copy( A ) ) 
B = np.matrix( np.eye( 4 ) )

print( 'A = ' )
print( A )

print( 'B = ' )
print( B )


# insert code here

GJ_Inverse_alt( A, B )

print( A )
print( B )

print( Aold * B )