D.3.2.5 orthogonalize

Procedure from library linalg.lib (see linalg_lib).

Usage:

orthogonalize(A); A = constant matrix

Return:

matrix, orthogonal basis of the column space of A

Example:

LIB "linalg.lib";
ring r=0,(x),lp;
matrix A[4][4]=5,6,12,4,7,3,2,6,12,1,1,2,6,4,2,10;
print(A);
print(orthogonalize(A));