samer@4: function B1=orthogonalise(B),
samer@4: % orthogonalise - Orthogonalise a basis matrix
samer@4: %
samer@4: % orthogonalise :: [[N,M]] -> [[N,M]].
samer@4: %
samer@4: % Works using SVD.
samer@4: 
samer@4: [U,S,V] = svd(B,0);
samer@4: B1 = U*V';
samer@4: 
samer@4: % alternative method, seems to be slower
samer@4: % B1 = B*real((B'*B)^(-0.5));
samer@4: