|  | LIB "linalg.lib";
ring r=0,(x,y),lp;
matrix A[3][3]=x,y,1,1,x2,y,x,6,0;
print(A);
==> x,y, 1,
==> 1,x2,y,
==> x,6, 0 
list Inv=inverse_B(A);
print(Inv[1]);
==> 6y,  -6,    x2-y2,
==> -xy, x,     xy-1, 
==> x3-6,-xy+6x,-x3+y 
print(Inv[2]);
==> x3-xy2+6xy-6
print(Inv[1]*A);
==> x3-xy2+6xy-6,0,           0,          
==> 0,           x3-xy2+6xy-6,0,          
==> 0,           0,           x3-xy2+6xy-6
 |