D.3.2.2 inverse_B

Procedure from library linalg.lib (see linalg_lib).

Usage:

inverse_B(A); A = square matrix

Return:

list Inv with
- Inv[1] = matrix I and
- Inv[2] = poly p
such that I*A = unitmat(n)*p;

Note:

p=1 if 1/det(A) is computable and p=det(A) if not;
the computation uses busadj.

Example:

LIB "linalg.lib";
ring r=0,(x,y),lp;
matrix A[3][3]=x,y,1,1,x2,y,x,6,0;
print(A);
list Inv=inverse_B(A);
print(Inv[1]);
print(Inv[2]);
print(Inv[1]*A);

See also: inverse; inverse_L.