| 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_L(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
|