<s> '''<STR_LIT>''' <EOL> import sympy <EOL> from sympy import Matrix , Symbol , symbols , zeros , ones , eye <EOL> def numbered_matrix ( nrows , ncols , symbol_prefix ) : <EOL> return Matrix ( nrows , ncols , lambda i , j : Symbol ( <EOL> symbol_prefix + '<STR_LIT>' % ( i + <NUM_LIT:1> , j + <NUM_LIT:1> ) ) ) <EOL> def numbered_vector ( nrows , symbol_prefix ) : <EOL> return Matrix ( nrows , <NUM_LIT:1> , lambda i , j : Symbol ( <EOL> symbol_prefix + '<STR_LIT>' % ( i + <NUM_LIT:1> ) ) ) <EOL> RS = numbered_matrix ( <NUM_LIT:3> , <NUM_LIT:3> , '<STR_LIT>' ) <EOL> mdc , mdr , rdc , rdr = symbols ( <EOL> '<STR_LIT>' , '<STR_LIT>' , '<STR_LIT>' , '<STR_LIT>' ) <EOL> md_adj = Matrix ( ( mdc - <NUM_LIT:1> , mdr - <NUM_LIT:1> , <NUM_LIT:0> ) ) / - <NUM_LIT:2> <EOL> rd_adj = Matrix ( ( rdc - <NUM_LIT:1> , rdr - <NUM_LIT:1> , <NUM_LIT:0> ) ) / - <NUM_LIT:2> <EOL> adj = - ( RS * md_adj ) + RS * rd_adj <EOL> adj . simplify ( ) <EOL> Q = RS [ : , : <NUM_LIT:2> ] * Matrix ( ( <EOL> ( mdc - rdc ) / <NUM_LIT:2> , <EOL> ( mdr - rdr ) / <NUM_LIT:2> ) ) <EOL> Q . simplify ( ) <EOL> assert adj == Q </s>
