|
 |
 |
 |
 |
 |
 |
 |
 |
 |
D.3.2.5 orthogonalize
Procedure from librarylinalg.lib (see linalg_lib).
- Usage:
- orthogonalize(A); A = matrix of constants
- Return:
- matrix, orthogonal basis of the column space of A
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); ==> 5, 6,12,4, ==> 7, 3,2, 6, ==> 12,1,1, 2, ==> 6, 4,2, 10 print(orthogonalize(A)); ==> 1,0, 0, 0, ==> 0,23,0, 0, ==> 0,0, 21,0, ==> 0,0, 0, 6 |